<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 28, 2013 at 3:00 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div class="im">Christophe Ortiz <<a href="mailto:christophe.ortiz@ciemat.es">christophe.ortiz@ciemat.es</a>> writes:<br>

> But just to understand what I'm doing, what does TSSetEquationType do<br>
> internally ?<br>
<br>
</div>It tells the method what is necessary to evaluate stage derivatives, and<br>
affects the starting method and (possibly) what can be carried over from<br>
the last step.<br>
<div class="im"><br></div></blockquote><div><br></div><div>Ok, I understand. For the moment it works well with this option.</div><div><br></div><div>I'm still doing some tests and I would like to know if it's possible to view a solution in log scale ?</div>
<div>To view the solution after diffusion I use</div><div> ierr = PetscViewerDrawOpen(PETSC_COMM_WORLD,NULL,"",0,0,300,300,&viewer);<br></div><div>VecView(X,viewer);</div><div><br></div><div>Since the values of X are large (~1e20) and vary logarithmically, I would like to view them in log scale. How can I do that ?</div>
<div><br></div><div>Also, is it possible to view the "real coordinates" in x-axis, instead of the node number ?</div><div><br></div><div>Christophe</div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div class="im">
> For the type of problem I'm trying to solve, is TS_EQ_DAE_IMPLICIT_INDEX1<br>
> ok ?<br>
<br>
</div>Yes.<br>
<div class="im"><br>
> What is the difference with TS_EQ_DAE_IMPLICIT_INDEX2,3 or _INDEXHI ?<br>
<br>
</div>I believe these are currently unused, but if we add support for<br>
computing consistent initial conditions, we'll have to know the index.<br>
</blockquote></div><br></div></div>