<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, May 15, 2017 at 5:03 PM, Mark Adams <span dir="ltr"><<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>></span> wrote:<br><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">I could use this fix for this global op counter that is getting trigger because I have an operator inside of another operator.<div></div></div></blockquote><div><br></div><div>Set the last arugment to PETSC_MAX_INT</div><div><br></div><div> <a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/SNES/SNESSetTolerances.html#SNESSetTolerances">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/SNES/SNESSetTolerances.html#SNESSetTolerances</a></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"><div dir="ltr"><div>Thanks,</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, May 4, 2017 at 9:09 AM, Mark Adams <span dir="ltr"><<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>></span> wrote:<br><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">OK, that makes sense, it fails when my velocity grid gets not tiny.<div><br></div><div>I can use tine velocity grids for now.</div></div><div class="gmail-m_4348075138553273734HOEnZb"><div class="gmail-m_4348075138553273734h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, May 2, 2017 at 11:18 AM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><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 class="gmail_extra"><div class="gmail_quote"><span>On Tue, May 2, 2017 at 10:10 AM, Mark Adams <span dir="ltr"><<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>></span> wrote:<br><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>/Users/markadams/Codes/petsc/a<wbr>rch-macosx-gnu-O/bin/mpiexec -n 1 ./vml -v_coord_cylinder -x_dm_refine 2 -v_dm_refine 2 -snes_rtol 1.e-6 -snes_stol 1.e-6 -ts_type cn -snes_fd -pc_type lu -ksp_type preonly -x_petscspace_order 1 -x_petscspace_poly_tensor -v_petscspace_order 1 -v_petscspace_poly_tensor -ts_dt .1 -ts_max_steps 10 -ts_final_time 1e10 -verbose 3 -num_species 1 -snes_monitor -masses 1,2,4 -thermal_temps 30,30,30 -domainv_lo -2,-2 -domainv_hi 2,2 -domainx_lo -12,-12 -domainx_hi 12,12 -E 0,0 -blobx_radius 2 -x_dm_view hdf5:x.h5 -x_vec_view hdf5:x.h5::append -v_dm_view hdf5:v.h5 -v_vec_view hdf5:v.h5::append -x_pre_dm_view hdf5:prex.h5 -x_pre_vec_view hdf5:prex.h5::append -snes_converged_reason -snes_linesearch_monitor -ts_adapt_monitor</div><div>main call SetupXDiscretization</div><div>main call SetInitialConditionDomain</div><div> VMLViewX DMGetOutputSequenceNumber=-1, cmd_str=-x_pre_vec_view</div><div> 0) species 0: charge density= -2.3940791757186e+00, z-momentum= 5.9851979392559e-01, energy= 3.2314073646197e-01, thermal-flux= 2.4419137539877e-01</div><div> 0) Normalized: charge density= -2.3940791757186e+00, z momentum= 5.9851979392559e-01, energy= 3.2314073646197e-01, thermal flux= 2.4419137539877e-01, local: 64 X cells, 81 X vertices</div><div> VMLViewX DMGetOutputSequenceNumber=0, cmd_str=(null)</div><div> VMLViewV DMGetOutputSequenceNumber=-1</div><div> 0 SNES Function norm 4.097052680599e+00 </div><div> 1 SNES Function norm 1.213148652908e-09 </div><div> Nonlinear solve did not converge due to DIVERGED_FUNCTION_COUNT iterations 1</div></div></blockquote><div><br></div></span><div>Neat! Mark, I think this has to do with you calling SNESEvaluateFunc() inside another one. We limit the number of function evaluations</div><div>to 10,000 by default, mostly to corral line searches. I think you hit this, and thus need to up the count.</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div><div class="gmail-m_4348075138553273734m_3634726183855441197h5"><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> TSAdapt none step 0 stage rejected t=0 + 1.000e-01, nonlinear solve failures 1 greater than current TS allowed</div><div>[0]PETSC ERROR: --------------------- Error Message ------------------------------<wbr>------------------------------<wbr>--</div><div>[0]PETSC ERROR: </div><div>[0]PETSC ERROR: TSStep has failed due to DIVERGED_NONLINEAR_SOLVE, increase -ts_max_snes_failures or make negative to attempt recovery</div><div>[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">http://www.mcs.anl.gov/petsc/d<wbr>ocumentation/faq.html</a> for trouble shooting.</div><div>[0]PETSC ERROR: Petsc Development GIT revision: v3.7.6-3659-g699918129c GIT Date: 2017-04-26 08:18:35 -0400</div><div>[0]PETSC ERROR: ./vml on a arch-macosx-gnu-O named MarksMac-5.local by markadams Tue May 2 11:04:02 2017</div><div>[0]PETSC ERROR: Configure options --with-cc=clang --with-cc++=clang++ COPTFLAGS="-O3 -g -mavx2" CXXOPTFLAGS="-O3 -g -mavx2" FOPTFLAGS="-O3 -g -mavx2" --download-mpich=1 --download-parmetis=1 --download-metis=1 --download-hypre=1 --download-ml=1 --download-triangle=1 --download-ctetgen=1 --download-p4est=1 --with-x=0 --download-superlu_dist --download-superlu --download-ctetgen --with-debugging=0 --download-hdf5=1 PETSC_ARCH=arch-macosx-gnu-O --download-chaco --with-viewfromoptions=1</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, May 1, 2017 at 10:25 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
and<br>
<br>
-snes_linesearch_monitor<br>
-ts_adapt_monitor<br>
<div class="gmail-m_4348075138553273734m_3634726183855441197m_-443444922006130246m_-4248779029230745088HOEnZb"><div class="gmail-m_4348075138553273734m_3634726183855441197m_-443444922006130246m_-4248779029230745088h5"><br>
<br>
> On May 1, 2017, at 7:51 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> wrote:<br>
><br>
> Run with -snes_converged_reason.<br>
><br>
> Matt<br>
><br>
> On Mon, May 1, 2017 at 7:14 PM, Mark Adams <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>> wrote:<br>
> I get this SNES failure and I don't understand what the problem is. The rtol is 1.e-6 and the first iteration reduces the residual by 9 orders of magnitude. Yet, TS is not satisfied. What is going on here?<br>
><br>
> mpiexec -n 1 ./vml -v_coord_cylinder -x_dm_refine 2 -v_dm_refine 2 -snes_rtol 1.e-6 -snes_stol 1.e-6 -ts_type cn -snes_fd -pc_type lu -ksp_type preonly -x_petscspace_order 1 -x_petscspace_poly_tensor -v_petscspace_order 1 -v_petscspace_poly_tensor -ts_dt .1 -ts_max_steps 10 -ts_final_time 1e10 -verbose 3 -num_species 1 -snes_monitor -masses 1,2,4 -thermal_temps 30,30,30 -domainv_lo -2,-2 -domainv_hi 2,2 -domainx_lo -12,-12 -domainx_hi 12,12 -E 0,0 -blobx_radius 2 -x_dm_view hdf5:x.h5 -x_vec_view hdf5:x.h5::append -v_dm_view hdf5:v.h5 -v_vec_view hdf5:v.h5::append -x_pre_dm_view hdf5:prex.h5 -x_pre_vec_view hdf5:prex.h5::append<br>
> ....<br>
><br>
> 0 SNES Function norm 4.097052680599e+00<br>
> 1 SNES Function norm 1.213148652908e-09<br>
> [0]PETSC ERROR: --------------------- Error Message ------------------------------<wbr>------------------------------<wbr>--<br>
> [0]PETSC ERROR:<br>
> [0]PETSC ERROR: TSStep has failed due to DIVERGED_NONLINEAR_SOLVE, increase -ts_max_snes_failures or make negative to attempt recovery<br>
> [0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" rel="noreferrer" target="_blank">http://www.mcs.anl.gov/petsc/d<wbr>ocumentation/faq.html</a> for trouble shooting.<br>
> [0]PETSC ERROR: Petsc Development GIT revision: v3.7.6-3659-g699918129c GIT Date: 2017-04-26 08:18:35 -0400<br>
> [0]PETSC ERROR: ./vml on a arch-macosx-gnu-O named MarksMac-5.local by markadams Mon May 1 19:21:32 2017<br>
> [0]PETSC ERROR: Configure options --with-cc=clang --with-cc++=clang++ COPTFLAGS="-O3 -g -mavx2" CXXOPTFLAGS="-O3 -g -mavx2" FOPTFLAGS="-O3 -g -mavx2" --download-mpich=1 --download-parmetis=1 --download-metis=1 --download-hypre=1 --download-ml=1 --download-triangle=1 --download-ctetgen=1 --download-p4est=1 --with-x=0 --download-superlu_dist --download-superlu --download-ctetgen --with-debugging=0 --download-hdf5=1 PETSC_ARCH=arch-macosx-gnu-O --download-chaco --with-viewfromoptions=1<br>
><br>
><br>
><br>
><br>
> --<br>
> 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<br>
<br>
</div></div></blockquote></div><br></div>
</blockquote></div></div></div><div><div class="gmail-m_4348075138553273734m_3634726183855441197h5"><br><br clear="all"><span class="gmail-HOEnZb"><font color="#888888"><div><br></div>-- <br><div class="gmail-m_4348075138553273734m_3634726183855441197m_-443444922006130246gmail_signature">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>
</font></span></div></div></div></div><span class="gmail-HOEnZb"><font color="#888888">
</font></span></blockquote></div><span class="gmail-HOEnZb"><font color="#888888"><br></font></span></div>
</div></div></blockquote></div><br></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">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></div>