[petsc-users] TS_SUNDIALS/TS_PSEUDO overshoots the final time

Fri Oct 31 11:49:10 CDT 2014

Hi Jed,

On Thu, Oct 30, 2014 at 9:22 PM, Jed Brown <jed at jedbrown.org> wrote:

> Gautam Bisht <gbisht at lbl.gov> writes:
>
> > All,
> >
> > While running src/ts/examples/tutorials/ex8.c (in PETSc 3.5.1) for TSTYPE
> > sundials and pseudo, I noticed that the solver overshoots the final time.
>
> Is this what you're looking for?
>
>
> http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/TS/TSSetExactFinalTime.html

Yes. But it appears only '-ts_exact_final_time INTERPOLATE' for '-ts_type
sundials' does not overshoot the final time.

./ex8 -ts_monitor -ts_type pseudo -ts_final_time 0.01  -ts_exact_final_time
STEPOVER
0 TS dt 0.001 time 0
1 TS dt 0.0011 time 0.0011
2 TS dt 0.00122331 time 0.00232331
3 TS dt 0.00136211 time 0.00368542
4 TS dt 0.00151873 time 0.00520415
5 TS dt 0.00169599 time 0.00690014
6 TS dt 0.00189725 time 0.00879739
7 TS dt 0.0021266 time 0.010924
steps 7 (0 rejected, 0 SNES fails), ftime 0.010924, nonlinits 7, linits 7
step    7 t=1.09239877e-02 h= 2.13e-03  |x|= 1.03e-01  |x_e|= 1.03e-01
 |x-x_e|= 7.93e-04

./ex8 -ts_monitor -ts_type pseudo -ts_final_time 0.01 -ts_exact_final_time
MATCHSTEP
0 TS dt 0.001 time 0
1 TS dt 0.0011 time 0.0011
2 TS dt 0.00122331 time 0.00232331
3 TS dt 0.00136211 time 0.00368542
4 TS dt 0.00151873 time 0.00520415
5 TS dt 0.00169599 time 0.00690014
6 TS dt 0.00189725 time 0.00879739
7 TS dt 0.0021266 time 0.010924
steps 7 (0 rejected, 0 SNES fails), ftime 0.010924, nonlinits 7, linits 7
step    7 t=1.09239877e-02 h= 2.13e-03  |x|= 1.03e-01  |x_e|= 1.03e-01
 |x-x_e|= 7.93e-04

./ex8 -ts_monitor -ts_type pseudo -ts_final_time 0.01 -ts_exact_final_time
INTERPOLATE
0 TS dt 0.001 time 0
1 TS dt 0.0011 time 0.0011
2 TS dt 0.00122331 time 0.00232331
3 TS dt 0.00136211 time 0.00368542
4 TS dt 0.00151873 time 0.00520415
5 TS dt 0.00169599 time 0.00690014
6 TS dt 0.00189725 time 0.00879739
[0]PETSC ERROR: --------------------- Error Message
--------------------------------------------------------------
[0]PETSC ERROR: No support for this operation for this object type
[0]PETSC ERROR: pseudo does not provide interpolation
[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for
trouble shooting.
[0]PETSC ERROR: Petsc Release Version 3.5.1, unknown
[0]PETSC ERROR: ./ex8 on a darwin-gnu-fort-debug named gmbook.dhcp.lbnl.us
by gbisht Fri Oct 31 09:44:36 2014
[0]PETSC ERROR: Configure options -with-hdf5-dir=/opt/local/
--with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate
--download-parmetis=yes --download-metis=yes --download-sundials=yes
--with-mpich=/opt/local --with-cc=/opt/local/bin/mpicc-mpich-gcc48
--with-cxx=/opt/local/bin/mpicxx-mpich-gcc48
--with-fc=/opt/local/bin/mpifort-mpich-gcc48
PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1
[0]PETSC ERROR: #1 TSInterpolate() line 2582 in
/Users/gbisht/Research/Models/PETSc/petsc3.5.1/src/ts/interface/ts.c
[0]PETSC ERROR: #2 TSSolve() line 2763 in
/Users/gbisht/Research/Models/PETSc/petsc3.5.1/src/ts/interface/ts.c
[0]PETSC ERROR: #3 main() line 436 in
/Users/gbisht/Research/Models/PETSc/petsc3.5.1/src/ts/examples/tutorials/ex8.c
[0]PETSC ERROR: ----------------End of Error Message -------send entire
error message to petsc-maint at mcs.anl.gov----------
application called MPI_Abort(MPI_COMM_WORLD, 56) - process 0

./ex8 -ts_monitor -ts_type sundials -ts_final_time 0.01
-ts_exact_final_time STEPOVER
0 TS dt 0.001 time 0
1 TS dt 0.0001 time 0.0001
2 TS dt 0.0001 time 0.0002
3 TS dt 0.0001 time 0.0003
4 TS dt 0.000199571 time 0.000499571
5 TS dt 0.000199571 time 0.000699141
6 TS dt 0.000199571 time 0.000898712
7 TS dt 0.000199571 time 0.00109828
8 TS dt 0.000404911 time 0.00150319
9 TS dt 0.000804269 time 0.00230746
10 TS dt 0.00141825 time 0.00372572
11 TS dt 0.00141825 time 0.00514397
12 TS dt 0.00141825 time 0.00656222
13 TS dt 0.00141825 time 0.00798047
14 TS dt 0.00141825 time 0.00939873
15 TS dt 0.00193698 time 0.0113357
steps 15 (0 rejected, 0 SNES fails), ftime 0.0113357, nonlinits 18, linits
18
step   15 t=1.13357082e-02 h= 1.94e-03  |x|= 1.07e-01  |x_e|= 1.07e-01
 |x-x_e|= 2.57e-06

./ex8 -ts_monitor -ts_type sundials -ts_final_time 0.01
-ts_exact_final_time MATCHSTEP
0 TS dt 0.001 time 0
1 TS dt 0.0001 time 0.0001
2 TS dt 0.0001 time 0.0002
3 TS dt 0.0001 time 0.0003
4 TS dt 0.000199571 time 0.000499571
5 TS dt 0.000199571 time 0.000699141
6 TS dt 0.000199571 time 0.000898712
7 TS dt 0.000199571 time 0.00109828
8 TS dt 0.000404911 time 0.00150319
9 TS dt 0.000804269 time 0.00230746
10 TS dt 0.00141825 time 0.00372572
11 TS dt 0.00141825 time 0.00514397
12 TS dt 0.00141825 time 0.00656222
13 TS dt 0.00141825 time 0.00798047
14 TS dt 0.00141825 time 0.00939873
15 TS dt 0.00193698 time 0.0113357
steps 15 (0 rejected, 0 SNES fails), ftime 0.0113357, nonlinits 18, linits
18
step   15 t=1.13357082e-02 h= 1.94e-03  |x|= 1.07e-01  |x_e|= 1.07e-01
 |x-x_e|= 2.57e-06

./ex8 -ts_monitor -ts_type sundials -ts_final_time 0.01
-ts_exact_final_time INTERPOLATE
0 TS dt 0.001 time 0
1 TS dt 0.0001 time 0.0001
2 TS dt 0.0001 time 0.0002
3 TS dt 0.0001 time 0.0003
4 TS dt 0.000199571 time 0.000499571
5 TS dt 0.000199571 time 0.000699141
6 TS dt 0.000199571 time 0.000898712
7 TS dt 0.000199571 time 0.00109828
8 TS dt 0.000404911 time 0.00150319
9 TS dt 0.000804269 time 0.00230746
10 TS dt 0.00141825 time 0.00372572
11 TS dt 0.00141825 time 0.00514397
12 TS dt 0.00141825 time 0.00656222
13 TS dt 0.00141825 time 0.00798047
14 TS dt 0.00141825 time 0.00939873
15 TS dt 0.00193698 time 0.01
steps 15 (0 rejected, 0 SNES fails), ftime 0.01, nonlinits 18, linits 18
step   15 t=1.00000000e-02 h= 1.94e-03  |x|= 9.52e-02  |x_e|= 9.52e-02
 |x-x_e|= 2.11e-06

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