1
1
openmpi/ompi/mpi/fortran/use-mpi-f08/mpi-f08-sizeof.F90
Jeff Squyres 253444c6d0 == Highlights ==
1. New mpifort wrapper compiler: you can utilize mpif.h, use mpi, and use mpi_f08 through this one wrapper compiler
 1. mpif77 and mpif90 still exist, but are sym links to mpifort and may be removed in a future release
 1. The mpi module has been re-implemented and is significantly "mo' bettah"
 1. The mpi_f08 module offers many, many improvements over mpif.h and the mpi module

This stuff is coming from a VERY long-lived mercurial branch (3 years!); it'll almost certainly take a few SVN commits and a bunch of testing before I get it correctly committed to the SVN trunk.

== More details ==

Craig Rasmussen and I have been working with the MPI-3 Fortran WG and Fortran J3 committees for a long, long time to make a prototype MPI-3 Fortran bindings implementation.  We think we're at a stable enough state to bring this stuff back to the trunk, with the goal of including it in OMPI v1.7.  

Special thanks go out to everyone who has been incredibly patient and helpful to us in this journey:

 * Rolf Rabenseifner/HLRS (mastermind/genius behind the entire MPI-3 Fortran effort)
 * The Fortran J3 committee
 * Tobias Burnus/gfortran
 * Tony !Goetz/Absoft
 * Terry !Donte/Oracle
 * ...and probably others whom I'm forgetting :-(

There's still opportunities for optimization in the mpi_f08 implementation, but by and large, it is as far along as it can be until Fortran compilers start implementing the new F08 dimension(..) syntax.

Note that gfortran is currently unsupported for the mpi_f08 module and the new mpi module.  gfortran users will a) fall back to the same mpi module implementation that is in OMPI v1.5.x, and b) not get the new mpi_f08 module.  The gfortran maintainers are actively working hard to add the necessary features to support both the new mpi_f08 module and the new mpi module implementations.  This will take some time.

As mentioned above, ompi/mpi/f77 and ompi/mpi/f90 no longer exist.  All the fortran bindings implementations have been collated under ompi/mpi/fortran; each implementation has its own subdirectory:

{{{
ompi/mpi/fortran/
  base/               - glue code
  mpif-h/             - what used to be ompi/mpi/f77
  use-mpi-tkr/        - what used to be ompi/mpi/f90
  use-mpi-ignore-tkr/ - new mpi module implementation
  use-mpi-f08/        - new mpi_f08 module implementation
}}}

There's also a prototype 6-function-MPI implementation under use-mpi-f08-desc that emulates the new F08 dimension(..) syntax that isn't fully available in Fortran compilers yet.  We did that to prove it to ourselves that it could be done once the compilers fully support it.  This directory/implementation will likely eventually replace the use-mpi-f08 version.

Other things that were done:

 * ompi_info grew a few new output fields to describe what level of Fortran support is included
 * Existing Fortran examples in examples/ were renamed; new mpi_f08 examples were added
 * The old Fortran MPI libraries were renamed:
   * libmpi_f77 -> libmpi_mpifh
   * libmpi_f90 -> libmpi_usempi
 * The configury for Fortran was consolidated and significantly slimmed down.  Note that the F77 env variable is now IGNORED for configure; you should only use FC. Example:
{{{
shell$ ./configure CC=icc CXX=icpc FC=ifort ...
}}}

All of this work was done in a Mercurial branch off the SVN trunk, and hosted at Bitbucket.  This branch has got to be one of OMPI's longest-running branches.  Its first commit was Tue Apr 07 23:01:46 2009 -0400 -- it's over 3 years old!  :-)  We think we've pulled in all relevant changes from the OMPI trunk (e.g., Fortran implementations of the new MPI-3 MPROBE stuff for mpif.h, use mpi, and use mpi_f08, and the recent Fujitsu Fortran patches).

I anticipate some instability when we bring this stuff into the trunk, simply because it touches a LOT of code in the MPI layer in the OMPI code base.  We'll try our best to make it as pain-free as possible, but please bear with us when it is committed.

This commit was SVN r26283.
2012-04-18 15:57:29 +00:00

421 строка
12 KiB
Fortran

! -*- f90 -*-
!
! Copyright (c) 2006-2011 Cisco Systems, Inc. All rights reserved.
! Copyright (c) 2009-2012 Los Alamos National Security, LLC.
! All rights reserved.
!
! $COPYRIGHT$
!
! Additional copyrights may follow
!
! $HEADER$
!
#include "ompi/mpi/fortran/configure-fortran-output.h"
MODULE mpi_f08_sizeof
use, intrinsic :: ISO_C_BINDING
interface MPI_Sizeof
#if OMPI_HAVE_FORTRAN_INTEGER1
module procedure MPI_Sizeof_integer_s_1, MPI_Sizeof_integer_a_1
#endif
#if OMPI_HAVE_FORTRAN_INTEGER2
module procedure MPI_Sizeof_integer_s_2, MPI_Sizeof_integer_a_2
#endif
#if OMPI_HAVE_FORTRAN_INTEGER4
module procedure MPI_Sizeof_integer_s_4, MPI_Sizeof_integer_a_4
#endif
#if OMPI_HAVE_FORTRAN_INTEGER8
module procedure MPI_Sizeof_integer_s_8, MPI_Sizeof_integer_a_8
#endif
#if OMPI_HAVE_FORTRAN_INTEGER16
module procedure MPI_Sizeof_integer_s_16, MPI_Sizeof_integer_a_16
#endif
#if OMPI_HAVE_FORTRAN_REAL2
module procedure MPI_Sizeof_real_s_2, MPI_Sizeof_real_a_2
#endif
#if OMPI_HAVE_FORTRAN_REAL4
module procedure MPI_Sizeof_real_s_4, MPI_Sizeof_real_a_4
#endif
#if OMPI_HAVE_FORTRAN_REAL8
module procedure MPI_Sizeof_real_s_8, MPI_Sizeof_real_a_8
#endif
#if OMPI_HAVE_FORTRAN_REAL16
module procedure MPI_Sizeof_real_s_16, MPI_Sizeof_real_a_16
#endif
#if OMPI_HAVE_FORTRAN_COMPLEX4
module procedure MPI_Sizeof_complex_s_4, MPI_Sizeof_complex_a_4
#endif
#if OMPI_HAVE_FORTRAN_COMPLEX8
module procedure MPI_Sizeof_complex_s_8, MPI_Sizeof_complex_a_8
#endif
#if OMPI_HAVE_FORTRAN_COMPLEX16
module procedure MPI_Sizeof_complex_s_16, MPI_Sizeof_complex_a_16
#endif
#if OMPI_HAVE_FORTRAN_COMPLEX32
module procedure MPI_Sizeof_complex_s_32, MPI_Sizeof_complex_a_32
#endif
end interface MPI_Sizeof
CONTAINS
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Sizes of integers
!
#if OMPI_HAVE_FORTRAN_INTEGER1
subroutine MPI_Sizeof_integer_s_1(x, size, ierror)
implicit none
INTEGER(KIND=OMPI_KIND_FORTRAN_INTEGER1), INTENT(IN) :: x
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_INTEGER1
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_integer_s_1
subroutine MPI_Sizeof_integer_a_1(x, size, ierror)
use, intrinsic :: ISO_C_BINDING
implicit none
INTEGER(KIND=OMPI_KIND_FORTRAN_INTEGER1), INTENT(IN) :: x(*)
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_INTEGER1
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_integer_a_1
#endif
#if OMPI_HAVE_FORTRAN_INTEGER2
subroutine MPI_Sizeof_integer_s_2(x, size, ierror)
implicit none
INTEGER(KIND=OMPI_KIND_FORTRAN_INTEGER2), INTENT(IN) :: x
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_INTEGER2
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_integer_s_2
subroutine MPI_Sizeof_integer_a_2(x, size, ierror)
use, intrinsic :: ISO_C_BINDING
implicit none
INTEGER(KIND=OMPI_KIND_FORTRAN_INTEGER2), INTENT(IN) :: x(*)
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_INTEGER2
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_integer_a_2
#endif
#if OMPI_HAVE_FORTRAN_INTEGER4
subroutine MPI_Sizeof_integer_s_4(x, size, ierror)
implicit none
INTEGER(KIND=OMPI_KIND_FORTRAN_INTEGER4), INTENT(IN) :: x
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_INTEGER4
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_integer_s_4
subroutine MPI_Sizeof_integer_a_4(x, size, ierror)
use, intrinsic :: ISO_C_BINDING
implicit none
INTEGER(KIND=OMPI_KIND_FORTRAN_INTEGER4), INTENT(IN) :: x(*)
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_INTEGER4
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_integer_a_4
#endif
#if OMPI_HAVE_FORTRAN_INTEGER8
subroutine MPI_Sizeof_integer_s_8(x, size, ierror)
implicit none
INTEGER(KIND=OMPI_KIND_FORTRAN_INTEGER8), INTENT(IN) :: x
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_INTEGER8
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_integer_s_8
subroutine MPI_Sizeof_integer_a_8(x, size, ierror)
use, intrinsic :: ISO_C_BINDING
implicit none
INTEGER(KIND=OMPI_KIND_FORTRAN_INTEGER8), INTENT(IN) :: x(*)
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_INTEGER8
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_integer_a_8
#endif
#if OMPI_HAVE_FORTRAN_INTEGER16
subroutine MPI_Sizeof_integer_s_16(x, size, ierror)
implicit none
INTEGER(KIND=OMPI_KIND_FORTRAN_INTEGER16), INTENT(IN) :: x
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_INTEGER16
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_integer_s_16
subroutine MPI_Sizeof_integer_a_16(x, size, ierror)
use, intrinsic :: ISO_C_BINDING
implicit none
INTEGER(KIND=OMPI_KIND_FORTRAN_INTEGER16), INTENT(IN) :: x(*)
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_INTEGER16
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_integer_a_16
#endif
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Sizes of reals
!
#if OMPI_HAVE_FORTRAN_REAL2
subroutine MPI_Sizeof_real_s_2(x, size, ierror)
implicit none
REAL(KIND=OMPI_KIND_FORTRAN_REAL2), INTENT(IN) :: x
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_REAL2
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_real_s_2
subroutine MPI_Sizeof_real_a_2(x, size, ierror)
use, intrinsic :: ISO_C_BINDING
implicit none
REAL(KIND=OMPI_KIND_FORTRAN_REAL2), INTENT(IN) :: x(*)
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_REAL2
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_real_a_2
#endif
#if OMPI_HAVE_FORTRAN_REAL4
subroutine MPI_Sizeof_real_s_4(x, size, ierror)
implicit none
REAL(KIND=OMPI_KIND_FORTRAN_REAL4), INTENT(IN) :: x
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_REAL4
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_real_s_4
subroutine MPI_Sizeof_real_a_4(x, size, ierror)
use, intrinsic :: ISO_C_BINDING
implicit none
REAL(KIND=OMPI_KIND_FORTRAN_REAL4), INTENT(IN) :: x(*)
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_REAL4
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_real_a_4
#endif
#if OMPI_HAVE_FORTRAN_REAL8
subroutine MPI_Sizeof_real_s_8(x, size, ierror)
implicit none
REAL(KIND=OMPI_KIND_FORTRAN_REAL8), INTENT(IN) :: x
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_REAL8
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_real_s_8
subroutine MPI_Sizeof_real_a_8(x, size, ierror)
use, intrinsic :: ISO_C_BINDING
implicit none
REAL(KIND=OMPI_KIND_FORTRAN_REAL8), INTENT(IN) :: x(*)
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_REAL8
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_real_a_8
#endif
#if OMPI_HAVE_FORTRAN_REAL16
subroutine MPI_Sizeof_real_s_16(x, size, ierror)
implicit none
REAL(KIND=OMPI_KIND_FORTRAN_REAL16), INTENT(IN) :: x
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_REAL16
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_real_s_16
subroutine MPI_Sizeof_real_a_16(x, size, ierror)
use, intrinsic :: ISO_C_BINDING
implicit none
REAL(KIND=OMPI_KIND_FORTRAN_REAL16), INTENT(IN) :: x(*)
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_REAL16
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_real_a_16
#endif
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Sizes of complex
!
#if OMPI_HAVE_FORTRAN_COMPLEX4
subroutine MPI_Sizeof_complex_s_4(x, size, ierror)
implicit none
COMPLEX(KIND=OMPI_KIND_FORTRAN_COMPLEX4), INTENT(IN) :: x
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_COMPLEX4
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_complex_s_4
subroutine MPI_Sizeof_complex_a_4(x, size, ierror)
use, intrinsic :: ISO_C_BINDING
implicit none
COMPLEX(KIND=OMPI_KIND_FORTRAN_COMPLEX4), INTENT(IN) :: x(*)
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_COMPLEX4
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_complex_a_4
#endif
#if OMPI_HAVE_FORTRAN_COMPLEX8
subroutine MPI_Sizeof_complex_s_8(x, size, ierror)
implicit none
COMPLEX(KIND=OMPI_KIND_FORTRAN_COMPLEX8), INTENT(IN) :: x
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_COMPLEX8
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_complex_s_8
subroutine MPI_Sizeof_complex_a_8(x, size, ierror)
use, intrinsic :: ISO_C_BINDING
implicit none
COMPLEX(KIND=OMPI_KIND_FORTRAN_COMPLEX8), INTENT(IN) :: x(*)
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_COMPLEX8
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_complex_a_8
#endif
#if OMPI_HAVE_FORTRAN_COMPLEX16
subroutine MPI_Sizeof_complex_s_16(x, size, ierror)
implicit none
COMPLEX(KIND=OMPI_KIND_FORTRAN_COMPLEX16), INTENT(IN) :: x
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_COMPLEX16
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_complex_s_16
subroutine MPI_Sizeof_complex_a_16(x, size, ierror)
use, intrinsic :: ISO_C_BINDING
implicit none
COMPLEX(KIND=OMPI_KIND_FORTRAN_COMPLEX16), INTENT(IN) :: x(*)
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_COMPLEX16
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_complex_a_16
#endif
#if OMPI_HAVE_FORTRAN_COMPLEX32
subroutine MPI_Sizeof_complex_s_32(x, size, ierror)
implicit none
COMPLEX(KIND=OMPI_KIND_FORTRAN_COMPLEX32), INTENT(IN) :: x
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_COMPLEX32
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_complex_s_32
subroutine MPI_Sizeof_complex_a_32(x, size, ierror)
use, intrinsic :: ISO_C_BINDING
implicit none
COMPLEX(KIND=OMPI_KIND_FORTRAN_COMPLEX32), INTENT(IN) :: x(*)
INTEGER, INTENT(OUT) :: size
INTEGER, OPTIONAL, INTENT(OUT) :: ierror
size = OMPI_SIZEOF_FORTRAN_COMPLEX32
if (present(ierror)) ierror = 0
end subroutine MPI_Sizeof_complex_a_32
#endif
END MODULE mpi_f08_sizeof