dolfin team mailing list archive

[Andrew Benson <question135129@xxxxxxxxxxxxxxxxxxxxx>]

New question #135129 on DOLFIN:
https://answers.launchpad.net/dolfin/+question/135129

I'm trying to modify the "stiff ODE" demos to use a multi-adaptive solver (I'd like to use this for a more interesting calculation and figured that understanding the demos would be a good way to begin). If I modify the line:

  ode.parameters["method"] = "dg";

in main.cpp to:

  ode.parameters["method"] = "mdg";

(and also modify the ODE function to have the correct form for a multi-adaptive solver) then I get an error:

Solving test problem number 9.
A mixed stiff/nonstiff test problem.
Solving ODE over the time interval [0.0, 30]
  Using multi-adaptive Newton solver.
  Using discrete tolerance tol = 1e-05.
  Time-stepping over the time interval [0, 30]
terminate called after throwing an instance of 'std::runtime_error'
  what():  *** Error: No init(..) function for preconditioning uBLASKrylovMatrix

Any ideas as to why there's no suitable init function?

Thanks,
Andrew.

-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.