[petsc-users] TSStep error

[Jin, Shuangshuang]

Thanks, Shri, that's a really good link. Very useful!

I have identified the error according to the instruction. The Jacobian was wrong.

Thanks,
Shuangshuang

From: petsc-users-bounces at mcs.anl.gov [mailto:petsc-users-bounces at mcs.anl.gov] On Behalf Of Shri
Sent: Thursday, March 28, 2013 8:51 AM
To: PETSc users list
Subject: Re: [petsc-users] TSStep error


On Mar 27, 2013, at 8:18 PM, Jin, Shuangshuang wrote:


Hi, I am using TS to solve a nonlinear DAE problem. I got an error below:

[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: ------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Development HG revision: 6e0ddc6e9b6d8a9d8eb4c0ede0105827a6b58dfb  HG Date: Mon Mar 11 22:54:30 2013 -0500
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: ./dynSim on a arch-complex named olympus.local by d3m956 Wed Mar 27 18:12:05 2013
[0]PETSC ERROR: Libraries linked from /pic/projects/ds/petsc-dev/arch-complex/lib
[0]PETSC ERROR: Configure run at Tue Mar 12 14:32:37 2013
[0]PETSC ERROR: Configure options --with-scalar-type=complex --with-clanguage=C++ PETSC_ARCH=arch-complex --with-fortran-kernels=generic
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: TSStep() line 2442 in src/ts/interface/ts.c
[0]PETSC ERROR: TSSolve() line 2553 in src/ts/interface/ts.c
[0]PETSC ERROR: simu() line 420 in "unknowndirectory/"simulation.C
[0]PETSC ERROR: runSimulation() line 83 in "unknowndirectory/"dynSim.h
Run simulation time: 0.0317168

What does this error imply for?
The error implies that the nonlinear solve (Newton's method), used at each time step, did not converge. This could be due to various reasons as described here
http://www.mcs.anl.gov/petsc/documentation/faq.html#newton


My ftime = 0.5, and the timestep is 0.001, and the solution method I'm using is TSBEULER
PetscReal ftime=0.5;
ierr = TSSetInitialTimeStep(ts, 0.0, 0.001); CHKERRQ(ierr);
ierr = TSSetType(ts, TSBEULER); CHKERRQ(ierr);

Is there any command option I should use when I run the code?

Thanks,
Shuangshuang

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130328/2ef9ce63/attachment.html>