dolfin team mailing list archive

[Robert Kirby <robert.c.kirby@xxxxxxxxx>]

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 &'

Not very revealing.  I've attached a horribly ugly source code for anyone
interested in doing a post-mortem on this.

Thanks,
Rob

Attachment: helmholtz.py
Description: Binary data