dolfin team mailing list archive

["Garth N. Wells" <gnw20@xxxxxxxxx>]

Robert Kirby wrote:
> Hi all,
> I've gotten a prototype working where I construct
> a uBLASSparseMatrix and get it to compute the same result as my
> matrix-free uBLASKrylovMatrix.  Both of these are constructed
> in PyDOLFIN.  However, I can feed the uBLASSparseMatrix
> I construct into a dolfin.uBLASKrylovSolver solve method,
> but I get a type error when I try to put my
> uBLASKrylovMatrix with mult and dim implemented in Python
> into the solver.
>
> In fact, I get a 
>
> TypeError: in method 'uBLASKrylovSolver_solve', argument 2 of type 
> 'dolfin::uBLASKrylovMatrix const &'

I get an error in 'assemble_matrix' with the attached program, so I'm 
not getting to the Krylov solver.

Try setting the preconditioner to 'none' when creating the Krylov solver.

Garth

> Not very revealing.  I've attached a horribly ugly source code for 
> anyone interested in doing a post-mortem on this.
>
> Thanks,
> Rob
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev