[petsc-users] SNES error
Matthew Knepley
knepley at gmail.com
Tue May 2 10:18:53 CDT 2017
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/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 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20170502/425b3bd1/attachment.html>
More information about the petsc-users
mailing list