dolfin team mailing list archive

[Kristen Kaasbjerg <cosby@xxxxxxxxx>]

Hi

I noticed that the LinearPDE class in python always uses
the direct LU solver (the solver type is set by a call to
get("PDE linear solver") which always returns "direct").
I tried to change it to "iterative" and didn't experience
any problems. Should I be careful using the iterative
solvers since the direct are always used by default ?

Kristen

PS: there is a typo in line 205 of assemble.py:

elif solve_type == "iterative":

should of course be

elif solver_type == "iterative":