On Wed, Apr 4, 2012 at 4:47 PM, Andrew Spott <span dir="ltr"><<a href="mailto:andrew.spott@gmail.com">andrew.spott@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I'm interested in moving a vector to a single processor (so I can output it in some custom way), but I can't find a simple way to do that.<br>
<br>
Is the easiest way to output the parallel vector to a file and then read it in to a single processor vector?<br></blockquote><div><br></div><div><a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Vec/VecScatterCreateToZero.html">http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Vec/VecScatterCreateToZero.html</a></div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Also, an unrelated question: If I have a vector that I'm using as the result of a matrix vector product, but it is initially empty, do I need to initialize it in any way? or can I just create it and run with it:<br>
<br>
… create x and A ...<br>
VecCreate(world, &b);<br></blockquote><div><br></div><div>This create is wrong, that Vec will be overwritten.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
MatGetVecs(A,&x,&b);<br></blockquote><div><br></div><div>You can jsut use these vectors.</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">
… do stuff with x and A ...<br>
<br>
MatMult(A,x,b);<br>
<br>
Thanks for the help,<br>
<br>
-Andrew</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<br>