Taylor McDonnell
Taylor McDonnell
The following no longer works. ```julia function f(out,du,u,p,t) out[1] = - 0.04u[1] + 1e4*u[2]*u[3] - du[1] out[2] = + 0.04u[1] - 3e7*u[2]^2 - 1e4*u[2]*u[3] - du[2] out[3] = u[1] +...
How much might be involved to implement the discrete adjoint analytically (rather than through AD)? I did a derivation, and it seems like a general implementation should be possible with...
Different amounts of computational noise in the value and partials of a dual number has the potential to corrupt even NaN-safe computations. Minimum working example: ```julia using ForwardDiff, LinearAlgebra #...
The following mass matrix DAE solvers fail to solve when used with callbacks that modify the current state: `Rodas4`, `Rodas42`, `Rodas4P`, `Rodas4P2`, `Rodas5`. MWE: ```julia # original problem using DifferentialEquations...
This package is currently compatible with ForwardDiff, but not ReverseDiff. ```julia using NLsolve, ForwardDiff, ReverseDiff function residual!(r, y, x) r[1] = (y[1] + x[1])*(y[2]^3-x[2])+x[3] r[2] = sin(y[2]*exp(y[1])-1)*x[4] end function solve(x)...
This package doesn't seem to load correctly in Julia 1.8. I think I might be able to access a useful stack trace if I were to try loading this package...
This package currently only allows rigid connections. A nice feature would be expand this to allow other types of joints. My thoughts for making this happen is to introduce a...
DAEs can be converted to ODEs using index reduction. It would be nice to try that for this package so that more solvers/features from DifferentialEquations can be used. The sensitivity...
While it is implied that dual numbers shouldn't show up in the (constant) parameters, it could still happen in practice. It might be worth adding a note to the documentation...