<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><br class="Apple-interchange-newline"><blockquote type="cite"><div class="gmail_quote">On Wed, Oct 26, 2011 at 12:45, Blaise Bourdin <span dir="ltr"><<a href="mailto:bourdin@lsu.edu">bourdin@lsu.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0.8ex; border-left-width: 1px; border-left-color: rgb(204, 204, 204); border-left-style: solid; padding-left: 1ex; position: static; z-index: auto; ">
<div id=":1p5">Any suggestion? FWIW, in this test, both the IJacobian matrix and RHS Vector are constant, but this will not be the case once in the future.<br>
<br>
Something else bugging me is that the Jacobian function does not seem to be ever called (even in ex22 from src/ts/examples/tutorial).</div></blockquote><div><br></div><div>How are you coming to that conclusion? The default method for ex22 is a method with an implicit part. If you don't set any type at all, then the default method is forward Euler (we should probably change this, maybe to a Rosenbrock-W), which does not evaluate a Jacobian, but that code sets the default type to ARKIMEX which does evaluate the IJacobian.</div></div></blockquote><div>That is correct. In both cases, I was setting a breakpoint at the wrong location.</div><br><blockquote type="cite"><div class="gmail_quote">
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div id=":1p5"> Is that expected? Is there a written description (or link to the literature) for each of the ODE solvers implemented in TS?</div>
</blockquote></div><br><div>The man page for each method should have a citation to a paper that introduced it or discussed it. Some of the methods in there were designed by Emil and have not been published yet.</div>
</blockquote>Thanks,</div><div>I had look at the manual and the man page for ARKIMEX but not the individual pages.</div><div><br></div><div>At this point we have overcome the original issue but are still having difficulties. We followed the lines of ex22 to setup our TS, with the most notable difference being that our DM is a DMMesh and not a DMDA. For some reason, the second time the Ifunction is called, the Udot argument is a Vec of NaN. We checked the TS in the debugger and dan't catch anything obviously wrong. We followed the flow of the code in the debugger, and the NaN get generated here:</div><div><br></div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">(gdb) </div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">713 ierr = VecAXPBYPCZ(ark->Ydot,-ark->shift,ark->shift,0,ark->Z,X);CHKERRQ(ierr); /* Ydot = shift*(X-Z) */</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">(gdb) where</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); position: static; z-index: auto; ">#0 SNESTSFormFunction_ARKIMEX (snes=0x18b3850, X=0x1951230, F=0x1919800, ts=0x1896910) at /opt/HPC/petsc-dev/src/ts/impls/arkimex/arkimex.c:713</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); position: static; z-index: auto; ">#1 0x00007f707689675b in SNESTSFormFunction (snes=0x18b3850, X=0x1951230, F=0x1919800, ctx=0x1896910) at /opt/HPC/petsc-dev/src/ts/interface/ts.c:2536</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">#2 0x00007f70767ab7f2 in SNESComputeFunction (snes=0x18b3850, x=0x1951230, y=0x1919800) at /opt/HPC/petsc-dev/src/snes/interface/snes.c:1314</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">#3 0x00007f70767e7c2a in SNESSolve_LS (snes=0x18b3850) at /opt/HPC/petsc-dev/src/snes/impls/ls/ls.c:609</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">#4 0x00007f70767bb791 in SNESSolve (snes=0x18b3850, b=0x197ac70, x=0x1951230) at /opt/HPC/petsc-dev/src/snes/interface/snes.c:2850</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">#5 0x00007f707684a5df in TSStep_ARKIMEX (ts=0x1896910) at /opt/HPC/petsc-dev/src/ts/impls/arkimex/arkimex.c:576</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">#6 0x00007f7076890411 in TSStep (ts=0x1896910) at /opt/HPC/petsc-dev/src/ts/interface/ts.c:1813</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">#7 0x00007f7076891860 in TSSolve (ts=0x1896910, x=0x15c16b0, ftime=0x7fff05be8a00) at /opt/HPC/petsc-dev/src/ts/interface/ts.c:1933</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">#8 0x00007f7076831a08 in tssolve_ (ts=0x7fff05be89f0, x=0x7fff05be8698, ftime=0x7fff05be8a00, __ierr=0x7fff05be7dec) at /opt/HPC/petsc-dev/src/ts/interface/ftn-auto/tsf.c:397</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">#9 0x000000000040d75c in solvetransient (appctx=...) at m_TransientHeat.F90:493</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">#10 0x000000000040be2e in transientheat () at TransientHeat.F90:98</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">#11 0x000000000047d5ea in main ()</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">#12 0x00007f7072ee8c4d 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 style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); "> rtld_fini=<value optimized out>, stack_end=0x7fff05be8b98) at libc-start.c:226</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">#13 0x00000000004061f9 in _start ()</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); "><br></div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">Do you have any suggestion of potential causes? At this point, we are running out of ideas.</div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); "><br></div><div style="font-family: arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255); ">Blaise</div><div> </div><div><br></div><div apple-content-edited="true"><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">-- </div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Department of Mathematics and Center for Computation & Technology</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Louisiana State University, Baton Rouge, LA 70803, USA</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Tel. +1 (225) 578 1612, Fax +1 (225) 578 4276 <a href="http://www.math.lsu.edu/~bourdin">http://www.math.lsu.edu/~bourdin</a></div><div><br class="khtml-block-placeholder"></div><div><br class="khtml-block-placeholder"></div><br class="Apple-interchange-newline"></div></span></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></span><br class="Apple-interchange-newline">
</div>
<br></body></html>