[petsc-dev] Problem with conjugate gradient solver

[Anush Krishnan]

Hi,

I was trying to compile and run the example file
$PETSC_DIR/src/ksp/pc/examples/tutorials/ex3.c, which demonstrates PETSc's
preconditioned conjugate gradient solver for a system of linear equations.
After compiling, running make runex3 is supposed to throw up the error
"Divergence because of indefinite preconditioner" and running make runex3_pd
is supposed to work and give the correct solution. But instead I get the
error "Other kind of divergence: this should not happen." for the latter. Is
there some bug in the code, or are there any other options that need to be
set?

Thanks,
Anush
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20100823/5e4dbb68/attachment.html>