[dune-fem] problem with diagonal precond and istl

Andreas Dedner a.s.dedner at warwick.ac.uk
Sun Mar 30 14:10:53 CEST 2014


Hi.
There is a problem when setting the parameter
stl.preconditioning.method: jacobi
(for example in fem-poisson). An exception is thrown

Error: NotImplemented
[applyToISTLBlockVector:/home/dedner/TRUNK/dune-fem/dune/fem/solver/diagonalpreconditioner.hh:46]:
preconditioning not possible for non-assembled operators

Q1: the problem seems to be that the ISTLMatrixObject is used to check
if it is an assembled
       operator but that class is not derived from AssembledOperator -
that class also
       formally wants DF and the ISTLMatrixObject only knows DFS.
Q2: why is the dune-fem implementation used instead of the istl internal
one.
       Interestingly enough I even have
       fem.preconditioning: false

Best
Andreas




More information about the dune-fem mailing list