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