On Sun, Sep 14, 2014 at 4:16 PM, Pierre Jolivet <[email protected]> wrote:
> Hello, > Could you please explain to me why the following example is not working > properly when <petsc.h> (from master, with OpenMPI 1.8.1) is included ? > > $ mpicxx in-place.cpp -I$PETSC_DIR/include > -I$PETSC_DIR/$PETSC_ARCH/include -L$PETSC_DIR/$PETSC_ARCH/lib -lpetsc > $ mpirun -np 2 ./a.out > Done with the scatter ! > 0 0 0 0 (this line should be filled with 0) > 1 1 1 1 (this line should be filled with 1) > Done with the gather ! > > $ mpicxx in-place.cpp -I$PETSC_DIR/include > -I$PETSC_DIR/$PETSC_ARCH/include -L$PETSC_DIR/$PETSC_ARCH/lib -lpetsc > -DPETSC_BUG > $ mpirun -np 2 ./a.out > [:3367] *** An error occurred in MPI_Type_size > [:3367] *** reported by process [4819779585,140733193388032] > [:3367] *** on communicator MPI_COMM_WORLD > [:3367] *** MPI_ERR_TYPE: invalid datatype > [:3367] *** MPI_ERRORS_ARE_FATAL (processes in this communicator will now > abort, > [:3367] *** and potentially your MPI job) > I just built this with MPICH and it runs fine: master:/PETSc3/petsc/petsc-pylith$ /PETSc3/petsc/petsc-pylith/arch-pylith-cxx-debug/bin/mpiexec -host localhost -n 2 /PETSc3/petsc/petsc-pylith/arch-pylith-cxx-debug/lib/in-place-obj/in-place Done with the scatter ! 0 0 0 0 (this line should be filled with 0) 1 1 1 1 (this line should be filled with 1) Done with the gather ! Will valgrind. Matt > Thank you for looking, > Pierre > > -- 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
