Fix for Intel compiler (suppress OpenMP for loop)

cp2k · Aug 23, 2023 · 73942e2 · 73942e2
1 parent 118b17c
commit 73942e2
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/src/pw/ps_implicit_methods.F b/src/pw/ps_implicit_methods.F
@@ -1279,6 +1279,11 @@ SUBROUTINE ps_implicit_prepare_blocks(pw_pool_orig, dct_pw_grid, green, &
 
          ALLOCATE (ps_implicit_env%idx_1dto3d(data_size))
          l = 1
+         ! Suppress OpenMP (at least the Intel compiler has an issue here)
+         ! An automatic OpenMP parallelization of this loop might be tricky
+         ! because of the incrementation of l
+!$OMP PARALLEL IF(.FALSE.)
+!$OMP DO
          DO k = lb3, ub3
             DO j = lb2, ub2
                DO i = lb1, ub1
@@ -1289,6 +1294,8 @@ SUBROUTINE ps_implicit_prepare_blocks(pw_pool_orig, dct_pw_grid, green, &
                END DO
             END DO
          END DO
+!$OMP END DO
+!$OMP END PARALLEL
 
          n_contacts = SIZE(ps_implicit_env%contacts)
          n_tiles_tot = 0