| Subroutine | Description | Category | Author |
|---|---|---|---|
| r1mach | Return floating point machine dependent constants. | R1 | Fox, P. A., (Bell Labs) |
| r1mpyq | Subsidiary to SNSQ and SNSQE | none | (UNKNOWN) |
| r1updt | Subsidiary to SNSQ and SNSQE | none | (UNKNOWN) |
| r9aimp | Evaluate the Airy modulus and phase. | C10D | Fullerton, W., (LANL) |
| r9atn1 | Evaluate ATAN(X) from first order relative accuracy so that ATAN(X) = X + X**3*R9ATN1(X). | C4A | Fullerton, W., (LANL) |
| r9chu | Evaluate for large Z Z**A * U(A,B,Z) where U is the logarithmic confluent hypergeometric function. | C11 | Fullerton, W., (LANL) |
| r9gmic | Compute the complementary incomplete Gamma function for A near a negative integer and for small X. | C7E | Fullerton, W., (LANL) |
| r9gmit | Compute Tricomi's incomplete Gamma function for small arguments. | C7E | Fullerton, W., (LANL) |
| r9knus | Compute Bessel functions EXP(X)*K-SUB-XNU(X) and EXP(X)* K-SUB- XNU+1(X) for 0.0 XNU 1.0. | C10B3 | Fullerton, W., (LANL) |
| r9lgic | Compute the log complementary incomplete Gamma function for large X and for A X. | C7E | Fullerton, W., (LANL) |
| r9lgit | Compute the logarithm of Tricomi's incomplete Gamma function with Perron's continued fraction for large X and A X. | C7E | Fullerton, W., (LANL) |
| r9lgmc | Compute the log Gamma correction factor so that LOG(GAMMA(X)) = LOG(SQRT(2*PI)) + (X-.5)*LOG(X) - X + R9LGMC(X). | C7E | Fullerton, W., (LANL) |
| r9ln2r | Evaluate LOG(1+X) from second order relative accuracy so that LOG(1+X) = X - X**2/2 + X**3*R9LN2R(X). | C4B | Fullerton, W., (LANL) |
| r9pak | Pack a base 2 exponent into a floating point number. | A6B | Fullerton, W., (LANL) |
| r9upak | Unpack a floating point number X so that X = Y*2**N. | A6B | Fullerton, W., (LANL) |
| radb2 | Calculate the fast Fourier transform of subvectors of length two. | none | Swarztrauber, P. N., (NCAR) |
| radb3 | Calculate the fast Fourier transform of subvectors of length three. | none | Swarztrauber, P. N., (NCAR) |
| radb4 | Calculate the fast Fourier transform of subvectors of length four. | none | Swarztrauber, P. N., (NCAR) |
| radb5 | Calculate the fast Fourier transform of subvectors of length five. | none | Swarztrauber, P. N., (NCAR) |
| radbg | Calculate the fast Fourier transform of subvectors of arbitrary length. | none | Swarztrauber, P. N., (NCAR) |
| radf2 | Calculate the fast Fourier transform of subvectors of length two. | none | Swarztrauber, P. N., (NCAR) |
| radf3 | Calculate the fast Fourier transform of subvectors of length three. | none | Swarztrauber, P. N., (NCAR) |
| radf4 | Calculate the fast Fourier transform of subvectors of length four. | none | Swarztrauber, P. N., (NCAR) |
| radf5 | Calculate the fast Fourier transform of subvectors of length five. | none | Swarztrauber, P. N., (NCAR) |
| radfg | Calculate the fast Fourier transform of subvectors of arbitrary length. | none | Swarztrauber, P. N., (NCAR) |
| rand | Generate a uniformly distributed random number. | L6A21 | Fullerton, W., (LANL) |
| ratqr | Compute the largest or smallest eigenvalues of a symmetric tridi agonal matrix using the rational QR method with Newton correc tion. | D4A5, | Smith, B. T., et al. |
| rc | Calculate an approximation to RC(X,Y) = Integral from zero to in finity of -1/2 -1 (1/2)(t+X) (t+Y) dt, where X is nonnegative and Y is positive. | C14 | Carlson, B. C. |
| rc3jj | Evaluate the 3j symbol f(L1) = ( L1 L2 L3) (-M2-M3 M2 M3) for all allowed values of L1, the other parameters being held fixed. | C19 | Gordon, R. G., Harvard University |
| rc3jm | Evaluate the 3j symbol g(M2) = (L1 L2 L3 ) (M1 M2 -M1-M2) for all allowed values of M2, the other parameters being held fixed. | C19 | Gordon, R. G., Harvard University |
| rc6j | Evaluate the 6j symbol h(L1) = {L1 L2 L3} {L4 L5 L6} for all al lowed values of L1, the other parameters being held fixed. | C19 | Gordon, R. G., Harvard University |
| rd | Compute the incomplete or complete elliptic integral of the 2nd kind. For X and Y nonnegative, X+Y and Z positive, RD(X,Y,Z) = Integral from zero to infinity of -1/2 -1/2 -3/2 (3/2)(t+X) (t+Y) (t+Z) dt. If X or Y is zero, the integral is complete. | C14 | Carlson, B. C. |
| rebak | Form the eigenvectors of a generalized symmetric eigensystem from the eigenvectors of derived matrix output from REDUC or REDUC2. | D4C4 | Smith, B. T., et al. |
| rebakb | Form the eigenvectors of a generalized symmetric eigensystem from the eigenvectors of derived matrix output from REDUC2. | D4C4 | Smith, B. T., et al. |
| reduc | Reduce a generalized symmetric eigenproblem to a standard symmet ric eigenproblem using Cholesky factorization. | D4C1C | Smith, B. T., et al. |
| reduc2 | Reduce a certain generalized symmetric eigenproblem to a standard symmetric eigenproblem using Cholesky factorization. | D4C1C | Smith, B. T., et al. |
| reort | Subsidiary to BVSUP | none | Watts, H. A., (SNLA) 910722 Updated AUTHOR section. (ALS) |
| rf | Compute the incomplete or complete elliptic integral of the 1st kind. For X, Y, and Z non-negative and at most one of them zero, RF(X,Y,Z) = Integral from zero to infinity of -1/2 -1/2 -1/2 (1/2)(t+X) (t+Y) (t+Z) dt. If X, Y or Z is zero, the integral is | C14 | Carlson, B. C. |
| rfftb | Compute the backward fast Fourier transform of a real coefficient array. | J1A1 | Swarztrauber, P. N., (NCAR) |
| rfftb1 | Compute the backward fast Fourier transform of a real coefficient array. | J1A1 | Swarztrauber, P. N., (NCAR) |
| rfftf | Compute the forward transform of a real, periodic sequence. | J1A1 | Swarztrauber, P. N., (NCAR) |
| rfftf1 | Compute the forward transform of a real, periodic sequence. | J1A1 | Swarztrauber, P. N., (NCAR) |
| rffti | Initialize a work array for RFFTF and RFFTB. | J1A1 | Swarztrauber, P. N., (NCAR) |
| rffti1 | Initialize a real and an integer work array for RFFTF1 and RFFTB1. | J1A1 | Swarztrauber, P. N., (NCAR) |
| rg | Compute the eigenvalues and, optionally, the eigenvectors of a real general matrix. | D4A2 | Smith, B. T., et al. |
| rgauss | Generate a normally distributed (Gaussian) random number. | L6A14 | Fullerton, W., (LANL) |
| rgg | Compute the eigenvalues and eigenvectors for a real generalized eigenproblem. | D4B2 | Smith, B. T., et al. |
| rj | Compute the incomplete or complete (X or Y or Z is zero) elliptic integral of the 3rd kind. For X, Y, and Z non- negative, at most one of them zero, and P positive, RJ(X,Y,Z,P) = Integral from ze ro to infinity of -1/2 -1/2 -1/2 -1 (3/2)(t+X) (t+Y) (t+Z) (t+P) | C14 | Carlson, B. C. |
| rkfab | Subsidiary to BVSUP | none | Watts, H. A., (SNLA) 910722 Updated AUTHOR section. (ALS) |
| rpqr79 | Find the zeros of a polynomial with real coefficients. | F1A1A | Vandevender, W. H., (SNLA) |
| rpzero | Find the zeros of a polynomial with real coefficients. | F1A1A | Kahaner, D. K., (NBS) |
| rs | Compute the eigenvalues and, optionally, the eigenvectors of a real symmetric matrix. | D4A1 | Smith, B. T., et al. |
| rsb | Compute the eigenvalues and, optionally, the eigenvectors of a symmetric band matrix. | D4A6 | Smith, B. T., et al. |
| rsco | Subsidiary to DEBDF | none | Watts, H. A., (SNLA) 910722 Updated AUTHOR section. (ALS) |
| rsg | Compute the eigenvalues and, optionally, the eigenvectors of a symmetric generalized eigenproblem. | D4B1 | Smith, B. T., et al. |
| rsgab | Compute the eigenvalues and, optionally, the eigenvectors of a symmetric generalized eigenproblem. | D4B1 | Smith, B. T., et al. |
| rsgba | Compute the eigenvalues and, optionally, the eigenvectors of a symmetric generalized eigenproblem. | D4B1 | Smith, B. T., et al. |
| rsp | Compute the eigenvalues and, optionally, the eigenvectors of a real symmetric matrix packed into a one dimensional array. | D4A1 | Smith, B. T., et al. |
| rst | Compute the eigenvalues and, optionally, the eigenvectors of a real symmetric tridiagonal matrix. | D4A5 | Smith, B. T., et al. |
| rt | Compute the eigenvalues and eigenvectors of a special real tridi agonal matrix. | D4A5 | Smith, B. T., et al. |
| runif | Generate a uniformly distributed random number. | L6A21 | Fullerton, W., (LANL) |
| rwupdt | Subsidiary to SNLS1 and SNLS1E | none | (UNKNOWN) |