Awesome! That is fantastic. :-)<div><br></div><div>Thanks Barry for your prompt response. </div><div>Mohamad</div><div><br><br><div class="gmail_quote">On Mon, Oct 10, 2011 at 5:45 PM, Barry Smith <span dir="ltr">&lt;<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>&gt;</span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div class="im"><br>
On Oct 10, 2011, at 7:41 PM, Mohamad M. Nasr-Azadani wrote:<br>
<br>
&gt; Thanks Barry,<br>
&gt;<br>
&gt; I am still not 100% sure if I can do this.<br>
&gt; Say I have save the global vector obtained by a DA (3D) that is shared amongst 16 processors. :.<br>
&gt; Can I load that data into a vector obtained from a DA (3D, same size obviously) that is shared on 1 processor?<br>
<br>
</div>  ABSOLUTELY. Or a DA on 2 processes etc.<br>
<font color="#888888"><br>
   Barry<br>
</font><div><div></div><div class="h5"><br>
&gt;<br>
&gt; Thanks,<br>
&gt; Best,<br>
&gt; Mohamad<br>
&gt;<br>
&gt;<br>
&gt; On Mon, Oct 10, 2011 at 5:36 PM, Barry Smith &lt;<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>&gt; wrote:<br>
&gt;<br>
&gt; On Oct 10, 2011, at 7:31 PM, Mohamad M. Nasr-Azadani wrote:<br>
&gt;<br>
&gt; &gt; Hi<br>
&gt; &gt; I was wondering if it would be possible if I write a global vector (associated with a certain 3D distributed array)  to file<br>
&gt; &gt; via:<br>
&gt; &gt;<br>
&gt; &gt;             ierr = PetscViewerBinaryOpen(PCW1, filename, FILE_MODE_WRITE, &amp;writer);<br>
&gt; &gt;             ierr = VecView(vec_data, writer);<br>
&gt; &gt;             ierr = PetscViewerDestroy(writer);<br>
&gt; &gt;<br>
&gt; &gt; And then load the data into a global vector which is not created using the same parallel layout?<br>
&gt; &gt; A simple example for this case would be to write the runtime data (parallel vector) to file and then just load the saved vector to do some simple SERIAL post processing.<br>
&gt;<br>
&gt;  YES.<br>
&gt;<br>
&gt;   The vector is saved to the file in the &quot;natural ordering&quot; that is starting with the logically 0,0,0 coordinate then increasing through the x axis, then the y axis then the z axis. To load back in in parallel you need to pass to VecLoad() a vector obtained with the appropriate DMCreateGlobalVector().<br>


&gt;<br>
&gt;<br>
&gt;  Barry<br>
&gt;<br>
&gt; &gt;<br>
&gt; &gt; Thanks in advance,<br>
&gt; &gt; Best<br>
&gt; &gt; Mohamad<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt;<br>
&gt;<br>
<br>
</div></div></blockquote></div><br></div>