<div class="gmail_quote">On Thu, Jun 23, 2011 at 14:53, John Fettig <span dir="ltr"><<a href="mailto:john.fettig@gmail.com">john.fettig@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":42">I was wrong on both counts, and I still don't understand how the<br>
column numbering for the off-diagonal block is calculated. It appears<br>
you are doing something more than just compressing the columns, you<br>
are also shifting the indices if an entire column is zero?</div></blockquote><div><br></div><div>Yeah, that's what we mean by "compressing the columns".</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":42"> Is there<br>
anyplace in the documentation that would help me understand this?<br>
Alternatively, where would I look in mpiaij.c?</div></blockquote></div><br><div>The MPIAIJ format works with blocks that are in SeqAIJ format. In most cases, it's not necessary to be familiar with the internal storage format. Are you just curious how it's implemented, or is there something you are trying to do that depends on the storage format?</div>