On Wed, Jun 22, 2011 at 8:55 AM, John Fettig <span dir="ltr"><<a href="mailto:john.fettig@gmail.com">john.fettig@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 having trouble understanding how the size of off-diagonal matrix<br>
is determined in the MPIAIJ format. For example, if I run on 2<br>
processors the global size is 4177x4177, but the local matrices are<br>
<br>
Rank 1:<br>
2089x2089 (diagonal block)<br>
2089x259 (off diagonal block)<br>
<br>
Rank2:<br>
2088x274 (off diagonal block)<br>
2088x2088 (diagonal block)<br>
<br>
Why isn't the off-diagonal block sized for all of the off-diagonal columns?</blockquote><div><br></div><div>It is compressed, so it only has columns that have a nonzero.</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;"><font color="#888888"><br>
John<br>
</font></blockquote></div><br><br clear="all"><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>