dolfin team mailing list archive

[dolfin@xxxxxxxxxx]

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