for all non-sp kinds, use max available precision

Author: perazz

src/stdlib_specialfunctions_gamma.fypp

@@ -219,7 +219,7 @@ contains
     #! gamma will use the next available more accurate KIND for the 
     #! internal more accurate solver.
     #:for i, k1, t1, i1 in IDX_CMPLX_KINDS_TYPES[:-1]
-    #:set k2 = CMPLX_KINDS[i + 1]
+    #:set k2 = CMPLX_KINDS[i + 1] if k1 == "sp" else CMPLX_KINDS[-1]
     #:set t2 = "real({})".format(k2)
     impure elemental function gamma_${t1[0]}$${k1}$(z) result(res)
         ${t1}$, intent(in) :: z
@@ -408,7 +408,7 @@ contains
     #! gamma will use the next available more accurate KIND for the 
     #! internal more accurate solver.
     #:for i, k1, t1, i1 in IDX_CMPLX_KINDS_TYPES[:-1]
-    #:set k2 = CMPLX_KINDS[i + 1]
+    #:set k2 = CMPLX_KINDS[i + 1] if k1 == "sp" else CMPLX_KINDS[-1]
     #:set t2 = "real({})".format(k2)
     impure elemental function l_gamma_${t1[0]}$${k1}$(z) result (res)
     !
@@ -549,7 +549,7 @@ contains
     #! gamma will use the next available more accurate KIND for the 
     #! internal more accurate solver.
     #:for i, k1, t1, i1 in IDX_REAL_KINDS_TYPES[:-1]
-    #:set k2 = REAL_KINDS[i + 1]
+    #:set k2 = REAL_KINDS[i + 1] if k1 == "sp" else REAL_KINDS[-1]
     #:set t2 = REAL_TYPES[i + 1]
     impure elemental function gpx_${t1[0]}$${k1}$(p, x) result(res)
     !

test/specialfunctions/test_specialfunctions_gamma.fypp

@@ -185,15 +185,14 @@ contains
 
         do i = 1, n
 
-            err = state_type(STDLIB_VALUE_ERROR,'Complex ${k1}$ failed: x=',x(i),' gamma=',gamma(x(i)), &
-                             'expected=',ans(i),' tol=',tol_${k1}$)
+            err = state_type(STDLIB_VALUE_ERROR,'Complex ${k1}$ failed: x=',x(i), &
+                             ' gamma=',gamma(x(i)), &
+                             'expected=',ans(i), &
+                             ' tol=',tol_${k1}$)
 
             call check(error, gamma(x(i)), ans(i), err%print(),&
                        thr = tol_${k1}$, rel = .true.)
 
-            !call check(error, gamma(x(i)), ans(i), "Complex kind ${k1}$ failed",&
-            !           thr = tol_${k1}$, rel = .true.)
-
         end do
 
         #:endif
@@ -257,8 +256,13 @@ contains
 
         do i = 1, n
 
-            call check(error, log_gamma(x(i)), ans(i), "Complex kind ${k1}$ "  &
-              //"failed", thr = tol_${k1}$, rel = .true.)
+            err = state_type(STDLIB_VALUE_ERROR,'Complex ${k1}$ failed: x=',x(i), &
+                             ' log(gamma)=',log_gamma(x(i)), &
+                             'expected=',ans(i), &
+                             ' tol=',tol_${k1}$)
+
+            call check(error, log_gamma(x(i)), ans(i), err%print(),&
+                       thr = tol_${k1}$, rel = .true.)
 
         end do