Skip to content

Commit

Permalink
Promote real to real_dp as it should be
Browse files Browse the repository at this point in the history
  • Loading branch information
Cyril Gandon committed Sep 27, 2024
1 parent c14c9d8 commit a65cc2f
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions src/stdlib_specialfunctions_gamma.fypp
Original file line number Diff line number Diff line change
Expand Up @@ -528,9 +528,9 @@ contains
! Log(n!)
!
${t1}$, intent(in) :: n
real :: res
real(dp) :: res
${t1}$, parameter :: zero = 0_${k1}$, one = 1_${k1}$, two = 2_${k1}$
real, parameter :: zero_k2 = 0.0
real(dp), parameter :: zero_dp = 0.0_dp

if(n < zero) call error_stop("Error(l_factorial): Logarithm of" &
//" factorial function argument must be non-negative")
Expand All @@ -539,11 +539,11 @@ contains

case (zero)

res = zero_k2
res = zero_dp

case (one)

res = zero_k2
res = zero_dp

case (two : )

Expand Down

0 comments on commit a65cc2f

Please sign in to comment.