[petsc-users] Using superlu_dist in a direct solve

[Sanjay Govindjee]

I wanted to use SuperLU Dist to perform a direct solve but seem to be 
encountering
a problem.  I was wonder if this is a know issue and if there is a 
solution for it.

The problem is easily observed using ex6.c in src/ksp/ksp/examples/tests.

Out of the box: make runex6 produces a residual error of O(1e-11), all 
is well.

I then changed the run to run on two processors and add the flag
-pc_factor_mat_solver_package spooles  this produces a residual error of 
O(1e-11), all is still well.

I then switch over to -pc_factor_mat_solver_package superlu_dist and the
residual error comes back as 22.6637!  Something seems very wrong.

My build is perfectly vanilla:

export PETSC_DIR=/Users/sg/petsc-3.3-p5/
export PETSC_ARCH=intel

./configure --with-cc=icc --with-fc=ifort  \ 
-download-{spooles,parmetis,superlu_dist,prometheus,mpich,ml,hypre,metis}

make PETSC_DIR=/Users/sg/petsc-3.3-p5/ PETSC_ARCH=intel all
make PETSC_DIR=/Users/sg/petsc-3.3-p5/ PETSC_ARCH=intel test

-sanjay
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20121223/bafbafc7/attachment.html>