Running squential and parallel linear solver in the same code?

[Ben Tay]

Hi,

In my MPI CFD code, I have to solve 2 linear eqns - momentum and poisson
eqn. There's some error in the code and I don't know which one causes the
error. May I know if it is possible to solve 1 eqn in sequential mode while
the other in MPI to pinpoint the error since my code works ok in a
sequential mode.

I tried to use this for one of the eqn:


call
MatCreateSeqAIJ(PETSC_COMM_SELF,size_x*size_y,size_x*size_y,9,PETSC_NULL_INTEGER,A_mat,ierr)

call VecCreateSeq(PETSC_COMM_SELF,size_x*size_y,b_rhs,ierr)

but there's error.

Thank you.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20070123/e321cad8/attachment.htm>