I'll check it. There are other routines need this fix. On my 'todo' list :-(<div><br></div><div>BTW, there are progress along hierarchical ksp work. I like to show it to you</div><div>and get your suggestions. Call me next week when you come to the lab.</div>
<div><br></div><div>Hong<br><br><div class="gmail_quote">On Fri, Mar 2, 2012 at 4:54 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
  Hong,<br>
<br>
  MatMerge_Seq() is not doing preallocation nor setting the flag that allows improper preallocation<br>
<br>
-*- mode: compilation; default-directory: "/Users/barrysmith/Src/petsc-dev/src/mat/examples/tests/" -*-<br>
Compilation started at Fri Mar  2 16:52:06<br>
<br>
make PETSC_DIR=/Users/barrysmith/Src/petsc-dev PETSC_ARCH=arch-gnu runex95<br>
[0]PETSC ERROR: --------------------- Error Message ------------------------------------<br>
[0]PETSC ERROR: Argument out of range!<br>
[0]PETSC ERROR: New nonzero at (0,0) caused a malloc!<br>
[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: Petsc Development HG revision: d35fa29f81218e6e5376faf87f3925c16b47cb76  HG Date: Thu Mar 01 19:13:43 2012 -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>
[0]PETSC ERROR: ./ex95 on a arch-gnu named <a href="http://anlextwls003-010.wl.anl-external.org" target="_blank">anlextwls003-010.wl.anl-external.org</a> by barrysmith Fri Mar  2 16:52:06 2012<br>
[0]PETSC ERROR: Libraries linked from /Users/barrysmith/Src/petsc-dev/arch-gnu/lib<br>
[0]PETSC ERROR: Configure run at Thu Mar  1 20:56:05 2012<br>
[0]PETSC ERROR: Configure options --download-blacs --download-fftw --download-hypre --download-metis --download-mpich --download-mumps --download-parmetis --download-ptscotch --download-scalapack --download-superlu --download-superlu_dist --download-yaml --with-shared-libraries PETSC_ARCH=arch-gnu<br>

[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: MatSetValues_MPIAIJ() line 506 in /Users/barrysmith/Src/petsc-dev/src/mat/impls/aij/mpi/mpiaij.c<br>
[0]PETSC ERROR: [1]PETSC ERROR: --------------------- Error Message ------------------------------------<br>
[1]PETSC ERROR: Argument out of range!<br>
[1]PETSC ERROR: New nonzero at (0,0) caused a malloc!<br>
[1]PETSC ERROR: ------------------------------------------------------------------------<br>
[1]PETSC ERROR: Petsc Development HG revision: d35fa29f81218e6e5376faf87f3925c16b47cb76  HG Date: Thu Mar 01 19:13:43 2012 -0600<br>
[1]PETSC ERROR: See docs/changes/index.html for recent updates.<br>
[1]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>
[1]PETSC ERROR: See docs/index.html for manual pages.<br>
[1]PETSC ERROR: ------------------------------------------------------------------------<br>
[1]PETSC ERROR: ./ex95 on a arch-gnu named <a href="http://anlextwls003-010.wl.anl-external.org" target="_blank">anlextwls003-010.wl.anl-external.org</a> by barrysmith Fri Mar  2 16:52:06 2012<br>
[1]PETSC ERROR: Libraries linked from /Users/barrysmith/Src/petsc-dev/arch-gnu/lib<br>
[1]PETSC ERROR: Configure run at Thu Mar  1 20:56:05 2012<br>
[1]PETSC ERROR: Configure options --download-blacs --download-fftw --download-hypre --download-metis --download-mpich --download-mumps --download-parmetis --download-ptscotch --download-scalapack --download-superlu --download-superlu_dist --download-yaml --with-shared-libraries PETSC_ARCH=arch-gnu<br>

[1]PETSC ERROR: ------------------------------------------------------------------------<br>
[1]PETSC ERROR: MatSetValues_MPIAIJ() line 506 in /Users/barrysmith/Src/petsc-dev/src/mat/impls/aij/mpi/mpiaij.c<br>
[1]PETSC ERROR: MatSetValues() line 1141 in /Users/barrysmith/Src/petsc-dev/src/mat/interface/matrix.c<br>
[1]PETSC ERROR: MatMerge_SeqsToMPINumeric() line 4548 in /Users/barrysmith/Src/petsc-dev/src/mat/impls/aij/mpi/mpiaij.c<br>
[1]PETSC ERROR: MatMerge_SeqsToMPI() line 4848 in /Users/barrysmith/Src/petsc-dev/src/mat/impls/aij/mpi/mpiaij.c<br>
[1]PETSC ERROR: main() line 41 in src/mat/examples/tests/ex95.c<br>
application called MPI_Abort(MPI_COMM_WORLD, 63) - process 1<br>
[cli_1]: aborting job:<br>
application called MPI_Abort(MPI_COMM_WORLD, 63) - process 1<br>
[2]PETSC ERROR: --------------------- Error Message ------------------------------------<br>
[2]PETSC ERROR: Argument out of range!<br>
[2]PETSC ERROR: New nonzero at (0,0) caused a malloc!<br>
[2]PETSC ERROR: ------------------------------------------------------------------------<br>
[2]PETSC ERROR: Petsc Development HG revision: d35fa29f81218e6e5376faf87f3925c16b47cb76  HG Date: Thu Mar 01 19:13:43 2012 -0600<br>
[2]PETSC ERROR: See docs/changes/index.html for recent updates.<br>
[2]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>
[2]PETSC ERROR: See docs/index.html for manual pages.<br>
[2]PETSC ERROR: ------------------------------------------------------------------------<br>
[2]PETSC ERROR: ./ex95 on a arch-gnu named <a href="http://anlextwls003-010.wl.anl-external.org" target="_blank">anlextwls003-010.wl.anl-external.org</a> by barrysmith Fri Mar  2 16:52:06 2012<br>
[2]PETSC ERROR: Libraries linked from /Users/barrysmith/Src/petsc-dev/arch-gnu/lib<br>
[2]PETSC ERROR: Configure run at Thu Mar  1 20:56:05 2012<br>
[2]PETSC ERROR: Configure options --download-blacs --download-fftw --download-hypre --download-metis --download-mpich --download-mumps --download-parmetis --download-ptscotch --download-scalapack --download-superlu --download-superlu_dist --download-yaml --with-shared-libraries PETSC_ARCH=arch-gnu<br>

[2]PETSC ERROR: ------------------------------------------------------------------------<br>
[2]PETSC ERROR: MatSetValues_MPIAIJ() line 506 in /Users/barrysmith/Src/petsc-dev/src/mat/impls/aij/mpi/mpiaij.c<br>
[2]PETSC ERROR: MatSetValues() line 1141 in /Users/barrysmith/Src/petsc-dev/src/mat/interface/matrix.c<br>
[2]PETSC ERROR: MatMerge_SeqsToMPINumeric() line 4548 in /Users/barrysmith/Src/petsc-dev/src/mat/impls/aij/mpi/mpiaij.c<br>
[2]PETSC ERROR: MatMerge_Seq<br>
Compilation finished at Fri Mar  2 16:52:06<br>
<br>
<br>
</blockquote></div><br></div>