On Tue, Nov 22, 2011 at 8:50 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 class="gmail_quote"><div class="im">On Tue, Nov 22, 2011 at 08:30, Andrej Mesaros <span dir="ltr"><<a href="mailto:andrej.mesaros@bc.edu" target="_blank">andrej.mesaros@bc.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>Given that I indeed call MatSetValues exclusively with row indices within the range determined by MatGetOwnershipRange, it should be impossible to generate entries on the wrong process, right? In such a case, could the corruption you mention be somehow due to the way I call other PETSc functions? Or is it at all possible that too small preallocation is making a problem?<br>
</div></blockquote><div><br></div></div><div>Try setting these options and running in debug mode.</div><div><br></div><div>MatSetOption(A,MAT_NO_OFF_PROC_ENTRIES,PETSC_TRUE);</div><div>MatSetOption(A,MAT_NEW_NONZERO_ALLOCATION_ERR,PETSC_TRUE);</div>
</div></blockquote><div><br></div><div>I believe -info gives info on communicated entries after MatAssemble.</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 class="gmail_quote"><div class="im"><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>
<br>
Also, what is the meaning of the memories in the report: "allocated",</div></blockquote><div><br></div></div><div>obtained with PetscMalloc()</div><div><br></div><div><a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Sys/PetscMallocGetCurrentUsage.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Sys/PetscMallocGetCurrentUsage.html</a></div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div> "used by process" </div></blockquote><div><br></div><div>resident set size returned by getrusage(), procfs, or similar</div>
<div><br></div><div><a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Sys/PetscMemoryGetCurrentUsage.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Sys/PetscMemoryGetCurrentUsage.html</a></div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>and "requested"?</div></blockquote><div><br></div><div>Amount you are trying to allocate now.</div>
</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>