>> which starts to >> look more like: Preconditioner.solve(Matrix, q, y). > > Or: p(Matrix); p.matvec(q, y). p(Matrix); p.matvec(q, y).
Or both... Preconditioner p(Matrix); p.matvec(q, y); ... (keep using Matrix and whatever you have done to it...) p.matvec(q2,y2); p.matvec(Matrix,q,y); ...(reset the matrix, do whatever you have to do...) Can't we all just get along?? -Ben ------------------------------------------------------------------------------ This SF.net email is sponsored by: SourcForge Community SourceForge wants to tell your story. http://p.sf.net/sfu/sf-spreadtheword _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
