Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion docs/Project.toml
Original file line number Diff line number Diff line change
Expand Up @@ -159,7 +159,7 @@ IterTools = "1"
Juniper = "0.9"
Lux = "1"
MLUtils = "0.4.4"
Manifolds = "0.10"
Manifolds = "0.11"
Manopt = "0.5"
ModelingToolkit = "10.23, 11"
NLPModels = "0.21, 0.22"
Expand Down
2 changes: 1 addition & 1 deletion lib/OptimizationSpeedMapping/Project.toml
Original file line number Diff line number Diff line change
Expand Up @@ -21,7 +21,7 @@ Pkg = "1"
Reexport = "1.2"
SafeTestsets = "0.1"
SciMLBase = "2.122.1, 3"
SpeedMapping = "0.3"
SpeedMapping = "0.4"
Test = "1.10"
julia = "1.10"

Expand Down
46 changes: 33 additions & 13 deletions lib/OptimizationSpeedMapping/src/OptimizationSpeedMapping.jl
Original file line number Diff line number Diff line change
Expand Up @@ -13,13 +13,22 @@ SciMLBase.allowscallback(::SpeedMappingOpt) = false
SciMLBase.has_init(opt::SpeedMappingOpt) = true
SciMLBase.requiresgradient(opt::SpeedMappingOpt) = true

function _retcode(opt_res)
status = opt_res.status
status === :first_order && return SciMLBase.ReturnCode.Success
status === :max_time && return SciMLBase.ReturnCode.MaxTime
status in (:max_iter, :max_eval) && return SciMLBase.ReturnCode.MaxIters
return SciMLBase.ReturnCode.Failure
end

function __map_optimizer_args(
cache::OptimizationBase.OptimizationCache, opt::SpeedMappingOpt;
cache::OptimizationBase.OptimizationCache,
opt::SpeedMappingOpt;
callback = nothing,
maxiters::Union{Number, Nothing} = nothing,
maxtime::Union{Number, Nothing} = nothing,
abstol::Union{Number, Nothing} = nothing,
reltol::Union{Number, Nothing} = nothing
reltol::Union{Number, Nothing} = nothing,
)

# add optimiser options from kwargs
Expand All @@ -37,14 +46,16 @@ function __map_optimizer_args(
if !isnothing(abstol)
@SciMLMessage(
lazy"common abstol is currently not used by $(opt)",
cache.verbose, :unsupported_kwargs
cache.verbose,
:unsupported_kwargs
)
end

if !isnothing(reltol)
@SciMLMessage(
lazy"common reltol is currently not used by $(opt)",
cache.verbose, :unsupported_kwargs
cache.verbose,
:unsupported_kwargs
)
end

Expand All @@ -66,26 +77,35 @@ function SciMLBase.__solve(cache::OptimizationCache{O}) where {O <: SpeedMapping
maxiters = OptimizationBase._check_and_convert_maxiters(cache.solver_args.maxiters)
maxtime = OptimizationBase._check_and_convert_maxtime(cache.solver_args.maxtime)
opt_args = __map_optimizer_args(
cache, cache.opt, maxiters = maxiters,
cache,
cache.opt,
maxiters = maxiters,
maxtime = maxtime,
abstol = cache.solver_args.abstol,
reltol = cache.solver_args.reltol; cache.solver_args...
reltol = cache.solver_args.reltol;
cache.solver_args...,
)

t0 = time()
opt_res = SpeedMapping.speedmapping(
cache.u0; f = _loss, (g!) = cache.f.grad,
cache.u0;
f = _loss,
(g!) = cache.f.grad,
lower = cache.lb,
upper = cache.ub, opt_args...
upper = cache.ub,
opt_args...,
)
t1 = time()
opt_ret = opt_res.converged ? SciMLBase.ReturnCode.Success :
SciMLBase.ReturnCode.Failure
opt_ret = _retcode(opt_res)
stats = OptimizationBase.OptimizationStats(; time = t1 - t0)
return SciMLBase.build_solution(
cache, cache.opt,
opt_res.minimizer, _loss(opt_res.minimizer);
original = opt_res, retcode = opt_ret, stats = stats
cache,
cache.opt,
opt_res.minimizer,
_loss(opt_res.minimizer);
original = opt_res,
retcode = opt_ret,
stats = stats,
)
end

Expand Down
17 changes: 16 additions & 1 deletion lib/OptimizationSpeedMapping/test/core_tests.jl
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
using OptimizationSpeedMapping, OptimizationBase, ForwardDiff
using OptimizationSpeedMapping, OptimizationBase, ForwardDiff, SciMLBase
using Test

@testset "OptimizationSpeedMapping.jl" begin
Expand All @@ -10,6 +10,8 @@ using Test
prob = OptimizationProblem(f, x0, _p)
sol = solve(prob, SpeedMappingOpt())
@test 10 * sol.objective < l1
@test SciMLBase.successful_retcode(sol)
@test hasproperty(sol.original, :status)

prob = OptimizationProblem(f, x0, _p; lb = [-1.0, -1.0], ub = [1.5, 1.5])
sol = solve(prob, SpeedMappingOpt())
Expand Down Expand Up @@ -38,4 +40,17 @@ using Test
sol = OptimizationBase.solve!(cache)
@test sol.u ≈ [2.0] atol = 1.0e-3
end

@testset "return code mapping" begin
@test OptimizationSpeedMapping._retcode((status = :first_order,)) ==
SciMLBase.ReturnCode.Success
@test OptimizationSpeedMapping._retcode((status = :max_time,)) ==
SciMLBase.ReturnCode.MaxTime
@test OptimizationSpeedMapping._retcode((status = :max_iter,)) ==
SciMLBase.ReturnCode.MaxIters
@test OptimizationSpeedMapping._retcode((status = :max_eval,)) ==
SciMLBase.ReturnCode.MaxIters
@test OptimizationSpeedMapping._retcode((status = :failure,)) ==
SciMLBase.ReturnCode.Failure
end
end
Loading