gamr.f

Go to the documentation of this file.

*DECK GAMR
      FUNCTION gamr (X)
C***BEGIN PROLOGUE  GAMR
C***PURPOSE  Compute the reciprocal of the Gamma function.
C***LIBRARY   SLATEC (FNLIB)
C***CATEGORY  C7A
C***TYPE      SINGLE PRECISION (GAMR-S, DGAMR-D, CGAMR-C)
C***KEYWORDS  FNLIB, RECIPROCAL GAMMA FUNCTION, SPECIAL FUNCTIONS
C***AUTHOR  Fullerton, W., (LANL)
C***DESCRIPTION
C
C GAMR is a single precision function that evaluates the reciprocal
C of the gamma function for single precision argument X.
C
C***REFERENCES  (NONE)
C***ROUTINES CALLED  ALGAMS, GAMMA, XERCLR, XGETF, XSETF
C***REVISION HISTORY  (YYMMDD)
C   770701  DATE WRITTEN
C   861211  REVISION DATE from Version 3.2
C   891214  Prologue converted to Version 4.0 format.  (BAB)
C   900727  Added EXTERNAL statement.  (WRB)
C***END PROLOGUE  GAMR
      EXTERNAL gamma
C***FIRST EXECUTABLE STATEMENT  GAMR
      gamr = 0.0
      IF (x.LE.0.0 .AND. aint(x).EQ.x) RETURN
C
      CALL xgetf (irold)
      CALL xsetf (1)
      IF (abs(x).GT.10.0) GO TO 10
      gamr = 1.0/gamma(x)
      CALL xerclr
      CALL xsetf (irold)
      RETURN
C
 10   CALL algams (x, alngx, sgngx)
      CALL xerclr
      CALL xsetf (irold)
      gamr = sgngx * exp(-alngx)
      RETURN
C
      END