Open MPI logo

Open MPI Development Mailing List Archives

  |   Home   |   Support   |   FAQ   |   all Development mailing list

From: Daniel Spångberg (Daniel.Spangberg_at_[hidden])
Date: 2007-06-12 17:15:26

Dear openmpi developers,

I have very recently joined this list to be able to submit this report, as
recommended in the FAQ, why I have not had time to familiarize myself with
the etiquette of this list, for which I apologize in advance. If you
require any additional information I am more than happy to provide you
with such information.

One of the users of our computer cluster have reported a problem with his
simulation code. We have been able to trace this down to a problem with
the MPI struct when absolute addressing is used for the members in the
struct, necessating the use of MPI_BOTTOM. The problem occurs with
MPI_BCAST, but works fine with MPI_SEND and MPI_RECV. I digged down into
the openmpi code and found out that the problem occurs only in fortran,
which seems to be because when MPI_BCAST is called, the OMPI_ADDR macro in
ompi/mpi/f77/constants.h is never evaluated which thus never turns the
fortran MPI_BOTTOM into a C MPI_BOTTOM. When MPI_SEND and MPI_RECV are
used the OMPI_ADDR macro is evaluated and no problem occurs.

An example of a problematic code (tested on two processes):

       program testme
       implicit none
       include 'mpif.h'
       integer ierr,rank,size
       integer btype(1),blen(1),param_batch
       real param
       integer(kind=mpi_address_kind) :: disp(1)
       integer status(MPI_STATUS_SIZE)


       call MPI_INIT(ierr)
       call MPI_GET_ADDRESS(param,disp(1),ierr)
       if ( write(*,*) 'MPI_GET_ADDRESS FAILED'
       call MPI_COMM_RANK(MPI_COMM_WORLD,rank,ierr)
       call MPI_COMM_SIZE(MPI_COMM_WORLD,size,ierr)
       call MPI_TYPE_CREATE_STRUCT(1,blen,disp,btype,param_batch,ierr)
       if ( then
         write(*,*) 'MPI_TYPE_CREATE_STRUCT FAILED'
       call MPI_TYPE_COMMIT(param_batch,ierr)
       if ( write(*,*) 'MPI_TYPE_COMMIT FAILED'
       if (rank.eq.0) param=100.
c if (rank.eq.0) then
c x param_batch,1,0,MPI_COMM_WORLD,ierr)
c else
c x param_batch,0,0,MPI_COMM_WORLD,status,ierr)
c endif
       call MPI_BCAST(MPI_BOTTOM,1,param_batch,0,MPI_COMM_WORLD,ierr)
       if ( write(*,*) 'MPI_BCAST FAILED'
       write(*,*) 'Rank:',rank,'Size:',size,'Param=',param
       call MPI_Finalize(ierr)

       end program testme

mpirun -np 2 testme_submitted_to_devel
[auchentoshan:21021] *** Process received signal ***
[auchentoshan:21021] Signal: Segmentation fault (11)
[auchentoshan:21021] Signal code: Address not mapped (1)
[auchentoshan:21021] Failing at address: 0x7fc04ffd7c
[auchentoshan:21021] [ 0] /lib64/tls/ [0x3d9660c430]
[auchentoshan:21021] [ 1] /lib64/tls/ [0x3d95571ec0]
[auchentoshan:21021] [ 2]
[auchentoshan:21021] [ 3]
[auchentoshan:21021] [ 4]
[auchentoshan:21021] [ 5]
[auchentoshan:21021] [ 6]
[auchentoshan:21021] [ 7]
/opt/openmpi-1.2.1-gcc4/lib/ [0x2a957d62ac]
[auchentoshan:21021] [ 8]
[auchentoshan:21021] [ 9] testme_submitted_to_devel(MAIN__+0x1f8)
[auchentoshan:21021] [10] testme_submitted_to_devel(main+0xe) [0x4011be]
[auchentoshan:21021] [11] /lib64/tls/
[auchentoshan:21021] [12] testme_submitted_to_devel [0x400dfa]
[auchentoshan:21021] *** End of error message ***
mpirun noticed that job rank 0 with PID 21021 on node auchentoshan exited
on signal 11 (Segmentation fault).
1 additional process aborted (not shown)

The openmpi version we have tested this on includes the latest, 1.2.2
version as well although the log message above is for version 1.2.1.
Distribution: Scientific Linux 4.4 (RHEL4 clone). The problem occurs on
both AMD64 and i386. The problem occurs both when using gcc/gfortran and
the portland group compilers version 7.0-4. Interconnect GBE.

daniels_at_auchentoshan:~ > gcc4 --version
gcc4 (GCC) 4.1.0 20060515 (Red Hat 4.1.0-18)
Copyright (C) 2006 Free Software Foundation, Inc.
This is free software; see the source for copying conditions. There is NO

daniels_at_auchentoshan:~ > gfortran --version
GNU Fortran 95 (GCC) 4.1.0 20060515 (Red Hat 4.1.0-18)
Copyright (C) 2006 Free Software Foundation, Inc.

GNU Fortran comes with NO WARRANTY, to the extent permitted by law.
You may redistribute copies of GNU Fortran
under the terms of the GNU General Public License.
For more information about these matters, see the file named COPYING

Working C code which does the same thing:
#include <stdio.h>
#include <stdlib.h>
#include <mpi.h>

int main(int argc, char **argv)
   int rank,size;
   float param=0.;
   MPI_Datatype btype[1],param_batch;
   MPI_Aint disp[1];
   int blen[1];

   if (MPI_Get_address(&param,&disp[0])!=MPI_SUCCESS)
     fprintf(stderr,"MPI_Get_address failed.\n");
   if (MPI_Type_create_struct(1,blen,disp,btype,&param_batch)!=MPI_SUCCESS)
     fprintf(stderr,"MPI_Type_Create_Struct failed.\n");
   if (MPI_Type_commit(&param_batch)!=MPI_SUCCESS)
     fprintf(stderr,"MPI_Type_Commit failed.\n");
   if (rank==0)
   if (MPI_Bcast(MPI_BOTTOM,1,param_batch,0,MPI_COMM_WORLD)!=MPI_SUCCESS)
     fprintf(stderr,"MPI_Bcast failed.\n");
   printf("Rank:%d, Size:%d, Param=%f\n",rank,size,param);
   return 0;

Best regards
Daniel Spångberg