<div dir="ltr"><div><div><div>Jed thank you for the response.<br><br></div>1) In my case, I want to output a transient diffusion problem (with DMPlex). I see that in ex11 there's a TSMonitorSet(). Is this (and its calling function) be the place I have something along the lines that you have suggested? That is, if TSMonitorSet() calls MonitorVTK(), then MonitorVTK() would have PetscViewerVTKOpen(...). Am I along the right track?<br><br></div>2) What I meant by parallel was if each rank would output its local part of the solution, kind of like the way it's done with FEniCS. That way, if I had an extremely large solution set (at least over 10 million DOFS) then I wouldn't be stuck with one giant ASCII file. Unless that's what binary formats like HDF5 are for?<br><br></div>Thanks,<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Mar 30, 2015 at 10:54 AM, 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"><span class="">Justin Chang <<a href="mailto:jchang27@uh.edu">jchang27@uh.edu</a>> writes:<br>
<br>
> Hi all,<br>
><br>
> I have a couple questions:<br>
><br>
> 1) Is there any documentation somewhere that talks specifically on how to<br>
> output using the *.vtu format? Specifically for transient problems. When I<br>
> look at TS ex11 it seems to me the user has to create several additional<br>
> routines, but I was wondering if there are any simpler<br>
> examples/instructions.<br>
<br>
</span>    ierr = PetscViewerVTKOpen(PETSC_COMM_WORLD,"your-file.vtu",FILE_MODE_WRITE,&viewer);CHKERRQ(ierr);<br>
    ierr = VecView(x,viewer);CHKERRQ(ierr);<br>
    ierr = PetscViewerDestroy(&viewer);CHKERRQ(ierr);<br>
<br>
... assuming x is a vector obtained from a DM capable of writing *.vtu.<br>
Note that the vtu format is not designed for time-series, so you'll need<br>
to write a separate file per time step.<br>
<span class=""><br>
> 2) Can vtu files be outputted in parallel?<br>
<br>
</span>Yes, we write one file in parallel (not one file per rank which is a<br>
scalability disaster).<br>
<br>
</blockquote></div><br></div><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div><div>Justin Chang<br></div>PhD Candidate, Civil Engineering - Computational Sciences<br></div>University of Houston, Department of Civil and Environmental Engineering<br></div>Houston, TX 77004<br>(512) 963-3262<br></div></div>