<div dir="ltr"><div dir="ltr">On Fri, Feb 3, 2023 at 3:03 PM Nicolas Barnafi <<a href="mailto:nbarnafi@cmm.uchile.cl">nbarnafi@cmm.uchile.cl</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">> There are a number of common errors:<br>
> <br>
> 1) Your PC has a prefix<br>
> <br>
> 2) You have not called KSPSetFromOptions() here<br>
> <br>
> Can you send the -ksp_view output?<br>
<br>
The PC at least has no prefix. I had to set ksp_rtol to 1 to get through <br>
the solution process, you will find both the petsc_rc and the ksp_view <br>
at the bottom of this message.<br>
<br>
Options are indeed being set from the options file, but there must be <br>
something missing at a certain level. Thanks for looking into this.<br></blockquote><div><br></div><div>Okay, the next step is to pass</div><div><br></div><div> -info</div><div><br></div><div>and send the output. This will tell us how the default splits were done. If that</div><div>is not conclusive, we will have to use the debugger.</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Best<br>
<br>
---- petsc_rc file<br>
<br>
-ksp_monitor<br>
-ksp_type gmres<br>
-ksp_view<br>
-mat_type aij<br>
-ksp_norm_type unpreconditioned<br>
-ksp_atol 1e-14<br>
-ksp_rtol 1<br>
-pc_type fieldsplit<br>
-pc_fieldsplit_type multiplicative<br>
<br>
---- ksp_view<br>
<br>
KSP Object: 1 MPI process<br>
type: gmres<br>
restart=500, using Classical (unmodified) Gram-Schmidt <br>
Orthogonalization with no iterative refinement<br>
happy breakdown tolerance 1e-30<br>
maximum iterations=10000, nonzero initial guess<br>
tolerances: relative=1., absolute=1e-14, divergence=10000.<br>
right preconditioning<br>
using UNPRECONDITIONED norm type for convergence test<br>
PC Object: 1 MPI process<br>
type: fieldsplit<br>
FieldSplit with MULTIPLICATIVE composition: total splits = 4<br>
Solver info for each split is in the following KSP objects:<br>
Split number 0 Defined by IS<br>
KSP Object: (fieldsplit_0_) 1 MPI process<br>
type: preonly<br>
maximum iterations=10000, initial guess is zero<br>
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.<br>
left preconditioning<br>
using DEFAULT norm type for convergence test<br>
PC Object: (fieldsplit_0_) 1 MPI process<br>
type: ilu<br>
PC has not been set up so information may be incomplete<br>
out-of-place factorization<br>
0 levels of fill<br>
tolerance for zero pivot 2.22045e-14<br>
matrix ordering: natural<br>
matrix solver type: petsc<br>
matrix not yet factored; no additional information available<br>
linear system matrix = precond matrix:<br>
Mat Object: (fieldsplit_0_) 1 MPI process<br>
type: seqaij<br>
rows=615, cols=615<br>
total: nonzeros=9213, allocated nonzeros=9213<br>
total number of mallocs used during MatSetValues calls=0<br>
not using I-node routines<br>
Split number 1 Defined by IS<br>
KSP Object: (fieldsplit_1_) 1 MPI process<br>
type: preonly<br>
maximum iterations=10000, initial guess is zero<br>
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.<br>
left preconditioning<br>
using DEFAULT norm type for convergence test<br>
PC Object: (fieldsplit_1_) 1 MPI process<br>
type: ilu<br>
PC has not been set up so information may be incomplete<br>
out-of-place factorization<br>
0 levels of fill<br>
tolerance for zero pivot 2.22045e-14<br>
matrix ordering: natural<br>
matrix solver type: petsc<br>
matrix not yet factored; no additional information available<br>
linear system matrix = precond matrix:<br>
Mat Object: (fieldsplit_1_) 1 MPI process<br>
type: seqaij<br>
rows=64, cols=64<br>
total: nonzeros=0, allocated nonzeros=0<br>
total number of mallocs used during MatSetValues calls=0<br>
using I-node routines: found 13 nodes, limit used is 5<br>
Split number 2 Defined by IS<br>
KSP Object: (fieldsplit_2_) 1 MPI process<br>
type: preonly<br>
maximum iterations=10000, initial guess is zero<br>
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.<br>
left preconditioning<br>
using DEFAULT norm type for convergence test<br>
PC Object: (fieldsplit_2_) 1 MPI process<br>
type: ilu<br>
PC has not been set up so information may be incomplete<br>
out-of-place factorization<br>
0 levels of fill<br>
tolerance for zero pivot 2.22045e-14<br>
matrix ordering: natural<br>
matrix solver type: petsc<br>
matrix not yet factored; no additional information available<br>
linear system matrix = precond matrix:<br>
Mat Object: (fieldsplit_2_) 1 MPI process<br>
type: seqaij<br>
rows=240, cols=240<br>
total: nonzeros=2140, allocated nonzeros=2140<br>
total number of mallocs used during MatSetValues calls=0<br>
not using I-node routines<br>
Split number 3 Defined by IS<br>
KSP Object: (fieldsplit_3_) 1 MPI process<br>
type: preonly<br>
maximum iterations=10000, initial guess is zero<br>
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.<br>
left preconditioning<br>
using DEFAULT norm type for convergence test<br>
PC Object: (fieldsplit_3_) 1 MPI process<br>
type: ilu<br>
PC has not been set up so information may be incomplete<br>
out-of-place factorization<br>
0 levels of fill<br>
tolerance for zero pivot 2.22045e-14<br>
matrix ordering: natural<br>
matrix solver type: petsc<br>
matrix not yet factored; no additional information available<br>
linear system matrix = precond matrix:<br>
Mat Object: (fieldsplit_3_) 1 MPI process<br>
type: seqaij<br>
rows=300, cols=300<br>
total: nonzeros=2292, allocated nonzeros=2292<br>
total number of mallocs used during MatSetValues calls=0<br>
not using I-node routines<br>
linear system matrix = precond matrix:<br>
Mat Object: 1 MPI process<br>
type: seqaij<br>
rows=1219, cols=1219<br>
total: nonzeros=26443, allocated nonzeros=26443<br>
total number of mallocs used during MatSetValues calls=0<br>
not using I-node routines<br>
solving time: 0.00449609<br>
iterations: 0<br>
estimated error: 25.4142<br>
<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>