[petsc-users] The ColPerm option for superlu

[Jed Brown]

On Thu, Jan 19, 2012 at 19:20, Xuefei (Rebecca) Yuan <xyuan at lbl.gov> wrote:

> Hello,
>
> I use ksp_view to get the information about the lu pc as follow:
>
> PC Object: 1 MPI processes
>  type: lu
>    LU: out-of-place factorization
>    tolerance for zero pivot 2.22045e-14
>    matrix ordering: nd
>    factor fill ratio given 0, needed 0
>      Factored matrix follows:
>        Matrix Object:         1 MPI processes
>          type: seqaij
>          rows=262144, cols=262144
>          package used to perform factorization: superlu
>          total: nonzeros=0, allocated nonzeros=0
>          total number of mallocs used during MatSetValues calls =0
>            SuperLU run parameters:
>              Equil: NO
>              ColPerm: 3
>              IterRefine: 0
>              SymmetricMode: NO
>              DiagPivotThresh: 1
>              PivotGrowth: NO
>              ConditionNumber: YES
>              RowPerm: 0
>              ReplaceTinyPivot: NO
>              PrintStat: YES
>              lwork: 0
>  linear system matrix = precond matrix:
>  Matrix Object:   1 MPI processes
>    type: seqaij
>    rows=262144, cols=262144
>    total: nonzeros=25969216, allocated nonzeros=25969216
>    total number of mallocs used during MatSetValues calls =0
>      using I-node routines: found 65536 nodes, limit used is 5
>
>
> In the SuperLU run parameters, how could I change ColPerm:3 to ColPerm has
> the value of 2?
>

$ ./ex2 -pc_type lu -pc_factor_mat_solver_package superlu -help |grep
colperm
  -mat_superlu_colperm <COLAMD> (choose one of) NATURAL MMD_ATA
MMD_AT_PLUS_A COLAMD (None)


>
> Is there any option that I can use to make such a change through PETSc?
>
> Thanks very much!
>
> Best regards,
>
> Rebecca
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/4e878a11/attachment.htm>