<div class="gmail_quote">On Thu, Jul 12, 2012 at 4:14 PM, Michael R Nowak <span dir="ltr"><<a href="mailto:nowakmr@umich.edu" target="_blank">nowakmr@umich.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":19i">I'm setting up a matrix with Slepc, and manually setting the preallocation<br>
for nonzero entries row-by-row using the MatMPISBAIJSetPreallocation<br>
function as follows:<br>
<br>
MatMPISBAIJSetPreallocation(LaplaceOperator,<br>
                            1,<br>
                            0,<br>
                            nonzeroEntries,<br>
                            0,<br>
                            PETSC_NULL);<br>
<br>
The problem with this is that the number of nonzero entries per row that I<br>
give Slepc is not the same as the number of nonzero entries that Slepc<br>
allocates space for, slepc always seems to allocate N more nonzeroes than<br>
needed,</div></blockquote><div><br></div><div>What do you mean by SLEPc allocating nonzeros?</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div id=":19i">
 with N being the number of rows in the matrix. This is a problem<br>
for me, since I would like to reduce the number of unneeded nonzero entries<br>
to exactly zero, in order to save memory. The number of nonzero entries<br>
that I give in nonzeroEntries is exactly what is used, so why are extra<br>
nonzeroes being allocated? Any help with this would be much appreciated.</div></blockquote></div><br><div>What evidence do you have that this is happening?</div><div><br></div>