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

[Hong Zhang]

Longmin :

I cannot reproduce it with petsc example:
petsc-dev/src/ksp/ksp/examples/tutorials>./ex2 -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
Norm of error < 1.e-12 iterations 1

Please use a debugger to check where it hangs. Then use valgrind to check
posible memory corruption.

Hong
>
> 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
>