GridapODEs.jl
GridapODEs.jl copied to clipboard
Time stepping for Gridap
Hi @fverdugo I have been taking a look at `DifferentialEquations.jl`, in order to create `ODESolver`s as wrappers of this library solvers. In our current version of the `GridapTimeStepper` package, I...
The issue you had was that you didn't utilize a sparse Jacobian solver.
```julia sol_lazy =solve(solver,op) # We need a forward solution that can be iterated in reverse sol = collect(sol_lazy) # Option 1 sol = CeckpointedSol(sol_lazy) # Option 2 sol_reversed = reverse(sol)...