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>
