On Mon, 16 Jul 2012, Blaise Bourdin wrote:

On Jul 16, 2012, at 10:11 AM, Jed Brown wrote:

...

On Mon, Jul 16, 2012 at 10:09 AM, Shri <[email protected]> wrote: Yes, i think so too. If no one objects then i can do that. Btw, MPI_COMM_SELF and MPI_COMM_WORLD are same in mpiuni/mpi.h #define MPI_COMM_WORLD 1 #define MPI_COMM_SELF MPI_COMM_WORLD

while the fortran version declares these communicators as different. mpiuni/mpif.h integer MPI_COMM_WORLD parameter (MPI_COMM_WORLD = 1) integer MPI_COMM_SELF parameter (MPI_COMM_SELF = 2)

WAT? The Fortran version did something right?

That's because it is tested by people who do stuff with it ;)

Actually for np=1 MPI_COMM_SELF==MPI_COMM_WORLD. mpi.h was changed to reflect this - but mpif.h was not. Satish ---------- balay@petsc:~/Petsc-historic/src/sys/src/mpiuni $ rcsdiff -r1.52 -r1.53 mpi.h =================================================================== RCS file: RCS/mpi.h,v retrieving revision 1.52 retrieving revision 1.53 diff -r1.52 -r1.53 1c1 < /* $Id: mpi.h,v 1.52 1997/07/28 21:41:02 balay Exp $ */ ---

/* $Id: mpi.h,v 1.53 1997/08/21 01:06:19 balay Exp $ */ 22c22 < #define MPI_COMM_SELF 2

#define MPI_COMM_SELF MPI_COMM_WORLD 427c427 < (MPIUNI_TMP = (void *) (comm), \

(*(newcomm) = (comm), \

429d428 < MPIUNI_TMP = (void *) (newcomm), \ balay@petsc:~/Petsc-historic/src/sys/src/mpiuni $