From 1e7a95267171457d76235364b1bcf2b7847fef8e Mon Sep 17 00:00:00 2001
From: Avraham Adler <avraham.adler@gmail.com>
Date: Sun, 9 Jun 2024 21:42:16 -0400
Subject: [PATCH] More SIMD declarations

---
 src/delaporte.f90 | 6 +++---
 src/utils.f90     | 2 +-
 2 files changed, 4 insertions(+), 4 deletions(-)

diff --git a/src/delaporte.f90 b/src/delaporte.f90
index 416fcc7..f7304e3 100644
--- a/src/delaporte.f90
+++ b/src/delaporte.f90
@@ -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
@@ -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
@@ -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
diff --git a/src/utils.f90 b/src/utils.f90
index fae8ba1..862c504 100644
--- a/src/utils.f90
+++ b/src/utils.f90
@@ -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