Change STAR to BOX as in the patch below. I'll push this to petsc-dev later.
diff --git a/src/ts/examples/tutorials/ex23.c b/src/ts/examples/tutorials/ex23.c --- a/src/ts/examples/tutorials/ex23.c +++ b/src/ts/examples/tutorials/ex23.c @@ -52,7 +52,7 @@ /* Get physics and time parameters */ ierr = GetParams(&user);CHKERRQ(ierr); /* Create a 2D DA with dof = 2 */ - ierr = DMDACreate2d(PETSC_COMM_WORLD,DMDA_BOUNDARY_NONE,DMDA_BOUNDARY_NONE,DMDA_STENCIL_STAR,-4,-4,PETSC_DECIDE,PETSC_DECIDE,2,1,PETSC_NULL,PETSC_NULL,&user.da);CHKERRQ(ierr); + ierr = DMDACreate2d(PETSC_COMM_WORLD,DMDA_BOUNDARY_NONE,DMDA_BOUNDARY_NONE,DMDA_STENCIL_BOX,-4,-4,PETSC_DECIDE,PETSC_DECIDE,2,1,PETSC_NULL,PETSC_NULL,&user.da);CHKERRQ(ierr); /* Set Element type (triangular) */ ierr = DMDASetElementType(user.da,DMDA_ELEMENT_P1);CHKERRQ(ierr); On Mon, Nov 26, 2012 at 10:28 AM, Milad Fatenejad <icksa1 at gmail.com> wrote: > Hello: > > I was trying to run example ts/examples/tutorial/ex23 and I seem to get an > error similar to the one I've shown below no matter what options I use. For > the error below, I simply ran ex23 with no arguments. I've tried running > with the released petsc version and the development version from bitbucket. > Any suggestions? > > Thanks > Milad > > [0]PETSC ERROR: --------------------- Error Message > ------------------------------------ > [0]PETSC ERROR: Argument out of range! > [0]PETSC ERROR: New nonzero at (2,9) caused a malloc! > [0]PETSC ERROR: > ------------------------------------------------------------------------ > [0]PETSC ERROR: Petsc Development HG revision: > f86d5620e977ae10ba0f45b74b5f06746e16ab61 HG Date: Mon Nov 19 18:09:55 2012 > +0100 > [0]PETSC ERROR: See docs/changes/index.html for recent updates. > [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. > [0]PETSC ERROR: See docs/index.html for manual pages. > [0]PETSC ERROR: > ------------------------------------------------------------------------ > [0]PETSC ERROR: ./ex23 on a arch-linu named animal5 by milad Mon Nov 26 > 10:12:42 2012 > [0]PETSC ERROR: Libraries linked from > /home/milad/petsc/petsc-dev-bb/arch-linux2-c-debug/lib > [0]PETSC ERROR: Configure run at Mon Nov 19 13:22:07 2012 > [0]PETSC ERROR: Configure options --with-moab=1 > --with-moab-dir=/home/milad/Downloads/moab/moab-opt-parallel/install > --with-hdf5=1 --with-hdf5-dir=/usr/local/hdf5-1.8.8-parallel > --with-netcdf=1 --with-netcdf-dir=/usr/local/netcdf-4.1.3-serial > --with-clanguage=C++ > [0]PETSC ERROR: > ------------------------------------------------------------------------ > [0]PETSC ERROR: MatSetValues_SeqAIJ() line 353 in > /home/milad/petsc/petsc-dev-bb/src/mat/impls/aij/seq/aij.c > [0]PETSC ERROR: MatSetValues() line 1080 in > /home/milad/petsc/petsc-dev-bb/src/mat/interface/matrix.c > [0]PETSC ERROR: MatSetValuesLocal() line 1941 in > /home/milad/petsc/petsc-dev-bb/src/mat/interface/matrix.c > [0]PETSC ERROR: SetUpMatrices() line 396 in > src/ts/examples/tutorials/ex23.c > [0]PETSC ERROR: main() line 76 in src/ts/examples/tutorials/ex23.c > application called MPI_Abort(MPI_COMM_WORLD, 63) - process 0 > > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20121127/3fd4568e/attachment.html>
