<div dir="ltr"><div dir="ltr"><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Feb 4, 2019 at 6:05 AM Aman Saxena via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr">I am trying to use PETSc's TS to implement implicit time-stepping, for 1d-Euler Equation discretized using 5th order WENO schemes. I have following queries related to using some of the snes options::</div><div dir="ltr"><br></div><div dir="ltr"><div dir="ltr">-- What does -ts_max_snes_failure means and what does the non-linear solver do once it diverges to ensure convergence?</div></div></div></div></blockquote><div><br></div><div>You want to look at the manual but I think this is the number of nonlinear solver failures that it will try to recover from before quitting.</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"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><br></div><div dir="ltr">-- If I do not specify “-snes_fd” or “-snes_mf”, does by default it uses the jacobian I have coded instead of building it using Finite-differences?</div></div></div></div></blockquote><div><br></div><div>test it but if you give it a Jacobian I would hope that it uses it! </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><br></div><div dir="ltr">-- As “-snes_fd ” replaces my routine to calculate the jacobian with petsc’s fd routine, does it stores the jacobian in the same Matrix that I pass to my hand-coded routine. It might lead to some inconsistency as my actual scheme involves stencils width equals 3 while I am approximating the Jacobian using first order scheme which requires stencils width equals 1. So I have used another “dmda” with sw = 1, to initialize the matrix using “DMCreateMatrix()” while fd needs sw = 3.</div></div></div></div></blockquote><div><br></div><div>This should be OK. Try it.</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"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><br></div><div dir="ltr">-- I am not able to use “-snes_type test” it gives me the following error </div><div dir="ltr">“[0]PETSC ERROR: Unknown type. Check for miss-spelling or missing package: <a href="http://www.mcs.anl.gov/petsc/documentation/installation.html#external" target="_blank">http://www.mcs.anl.gov/petsc/documentation/installation.html#external</a></div><div dir="ltr"> [0]PETSC ERROR: Unable to find requested SNES type test”</div><div><br></div></div></div></div></blockquote><div><br></div><div>This has changed. Use -snes_test_jacobian</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"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div></div><div>-- What exactly "PCMAT" do? Do I have to supply matrix other than in RHSJacobian for PCMAT?</div></div></div></div></blockquote><div><br></div><div>I'm not sure where you picked this up you don't need to worry about it:</div><div><br></div><div><a href="https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCMAT.html">https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCMAT.html</a><br></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"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><br></div><div>Thanks </div></div></div></div>
</blockquote></div></div></div>