[petsc-users] SuperLU_dist computes rubbish

[Jan Blechta]

Using attached petsc4py code, matrix and right-hand side, SuperLU_dist
returns totally wrong solution for mixed Laplacian:

$ tar -xzf report.tar.gz
$ python test-solve.py -pc_factor_mat_solver_package mumps -ksp_final_residual
KSP final norm of residual 3.81865e-15
$ python test-solve.py -pc_factor_mat_solver_package umfpack -ksp_final_residual
KSP final norm of residual 3.68546e-14
$ python test-solve.py -pc_factor_mat_solver_package superlu_dist -ksp_final_residual
KSP final norm of residual 1827.72

Moreover final residual is random when run using mpirun -np 3. Maybe
a memory corruption issue? This is reproducible using PETSc 3.6.2 (and
SuperLU_dist configured by PETSc) and much older, see
http://fenicsproject.org/pipermail/fenics-support/2014-March/000439.html
but has never been reported upstream.

The code for assembling the matrix and rhs using FEniCS is also
included for the sake of completeness.

Jan
-------------- next part --------------
A non-text attachment was scrubbed...
Name: report.tar.gz
Type: application/gzip
Size: 54716 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20151115/11fdbfdc/attachment-0001.bin>