Fixed by Benjamin, closing bug.
** Changed in: dolfin
Status: Confirmed => Fix Committed
--
ODE Courtemanche demo does not converge
https://bugs.launchpad.net/bugs/427473
You received this bug notification because you are subscribed to DOLFIN.
Status in DOLFIN: Fix Committed
Bug description:
bkeh...@bkehlet-laptop:~/software/dolfin/demo/ode/courtemanche/cpp$ ./demo
Creating ODE of size 21.
Solving ODE over the time interval [0.0, 300]
Initializing continuous Galerkin method cG(1).
Using mono-adaptive Newton solver.
Using discrete tolerance tol = 1e-08.
Using uBLAS Krylov solver with no preconditioning.
Using BiCGStab Krylov solver for matrix Jacobian.
Time-stepping over the time interval [0, 300]
Warning: consider implementing Jacobian ODE::J() to improve efficiency.
*** Warning: Time slab system seems to be diverging.
*** Warning: ODE solver did not converge at t = 0.25
ODE solution computed in 0.010 seconds.
Total number of (macro) time steps: 0
Number of rejected time steps: 0
Function evaluations: 20
Potential at end time: 0.000000
_______________________________________________
DOLFIN-dev mailing list
[email protected]
http://www.fenics.org/mailman/listinfo/dolfin-dev
