[petsc-dev] Including petsc.h breaks user code

[Matthew Knepley]

On Sun, Sep 14, 2014 at 4:16 PM, Pierre Jolivet <jolivet at ann.jussieu.fr>
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20140914/eec3938f/attachment.html>