<div class="gmail_quote">On Thu, Jan 19, 2012 at 19:20, Xuefei (Rebecca) Yuan <span dir="ltr"><<a href="mailto:xyuan@lbl.gov">xyuan@lbl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hello,<br>
<br>
I use ksp_view to get the information about the lu pc as follow:<br>
<br>
PC Object: 1 MPI processes<br>
type: lu<br>
LU: out-of-place factorization<br>
tolerance for zero pivot 2.22045e-14<br>
matrix ordering: nd<br>
factor fill ratio given 0, needed 0<br>
Factored matrix follows:<br>
Matrix Object: 1 MPI processes<br>
type: seqaij<br>
rows=262144, cols=262144<br>
package used to perform factorization: superlu<br>
total: nonzeros=0, allocated nonzeros=0<br>
total number of mallocs used during MatSetValues calls =0<br>
SuperLU run parameters:<br>
Equil: NO<br>
ColPerm: 3<br>
IterRefine: 0<br>
SymmetricMode: NO<br>
DiagPivotThresh: 1<br>
PivotGrowth: NO<br>
ConditionNumber: YES<br>
RowPerm: 0<br>
ReplaceTinyPivot: NO<br>
PrintStat: YES<br>
lwork: 0<br>
linear system matrix = precond matrix:<br>
Matrix Object: 1 MPI processes<br>
type: seqaij<br>
rows=262144, cols=262144<br>
total: nonzeros=25969216, allocated nonzeros=25969216<br>
total number of mallocs used during MatSetValues calls =0<br>
using I-node routines: found 65536 nodes, limit used is 5<br>
<br>
<br>
In the SuperLU run parameters, how could I change ColPerm:3 to ColPerm has the value of 2?<br></blockquote><div><br></div><div><div>$ ./ex2 -pc_type lu -pc_factor_mat_solver_package superlu -help |grep colperm</div><div>
-mat_superlu_colperm <COLAMD> (choose one of) NATURAL MMD_ATA MMD_AT_PLUS_A COLAMD (None)</div></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
Is there any option that I can use to make such a change through PETSc?<br>
<br>
Thanks very much!<br>
<br>
Best regards,<br>
<font color="#888888"><br>
Rebecca<br>
<br>
<br>
</font></blockquote></div><br>