Ignore:
Timestamp:
Sep 25, 2019 11:26:18 AM (21 months ago)
Author:
knoop
Message:

Added missing OpenMP directives within "transpose" and "tridia_solver"

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/transpose.f90

    r4182 r4236  
    2525! -----------------
    2626! $Id$
     27! Added missing OpenMP directives
     28!
     29! 4182 2019-08-22 15:20:23Z scharf
    2730! Corrected "Former revisions" section
    2831!
     
    162165!--    Reorder transposed array
    163166!$OMP  PARALLEL PRIVATE ( i, j, k, l, ys )
    164 !$OMP  DO
    165167       DO  l = 0, pdims(2) - 1
    166168          ys = 0 + l * ( nyn_x - nys_x + 1 )
     
    169171          !$ACC PRESENT(f_out, work)
    170172#endif
     173          !$OMP DO
    171174          DO  i = nxl_y, nxr_y
    172175             DO  k = nzb_y, nzt_y
     
    176179             ENDDO
    177180          ENDDO
     181          !$OMP END DO NOWAIT
    178182       ENDDO
    179183!$OMP  END PARALLEL
     
    298302!--    Reorder input array for transposition
    299303!$OMP  PARALLEL PRIVATE ( i, j, k, l, xs )
    300 !$OMP  DO
    301304       DO  l = 0, pdims(1) - 1
    302305          xs = 0 + l * nnx
     
    305308          !$ACC PRESENT(work, f_in)
    306309#endif
     310          !$OMP DO
    307311          DO  k = nzb_x, nzt_x
    308312             DO  i = xs, xs + nnx - 1
     
    312316             ENDDO
    313317          ENDDO
     318          !$OMP END DO NOWAIT
    314319       ENDDO
    315320!$OMP  END PARALLEL
     
    460465!--    Reorder input array for transposition
    461466!$OMP  PARALLEL PRIVATE ( i, j, k, l, ys )
    462 !$OMP  DO
    463467       DO  l = 0, pdims(2) - 1
    464468          ys = 0 + l * ( nyn_x - nys_x + 1 )
     
    467471          !$ACC PRESENT(work, f_in)
    468472#endif
     473          !$OMP DO
    469474          DO  i = nxl_y, nxr_y
    470475             DO  k = nzb_y, nzt_y
     
    474479             ENDDO
    475480          ENDDO
     481          !$OMP END DO NOWAIT
    476482       ENDDO
    477483!$OMP  END PARALLEL
     
    749755!--    Reorder transposed array
    750756!$OMP  PARALLEL PRIVATE ( i, j, k, l, zs )
    751 !$OMP  DO
    752757       DO  l = 0, pdims(1) - 1
    753758          zs = 1 + l * ( nzt_y - nzb_y + 1 )
     
    756761          !$ACC PRESENT(f_out, work)
    757762#endif
     763          !$OMP DO
    758764          DO  j = nys_z, nyn_z
    759765             DO  k = zs, zs + nzt_y - nzb_y
     
    763769             ENDDO
    764770          ENDDO
     771          !$OMP END DO NOWAIT
    765772       ENDDO
    766773!$OMP  END PARALLEL
     
    909916!--    Reorder transposed array
    910917!$OMP  PARALLEL PRIVATE ( i, j, k, l, xs )
    911 !$OMP  DO
    912918       DO  l = 0, pdims(1) - 1
    913919          xs = 0 + l * nnx
     
    916922          !$ACC PRESENT(f_out, work)
    917923#endif
     924          !$OMP DO
    918925          DO  k = nzb_x, nzt_x
    919926             DO  i = xs, xs + nnx - 1
     
    923930             ENDDO
    924931          ENDDO
     932          !$OMP END DO NOWAIT
    925933       ENDDO
    926934!$OMP  END PARALLEL
     
    10281036!--    Reorder input array for transposition
    10291037!$OMP  PARALLEL PRIVATE ( i, j, k, l, zs )
    1030 !$OMP  DO
    10311038       DO  l = 0, pdims(1) - 1
    10321039          zs = 1 + l * ( nzt_y - nzb_y + 1 )
     
    10351042          !$ACC PRESENT(work, f_in)
    10361043#endif
     1044          !$OMP DO
    10371045          DO  j = nys_z, nyn_z
    10381046             DO  k = zs, zs + nzt_y - nzb_y
     
    10421050             ENDDO
    10431051          ENDDO
     1052          !$OMP END DO NOWAIT
    10441053       ENDDO
    10451054!$OMP  END PARALLEL
Note: See TracChangeset for help on using the changeset viewer.