multiple rhs

[David Fuentes]

Hello,

I have a sparse matrix, A, with which I want to solve multiple right hand sides
with a direct solver.  Is this the correct call sequence ?


   MatGetFactor(A,MAT_SOLVER_PETSC,MAT_FACTOR_LU,&Afact);
   IS isrow,iscol;
   MatGetOrdering(A,MATORDERING_ND,&isrow,&iscol);
   MatLUFactorSymbolic(Afact,A,isrow,iscol,&info);
   MatLUFactorNumeric(Afact,A,&info);
   MatMatSolve(Afact,B,X);


my solve keeps running out of memory

"[0]PETSC ERROR: Memory requested xxx!"


is this in bytes? I can't tell if the problem I'm trying to solve
is too large form my machine or if I just have bug in the call sequence.




thank you,
David Fuentes