Mercurial > octave-nkf
diff libcruft/slatec-fn/algams.f @ 7789:82be108cc558
First attempt at single precision tyeps
* * *
corrections to qrupdate single precision routines
* * *
prefer demotion to single over promotion to double
* * *
Add single precision support to log2 function
* * *
Trivial PROJECT file update
* * *
Cache optimized hermitian/transpose methods
* * *
Add tests for tranpose/hermitian and ChangeLog entry for new transpose code
author | David Bateman <dbateman@free.fr> |
---|---|
date | Sun, 27 Apr 2008 22:34:17 +0200 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/libcruft/slatec-fn/algams.f Sun Apr 27 22:34:17 2008 +0200 @@ -0,0 +1,38 @@ +*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