cgamr.f 1.0 KB

123456789101112131415161718192021222324252627282930313233343536
  1. *DECK CGAMR
  2. COMPLEX FUNCTION CGAMR (Z)
  3. C***BEGIN PROLOGUE CGAMR
  4. C***PURPOSE Compute the reciprocal of the Gamma function.
  5. C***LIBRARY SLATEC (FNLIB)
  6. C***CATEGORY C7A
  7. C***TYPE COMPLEX (GAMR-S, DGAMR-D, CGAMR-C)
  8. C***KEYWORDS FNLIB, RECIPROCAL GAMMA FUNCTION, SPECIAL FUNCTIONS
  9. C***AUTHOR Fullerton, W., (LANL)
  10. C***DESCRIPTION
  11. C
  12. C CGAMR(Z) calculates the reciprocal gamma function for COMPLEX
  13. C argument Z. This is a preliminary version that is not accurate.
  14. C
  15. C***REFERENCES (NONE)
  16. C***ROUTINES CALLED CLNGAM, XERCLR, XGETF, XSETF
  17. C***REVISION HISTORY (YYMMDD)
  18. C 770701 DATE WRITTEN
  19. C 861211 REVISION DATE from Version 3.2
  20. C 891214 Prologue converted to Version 4.0 format. (BAB)
  21. C***END PROLOGUE CGAMR
  22. COMPLEX Z, CLNGAM
  23. C***FIRST EXECUTABLE STATEMENT CGAMR
  24. CGAMR = (0.0, 0.0)
  25. X = REAL (Z)
  26. IF (X.LE.0.0 .AND. AINT(X).EQ.X .AND. AIMAG(Z).EQ.0.0) RETURN
  27. C
  28. CALL XGETF (IROLD)
  29. CALL XSETF (1)
  30. CGAMR = CLNGAM(Z)
  31. CALL XERCLR
  32. CALL XSETF (IROLD)
  33. CGAMR = EXP (-CGAMR)
  34. C
  35. RETURN
  36. END