dolfin team mailing list archive
-
dolfin team
-
Mailing list archive
-
Message #15313
[Bug 427473] Re: ODE Courtemanche demo does not converge
** Changed in: dolfin
Status: New => Confirmed
** Changed in: dolfin
Importance: Undecided => Medium
--
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: Confirmed
Bug description:
bkehlet@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