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