<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Dec 19, 2013 at 5:20 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">My original matrix is serial. I then want to distribute this across MPI processes for use with Arpack. MatGetSubMatrix wouldn't change communicators right?<br>
</blockquote><div><br></div><div>When you create the matrix, use PETSC_COMM_WORLD and 0 sizes on the other procs.</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">

----- Original Message -----<br>
From: "Matthew Knepley" <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>><br>
To: "James A Charles" <<a href="mailto:charlesj@purdue.edu">charlesj@purdue.edu</a>><br>
Cc: <a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a><br>
Sent: Thursday, December 19, 2013 5:58:25 PM<br>
Subject: Re: [petsc-users] Petsc Matrix Redistribution<br>
<br>
<br>
<br>
<br>
On Thu, Dec 19, 2013 at 4:44 PM, James A Charles < <a href="mailto:charlesj@purdue.edu">charlesj@purdue.edu</a> > wrote:<br>
<br>
<br>
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>

<br>
<br>
<br>
The easiet way is probably <a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatGetSubMatrix.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatGetSubMatrix.html</a> , but there is also<br>

<a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatPermute.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatPermute.html</a> .<br>
<br>
<br>
Matt<br>
<br>
<br>
Thanks,<br>
James Charles<br>
<span class="HOEnZb"><font color="#888888"><br>
<br>
<br>
<br>
--<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>
</font></span></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>