We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 069501e commit c0ce8dcCopy full SHA for c0ce8dc
src/stdlib_specialfunctions_gamma.fypp
@@ -671,29 +671,6 @@ contains
671
672
end do
673
674
- else !Algorithm 2 in the reference
675
-
676
- a = - x
677
- c = one / a
678
- d = p - one
679
- b = c * (a - d)
680
- n = 1
681
682
- do
683
684
- c = d * (d - one) / (a * a)
685
- d = d - 2
686
- y = c * (a - d)
687
- b = b + y
688
- n = n + 1
689
690
- if(n > int((p - 2) / 2) .or. y < b * tol_${k2}$) exit
691
692
- end do
693
694
- if(y >= b * tol_${k2}$ .and. mod(int(p), 2) /= 0) b = b + d * c / a
695
696
- g = ((-1) ** p * exp(-a + log_gamma(p) - (p - 1) * log(a)) + b) / a
697
end if
698
699
res = g
0 commit comments