[petsc-users] SNES error

Matthew Knepley knepley at gmail.com
Mon May 15 17:12:20 CDT 2017 

On Mon, May 15, 2017 at 5:03 PM, Mark Adams <mfadams at lbl.gov> wrote:

> I could use this fix for this global op counter that is getting trigger
> because I have an operator inside of another operator.
>

Set the last arugment to PETSC_MAX_INT


http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/SNES/SNESSetTolerances.html#SNESSetTolerances

    Matt


> Thanks,
>
> On Thu, May 4, 2017 at 9:09 AM, Mark Adams <mfadams at lbl.gov> wrote:
>
>> OK, that makes sense, it fails when my velocity grid gets not tiny.
>>
>> I can use tine velocity grids for now.
>>
>> On Tue, May 2, 2017 at 11:18 AM, Matthew Knepley <knepley at gmail.com>
>> wrote:
>>
>>> On Tue, May 2, 2017 at 10:10 AM, Mark Adams <mfadams at lbl.gov> wrote:
>>>
>>>> /Users/markadams/Codes/petsc/arch-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
>>>> main call SetupXDiscretization
>>>> main call SetInitialConditionDomain
>>>>                 VMLViewX DMGetOutputSequenceNumber=-1,
>>>> cmd_str=-x_pre_vec_view
>>>>   0) species 0: charge density= -2.3940791757186e+00, z-momentum=
>>>>  5.9851979392559e-01, energy=  3.2314073646197e-01, thermal-flux=
>>>>  2.4419137539877e-01
>>>>           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
>>>>                 VMLViewX DMGetOutputSequenceNumber=0, cmd_str=(null)
>>>>         VMLViewV DMGetOutputSequenceNumber=-1
>>>>     0 SNES Function norm 4.097052680599e+00
>>>>     1 SNES Function norm 1.213148652908e-09
>>>>   Nonlinear solve did not converge due to DIVERGED_FUNCTION_COUNT
>>>> iterations 1
>>>>
>>>
>>> Neat! Mark, I think this has to do with you calling SNESEvaluateFunc()
>>> inside another one. We limit the number of function evaluations
>>> to 10,000 by default, mostly to corral line searches. I think you hit
>>> this, and thus need to up the count.
>>>
>>>   Thanks,
>>>
>>>      Matt
>>>
>>>
>>>>       TSAdapt none step   0 stage rejected t=0          + 1.000e-01,
>>>> nonlinear solve failures 1 greater than current TS allowed
>>>> [0]PETSC ERROR: --------------------- Error Message
>>>> --------------------------------------------------------------
>>>> [0]PETSC ERROR:
>>>> [0]PETSC ERROR: TSStep has failed due to DIVERGED_NONLINEAR_SOLVE,
>>>> increase -ts_max_snes_failures or make negative to attempt recovery
>>>> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html
>>>> for trouble shooting.
>>>> [0]PETSC ERROR: Petsc Development GIT revision: v3.7.6-3659-g699918129c
>>>>  GIT Date: 2017-04-26 08:18:35 -0400
>>>> [0]PETSC ERROR: ./vml on a arch-macosx-gnu-O named MarksMac-5.local by
>>>> markadams Tue May  2 11:04:02 2017
>>>> [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
>>>>
>>>>
>>>> On Mon, May 1, 2017 at 10:25 PM, Barry Smith <bsmith at mcs.anl.gov>
>>>> wrote:
>>>>
>>>>>
>>>>>   and
>>>>>
>>>>>   -snes_linesearch_monitor
>>>>>   -ts_adapt_monitor
>>>>>
>>>>>
>>>>> > On May 1, 2017, at 7:51 PM, Matthew Knepley <knepley at gmail.com>
>>>>> wrote:
>>>>> >
>>>>> > Run with -snes_converged_reason.
>>>>> >
>>>>> >    Matt
>>>>> >
>>>>> > On Mon, May 1, 2017 at 7:14 PM, Mark Adams <mfadams at lbl.gov> wrote:
>>>>> > 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?
>>>>> >
>>>>> > 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
>>>>> >  ....
>>>>> >
>>>>> >    0 SNES Function norm 4.097052680599e+00
>>>>> >     1 SNES Function norm 1.213148652908e-09
>>>>> > [0]PETSC ERROR: --------------------- Error Message
>>>>> --------------------------------------------------------------
>>>>> > [0]PETSC ERROR:
>>>>> > [0]PETSC ERROR: TSStep has failed due to DIVERGED_NONLINEAR_SOLVE,
>>>>> increase -ts_max_snes_failures or make negative to attempt recovery
>>>>> > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/d
>>>>> ocumentation/faq.html for trouble shooting.
>>>>> > [0]PETSC ERROR: Petsc Development GIT revision:
>>>>> v3.7.6-3659-g699918129c  GIT Date: 2017-04-26 08:18:35 -0400
>>>>> > [0]PETSC ERROR: ./vml on a arch-macosx-gnu-O named MarksMac-5.local
>>>>> by markadams Mon May  1 19:21:32 2017
>>>>> > [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
>>>>> >
>>>>> >
>>>>> >
>>>>> >
>>>>> > --
>>>>> > What most experimenters take for granted before they begin their
>>>>> experiments is infinitely more interesting than any results to which their
>>>>> experiments lead.
>>>>> > -- Norbert Wiener
>>>>>
>>>>>
>>>>
>>>
>>>
>>> --
>>> What most experimenters take for granted before they begin their
>>> experiments is infinitely more interesting than any results to which their
>>> experiments lead.
>>> -- Norbert Wiener
>>>
>>
>>
>


-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20170515/103e814c/attachment-0001.html>

More information about the petsc-users mailing list