dolfin team mailing list archive
-
dolfin team
-
Mailing list archive
-
Message #00280
Update computation of mono-adaptive Jacobian to handle implicit ODEs. [...]
Commit from logg (2005-02-02 16:47 CET)
----------------
Update computation of mono-adaptive Jacobian to handle implicit ODEs.
At this point, the following solvers are implemented:
1. multi-adaptive + fixed point
2. multi-adaptive + Newton (preconditioner missing)
3. mono-adaptive + fixed point
4. mono-adaptive + Newton
All solvers are implemented for general order cG(q) or dG(q).
Solvers 1 and 2 can solve u' = f.
Solvers 3 and 4 can solve u' = f and Mu' = f.
Remains:
- Testing
- Optimizations (we might benefit from using PETSc vector operations)
- Benchmarks
- Implement preconditioner for case 2
dolfin src/demo/solvers/ode/test/main.cpp 1.47
dolfin src/kernel/la/NewGMRES.cpp 1.8
dolfin src/kernel/ode/MonoAdaptiveJacobian.cpp 1.5
dolfin src/kernel/ode/MonoAdaptiveNewtonSolver.cpp 1.9
dolfin src/kernel/ode/MultiAdaptiveNewtonSolver.cpp 1.5
dolfin src/kernel/ode/dolfin/MonoAdaptiveJacobian.h 1.4