I need to benchmark the Krylov solver against UMFPACK for analysis
purposes, in order to do this I just changed the declaration:
uBlasLUsolver solver
to:
uBlasKrylovSolver solver;
and then I just use solve() likewise I do with the LU solver. I tried with
a very simple diagonal matrix and the Krylov solver does not converge. I
tried both with BiCGStab and GMRES. It's not a problem due to the
tolerance because the result is clearly wrong (diverging). I'm also not
able to use:
set("Krylov monitor convergence", true);
in the main file because no informations are displayed more than the
non-converging warning. Should I change something else?
Alessio
_______________________________________________
DOLFIN-dev mailing list
[email protected]
http://www.fenics.org/mailman/listinfo/dolfin-dev
