I've attached a small sample program that demonstrates the problem. You can toggle working/non-working behaviour by toggling commenting on line 27.
I've tried to open a bug report, but the system isn't letting me register for Trac: Trac detected an internal error: KeyError: 'recaptcha_challenge_field' Shane On Friday, November 30, 2012 10:35:57 AM users-requ...@open-mpi.org wrote: All, I have a Fortran code that works quite well with OpenMPI 1.4.3 where I create a handle using: call MPI_TYPE_CREATE_F90_INTEGER(9, COMM_INT4, ierror) and then do a reduction with: call MPI_ALLREDUCE(send_buffer, buffer, count, COMM_INT4, MPI_SUM, communicator, ierror) However, this fails with OpenMPI 1.6.3 stating: An error occurred in MPI_Allreduce: the reduction operation MPI_SUM is not defined for non-intrinsic datatypes (attempted with datatype named "Dup MPI_INT") The MPI standard states that MPI_SUM works with Fortran integers, and Fortran integers are defined in Section 5.9.2 as: MPI_INTEGER, MPI_AINT, MPI_OFFSET, *and handles returned from MPI_TYPE_CREATE_F90_INTEGER*, and if available: MPI_INTEGER1, MPI_INTEGER2, MPI_INTEGER4, MPI_INTEGER8, MPI_INTEGER16 (emphasis mine) However the manual page on MPI_Reduce for OpenMPI only states that Fortran integers are: Fortran integer: MPI_INTEGER Does OpenMPI not support using MPI_SUM with Fortran integer handles? I'm hoping that it's just an oversight as it used to work... System: Fedora 17 Compilers: GCC 4.7.2 OpenMPI's tested: 1.4.3 (works), 1.5.4 (doesn't work), 1.6.3 (doesn't work) Thanks for any insight, Shane
program mpi_sum_test implicit none include 'mpif.h' integer, parameter :: my_int_kind = selected_int_kind(9) integer(kind=my_int_kind) :: my_ints integer(kind=my_int_kind) :: my_ints_sum integer :: test_fortran_integer_handle integer :: test_fortran_integer_size integer :: ierror integer :: node_num integer :: num_nodes call mpi_init(ierror) call mpi_comm_rank(mpi_comm_world, node_num, ierror) call mpi_comm_size(mpi_comm_world, num_nodes, ierror) call mpi_type_create_f90_integer(9, test_fortran_integer_handle, ierror) test_fortran_integer_size = sizeof(my_ints) if (node_num == 0) write (*,*) 'integer handle created as: ', test_fortran_integer_handle, ' mpi_int is ', mpi_int ! *** An error occurred in MPI_Reduce: the reduction operation MPI_SUM is not defined for non-intrinsic datatypes ! (attempted with datatype named "Dup MPI_INT") ! Uncomment next line to make this program work (since MPI reports my handle is a duplicate of MPI_INT we can just use that): !test_fortran_integer_handle = mpi_int my_ints = int(node_num+1, kind=my_int_kind) call mpi_reduce(my_ints, my_ints_sum, 1, test_fortran_integer_handle, mpi_sum, 0, mpi_comm_world, ierror) if (node_num == 0) write (*,*) 'Sum = ', my_ints_sum call mpi_finalize(ierror) end program mpi_sum_test