<div dir="ltr"><div dir="ltr"><div dir="ltr">Good afternoon petsc team, <div><br></div><div>I have a few questions regarding the matrixfree in the snes context.</div><div><br></div><div>I am currently using the matrix-free method with my own preconditioner to solve a TS problem. I currently trigger this with the -snes_mf_operator option.</div><div><br></div><div>So far the default values have been working great, but I ran into trouble during a convergence study when I made the time step very small, which I suspect is probably happening because it's dropping below the matrix-vector product tolerance.</div><div><br></div><div>I wanted to explore the parameters of the matrix-vector product but using the -snes_view option I notice that the following full message comes out:</div><div><br></div><div>-----------------------------------------------------------------------------------</div><div><div>SNES Object: 96 MPI processes</div><div>  type: newtonls</div><div>  maximum iterations=20, maximum function evaluations=10000</div><div>  tolerances: relative=1e-06, absolute=1e-06, solution=1e-15</div><div>  total number of linear solver iterations=14</div><div>  total number of function evaluations=35</div><div>  norm schedule ALWAYS</div><div>  Jacobian is never rebuilt</div><div>  Jacobian is applied matrix-free with differencing</div><div>  Jacobian is applied matrix-free with differencing, no explict Jacobian</div><div>  SNESLineSearch Object: 96 MPI processes</div><div>    type: l2</div><div>    maxstep=1.000000e+08, minlambda=1.000000e-12</div><div>    tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08</div><div>    maximum iterations=1</div><div>  KSP Object: 96 MPI processes</div><div>    type: gmres</div><div>      restart=100, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement</div><div>      happy breakdown tolerance 1e-30</div><div>    maximum iterations=400, initial guess is zero</div><div>    tolerances:  relative=0.01, absolute=1e-50, divergence=10000.</div><div>    left preconditioning</div><div>    using PRECONDITIONED norm type for convergence test</div><div>  PC Object: 96 MPI processes</div><div>    type: shell</div><div>      MyPreconditioner</div><div>    linear system matrix followed by preconditioner matrix:</div><div>    Mat Object: 96 MPI processes</div><div>      type: mffd</div><div>      rows=8400000, cols=8400000</div><div>        Matrix-free approximation:</div><div>          err=1.49012e-08 (relative error in function evaluation)</div><div>          Using wp compute h routine</div><div>              Does not compute normU</div><div>    Mat Object: 96 MPI processes</div><div>      type: mpiaij</div><div>      rows=8400000, cols=8400000, bs=4</div><div>      total: nonzeros=2041029504, allocated nonzeros=2041029504</div><div>      total number of mallocs used during MatSetValues calls=0</div></div><div>--------------------------------------------------------------------------------</div><div><br></div><div>I have been focusing specifically on the part:</div><div><br></div><div><div>        Matrix-free approximation:</div><div>          err=1.49012e-08 (relative error in function evaluation)</div><div>          Using wp compute h routine</div><div>              Does not compute normU</div></div><div><br></div><div>I notice that this part of the message does not change even when I specify "-snes_mf_type default" (which is the alternative to wp), so I am a bit confused on how I interact to the -snes_mf_type or if I am interpreting the message wrong. Could you clarify this?</div><div><br></div><div>Also let me know if you have any general suggestions on tightening the tolerance of the matricxvector products.</div><div><br></div><div>Thank you,</div><div><br></div><div>-Alfredo  </div><div><br></div><div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><font face="arial, sans-serif">Alfredo Duarte</font><div><font face="arial, sans-serif">Graduate Research Assistant</font></div><div><font face="arial, sans-serif">The University of Texas at Austin</font></div></div></div></div></div></div></div></div></div>