VTK output from DA vectors

Matthew Knepley knepley at gmail.com
Sun Jul 12 16:48:25 CDT 2009


1) I believe the standard VTK viewer (ASCII Viewer with VTK format) gets the
coordinates
    right. Can you verify this?

2) You should be able to split the 4- component field into 4 fields using
the Split filter in VTK
     or whatever viewer you use (I do this in Mayavi2, but Paraview also
works).

  Matt

On Sun, Jul 12, 2009 at 10:34 AM, Saswata Hier-Majumder <saswata at umd.edu>wrote:

> Hi,
> I would like to generate a vtk output from   a multicomponent problem. In
> the vtk file, I would like the DA coordinates as well as all 4 components
> stored separately as scalar point data.
>
> I have been using the  VecView_VTK routine from
> /ksp/ksp/ksp/examples/tutorials/ex29.c. But the vtk output seems to contain
> only the local coordinates (may be because DAGetCoordinates is not
> collective?).  Is there a way to fix this?
>
> Also, using the same routine, all components of the solution corresponding
> to a node are dumped together. Is there a way to extract each component
> separately and ouput them separately as scalar point data?
>
> Thanks
>
> --
> www.geol.umd.edu/~saswata <http://www.geol.umd.edu/%7Esaswata>
>
>


-- 
What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090712/fa1f04a0/attachment.htm>


More information about the petsc-users mailing list