Mercurial > octave-nkf
view libcruft/slatec-fn/algams.f @ 14026:3781981be535 ss-3-5-90
snapshot 3.5.90
* configure.ac (AC_INIT): Version is now 3.5.90.
(OCTAVE_API_VERSION_NUMBER): Now 46.
(OCTAVE_RELEASE_DATE): Now 2011-12-11.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Sun, 11 Dec 2011 23:18:31 -0500 |
parents | 82be108cc558 |
children |
line wrap: on
line source
*DECK ALGAMS SUBROUTINE ALGAMS (X, ALGAM, SGNGAM) C***BEGIN PROLOGUE ALGAMS C***PURPOSE Compute the logarithm of the absolute value of the Gamma C function. C***LIBRARY SLATEC (FNLIB) C***CATEGORY C7A C***TYPE SINGLE PRECISION (ALGAMS-S, DLGAMS-D) C***KEYWORDS ABSOLUTE VALUE OF THE LOGARITHM OF THE GAMMA FUNCTION, C FNLIB, SPECIAL FUNCTIONS C***AUTHOR Fullerton, W., (LANL) C***DESCRIPTION C C Evaluates the logarithm of the absolute value of the gamma C function. C X - input argument C ALGAM - result C SGNGAM - is set to the sign of GAMMA(X) and will C be returned at +1.0 or -1.0. C C***REFERENCES (NONE) C***ROUTINES CALLED ALNGAM C***REVISION HISTORY (YYMMDD) C 770701 DATE WRITTEN C 890531 Changed all specific intrinsics to generic. (WRB) C 890531 REVISION DATE from Version 3.2 C 891214 Prologue converted to Version 4.0 format. (BAB) C***END PROLOGUE ALGAMS C***FIRST EXECUTABLE STATEMENT ALGAMS ALGAM = ALNGAM(X) SGNGAM = 1.0 IF (X.GT.0.0) RETURN C INT = MOD (-AINT(X), 2.0) + 0.1 IF (INT.EQ.0) SGNGAM = -1.0 C RETURN END