VTK output from DA vectors

[Saswata Hier-Majumder]

1) I did. It returns correct values for the x nodes, but fails to do so 
for the y nodes. Here's a sample output generated from the driven cavity 
problem. ( I suppressed printing the solutions for brevity).

2) Thanks, I'll try that.

Matthew Knepley wrote:
> 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>
>>
>>
>>     
>
>
>   

-- 
www.geol.umd.edu/~saswata

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: solution.vtk
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090712/02cab317/attachment.diff>