multiple rhs

[Hong Zhang]

David,

I do not see any problem with the calling sequence.

The memory is determined in MatLUFactorSymbolic().
Does your code crashes within MatLUFactorSymbolic()?
Please send us complete error message.

Hong

On Wed, 11 Mar 2009, David Fuentes wrote:

>
> 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
>