Open MPI User's Mailing List Archives

Subject: [OMPI users] memory allocation - Fortran
From: Martin Siegert (siegert_at_[hidden])
Date: 2010-11-12 16:34:28

• Next message: George Bosilca: "Re: [OMPI users] memcpy overlap in ompi_ddt_copy_content_same_ddt and glibc 2.12"
• Previous message: Jeff Squyres: "Re: [OMPI users] Compile problems with 1.3.2"

Hi,

section 4.11

http://www.mpi-forum.org/docs/mpi-20-html/node54.htm

describes how memory allocation can be done in a Fortran program:

REAL A
POINTER (P, A(100,100)) ! no memory is allocated
CALL MPI_ALLOC_MEM(4*100*100, MPI_INFO_NULL, P, IERR)
! memory is allocated
...
A(3,5) = 2.71;
...
CALL MPI_FREE_MEM(A, IERR) ! memory is freed

However, this uses non-standard Cray pointers. I am wondering how this
can be coded using the standard iso_c_binding module - I do not have
much experience using that module, but apparently it provides C
pointers:

use, intrinsic :: iso_c_binding, only: c_float, c_ptr, c_f_pointer
real (c_float), pointer :: a(:,:)
type (c_ptr) :: p

call MPI_Alloc_mem(4*100*100, MPI_INFO_NULL, p, ierr)
call c_f_pointer(p, a, (/100, 100/))
...
a(3,5) = 2.71
...
call MPI_Free_mem(a, ierr)

Is that going to work? Or is there a different way of recoding the
example?

Thanks in advance!

Cheers,
Martin

-- 
Martin Siegert
Head, Research Computing
WestGrid/ComputeCanada Site Lead
IT Services                                phone: 778 782-4691
Simon Fraser University                    fax:   778 782-4242
Burnaby, British Columbia                  email: siegert_at_[hidden]
Canada  V5A 1S6

• Next message: George Bosilca: "Re: [OMPI users] memcpy overlap in ompi_ddt_copy_content_same_ddt and glibc 2.12"
• Previous message: Jeff Squyres: "Re: [OMPI users] Compile problems with 1.3.2"