code:

SUBROUTINE TEST()
...
INTEGER M, N
REAL,ALLOCATABLE:: TMP(:,:)
...
ALLOCATE(TMP(M,N))
....
DEALLOCATE(TMP)
...
END
 

#pragma omp parallel for
for( int i = 0; i < 10000000; ++i ) {
 //...
}

В ответ на:

---

с фортраном - нет, к сожалению
  

---

Quote:

---

Не очень хочется тыкаться, хочется понимать, что делаешь.

---

code:
MODULE MOD_OMP

PUBLIC :: TEST

CONTAINS

SUBROUTINE TEST(M, N)
IMPLICIT NONE

INTEGER, INTENT(IN) :: M, N
INTEGER :: I, J
REAL, DIMENSION(:,:), ALLOCATABLE :: TMP

ALLOCATE(TMP(M,N))

DO J = 1,N
   DO I = 1,M
      TMP(I,J) = I+J
   ENDDO
ENDDO

DEALLOCATE(TMP)

END SUBROUTINE TEST

END MODULE MOD_OMP

PROGRAM MAIN
USE MOD_OMP
IMPLICIT NONE

INTEGER :: I

!$OMP PARALLEL DO &
!$OMP DEFAULT(NONE) &
!$OMP PRIVATE(I)

DO I = 1,1000
  CALL TEST(1000,500)
ENDDO
!$OMP END PARALLEL DO

END PROGRAM

В ответ на:

---

так и есть: вместо mpi_init() надо вызывать MPI_INIT_THREAD(MPI_THREAD_FUNNELED, PROVIDED, IERROR) и по результату в PROVIDED определять, разрешено ли данному треду использовать MPI вызовы.
  

---

Quote:

---

А у тебя в программе количество тредов задается в переменной окружения?

---

Quote:

---

Проблема была в компиляторе - mpixlf77_r заменил на mpixlf90_r

---