[petsc-dev] Problem mixing sieve and TS
Paul Sicsic
paul.sicsic at polytechnique.edu
Thu Oct 27 00:09:46 CDT 2011
- The dot variable is still Nan when trying to compute the IFunction.
- These options seem to still be giving a shift that is infinity when
computing the IJacobian. (In step 0 here)
> Something must have computed a step size of 0 (or we're doing the wrong
> thing with a method that has explicit stages; we are changing some things
> with that). Of the recent methods that we have been playing with, I would
> recommend starting with "-ts_type rosw -ts_rosw_type ra3pw" or
>
[0]PETSC ERROR: MatAXPY() line 30 in /opt/HPC/petsc-dev/src/mat/utils/axpy.c
Invalid Pointer to Object: Parameter # 3
Program received signal SIGABRT, Aborted.
0x00007f5b55c94a75 in *__GI_raise (sig=<value optimized out>) at
../nptl/sysdeps/unix/sysv/linux/raise.c:64
64 ../nptl/sysdeps/unix/sysv/linux/raise.c: No such file or directory.
in ../nptl/sysdeps/unix/sysv/linux/raise.c
(gdb) where
#0 0x00007f5b55c94a75 in *__GI_raise (sig=<value optimized out>) at
../nptl/sysdeps/unix/sysv/linux/raise.c:64
#1 0x00007f5b55c985c0 in *__GI_abort () at abort.c:92
#2 0x00007f5b589e52a3 in PetscAbortErrorHandler (comm=1140850689, line=30,
fun=0x7f5b597a0e5d "MatAXPY", file=0x7f5b597a0b30
"/opt/HPC/petsc-dev/src/mat/utils/axpy.c",
dir=0x7f5b597a0b2d "", n=64, p=PETSC_ERROR_INITIAL, mess=0x7fffd7a3f2a0
"Invalid Pointer to Object: Parameter # 3", ctx=0x0) at
/opt/HPC/petsc-dev/src/sys/error/errabort.c:62
#3 0x00007f5b589dfd59 in PetscError (comm=1140850689, line=30,
func=0x7f5b597a0e5d "MatAXPY", file=0x7f5b597a0b30
"/opt/HPC/petsc-dev/src/mat/utils/axpy.c", dir=0x7f5b597a0b2d "", n=64,
p=PETSC_ERROR_INITIAL, mess=0x7f5b597a0b78 "Invalid Pointer to Object:
Parameter # %d") at /opt/HPC/petsc-dev/src/sys/error/err.c:356
#4 0x00007f5b5906d7b6 in MatAXPY (Y=0x26f3500, a=inf, X=0x282583000007f5b,
str=DIFFERENT_NONZERO_PATTERN) at /opt/HPC/petsc-dev/src/mat/utils/axpy.c:30
#5 0x00007f5b59092242 in mataxpy_ (Y=0x28aee00, a=0x7fffd7a3ffe8,
X=0x7fffd7a40acc, str=0x7fffd7a40ac4, __ierr=0x7fffd7a42af0) at
/opt/HPC/petsc-dev/src/mat/utils/ftn-auto/axpyf.c:65
#6 0x000000000040daed in ijacobianpoisson (dummyts=...,
t=2.4046180743743891e-57, u=..., udot=..., a=inf, jac=..., prejac=...,
appctx=..., ierr=0) at m_TransientHeat.F90:452
#7 0x00007f5b595c911d in ourijacobian (ts=0x2891cb0,
d=2.4046180743743891e-57, x=0x2958040, xdot=0x2954a00, shift=inf,
m=0x28aee00, p=0x28aee08, type=0x7fffd7a403d4, ctx=0x7fffd7a42af0)
at /opt/HPC/petsc-dev/src/ts/interface/ftn-custom/ztsf.c:82
#8 0x00007f5b5961a2f9 in TSComputeIJacobian (ts=0x2891cb0,
t=2.4046180743743891e-57, X=0x2958040, Xdot=0x2954a00, shift=inf,
A=0x28aee00, B=0x28aee08, flg=0x7fffd7a403d4, imex=PETSC_TRUE)
at /opt/HPC/petsc-dev/src/ts/interface/ts.c:459
#9 0x00007f5b595ee0d7 in SNESTSFormJacobian_RosW (snes=0x28aebf0,
X=0x2949830, A=0x28aee00, B=0x28aee08, str=0x7fffd7a403d4, ts=0x2891cb0) at
/opt/HPC/petsc-dev/src/ts/impls/rosw/rosw.c:781
#10 0x00007f5b5962e35c in SNESTSFormJacobian (snes=0x28aebf0, X=0x2949830,
A=0x28aee00, B=0x28aee08, flag=0x7fffd7a403d4, ctx=0x2891cb0) at
/opt/HPC/petsc-dev/src/ts/interface/ts.c:2578
#11 0x00007f5b59543e8b in SNESComputeJacobian (snes=0x28aebf0, X=0x2949830,
A=0x28aee00, B=0x28aee08, flg=0x7fffd7a403d4) at
/opt/HPC/petsc-dev/src/snes/interface/snes.c:1415
#12 0x00007f5b595710cc in SNESSolve_KSPONLY (snes=0x28aebf0) at
/opt/HPC/petsc-dev/src/snes/impls/ksponly/ksponly.c:32
#13 0x00007f5b59552791 in SNESSolve (snes=0x28aebf0, b=0x0, x=0x2949830) at
/opt/HPC/petsc-dev/src/snes/interface/snes.c:2850
#14 0x00007f5b595ec791 in TSStep_RosW (ts=0x2891cb0) at
/opt/HPC/petsc-dev/src/ts/impls/rosw/rosw.c:671
#15 0x00007f5b59627411 in TSStep (ts=0x2891cb0) at
/opt/HPC/petsc-dev/src/ts/interface/ts.c:1813
#16 0x00007f5b59628860 in TSSolve (ts=0x2891cb0, x=0x25bbb50,
ftime=0x7fffd7a433f0) at /opt/HPC/petsc-dev/src/ts/interface/ts.c:1933
#17 0x00007f5b595c8a08 in tssolve_ (ts=0x7fffd7a433e0, x=0x7fffd7a43088,
ftime=0x7fffd7a433f0, __ierr=0x7fffd7a427dc) at
/opt/HPC/petsc-dev/src/ts/interface/ftn-auto/tsf.c:397
#18 0x000000000040d75c in solvetransient (appctx=...) at
m_TransientHeat.F90:493
#19 0x000000000040be2e in transientheat () at TransientHeat.F90:98
#20 0x000000000047d5ea in main ()
#21 0x00007f5b55c7fc4d in __libc_start_main (main=<value optimized out>,
argc=<value optimized out>, ubp_av=<value optimized out>, init=<value
optimized out>, fini=<value optimized out>,
rtld_fini=<value optimized out>, stack_end=0x7fffd7a43588) at
libc-start.c:226
#22 0x00000000004061f9 in _start ()
> "-ts_type rosw -ts_rosw_type ra34pw2", both of which have nice stability
> properties, good error estimates, and are more permissive than the ARKs
> about how you write your equations.
>
[0]PETSC ERROR: MatAXPY() line 30 in /opt/HPC/petsc-dev/src/mat/utils/axpy.c
Invalid Pointer to Object: Parameter # 3
Program received signal SIGABRT, Aborted.
0x00007f8fe1adba75 in *__GI_raise (sig=<value optimized out>) at
../nptl/sysdeps/unix/sysv/linux/raise.c:64
64 ../nptl/sysdeps/unix/sysv/linux/raise.c: No such file or directory.
in ../nptl/sysdeps/unix/sysv/linux/raise.c
(gdb) where
#0 0x00007f8fe1adba75 in *__GI_raise (sig=<value optimized out>) at
../nptl/sysdeps/unix/sysv/linux/raise.c:64
#1 0x00007f8fe1adf5c0 in *__GI_abort () at abort.c:92
#2 0x00007f8fe482c2a3 in PetscAbortErrorHandler (comm=1140850689, line=30,
fun=0x7f8fe55e7e5d "MatAXPY", file=0x7f8fe55e7b30
"/opt/HPC/petsc-dev/src/mat/utils/axpy.c",
dir=0x7f8fe55e7b2d "", n=64, p=PETSC_ERROR_INITIAL, mess=0x7fff2d736490
"Invalid Pointer to Object: Parameter # 3", ctx=0x0) at
/opt/HPC/petsc-dev/src/sys/error/errabort.c:62
#3 0x00007f8fe4826d59 in PetscError (comm=1140850689, line=30,
func=0x7f8fe55e7e5d "MatAXPY", file=0x7f8fe55e7b30
"/opt/HPC/petsc-dev/src/mat/utils/axpy.c", dir=0x7f8fe55e7b2d "", n=64,
p=PETSC_ERROR_INITIAL, mess=0x7f8fe55e7b78 "Invalid Pointer to Object:
Parameter # %d") at /opt/HPC/petsc-dev/src/sys/error/err.c:356
#4 0x00007f8fe4eb47b6 in MatAXPY (Y=0x22ab500, a=inf, X=0x23dd83000007f8f,
str=DIFFERENT_NONZERO_PATTERN) at /opt/HPC/petsc-dev/src/mat/utils/axpy.c:30
#5 0x00007f8fe4ed9242 in mataxpy_ (Y=0x2466e00, a=0x7fff2d7371d8,
X=0x7fff2d737cbc, str=0x7fff2d737cb4, __ierr=0x7fff2d739ce0) at
/opt/HPC/petsc-dev/src/mat/utils/ftn-auto/axpyf.c:65
#6 0x000000000040daed in ijacobianpoisson (dummyts=...,
t=2.4046180743743891e-57, u=..., udot=..., a=inf, jac=..., prejac=...,
appctx=..., ierr=0) at m_TransientHeat.F90:452
#7 0x00007f8fe541011d in ourijacobian (ts=0x2449cb0,
d=2.4046180743743891e-57, x=0x2513680, xdot=0x2510040, shift=inf,
m=0x2466e00, p=0x2466e08, type=0x7fff2d7375c4, ctx=0x7fff2d739ce0)
at /opt/HPC/petsc-dev/src/ts/interface/ftn-custom/ztsf.c:82
#8 0x00007f8fe54612f9 in TSComputeIJacobian (ts=0x2449cb0,
t=2.4046180743743891e-57, X=0x2513680, Xdot=0x2510040, shift=inf,
A=0x2466e00, B=0x2466e08, flg=0x7fff2d7375c4, imex=PETSC_TRUE)
at /opt/HPC/petsc-dev/src/ts/interface/ts.c:459
#9 0x00007f8fe54350d7 in SNESTSFormJacobian_RosW (snes=0x2466bf0,
X=0x2501830, A=0x2466e00, B=0x2466e08, str=0x7fff2d7375c4, ts=0x2449cb0) at
/opt/HPC/petsc-dev/src/ts/impls/rosw/rosw.c:781
#10 0x00007f8fe547535c in SNESTSFormJacobian (snes=0x2466bf0, X=0x2501830,
A=0x2466e00, B=0x2466e08, flag=0x7fff2d7375c4, ctx=0x2449cb0) at
/opt/HPC/petsc-dev/src/ts/interface/ts.c:2578
#11 0x00007f8fe538ae8b in SNESComputeJacobian (snes=0x2466bf0, X=0x2501830,
A=0x2466e00, B=0x2466e08, flg=0x7fff2d7375c4) at
/opt/HPC/petsc-dev/src/snes/interface/snes.c:1415
#12 0x00007f8fe53b80cc in SNESSolve_KSPONLY (snes=0x2466bf0) at
/opt/HPC/petsc-dev/src/snes/impls/ksponly/ksponly.c:32
#13 0x00007f8fe5399791 in SNESSolve (snes=0x2466bf0, b=0x0, x=0x2501830) at
/opt/HPC/petsc-dev/src/snes/interface/snes.c:2850
#14 0x00007f8fe5433791 in TSStep_RosW (ts=0x2449cb0) at
/opt/HPC/petsc-dev/src/ts/impls/rosw/rosw.c:671
#15 0x00007f8fe546e411 in TSStep (ts=0x2449cb0) at
/opt/HPC/petsc-dev/src/ts/interface/ts.c:1813
#16 0x00007f8fe546f860 in TSSolve (ts=0x2449cb0, x=0x2173b50,
ftime=0x7fff2d73a5e0) at /opt/HPC/petsc-dev/src/ts/interface/ts.c:1933
#17 0x00007f8fe540fa08 in tssolve_ (ts=0x7fff2d73a5d0, x=0x7fff2d73a278,
ftime=0x7fff2d73a5e0, __ierr=0x7fff2d7399cc) at
/opt/HPC/petsc-dev/src/ts/interface/ftn-auto/tsf.c:397
#18 0x000000000040d75c in solvetransient (appctx=...) at
m_TransientHeat.F90:493
#19 0x000000000040be2e in transientheat () at TransientHeat.F90:98
#20 0x000000000047d5ea in main ()
#21 0x00007f8fe1ac6c4d in __libc_start_main (main=<value optimized out>,
argc=<value optimized out>, ubp_av=<value optimized out>, init=<value
optimized out>, fini=<value optimized out>,
rtld_fini=<value optimized out>, stack_end=0x7fff2d73a778) at
libc-start.c:226
#22 0x00000000004061f9 in _start ()
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20111027/07aad319/attachment.html>
More information about the petsc-dev
mailing list