[petsc-users] TSGL failure

[Mike Sussman]

Folks,

The tutorial example in TS, ex2f.F, is designed to use the theta method.
If, however, I use the runtime option "-ts_type gl", I get a segfault.
I also get a segfault with the line
CALL TSSetType(ts,TSGL,ierr);
in the code.

Is this supposed to happen?  What am I missing?

I am running 3.2p5.

The message says:
% ex2f -ts_type gl -ts_gl_type irks -ts_gl_rtol 1.e-5 -ts_gl_atol 1.e-10
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation,
probably memory access out of range
[0]PETSC ERROR: Try option -start_in_debugger or
-on_error_attach_debugger
[0]PETSC ERROR: or see
http://www.mcs.anl.gov/petsc/petsc-as/documentation/faq.html#valgrind[0]PETSC ERROR: or try http://valgrind.org on GNU/linux and Apple Mac OS X to find memory corruption errors
[0]PETSC ERROR: likely location of problem given in stack below
[0]PETSC ERROR: ---------------------  Stack Frames
------------------------------------
[0]PETSC ERROR: Note: The EXACT line numbers in the stack are not
available,
[0]PETSC ERROR:       INSTEAD the line number of the start of the
function
[0]PETSC ERROR:       is given.
[0]PETSC ERROR: [0] TSGLAdaptChoose line
232 /home/mike/Programming/petsc-3.2-p5/src/ts/impls/implicit/gl/gladapt.c
[0]PETSC ERROR: [0] TSGLChooseNextScheme line
774 /home/mike/Programming/petsc-3.2-p5/src/ts/impls/implicit/gl/gl.c
[0]PETSC ERROR: [0] TSSolve_GL line
822 /home/mike/Programming/petsc-3.2-p5/src/ts/impls/implicit/gl/gl.c
[0]PETSC ERROR: [0] TSSolve line
1824 /home/mike/Programming/petsc-3.2-p5/src/ts/interface/ts.c
[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: Signal received!
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.2.0, Patch 5, Sat Oct 29
13:45:54 CDT 2011 
[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: ex2f on a linux-gnu named ozhp by mike Mon Nov 28
19:38:30 2011
[0]PETSC ERROR: Libraries linked
from /home/mike/Programming/petsc-3.2-p5/linux-gnu-c-opt/lib
[0]PETSC ERROR: Configure run at Sat Nov 26 21:07:02 2011
[0]PETSC ERROR: Configure options --with-cc=mpicc --with-fc=mpif90
--with-cxx=mpicxx --download-hypre=1
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: User provided function() line 0 in unknown directory
unknown file
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD 
with errorcode 59.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on


_______
Mike Sussman