Ambient Forcing v1.0.0
Complete Refactor of the Ambient Forcing package. This version is significantly faster than the initial release.
This version replaces ForwardDiff.jl with Symbolics.jl as the standard solution for obtaining Jacobians of the constraint manifold.
Furthermore, the projection operator is now based directly on the Jacobian, instead of the more costly computation of an ortho-normal basis for the nullspace, see further: https://math.stackexchange.com/questions/1318637/projection-of-a-vector-onto-the-null-space-of-a-matrix.
The previous function ambient_forcing is split into two: ambient_forcing_problem creates an ODEProblem for a given system with constraints and can be solved with different initial conditions and random vectors, saving the time for precompilation. The new ambient_forcing is just a convenience wrapper to the solve call.