MPI_finalize 是否释放内存?
Does MPI_finalize release memory?
在下面的代码中,我有一个在 MPI 中使用的数组 b。据我了解,每个处理器甚至在调用 MPI_INIT
之前都会获得 b 的副本。但是我们调用 MPI_FINALIZE
之后会发生什么?那块内存是否仍然可供每个处理器使用?
以类似的方式,如果将 b 声明为指针,它会在 MPI_INIT
-MPI_FINALIZE
内部分配,但它是没有释放?完成 MPI 后该内存是否仍然可用?
program main
use mpi
implicit none
integer myr, numpr, ier
integer b(1000)
call MPI_INIT(ier)
call MPI_COMM_RANK(MPI_COMM_WORLD, myr, ier)
call MPI_COMM_SIZE(MPI_COMM_WORLD, numpr, ier)
if (myr .eq. 0) then
!initialize b array
endif
call MPI_BCAST(b, 100, MPI_INTEGER, 0, MPI_COMM_WORLD, ier)
call MPI_FINALIZE(ier)
!do more calculations with b
end
如果您想象自己编写的代码没有任何 MPI 内容,您会看到每个处理器都以大小为 1000 的 B 数组开头,因为您是这样声明它的:
integer b(1000)
MPI_Init 和 MPI_Finalise 均未参与分配或取消分配任何此内存。
同样,您可以在 运行 时间 (C) 分配一个数组,它会一直存在,直到您显式取消分配它:
PROGRAM main
use mpi
implicit none
integer myr, numpr, ier
integer b(1000)
INTEGER, ALLOCATABLE :: C(:)
call MPI_INIT(ier)
call MPI_COMM_RANK(MPI_COMM_WORLD, myr, ier)
call MPI_COMM_SIZE(MPI_COMM_WORLD, numpr, ier)
ALLOCATE(C(1000))
if (myr .eq. 0) then
b = 100 ! Set all values to 100
c = 99 ! Ditto 99
endif
call MPI_BCAST(b, 1000, MPI_INTEGER, 0, MPI_COMM_WORLD, ier)
call MPI_BCAST(c, 1000, MPI_INTEGER, 0, MPI_COMM_WORLD, ier)
call MPI_FINALIZE(ier)
PRINT *, myr, B(200)
PRINT *, myr, C(200)
DEALLOCATE(C)
END PROGRAM main
产生输出:
1 100
1 99
0 100
0 99
另外,请注意您的初始代码中有错字(我认为)。您只发送 B 的前 100 个成员(大小为 1000)。
在下面的代码中,我有一个在 MPI 中使用的数组 b。据我了解,每个处理器甚至在调用 MPI_INIT
之前都会获得 b 的副本。但是我们调用 MPI_FINALIZE
之后会发生什么?那块内存是否仍然可供每个处理器使用?
以类似的方式,如果将 b 声明为指针,它会在 MPI_INIT
-MPI_FINALIZE
内部分配,但它是没有释放?完成 MPI 后该内存是否仍然可用?
program main
use mpi
implicit none
integer myr, numpr, ier
integer b(1000)
call MPI_INIT(ier)
call MPI_COMM_RANK(MPI_COMM_WORLD, myr, ier)
call MPI_COMM_SIZE(MPI_COMM_WORLD, numpr, ier)
if (myr .eq. 0) then
!initialize b array
endif
call MPI_BCAST(b, 100, MPI_INTEGER, 0, MPI_COMM_WORLD, ier)
call MPI_FINALIZE(ier)
!do more calculations with b
end
如果您想象自己编写的代码没有任何 MPI 内容,您会看到每个处理器都以大小为 1000 的 B 数组开头,因为您是这样声明它的:
integer b(1000)
MPI_Init 和 MPI_Finalise 均未参与分配或取消分配任何此内存。
同样,您可以在 运行 时间 (C) 分配一个数组,它会一直存在,直到您显式取消分配它:
PROGRAM main
use mpi
implicit none
integer myr, numpr, ier
integer b(1000)
INTEGER, ALLOCATABLE :: C(:)
call MPI_INIT(ier)
call MPI_COMM_RANK(MPI_COMM_WORLD, myr, ier)
call MPI_COMM_SIZE(MPI_COMM_WORLD, numpr, ier)
ALLOCATE(C(1000))
if (myr .eq. 0) then
b = 100 ! Set all values to 100
c = 99 ! Ditto 99
endif
call MPI_BCAST(b, 1000, MPI_INTEGER, 0, MPI_COMM_WORLD, ier)
call MPI_BCAST(c, 1000, MPI_INTEGER, 0, MPI_COMM_WORLD, ier)
call MPI_FINALIZE(ier)
PRINT *, myr, B(200)
PRINT *, myr, C(200)
DEALLOCATE(C)
END PROGRAM main
产生输出:
1 100
1 99
0 100
0 99
另外,请注意您的初始代码中有错字(我认为)。您只发送 B 的前 100 个成员(大小为 1000)。