<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>Dmitry,</div><div><br></div><div>finally, I am able to make a simple program and reproduce the error. Since I had moved to</div><div>petsc-dev (before I used petsc-3.2), the error is a little different.</div><div><br></div><div>Download the program attached and run the run.sh. You will see the output from PCASM.</div><div>There is a submatrix (submat[1] on rank 0) which has a row like</div><div><br></div><div> row 1: (0, -1) (1, 0 + 0.4 i) (4, -1) (2, -1) </div><div><br></div><div>Note that the column 2 is placed in the last and the column numbers are not in the normal order.</div><div>Following that, we can also see the error message like</div><div><br></div><div> [0]PETSC ERROR: New nonzero at (1,2) caused a malloc!</div><div><br></div><div>You can also try with GASM, use</div><div><br></div><div> mpirun -np 2 ./gasm_test -dd_type gasm -n 4 -px 1 -py 2 -sx 2 -sy 1 </div><div><br></div><div>which partitions the domain to 1 * 2 and mapped to the two processors, and for each processor</div><div>we further partition its local domain to 2 * 1 subdomains. </div><div><br></div><div>Then, you will see that the row 1 of submat[1] on rank 1 becomes</div><div><br></div><div> row 1: (0, -1) (1, 0 + 0.4 i) (3, -1) (6, -1) </div><div><br></div><div>which is totally wrong because from the IS for this overlapping subdomain and the big matrix A,</div><div>row 1 should have non-zeros like the output from PCASM.</div><div><br></div><div>I guess the problem is due to that we set PCASMSetSortIndices to FALSE and something goes wrong</div><div>in asm.c and gasm.c.</div><div><br></div><div>Thanks!</div><div>Hui</div><div><br></div><div><br></div><div></div></body></html>