[petsc-users] superlu ILUT preconditioner break

[Gong Ding]

Hi,
I compiled latest petsc-dev with superlu and try the ILUT preconditioner.
However, it raising segmentation violation error.

The valgrind reported:

Poisson Solver init...
Using GMRES linear solver...
Using SuperLU ILUT preconditioner...
 its 0   |residual|_2 = 6.216736e+00  linear iter = 0
==18820== Invalid read of size 4
==18820==    at 0x269AF7F: dgsisx (dgsisx.c:692)
==18820==    by 0x163E5D6: MatSolve_SuperLU_Private (superlu.c:321)
==18820==    by 0x163F162: MatSolve_SuperLU (superlu.c:366)
==18820==    by 0x151866A: MatSolve (matrix.c:3055)
==18820==    by 0x1B56197: PCApply_ILU (ilu.c:270)
==18820==    by 0x1F95DBF: PCApply (precon.c:384)
==18820==    by 0x1FBFE83: KSPInitialResidual (itres.c:64)
==18820==    by 0x1C6CF31: KSPSolve_GMRES (gmres.c:243)
==18820==    by 0x1BEC8B5: KSPSolve (itfunc.c:426)
==18820==    by 0x1D2433F: SNES_KSPSolve (snes.c:3060)
==18820==    by 0x1D3DE74: SNESSolve_LS (ls.c:190)
==18820==    by 0x1D1AA5B: SNESSolve (snes.c:2360)
==18820==  Address 0x0 is not stack'd, malloc'd or (recently) free'd
==18820==
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger
[0]PETSC ERROR: or see http://www.mcs.anl.gov/petsc/petsc-as/documentation/faq.html#valgrind[0]PETSC ERROR: or try http://valgrind.org on GNU/linux and Apple Mac OS X to find memory corruption errors
[0]PETSC ERROR: likely location of problem given in stack below
[0]PETSC ERROR: ---------------------  Stack Frames ------------------------------------
[0]PETSC ERROR: Note: The EXACT line numbers in the stack are not available,
[0]PETSC ERROR:       INSTEAD the line number of the start of the function
[0]PETSC ERROR:       is given.
[0]PETSC ERROR: [0] MatSolve_SuperLU_Private line 275 src/mat/impls/aij/seq/superlu/superlu.c
[0]PETSC ERROR: [0] MatSolve_SuperLU line 364 src/mat/impls/aij/seq/superlu/superlu.c
[0]PETSC ERROR: [0] MatSolve line 3038 src/mat/interface/matrix.c
[0]PETSC ERROR: [0] PCApply_ILU line 269 src/ksp/pc/impls/factor/ilu/ilu.c
[0]PETSC ERROR: [0] PCApply line 374 src/ksp/pc/interface/precon.c
[0]PETSC ERROR: [0] KSPInitialResidual line 45 src/ksp/ksp/interface/itres.c
[0]PETSC ERROR: [0] KSPSolve_GMRES line 224 src/ksp/ksp/impls/gmres/gmres.c
[0]PETSC ERROR: [0] KSPSolve line 332 src/ksp/ksp/interface/itfunc.c
[0]PETSC ERROR: [0] SNES_KSPSolve line 3058 src/snes/interface/snes.c
[0]PETSC ERROR: [0] SNESSolve_LS line 142 src/snes/impls/ls/ls.c
[0]PETSC ERROR: [0] SNESSolve line 2334 src/snes/interface/snes.c
--------------------- Error Message ------------------------------------
Fatal Error: at line 0 in  unknown file
------------------------------------------------------------------------
    
The PC configureation is listed as follows:
ierr = PCSetType (pc, (char*) PCILU);     genius_assert(!ierr);
ierr = PCFactorSetMatSolverPackage (pc, "superlu"); genius_assert(!ierr);
ierr = PetscOptionsSetValue("-mat_superlu_ilu_droptol","1e-7"); genius_assert(!ierr);
ierr = PetscOptionsSetValue("-mat_superlu_ilu_filltol","1e-2"); genius_assert(!ierr);
ierr = PetscOptionsSetValue("-mat_superlu_ilu_fillfactor","200"); genius_assert(!ierr);


Or I tested in a much simple code (attached) with command line argument
/petsc_solver -fmat poisson.mat -fvec poisson.rhs -pc_type ilu -pc_factor_mat_solver_package superlu  -ksp_monitor
which has the same problem.

Is it the problem of superlu itself or petsc interface?

-------------- next part --------------
A non-text attachment was scrubbed...
Name: poisson.mat
Type: application/octet-stream
Size: 58656 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110327/e762c3ac/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: poisson.rhs
Type: application/octet-stream
Size: 5640 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110327/e762c3ac/attachment-0003.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: petsc_solver.c
Type: text/x-csrc
Size: 2298 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110327/e762c3ac/attachment-0001.c>