<div dir="ltr"><div dir="ltr">On Thu, Apr 4, 2019 at 12:43 PM Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> writes:<br>
<br>
> On Thu, Apr 4, 2019 at 11:52 AM Jed Brown <<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>> wrote:<br>
><br>
>> Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> writes:<br>
>><br>
>> > On Thu, Apr 4, 2019 at 12:26 AM Jed Brown via petsc-dev <<br>
>> > <a href="mailto:petsc-dev@mcs.anl.gov" target="_blank">petsc-dev@mcs.anl.gov</a>> wrote:<br>
>> ><br>
>> >> Patrick Sanan via petsc-dev <<a href="mailto:petsc-dev@mcs.anl.gov" target="_blank">petsc-dev@mcs.anl.gov</a>> writes:<br>
>> >><br>
>> >> > I was going to say something similar to Matt; to me it is worth the<br>
>> 1.5x<br>
>> >> > redundancy in the data if I can just view the output file with "doing<br>
>> >> > anything special". I have also done something similar to what you're<br>
>> >> > describing with save states, and it works fine; weirdly it seems like<br>
>> >> much<br>
>> >> > higher overhead than it is, though.<br>
>> >><br>
>> >> That overhead for 2D fields isn't a big deal for me, but how would you<br>
>> >> identify that three components should be interpreted as vectors?  What<br>
>> >> if they are three chemical concentrations?  Plotting as vectors would<br>
>> >> just be confusing in that case.<br>
>> >><br>
>> ><br>
>> > In PyLith, we annotate solution fields with transformation properties<br>
>> like<br>
>> > this.<br>
>><br>
>> How do you "annotate" in a VTK file?<br>
>><br>
><br>
> My intention was to have them annotated in the library, and then we<br>
> only generate VTK vectors for true vectors. For other fields we make a<br>
> bunch of scalars.<br>
<br>
This could map to your fields/components concept in Plex, but we don't<br>
presently have that distinction in DMDA.<br>
</blockquote></div><br clear="all"><div>Yep, I am aware. I mentioned the PyLith methodology because I think its a possible</div><div>avenue for the future.</div><div><br></div><div>   Matt</div><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>