[petsc-users] "-mat_superlu_colperm MMD_AT_PLUS_A" causes the program to hang

[Longmin RAN]

Dear all,


I'm using superlu within petsc to solve systems with symmetric sparse
matrix. In superlu manual I read that MMD_AT_PLUS_A column
permutation, together with little diagonal pivot threshold, should be
used for symmetric mode. So I launch my program with the following
options:
  -ksp_type preonly
  -pc_type lu
  -pc_factor_mat_solver_package superlu
  -mat_superlu_diagpivotthresh 0.001
  -mat_superlu_symmetricmode TRUE
  -mat_superlu_colperm MMD_AT_PLUS_A

It seems that "-mat_superlu_colperm MMD_AT_PLUS_A" causes the program
to hang: when I deleted this option, my calculation is executed
correctly. But It's always interesting to be able to use the column
permutation option. Do you guys have any ideas ?


Cheers,

Longmin