Changeset 4236


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

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

Location:
palm/trunk/SOURCE
Files:
2 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
  • palm/trunk/SOURCE/tridia_solver_mod.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!
     
    249252#endif
    250253
     254          !$OMP PARALLEL PRIVATE(i,j,k)
     255
    251256!
    252257!--       Forward substitution
     
    258263             !$ACC LOOP COLLAPSE(2)
    259264#endif
     265             !$OMP DO
    260266             DO  j = nys_z, nyn_z
    261267                DO  i = nxl_z, nxr_z
     
    286292             !$ACC LOOP COLLAPSE(2)
    287293#endif
     294             !$OMP DO
    288295             DO  j = nys_z, nyn_z
    289296                DO  i = nxl_z, nxr_z
     
    301308          !$ACC END PARALLEL
    302309#endif
     310
     311          !$OMP END PARALLEL
    303312
    304313!
Note: See TracChangeset for help on using the changeset viewer.