Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

gMAM HeymannVandenEijnden method broken #26

Open
oameye opened this issue Feb 23, 2024 · 1 comment
Open

gMAM HeymannVandenEijnden method broken #26

oameye opened this issue Feb 23, 2024 · 1 comment
Labels
bug Something isn't working enhancement New feature or request help wanted Extra attention is needed

Comments

@oameye
Copy link
Member

oameye commented Feb 23, 2024

the following code

using CriticalTransitions, Test

function meier_stein(u, p, t) # out-of-place
    x, y = u
    dx = x - x^3 - 10 * x * y^2
    dy = -(1 + x^2) * y
    [dx, dy]
end
σ = 0.25
sys = StochSystem(meier_stein, [], zeros(2), σ, idfunc, nothing, I(2), "WhiteGauss")

# initial path: parabola
xx = range(-1.0, 1.0, length = 30)
yy = 0.3 .* (-xx .^ 2 .+ 1)
init = Matrix([xx yy]')

method = CriticalTransitions.LBFGS()
gm = geometric_min_action_method(sys, init, maxiter = 100, method=method)
method = "HeymannVandenEijnden"
gm = geometric_min_action_method(sys, init, maxiter = 100, method=method)

gives

InexactError: Int64(4.384056944159697)
@oameye oameye changed the title gMAM InexactError gMAM HeymannVandenEijnden method broken Feb 23, 2024
@reykboerner
Copy link
Collaborator

Yes, sorry that the HeymannVandenEijnden option of gMAM is currently not working. I hope to fix it at some point but for now Optim.jl methods such as the default LBFGS() should be used instead.

@reykboerner reykboerner added bug Something isn't working enhancement New feature or request help wanted Extra attention is needed labels Mar 13, 2024
@oameye oameye added this to the Release version 1.0 milestone Mar 27, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working enhancement New feature or request help wanted Extra attention is needed
Projects
None yet
Development

No branches or pull requests

2 participants