fix(factorial): result should be of the default kind, not double default kind

Cyril Gandon · Cyril Gandon · commit cb6f989dfdc8 · 2024-09-23T09:44:43.000+02:00
Fixes #875
diff --git a/src/stdlib_specialfunctions_gamma.fypp b/src/stdlib_specialfunctions_gamma.fypp
@@ -547,7 +547,7 @@ contains
 
         case (two : )
 
-            res = l_gamma(n + 1, 1.0D0)
+            res = l_gamma(n + 1, 1.0)
 
         end select
     end function l_factorial_${t1[0]}$${k1}$
