[petsc-users] Bug in boundary conditions (ex25.c) ?

[Christophe Ortiz]

On Mon, Oct 28, 2013 at 3:00 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:

> Christophe Ortiz <christophe.ortiz at ciemat.es> writes:
> > But just to understand what I'm doing, what does TSSetEquationType do
> > internally ?
>
> It tells the method what is necessary to evaluate stage derivatives, and
> affects the starting method and (possibly) what can be carried over from
> the last step.
>
>
Ok, I understand. For the moment it works well with this option.

I'm still doing some tests and I would like to know if it's possible to
view a solution in log scale ?
To view the solution after diffusion I use
 ierr = PetscViewerDrawOpen(PETSC_COMM_WORLD,NULL,"",0,0,300,300,&viewer);
VecView(X,viewer);

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 ?

Also, is it possible to view the "real coordinates" in x-axis, instead of
the node number ?

Christophe



> > For the type of problem I'm trying to solve, is TS_EQ_DAE_IMPLICIT_INDEX1
> > ok ?
>
> Yes.
>
> > What is the difference with TS_EQ_DAE_IMPLICIT_INDEX2,3 or _INDEXHI ?
>
> I believe these are currently unused, but if we add support for
> computing consistent initial conditions, we'll have to know the index.
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20131029/df296c5c/attachment.html>