[Fortran] 纯文本查看 复制代码
module sx_type
type Mdata
integer :: nx(3)
real,allocatable :: sx(:)
end type
contains
subroutine deall(tmp)
type(Mdata) :: tmp
if(allocated(tmp%sx)) deallocate(tmp%sx)
end subroutine deall
end module
!----------------program-----------
program testopenmp
use sx_type
use omp_lib
implicit none
integer ix,iz
type(Mdata) :: test
call omp_set_num_threads(2)
test%nx(1)=4
test%nx(2)=5
test%nx(3)=6
!$omp parallel do private(test,ix,iz)
do ix=1,3
allocate(test%sx(test%nx(ix))) ! allocate test->sx on threads
test%sx=2.0
do iz=1,4
test%sx(iz)=iz+omp_get_thread_num()
print *,'thread ',omp_get_thread_num(),'sx=',test%sx(iz)
end do
call deall(test) !deallocate test
end do
!$omp end parallel do
end program