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
report.tar.gz
Description: application/gzip
