<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Dec 19, 2013 at 4:44 PM, James A Charles <span dir="ltr"><<a href="mailto:charlesj@purdue.edu" target="_blank">charlesj@purdue.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Hello,<br>
<br>
I want to redistribute a matrix across MPI processes (block tridiagonal of around size 20,000) for use with Arpack so that I can solve for eignepairs in parallel. Is this possible and if so what is the best way to do this?<br>
</blockquote><div><br></div><div>The easiet way is probably <a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatGetSubMatrix.html">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatGetSubMatrix.html</a>, but there is also</div>
<div><a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatPermute.html">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatPermute.html</a>.</div><div><br></div><div>    Matt</div><div>
 </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
Thanks,<br>
James Charles<br>
</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
</div></div>