<div dir="ltr"><div style>Hmm, this debugging check was added by Mark here.</div><div><br></div><a href="https://bitbucket.org/petsc/petsc-dev/commits/fb7e4343253b73d3fb40857991a8c6d5d0c7c56b#Lsrc/mat/impls/aij/mpi/mmaij.cT126">https://bitbucket.org/petsc/petsc-dev/commits/fb7e4343253b73d3fb40857991a8c6d5d0c7c56b#Lsrc/mat/impls/aij/mpi/mmaij.cT126</a><div>
<br></div><div>Barry touched the code a bit earlier to add this case:</div><div><br></div><div><a href="https://bitbucket.org/petsc/petsc-dev/commits/c997baff1558fd165412be034e812bb69f90972a">https://bitbucket.org/petsc/petsc-dev/commits/c997baff1558fd165412be034e812bb69f90972a</a></div>
<div><div class="gmail_extra"><br></div><div class="gmail_extra" style>Prior to that, Mark had reversed the meaning of that bit of code:</div><div class="gmail_extra" style><br></div><div class="gmail_extra" style><a href="https://bitbucket.org/petsc/petsc-dev/commits/a85625d104df12863d7d1735f9fdbe9b281ef7d2#chg-src/mat/impls/aij/mpi/mmaij.c">https://bitbucket.org/petsc/petsc-dev/commits/a85625d104df12863d7d1735f9fdbe9b281ef7d2#chg-src/mat/impls/aij/mpi/mmaij.c</a><br>
</div><div class="gmail_extra"><br></div></div><div><div class="gmail_extra"><div class="gmail_extra" style>Even in Barry's original code, I don't understand why we can expect this check to be producing a collective result. Surely there can be some ranks that see everything evenly divide and others that do not.</div>
<div class="gmail_extra" style><br></div><div class="gmail_extra" style>It seems to me that we should either (a) explicitly disallow matrices with a block size to be built when the blocks are not filled or (b) when the user requests blocks, pad blocks so that we can always use an ISBlock.</div>
</div></div><div><div class="gmail_extra"><div><br></div><br><div class="gmail_quote">On Mon, Feb 18, 2013 at 2:35 PM, Garth N. Wells <span dir="ltr"><<a href="mailto:gnw20@cam.ac.uk" target="_blank">gnw20@cam.ac.uk</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">I have some code using ML as a preconditioner, and  it works fines<br>

with PETSc 3.3 but crashes with PETSc dev. I'm not setting any of the<br>
ML parameters.  Hopefully someone will recognise the error (below) and<br>
save me having to bisect the changesets, which is hard because changes<br>
in PETSc dev make it hard for me to compile against earlier PETSc dev<br>
revisions.<br>
<br>
Garth<br>
<br>
<br>
[0]PETSC ERROR: --------------------- Error Message<br>
------------------------------------<br>
[0]PETSC ERROR: Petsc has generated inconsistent data!<br>
[0]PETSC ERROR: Use of blocked not consistant (I am usning blocked)!<br>
[0]PETSC ERROR:<br>
------------------------------------------------------------------------<br>
[0]PETSC ERROR: Petsc Development HG revision:<br>
96a09d6a701823d69aede247ed949c1dd4e7c8ff  HG Date: Mon Feb 18 11:08:29<br>
2013 -0600<br>
[0]PETSC ERROR: See docs/changes/index.html for recent updates.<br>
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>
[0]PETSC ERROR: See docs/index.html for manual pages.<br>
[0]PETSC ERROR:<br>
------------------------------------------------------------------------<br>
[0]PETSC ERROR: Unknown Name on a arch-linux2-c-debug named .......<br>
Mon Feb 18 19:25:17 2013<br>
[0]PETSC ERROR: Libraries linked from /home/garth/local/gcc/petsc-dev/lib<br>
[0]PETSC ERROR:<br>
[3]PETSC ERROR: See docs/index.html for manual pages.<br>
</blockquote></div><br></div></div></div>