<div class="gmail_quote"><div><ul><li>The dot variable is still Nan when trying to compute the IFunction. </li><li>These options seem to still be giving a shift that is infinity when computing the IJacobian. (In step 0 here)</li>
</ul></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div></div><div>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 </div>
</blockquote><div><br></div><div><div>[0]PETSC ERROR: MatAXPY() line 30 in /opt/HPC/petsc-dev/src/mat/utils/axpy.c Invalid Pointer to Object: Parameter # 3</div><div><br></div><div>Program received signal SIGABRT, Aborted.</div>
<div>0x00007f5b55c94a75 in *__GI_raise (sig=<value optimized out>) at ../nptl/sysdeps/unix/sysv/linux/raise.c:64</div><div>64 ../nptl/sysdeps/unix/sysv/linux/raise.c: No such file or directory.</div><div> in ../nptl/sysdeps/unix/sysv/linux/raise.c</div>
<div>(gdb) where</div><div>#0 0x00007f5b55c94a75 in *__GI_raise (sig=<value optimized out>) at ../nptl/sysdeps/unix/sysv/linux/raise.c:64</div><div>#1 0x00007f5b55c985c0 in *__GI_abort () at abort.c:92</div><div>
#2 0x00007f5b589e52a3 in PetscAbortErrorHandler (comm=1140850689, line=30, fun=0x7f5b597a0e5d "MatAXPY", file=0x7f5b597a0b30 "/opt/HPC/petsc-dev/src/mat/utils/axpy.c", </div>
<div> 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</div><div>#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, </div>
<div> p=PETSC_ERROR_INITIAL, mess=0x7f5b597a0b78 "Invalid Pointer to Object: Parameter # %d") at /opt/HPC/petsc-dev/src/sys/error/err.c:356</div><div>#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</div>
<div>#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</div><div>#6 0x000000000040daed in ijacobianpoisson (dummyts=..., t=2.4046180743743891e-57, u=..., udot=..., a=inf, jac=..., prejac=..., appctx=..., ierr=0) at m_TransientHeat.F90:452</div>
<div>#7 0x00007f5b595c911d in ourijacobian (ts=0x2891cb0, d=2.4046180743743891e-57, x=0x2958040, xdot=0x2954a00, shift=inf, m=0x28aee00, p=0x28aee08, type=0x7fffd7a403d4, ctx=0x7fffd7a42af0)</div><div> at /opt/HPC/petsc-dev/src/ts/interface/ftn-custom/ztsf.c:82</div>
<div>#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)</div><div> at /opt/HPC/petsc-dev/src/ts/interface/ts.c:459</div>
<div>#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</div><div>#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</div>
<div>#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</div><div>#12 0x00007f5b595710cc in SNESSolve_KSPONLY (snes=0x28aebf0) at /opt/HPC/petsc-dev/src/snes/impls/ksponly/ksponly.c:32</div>
<div>#13 0x00007f5b59552791 in SNESSolve (snes=0x28aebf0, b=0x0, x=0x2949830) at /opt/HPC/petsc-dev/src/snes/interface/snes.c:2850</div><div>#14 0x00007f5b595ec791 in TSStep_RosW (ts=0x2891cb0) at /opt/HPC/petsc-dev/src/ts/impls/rosw/rosw.c:671</div>
<div>#15 0x00007f5b59627411 in TSStep (ts=0x2891cb0) at /opt/HPC/petsc-dev/src/ts/interface/ts.c:1813</div><div>#16 0x00007f5b59628860 in TSSolve (ts=0x2891cb0, x=0x25bbb50, ftime=0x7fffd7a433f0) at /opt/HPC/petsc-dev/src/ts/interface/ts.c:1933</div>
<div>#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</div><div>#18 0x000000000040d75c in solvetransient (appctx=...) at m_TransientHeat.F90:493</div>
<div>#19 0x000000000040be2e in transientheat () at TransientHeat.F90:98</div><div>#20 0x000000000047d5ea in main ()</div><div>#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>, </div>
<div> rtld_fini=<value optimized out>, stack_end=0x7fffd7a43588) at libc-start.c:226</div><div>#22 0x00000000004061f9 in _start ()</div></div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>"-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.</div>
</blockquote></div><br><div>[0]PETSC ERROR: MatAXPY() line 30 in /opt/HPC/petsc-dev/src/mat/utils/axpy.c Invalid Pointer to Object: Parameter # 3</div><div><br></div><div>Program received signal SIGABRT, Aborted.</div><div>
0x00007f8fe1adba75 in *__GI_raise (sig=<value optimized out>) at ../nptl/sysdeps/unix/sysv/linux/raise.c:64</div><div>64 ../nptl/sysdeps/unix/sysv/linux/raise.c: No such file or directory.</div><div> in ../nptl/sysdeps/unix/sysv/linux/raise.c</div>
<div>(gdb) where</div><div>#0 0x00007f8fe1adba75 in *__GI_raise (sig=<value optimized out>) at ../nptl/sysdeps/unix/sysv/linux/raise.c:64</div><div>#1 0x00007f8fe1adf5c0 in *__GI_abort () at abort.c:92</div><div>
#2 0x00007f8fe482c2a3 in PetscAbortErrorHandler (comm=1140850689, line=30, fun=0x7f8fe55e7e5d "MatAXPY", file=0x7f8fe55e7b30 "/opt/HPC/petsc-dev/src/mat/utils/axpy.c", </div>
<div> 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</div><div>#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, </div>
<div> p=PETSC_ERROR_INITIAL, mess=0x7f8fe55e7b78 "Invalid Pointer to Object: Parameter # %d") at /opt/HPC/petsc-dev/src/sys/error/err.c:356</div><div>#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</div>
<div>#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</div><div>#6 0x000000000040daed in ijacobianpoisson (dummyts=..., t=2.4046180743743891e-57, u=..., udot=..., a=inf, jac=..., prejac=..., appctx=..., ierr=0) at m_TransientHeat.F90:452</div>
<div>#7 0x00007f8fe541011d in ourijacobian (ts=0x2449cb0, d=2.4046180743743891e-57, x=0x2513680, xdot=0x2510040, shift=inf, m=0x2466e00, p=0x2466e08, type=0x7fff2d7375c4, ctx=0x7fff2d739ce0)</div><div> at /opt/HPC/petsc-dev/src/ts/interface/ftn-custom/ztsf.c:82</div>
<div>#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)</div><div> at /opt/HPC/petsc-dev/src/ts/interface/ts.c:459</div>
<div>#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</div><div>#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</div>
<div>#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</div><div>#12 0x00007f8fe53b80cc in SNESSolve_KSPONLY (snes=0x2466bf0) at /opt/HPC/petsc-dev/src/snes/impls/ksponly/ksponly.c:32</div>
<div>#13 0x00007f8fe5399791 in SNESSolve (snes=0x2466bf0, b=0x0, x=0x2501830) at /opt/HPC/petsc-dev/src/snes/interface/snes.c:2850</div><div>#14 0x00007f8fe5433791 in TSStep_RosW (ts=0x2449cb0) at /opt/HPC/petsc-dev/src/ts/impls/rosw/rosw.c:671</div>
<div>#15 0x00007f8fe546e411 in TSStep (ts=0x2449cb0) at /opt/HPC/petsc-dev/src/ts/interface/ts.c:1813</div><div>#16 0x00007f8fe546f860 in TSSolve (ts=0x2449cb0, x=0x2173b50, ftime=0x7fff2d73a5e0) at /opt/HPC/petsc-dev/src/ts/interface/ts.c:1933</div>
<div>#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</div><div>#18 0x000000000040d75c in solvetransient (appctx=...) at m_TransientHeat.F90:493</div>
<div>#19 0x000000000040be2e in transientheat () at TransientHeat.F90:98</div><div>#20 0x000000000047d5ea in main ()</div><div>#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>, </div>
<div> rtld_fini=<value optimized out>, stack_end=0x7fff2d73a778) at libc-start.c:226</div><div>#22 0x00000000004061f9 in _start ()</div><div><br></div>