Hello,
This was in fact our own bug and not an error in PETSc. A misplaced call
to DMSetUp was the culprit.
I'm sorry to have bothered you with this.
We are going to update to the latest release, too.
Thanks,
Bastian
On 13.08.20 14:53, Matthew Knepley wrote:
On Thu, Aug 13, 2020 at 8:49 AM Bastian Löhrer
<[email protected] <mailto:[email protected]>>
wrote:
Dear PETSc people,
in PETSc 3.3
call DMView( dm, PETSC_VIEWER_STDOUT_WORLD, ierr)
printed out the ownership ranges like so:
Processor [0] M 32 N 34 P 32 m 1 n 2 p 2 w 1 s 1
X range of indices: 0 32, Y range of indices: 0 17, Z range of
indices: 0 16
Processor [1] M 32 N 34 P 32 m 1 n 2 p 2 w 1 s 1
X range of indices: 0 32, Y range of indices: 17 34, Z range
of indices: 0 16
Processor [2] M 32 N 34 P 32 m 1 n 2 p 2 w 1 s 1
X range of indices: 0 32, Y range of indices: 0 17, Z range of
indices: 16 32
Processor [3] M 32 N 34 P 32 m 1 n 2 p 2 w 1 s 1
X range of indices: 0 32, Y range of indices: 17 34, Z range
of indices: 16 32
In PETSc 3.8.4 (and later?) the same function call only prints out:
DM Object: 4 MPI processes
type: da
Does the feature to print out the ownership ranges still exist?
I am unable to find it.
Certainly the latest release prints what you expect:
knepley/feature-plex-stokes-tutorial
$:/PETSc3/petsc/petsc-dev/src/snes/tutorials$ make ex5
/PETSc3/petsc/apple/bin/mpicc -Wl,-multiply_defined,suppress
-Wl,-multiply_defined -Wl,suppress -Wl,-commons,use_dylibs
-Wl,-search_paths_first -Wl,-no_compact_unwind-Wall -Wwrite-strings
-Wno-strict-aliasing -Wno-unknown-pragmas -fstack-protector
-fno-stack-check -Qunused-arguments -fvisibility=hidden -g3 -Wall
-Wwrite-strings -Wno-strict-aliasing -Wno-unknown-pragmas
-fstack-protector -fno-stack-check -Qunused-arguments
-fvisibility=hidden -g3 -I/PETSc3/petsc/petsc-dev/include
-I/PETSc3/petsc/petsc-dev/arch-master-debug/include -I/opt/X11/include
-I/PETSc3/petsc/apple/include
-I/PETSc3/petsc/petsc-dev/arch-master-debug/include/eigen3ex5.c-Wl,-rpath,/PETSc3/petsc/petsc-dev/arch-master-debug/lib
-L/PETSc3/petsc/petsc-dev/arch-master-debug/lib
-Wl,-rpath,/PETSc3/petsc/petsc-dev/arch-master-debug/lib
-L/PETSc3/petsc/petsc-dev/arch-master-debug/lib
-Wl,-rpath,/opt/X11/lib -L/opt/X11/lib
-Wl,-rpath,/PETSc3/petsc/apple/lib -L/PETSc3/petsc/apple/lib
-Wl,-rpath,/usr/local/lib/gcc/x86_64-apple-darwin19/9.2.0
-L/usr/local/lib/gcc/x86_64-apple-darwin19/9.2.0
-Wl,-rpath,/usr/local/lib -L/usr/local/lib -lpetsc -lcmumps -ldmumps
-lsmumps -lzmumps -lmumps_common -lpord -lscalapack -lumfpack -lklu
-lcholmod -lbtf -lccolamd -lcolamd -lcamd -lamd -lsuitesparseconfig
-lsuperlu_dist -lml -lfftw3_mpi -lfftw3 -lp4est -lsc -llapack -lblas
-legadslite -ltriangle -lX11 -lexodus -lnetcdf -lpnetcdf
-lhdf5hl_fortran -lhdf5_fortran -lhdf5_hl -lhdf5 -lchaco -lparmetis
-lmetis -lz -lctetgen -lc++ -ldl -lmpifort -lmpi -lpmpi -lgfortran
-lquadmath -lm -lc++ -ldl -o ex5
knepley/feature-plex-stokes-tutorial
$:/PETSc3/petsc/petsc-dev/src/snes/tutorials$ ./ex5 -dm_view
DM Object: 1 MPI processes
type: da
Processor [0] M 4 N 4 m 1 n 1 w 1 s 1
X range of indices: 0 4, Y range of indices: 0 4
DM Object: 1 MPI processes
type: da
Processor [0] M 4 N 4 m 1 n 1 w 2 s 1
X range of indices: 0 4, Y range of indices: 0 4
knepley/feature-plex-stokes-tutorial
$:/PETSc3/petsc/petsc-dev/src/snes/tutorials$ $MPIEXEC -np 4 ./ex5
-dm_view
DM Object: 4 MPI processes
type: da
Processor [0] M 4 N 4 m 2 n 2 w 1 s 1
X range of indices: 0 2, Y range of indices: 0 2
Processor [1] M 4 N 4 m 2 n 2 w 1 s 1
X range of indices: 2 4, Y range of indices: 0 2
Processor [2] M 4 N 4 m 2 n 2 w 1 s 1
X range of indices: 0 2, Y range of indices: 2 4
Processor [3] M 4 N 4 m 2 n 2 w 1 s 1
X range of indices: 2 4, Y range of indices: 2 4
DM Object: 4 MPI processes
type: da
Processor [0] M 4 N 4 m 2 n 2 w 2 s 1
X range of indices: 0 2, Y range of indices: 0 2
Processor [1] M 4 N 4 m 2 n 2 w 2 s 1
X range of indices: 2 4, Y range of indices: 0 2
Processor [2] M 4 N 4 m 2 n 2 w 2 s 1
X range of indices: 0 2, Y range of indices: 2 4
Processor [3] M 4 N 4 m 2 n 2 w 2 s 1
X range of indices: 2 4, Y range of indices: 2 4
We can try and go back to debug 3.8.4, but that is a long time ago.
Can you use the latest release?
Thanks,
Matt
Best,
Bastian
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which
their experiments lead.
-- Norbert Wiener
https://www.cse.buffalo.edu/~knepley/
<http://www.cse.buffalo.edu/~knepley/>
