[petsc-users] MatMatSolve in sequential call

Alexander Grayver agrayver at gfz-potsdam.de
Mon Aug 13 10:29:06 CDT 2012


On 13.08.2012 17:23, Hong Zhang wrote:
> Alexander :
> Hmm, isrow is a matrix ordering, likely a bug :-(
> Can you send us a stand-alone code that reproduces this error?

Hong,

I can try to pull it out, but this may be tricky.

> Do you use an external package for MatMatSolve()?

Yes, MUMPS.
I use ksp with preonly.

Originally I thought it will be faster to solve multiple problems using 
MatMatSolve.
This turned out to be not true however, but the code is still like that.

-- 
Regards,
Alexander



More information about the petsc-users mailing list