On Tue, Sep 20, 2011 at 2:44 AM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</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;">
<div>I'm so tired of calling a different preallocation variant for each matrix format.</div><div><br></div>MatMPIBAIJSetPreallocation() has enough information to preallocate {Seq,MPI}{AIJ,BAIJ}. Should we make this one call enough?</blockquote>
<div><br></div><div>I disagree with this. I think the correct reduction in complexity is for the DM to call all relevant versions.</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>Downsides:</div><div><br></div><div>If the user wants to allocate some special pattern within blocks, this would fill it all in, so they would have to check the types and call only the version that they really want.</div>

<div><br></div><div>Old code calls all these functions which could mean a lot of reallocation, they can change the code to only call this one function.</div>
</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<br>