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":
_______________________________________________
DOLFIN-dev mailing list
[email protected]
http://www.fenics.org/mailman/listinfo/dolfin-dev
