tomaklutfu
tomaklutfu
/ odedual.jl
Last active
June 24, 2020 19:52
Partials of states have different values for non-stiff ODE solvers with a callback at the event location
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| using OrdinaryDiffEq, ForwardDiff, LinearAlgebra | |
| y0val = 0.6 | |
| vy0val = -0.1 | |
| gyval = 0.005 | |
| # val d/dy0 d/dvy0 d/dgy | |
| y0 = ForwardDiff.Dual( y0val, 1.0, 0.0, 0.0 ) | |
| vy0= ForwardDiff.Dual(vy0val, 0.0, 1.0, 0.0 ) | |
| gy = ForwardDiff.Dual( gyval, 0.0, 0.0, 1.0 ) | |
| u0dual = [y0; vy0] |