Skip to content

Commit

Permalink
More SIMD declarations
Browse files Browse the repository at this point in the history
  • Loading branch information
aadler committed Jun 10, 2024
1 parent 4a221ab commit 1e7a952
Show file tree
Hide file tree
Showing 2 changed files with 4 additions and 4 deletions.
6 changes: 3 additions & 3 deletions src/delaporte.f90
Original file line number Diff line number Diff line change
Expand Up @@ -107,7 +107,7 @@ module delaporte
!-------------------------------------------------------------------------------

pure elemental function ddelap_f_s(x, alpha, beta, lambda) result(pmf)
!$omp declare simd (ddelap_f_s) notinbranch
!$omp declare simd(ddelap_f_s) linear(ref(x)) notinbranch

real(kind = c_double), intent(in) :: x, alpha, beta, lambda
real(kind = c_double) :: pmf, ii, kk
Expand Down Expand Up @@ -180,7 +180,7 @@ end subroutine ddelap_f
!-------------------------------------------------------------------------------

pure elemental function pdelap_f_s(q, alpha, beta, lambda) result(cdf)
!$omp declare simd (pdelap_f_s) notinbranch
!$omp declare simd(pdelap_f_s) linear(ref(q)) notinbranch

real(kind = c_double), intent(in) :: q, alpha, beta, lambda
real(kind = c_double) :: cdf
Expand Down Expand Up @@ -274,7 +274,7 @@ end subroutine pdelap_f
!-------------------------------------------------------------------------------

pure elemental function qdelap_f_s(p, alpha, beta, lambda) result(value)
!$omp declare simd (qdelap_f_s) notinbranch
!$omp declare simd(qdelap_f_s) linear(ref(p)) notinbranch

real(kind = c_double), intent(in) :: p, alpha, beta, lambda
real(kind = c_double) :: testcdf, value
Expand Down
2 changes: 1 addition & 1 deletion src/utils.f90
Original file line number Diff line number Diff line change
Expand Up @@ -135,7 +135,7 @@ end subroutine sOMPT_f
!-------------------------------------------------------------------------------

pure elemental function imk(i, k) result(j)
!$omp declare simd(imk) linear(uval(i):1) uniform(k) notinbranch
!$omp declare simd(imk) linear(uval(i)) uniform(k) notinbranch

integer(kind = c_int), intent(in) :: i, k
integer(kind = c_int) :: j
Expand Down

0 comments on commit 1e7a952

Please sign in to comment.