<div dir="ltr"><div class="gmail_quote"><div dir="ltr">On Mon, Aug 6, 2018 at 5:06 PM Mike Wick <<a href="mailto:michael.wick.1980@gmail.com">michael.wick.1980@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Hi PETSc team:</div><div><br></div><div>I have a question regarding the compatibility between the nest matrix and the vec object. My understanding is that one should create and use nest vectors when the matrix is in the nest form.</div><div><br></div><div>In the tutorial example <br></div><div><a href="http://www.mcs.anl.gov/petsc/petsc-current/src/snes/examples/tutorials/ex70.c.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-current/src/snes/examples/tutorials/ex70.c.html</a></div><div><br></div><div>I noticed that the vectors are created as regular MPI vectors, with the local partition left for PETSC_DECIDE, in lines 171-173. I wander if the underlying structure of the vectors are compatible with the nest matrix. If not, will that create some unnecessary communication work?</div><div><br></div><div>Look forward to your suggestions.</div></div></blockquote><div><br></div><div>Nest matrices work with any vector.</div><div><br></div><div>  Thanks,</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 dir="ltr"><div>Best,</div><div><br></div><div>Mike<br></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><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.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div>