[petsc-users] Segmentation Fault in ts tutorial ex23

[Jed Brown]

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>