Hi,

Several months ago, I posted message "superlu ILUT preconditioner break" with 
following link
http://lists.mcs.anl.gov/pipermail/petsc-users/2011-March/008355.html

Today I studied this problem carefully to see why it crashes.

The segment fault is arising at line 692 of superlu's dgsisx.c, because the 
array "perm" is 0.
This array which contains the permutation calculated from MC64 procedure, can 
only be initialized when nofact flag is true. However, petsc will call dgsisx 
twice, with nofact value 1 and 0, respectively.  In the second calling, the 
perm array is not initialized but be used. That's why the code crash.

Hope this bug can be fixed without too many efforts.

Gong Ding





Reply via email to