specialfunctions: generalize gamma precision

example/specialfunctions_gamma/example_gamma.f90

@@ -1,21 +1,19 @@
 program example_gamma
-  use stdlib_kinds, only: dp, int64
+  use stdlib_kinds, only: sp, dp, int64
   use stdlib_specialfunctions_gamma, only: gamma
   implicit none
 
   integer :: i
   integer(int64) :: n
   real :: x
   real(dp) :: y
-  complex :: z
-  complex(dp) :: z1
+  complex(sp) :: z
 
   i = 10
   n = 15_int64
   x = 2.5
   y = 4.3_dp
   z = (2.3, 0.6)
-  z1 = (-4.2_dp, 3.1_dp)
 
   print *, gamma(i)              !integer gives exact result
 ! 362880
@@ -32,6 +30,4 @@ program example_gamma
   print *, gamma(z)
 ! (0.988054395, 0.383354813)
 
-  print *, gamma(z1)
-! (-2.78916032990983999E-005, 9.83164600163221218E-006)
 end program example_gamma

example/specialfunctions_gamma/example_log_gamma.f90

@@ -1,19 +1,18 @@
 program example_log_gamma
-  use stdlib_kinds, only: dp
+  use stdlib_kinds, only: sp, dp
   use stdlib_specialfunctions_gamma, only: log_gamma
   implicit none
 
   integer :: i
   real :: x
   real(dp) :: y
-  complex :: z
-  complex(dp) :: z1
+  complex(sp) :: z
 
   i = 10
   x = 8.76
   y = x
   z = (5.345, -3.467)
-  z1 = z
+
   print *, log_gamma(i)     !default single precision output
 !12.8018274
 
@@ -29,7 +28,4 @@ program example_log_gamma
 
 !(2.56165648, -5.73382425)
 
-  print *, log_gamma(z1)
-
-!(2.5616575105114614, -5.7338247782852498)
 end program example_log_gamma