On Wed, Mar 5, 2014 at 2:21 PM, Jed Brown <[email protected]> wrote: > Xiangdong <[email protected]> writes: > > > I got the following error messages on each processor when use mpiio: > > > > [5]PETSC ERROR: DMDAArrayMPIIO() line 532 in > > MyLocal/petsc/petsc-3.4.3/src/dm/impls/da/gr2.c > > Looks like this line is failing. Please either use a debugger or create > a representative test case so we can debug. > > ierr = > MPI_Type_create_subarray(dd->dim+1,gsizes,lsizes,lstarts,MPI_ORDER_FORTRAN,MPIU_SCALAR,&view);CHKERRQ(ierr); > > > > Does this work? > > petsc/src/ksp/ksp/examples/tutorials$ mpiexec -n 4 ./ex45 -da_refine 3 > -ksp_monitor -pc_type mg -ksp_view_solution binary:foo -viewer_binary_mpiio >
This works fine. However, when I add the option -da_grid_x 8, the same error messages pop up. It is strange that -da_grid_x 5 works fine. Moreover, if we set M=N=P=8, it also works fine. Thank you. Xiangdong > > > [5]PETSC ERROR: VecView_MPI_DA() line 605 in > > MyLocal/petsc/petsc-3.4.3/src/dm/impls/da/gr2.c > > [5]PETSC ERROR: VecView() line 717 in > > MyLocal/petsc/petsc-3.4.3/src/vec/vec/interface/vector.c > > > > The program does not crash, but skip writing the binary output. > > That's because you weren't checking error codes. >
