<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Jan 28, 2014 at 9:41 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im">Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> writes:<br>
> In ex62:<br>
><br>
>   if (user.runType == RUN_FULL) {<br>
>     PetscViewer viewer;<br>
>     Vec         uLocal;<br>
>     const char *name;<br>
><br>
>     ierr = PetscViewerCreate(PETSC_COMM_WORLD, &viewer);CHKERRQ(ierr);<br>
>     ierr = PetscViewerSetType(viewer, PETSCVIEWERVTK);CHKERRQ(ierr);<br>
>     ierr = PetscViewerSetFormat(viewer,<br>
> PETSC_VIEWER_ASCII_VTK);CHKERRQ(ierr);<br>
>     ierr = PetscViewerFileSetName(viewer, "ex62_sol.vtk");CHKERRQ(ierr);<br>
<br>
</div>Please don't do this.  Just set the type to PETSCVIEWERVTK, set the file<br>
name to some *.vtu; the VTK viewer will infer output format from<br>
extension and write the VTU XML format with binary-appended data in a<br>
fairly fast and memory-scalable way.  Then<br>
<br>
  VecView(u,viewer);</blockquote><div><br></div><div>I think Jed is wrong here. You need to be using local vectors since they have</div><div>BC and constraints in them. Please explain how your VTU format will magically</div>
<div>put them in.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
>     ierr = DMGetLocalVector(dm, &uLocal);CHKERRQ(ierr);<br>
>     ierr = PetscObjectGetName((PetscObject) u, &name);CHKERRQ(ierr);<br>
>     ierr = PetscObjectSetName((PetscObject) uLocal, name);CHKERRQ(ierr);<br>
>     ierr = DMGlobalToLocalBegin(dm, u, INSERT_VALUES, uLocal);CHKERRQ(ierr);<br>
>     ierr = DMGlobalToLocalEnd(dm, u, INSERT_VALUES, uLocal);CHKERRQ(ierr);<br>
>     ierr = VecView(uLocal, viewer);CHKERRQ(ierr);<br>
>     ierr = DMRestoreLocalVector(dm, &uLocal);CHKERRQ(ierr);<br>
<br>
</div>None of the above, then close the file:<br>
<br>
>     ierr = PetscViewerDestroy(&viewer);CHKERRQ(ierr);<br>
>   }<br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>