Source Codes in
Fortran77

Over the years, I have collected, modified, adapted, adopted or created a number of software packages in FORTRAN77. You might be able to use one of these libraries, or a routine or two from a library.

Available software includes:

1. anyplt, a library which was intended as a simple interface to a variety of graphics packages.
2. args a program which reports the command line arguments of a FORTRAN77 program;
3. asa005 a library which evaluates the lower tail of the noncentral Student's T distribution; this is Applied Statistics Algorithm 5;
4. asa006 a library which computes the Cholesky factor of a positive definite symmetric matrix; this is Applied Statistics Algorithm 6;
5. asa007 a library which computes the inverse of a positive definite symmetric matrix; this is Applied Statistics Algorithm 7;
6. asa032 a library which evaluates the incomplete Gamma function; this is Applied Statistics Algorithm 32;
7. asa047 a library which minimizes a scalar function of several variables using the Nelder-Mead algorithm; this is Applied Statistics Algorithm 47;
8. asa058 a library which carries out clustering of data; this is Applied Statistics Algorithm 58;
9. asa063 a library which evaluates the incomplete Beta function; this is Applied Statistics Algorithm 63;
10. asa066 a library which computes the cumulative density function of the standard normal distribution; this is Applied Statistics Algorithm 66;
11. asa076 a library which evaluates Owen's T function; this is Applied Statistics Algorithm 76;
12. asa091, a library which computes the percentage points of the Chi-square distribution; this is Applied Statistics Algorithm 91;
13. asa103, a library which evaluates the digamma or psi function; this is Applied Statistics Algorithm 103;
14. asa109 a library which inverts the incomplete Beta function; this is Applied Statistics Algorithm 109;
15. asa111 a library which evaluates the percentage points of the normal distribution; this is Applied Statistics Algorithm 111;
16. asa113 a library which clusters data using the Banfield and Bassill transfer and swap algorithm; this is Applied Statistics Algorithm 113;
17. asa121, a library which evaluates the trigamma function; this is Applied Statistics Algorithm 121;
18. asa136 a library which clusters data using the Hartigan and Wong K-means algorithm; this is Applied Statistics Algorithm 136;
19. asa144 a library which randomly generates tables with given row and column sums; this is Applied Statistics Algorithm 144;
20. asa147 a library which evaluates the incomplete Gamma function; this is Applied Statistics Algorithm 147;
21. asa152, a library which evaluates point and cumulative probabilities associated with the hypergeometric distribution; this is Applied Statistics Algorithm 152;
22. asa159 a library which generates a table with given row and column sums; this is Applied Statistics Algorithm 159;
23. asa189, a library which estimates the parameters of a beta binomial distribution based on a sample of values generated by the distribution, by D Smith; this is Applied Statistics Algorithm 189;
24. asa205 a library which generates tables with given row and column sums; this is Applied Statistics Algorithm 205;
25. asa226 a library which evaluates the noncentral incomplete Beta function; this is Applied Statistics Algorithm 226;
26. asa239 a library which evaluates the incomplete Gamma function; this is Applied Statistics Algorithm 239;
27. asa241, a library which computes the inverse of the normal CDF; this is Applied Statistics Algorithm 241;
28. asa243 a library which evaluates the lower tail of the noncentral Student's T distribution; this is Applied Statistics Algorithm 243;
29. asa245 a library which evaluates the logarithm of the Gamma function; this is Applied Statistics Algorithm 245;
30. asa266, a library which evaluates various properties of the Dirichlet probability density function; this is Applied Statistics Algorithm 266;
31. asa299 a library which computes the lattice points (integer coordinates) in an M-dimensional simplex; this is Applied Statistics Algorithm 299;
32. asa310 a library which evaluates the noncentral incomplete Beta function; this is Applied Statistics Algorithm 310;
33. beta_nc, a library which evaluates the CDF of the noncentral Beta distribution;
34. betis a program which applies the boundary element method to solve Laplace's equation in a 2D region.
35. blas1 a library which constitutes the Level 1 Basic Linear Algebra Subprograms, (vector/vector);
36. blas1_c a library which constitutes the Level 1 Basic Linear Algebra Subprograms, (vector/vector) using single precision complex arithmetic;
37. blas1_d a library which constitutes the Level 1 Basic Linear Algebra Subprograms, (vector/vector) using double precision real arithmetic;
38. blas1_s a library which constitutes the Level 1 Basic Linear Algebra Subprograms, (vector/vector) using single precision real arithmetic;
39. blas1_z a library which constitutes the Level 1 Basic Linear Algebra Subprograms, (vector/vector) using double precision complex arithmetic;
40. blas2 a library which constitutes the Level 2 Basic Linear Algebra Subprograms, (matrix/vector);
41. blas3 a library which constitutes the Level 3 Basic Linear Algebra Subprograms, matrix/matrix);
42. brent a library which contains Richard Brent's routines for finding zeroes or minima of functions, without the use of derivative information.
43. calcomp a library which emulates a few CALCOMP plotting calls;
44. channel a program which uses a finite element model of time independent viscous incompressible flow in a 2D channel;
45. chrpak a library which can be used to analyze and manipulate character strings;
46. clawpack a library which numerically solves hyperbolic differential equations using a wave propagation approach, by Randy LeVeque;
47. clebsch_gordan a library which evaluates Clebsch-Gordan coefficients involved in angular momentum coupling in quantum mechanics;
48. cnf_io is a library which reads or writes CNF data files, used to store boolean formulas in conjunctive normal form.
49. cwg_ode a library which contains three ODE solvers by William Gear.
50. dislin examples which illustrate the use of a scientific plotting package;
51. distances, a program which reads a multiple alighment file (perhaps preprocessed by READ_ALIGN) and creates a distance matrix suitable for a principal component analysis by the SPACER program.
52. dqed a library which solves bounded and constrained least squares problems and systems of nonlinear equations;
53. driv a library which solves real or complex systems of ordinary differential equations;
54. eispack, a library which carries out eigenvalue computations; superseded by LAPACK;
55. f77 examples which illustrate features of FORTRAN77;
56. f77_intrinsics examples which illustrate the use of some of the intrinsic functions included with the F77 standard language;
57. fd1d, a program which implements a finite difference algorithm for a 1D predator-prey system;
58. fem1d a program which applies the finite element method, with piecewise linear basis functions, to a linear two point boundary value problem in 1D;
59. fem1d_adaptive a program which uses an adaptive mesh when applying the finite element method, with piecewise linear basis functions, to a linear two point boundary value problem in 1D;
60. fem1d_nonlinear a program which applies the finite element method, with piecewise linear basis functions, to a nonlinear two point boundary value problem in 1D;
61. fem1d_pmethod a program which applies the p-method version of the finite element method to a linear two point boundary value problem in 1D;
62. fem2d_poisson a program which implements the finite element method for Poisson's equation on a triangulated square in 2D;
63. fft_open_mp is a program which demonstrates the computation of a Fast Fourier Transform in parallel, using OpenMP.
64. fftpack5 a library which contains version 5 of the FFTPACK Fast Fourier Transform package;
65. fftw3 examples which illustrate the use of the FFTW3 Fast Fourier Transform package;
66. filum, a library which performs various operations on files;
67. fishpack a library which solves the Poisson, Laplace, or Helmholtz equation in 2D or 3D on a variety of geometries.
68. g77_intrinsics examples which illustrate the use of some of the "extra" intrinsic functions included with the G77 compiler;
69. geometry a library which carries out a number of 2D/3D geometric calculations;
70. gfortran_intrinsics, a program which illustrates the use of some of the intrinsic functions included by the GFORTRAN compiler, above and beyond those specified in the FORTRAN language standards.
71. gprof, examples which illustrate the use of the GPROF program performance monitor;
72. hb_to_st a program which reads a sparse matrix in a Harwell-Boeing file and writes an equivalent Sparse Triplet file.
73. hcell a program which solves the time-dependent 2D incompressible Navier Stokes equations in an H-shaped region.
74. hcell_steady a program which solves the time-independent or "steady" 2D incompressible Navier Stokes equations in an H-shaped region.
75. heat_mpi, a program which solves the 1D time dependent heat equation using MPI.
76. inout a program which solves the time-dependent 2D incompressible Navier Stokes equations in a rectangular region with inflow and outflow.
77. kmedian a program which solves the K-Median problem;
78. lamp a library which solves linear assignment and matching problems.
79. lammps a program which performs a classical molecular dynamics simulation.
80. lapack a library which solves linear systems and eigenvalue problems in complex, double precision, single precision or double precision complex arithmetic, for a wide variety of matrix types and storage formats; a more modern library than LINPACK.
81. lapack_d examples which illustrate the use of LAPACK routines for linear algebra problems involving double precision real arithmetic.
82. linpack a library which solves linear systems and eigenvalue problems in complex, double precision, single precision or double precision complex arithmetic, for a wide variety of matrix types and storage formats; superseded by LAPACK.
83. linpack_bench, a program which carries out the LINPACK benchmark;
84. linpack_c a library which solves linear systems using single precision complex arithmetic;
85. linpack_d a library which solves linear systems using double precision real arithmetic;
86. linpack_s a library which solves linear systems using single precision real arithmetic;
87. linpack_z a library which solves linear systems using double precision complex arithmetic;
88. machar a library which computes the values of various machine characteristic constants;
89. machine a library which returns tabulated values of machine constants;
90. makefiles examples which illustrate the use of MAKEFILES for maintaining a software project;
91. matmul, a program which carries out a matrix multiplication benchmark;
92. mdbnch a program which is a molecular dynamics benchmark;
93. mixed examples which illustrate the use of mixed language programming in which the main program is written in FORTRAN77;
94. mpi examples which illustrate the use of the MPI library for the parallel execution of programs;
95. mpi_condor examples which illustrate the use of CONDOR scripts for submitting MPI jobs to the FSU SCS cluster.
96. mpi_eclipse examples which illustrate the use of LoadLeveler scripts for submitting MPI jobs to the ECLIPSE and TERAGOLD IBM SP systems, formerly available at FSU.
97. mpi_stubs a library which may be used if a parallel program is to be compiled, loaded and run in a serial machine.
98. mpi_sysx examples which illustrate the use of PBS scripts for submitting MPI jobs to System X.
99. mxv, a program which compares the performance of (DO I, DO J) loops and (DO J, DO I ) loops for computing the product of an MxN matrix A and an N vector X.
100. mxv_open_mp, a program which measures the performance of the matrix multiplication problem y=A*x, with and without parallelization by OpenMP.
101. nag examples which illustrate the use of the NAG library;
102. ncarg the NCAR graphics library;
103. netode a program which models the time-dependent behavior of a two-phase fluid in an abstract network of nodes and links.
104. nintlib a library which applies approximate integration (quadrature) in M dimensions;
105. nl2sol a library which implements an adaptive nonlinear least squares algorithm;
106. normal a library which generates normally distributed pseudorandom numbers.
107. ode, a library which handles ordinary differential equations, by Shampine and Gordon;
108. odepack a library which contains nine ODE solvers, including LSODE, LSODES, LSODA, LSODAR, LSODPK, LSODKR, LSODI, LSOIBT, and LSODIS, by Alan Hindmarsh.
109. open_mp examples which illustrate the use of the use of the OpenMP interface for parallel processing on shared memory systems;
110. open_mp_sgi examples which illustrate the use of the OpenMP programs executed on one of Virginia Tech's SGI systems.
111. open_mp_stubs a library which is a dummy implementation of the OpenMP routines, and can be used to compile, load, and run a program using OpenMP on a system which does not have OpenMP installed.
112. owens a library which evaluates Owen's T function;
113. petsc examples which illustrate the use of the Portable Extensible Toolkit for Scientific Computations;
114. pitcon66 a library which seeks to produce a sequence of points that satisfy a set of nonlinear equations with one degree of freedom; this is version 6.6 of ACM TOMS algorithm 596.
115. pltmg_double a library which carries out Piecewise Linear Multigrid finite element calculations (double precision version), by Randy Bank;
116. pltmg_single a library which carries out Piecewise Linear Multigrid finite element calculations (single precision version), by Randy Bank;
117. polpak a library which evaluates special functions, recursive polynomials, and so on.
118. power_method a library which carries out the power method for finding a dominant eigenvalue and its eigenvector.
119. pppack a library which computes piecewise polynomial interpolants and approximants, in particular, cubic splines;
120. praxis, a library which implements Richard Brent's method for minimizing a scalar function of several variables, without using derivatives;
121. quad_mpi, a program which applies a quadrature rule to estimate an integral, using MPI;
122. quad_open_mp, a program which applies a quadrature rule to estimate an integral, using OpenMP;
123. quadpack a library which approximates integrals of various types over finite, semi-infinite and infinite one dimensional intervals;
124. random_mpi, a program which demonstrates one way to generate the same sequence of random numbers for both sequential execution and parallel execution under MPI.
125. ranlib a library which produces random samples from a variety of probabilistic distributions, by Brown and Lovato.
126. read_align, a program which extracts data from a multiple alignment file for sequences and writes the data to a file suitable for input to DISTANCES.
127. rkf45, a library which solves a system of ordinary differential equations by the Runge-Kutta-Fehlberg method.
128. satisfiability is a program which demonstrates, for a particular circuit, an exhaustive search for solutions of the circuit satisfiability problem.
129. satisfiability_mpi is a program which demonstrates, for a particular circuit, an exhaustive search for solutions of the circuit satisfiability problem, using MPI to perform the calculation in parallel.
130. satisfiability_open_mp is a program which demonstrates, for a particular circuit, an exhaustive search for solutions of the circuit satisfiability problem, using OpenMP for parallel execution.
131. scpack a library which performs the numerical computation of Schwarz-Christoffel conformal mappings, by Nick Trefethen.
132. select, a library which carries out Nijenhuis and Wilf's generalized combinatorial selection algorithm;
133. serba a program which applies the boundary element method to solve the elasticity equation in a 2D region.
134. sgefa_open_mp, a program which compares a standard linear algebra solver against a revised version which can be run in parallel with OpenMP.
135. simpack a library which estimates an integral over a region that is the union of multidimensional simplexes.
136. spacer, a program which reads a file containing a distance matrix and performs a principal component analysis.
137. spacer_data_convert, a program which reads a text file of distances, and reformats it for use with SPACER.
138. sparsepak, a library which forms an obsolete version of the Waterloo Sparse Matrix Package;
139. specfun a library which evaluates certain special functions.
140. starpac a library which carries out statistical data analysis.
141. subpak a library which implements a wide variety of operations on real and integer vectors and matrices;
142. subset a library which can rank, unrank, and generate random subsets, combinations, permutations, and so on;
143. table_io, a library which reads and writes a set of data as a TABLE file.
144. tcell a program which solves the time-dependent 2D incompressible Navier Stokes equations in a T-shaped region.
145. tcell_mass a program which creates and saves the mass matrix associated with a finite element solution of the time dependent 2D incompressible Navier Stokes equations in the "TCELL" region.
146. templates a library which forms basic implementations of iterative methods for solving large sparse linear systems;
147. test_eigen, a library which implements test matrices for eigenvalue analysis.
148. test_values a library which stores a small set of values for a variety of mathematical functions.
149. testpack a library which contains a set of test integrand functions by Alan Genz.
150. timer examples which illustrate how executing code may be timed.
151. timestamp a library which prints out the YMDHMS date as a timestamp.
152. toms178 a library which seeks the minimizer of a scalar function of several variables using the Hooke-Jeeves method;
     this is a FORTRAN77 version of ACM TOMS algorithm 178.
153. toms179 a library which calculates the incomplete Beta ratio;
     this is a FORTRAN77 version of ACM TOMS algorithm 179.
154. toms291 a library which approximates the logarithm of the Gamma function;
     this is a FORTRAN77 version of ACM TOMS algorithm 291.
155. toms322 a library which evaluates the normal, Student's T, F and chi-square cumulative density functions;
     this is a FORTRAN77 version of ACM TOMS algorithm 322.
156. toms332 a library which evaluates Jacobi polynomials;
     this is a FORTRAN77 version of ACM TOMS algorithm 332.
157. toms343 a library which computes the eigenvalues and eigenvectors of a general real matrix;
     this is a FORTRAN77 version of ACM TOMS algorithm 343.
158. toms344, a library which calculates the Student T distribution;
     this is a FORTRAN77 version of ACM TOMS algorithm 344.
159. toms347 a library which sorts an integer vector;
     this is a FORTRAN77 version of ACM TOMS algorithm 347.
160. toms351 a library which estimates an integral using Romberg integration;
     this is a FORTRAN77 version of ACM TOMS algorithm 351.
161. toms352 a library which computes the characteristic values and associated solutions of Mattieu's differential equation;
     this is a FORTRAN77 version of ACM TOMS algorithm 352.
162. toms353 a library which estimates an integral involving a cosine or sine factor using Filon quadrature;
     this is a FORTRAN77 version of ACM TOMS algorithm 353.
163. toms358 a library which computes the singular value decomposition of a complex matrix;
     this is a FORTRAN77 version of ACM TOMS algorithm 358.
164. toms359 a library which computes the factorial analysis of variance;
     this is a FORTRAN77 version of ACM TOMS algorithm 359.
165. toms365 a library which finds a root of an analytic complex function by the downhill method;
     this is a FORTRAN77 version of ACM TOMS algorithm 365.
166. toms370 a library which produces random numbers sampled from an arbitrary PDF;
     this is a FORTRAN77 version of ACM TOMS algorithm 370.
167. toms379 a library which approximates the integral of a function;
     this is a FORTRAN77 version of ACM TOMS algorithm 379.
168. toms384 a library which computes the eigenvalues and eigenvectors of a symmetric matrix;
     this is a FORTRAN77 version of ACM TOMS algorithm 384.
169. toms385 a library which evaluates the exponential integral function;
     this is a FORTRAN77 version of ACM TOMS algorithm 385.
170. toms392 a library which approximates the time evolution of a hyperbolic PDE;
     this is a FORTRAN77 version of ACM TOMS algorithm 392.
171. toms418 a library which implements Filon quadrature of integrals involve a factor of COS(W*X) or SIN(W*X);
     this is ACM TOMS algorithm 418.
172. toms423 a library which implements Gaussian elimination, revised versions of DECOMP and SOLVE;
     this is ACM TOMS algorithm 423.
173. toms424 a library which implements Clenshaw-Curtis integration;
     this is ACM TOMS algorithm 424.
174. toms425 a library which generates a sequence of vectors from a multivariate normal distribution with given covariance;
     this is ACM TOMS algorithm 425.
175. toms427 a library which estimates a Fourier cosine integral;
     this is ACM TOMS algorithm 427.
176. toms429 a library which reports information about the estimated location of roots of a polynomial;
     this is ACM TOMS algorithm 429.
177. toms431 a library which solves quadratic and linear programming problems;
     this is ACM TOMS algorithm 431.
178. toms434 a library which computes the probability of a given RxC contingency table;
     this is ACM TOMS algorithm 434.
179. toms435 a library which computes the modified incomplete gamma function;
     this is ACM TOMS algorithm 435.
180. toms436 a library which carries out product-type trapezoidal integration;
     this is ACM TOMS algorithm 436.
181. toms437 a library which carries out product-type Simpson's integration;
     this is ACM TOMS algorithm 437.
182. toms438 a library which carries out product-type two point Gauss-Legendre-Simpson integration;
     this is ACM TOMS algorithm 438.
183. toms439 a library which carries out product-type three point Gauss-Legendre-Simpson integration;
     this is ACM TOMS algorithm 439.
184. toms441 a library which samples the dipole distribution;
     this is ACM TOMS algorithm 441.
185. toms443 a library which evaluates Lambert's W function;
     this is ACM TOMS algorithm 443.
186. toms446 a library which manipulates Chebyshev series;
     this is ACM TOMS algorithm 446.
187. toms448 a library which counts the number of partitions of an integer into a given set of integers;
     this is ACM TOMS algorithm 448.
188. toms450 a library which minimizes a scalar function of N variables using the method of Rosenbrock;
     this is ACM TOMS algorithm 450.
189. toms451 a library which determines quantiles of the Chi-Square distribution;
     this is ACM TOMS algorithm 451.
190. toms452 a library which generates a sequence of combinations of M things out of N;
     this is ACM TOMS algorithm 452.
191. toms453 a library which determines the abscissas and weights for a Gaussian quadrature formula for Bromwich's integral;
     this is ACM TOMS algorithm 453.
192. toms454 a library which implements the complex method of constrained optimization;
     this is ACM TOMS algorithm 454.
193. toms456 a library which implements the routing algorithm;
     this is ACM TOMS algorithm 456.
194. toms460 a library which determines the optimal parameters for an alternating direction implicit iteration;
     this is ACM TOMS algorithm 460.
195. toms461 a library which applies cubic spline solutions to a class of functional differential equations;
     this is ACM TOMS algorithm 461.
196. toms462 a library which evaluates the bivariate normal distribution;
     this is ACM TOMS algorithm 462.
197. toms463 a library which chooses "nice" ranges for a graph of computer-generated data;
     this is ACM TOMS algorithm 463.
198. toms467 a library which transposes a matrix in place;
     this is ACM TOMS algorithm 467.
199. toms468 a library which carries out automatic numerical integration;
     this is ACM TOMS algorithm 468.
200. toms470 a library which solves an "almost tridiagonal" system of linear equations;
     this is ACM TOMS algorithm 470.
201. toms493 a library which seeks the roots of a real polynomial;
     this library is commonly called RPOLY;
     this is ACM TOMS algorithm 493.
202. toms526 a library which interpolates scattered bivariate data, by Hiroshi Akima
     this is ACM TOMS algorithm 526.
203. toms550 a library which evaluates certain measures of a solid polyhedron, including surface area, centroid, volume, weight, moments and products of inertia;
     this is ACM TOMS algorithm 550.
204. toms552 a library which solves an overdetermined constrained linear system;
     this is ACM TOMS algorithm 552.
205. toms571 a library which evaluates or inverts the von Mises and Fisher distributions of directions on a circle or sphere;
     this is ACM TOMS algorithm 571.
206. toms596 a library which seeks to produce a sequence of points that satisfy a set of nonlinear equations with one degree of freedom; this library is commonly called PITCON;
     this is ACM TOMS algorithm 596.
207. toms611 a library which performs unconstrained optimization of a scalar function;
     this library is commonly called UNCMIN;
     this is ACM TOMS algorithm 611.
208. toms612 a library which approximates the integral of a function over a a triangle;
     the primary routine is known as TRIEX;
     this is ACM TOMS algorithm 612.
209. toms626 a library which creates a contour plot of data associated with a triangulated set of points;
     this is ACM TOMS algorithm 626.
210. toms647 a library which evaluates the Faure, Halton, and Sobol quasirandom sequences,
     this is ACM TOMS algorithm 647;
211. toms648_nsd a library which defines non-stiff test problems for ODE solvers, using double precision arithmetic;
     this is one part of ACM TOMS algorithm 648.
212. toms648_nss a library which defines non-stiff test problems for ODE solvers, using single precision arithmetic;
     this is one part of ACM TOMS algorithm 648.
213. toms655 a library which computes the weights for interpolatory quadrature rule;
     this library is commonly called IQPACK;
     this is ACM TOMS algorithm 655.
214. toms659 a library which computes elements of the Sobol quasirandom sequence;
     this is ACM TOMS algorithm 659.
215. toms664 a library which uses Gauss elimination to solve linear systems involving large banded matrices. Disk storage is used when there is not enough memory for the working data;
     this is ACM TOMS algorithm 664.
216. toms706 a library which estimates the integral of a function over a triangulated region;
     this is ACM TOMS algorithm 706.
217. toms708 a library which evaluates the incomplete beta function;
     this is ACM TOMS algorithm 708.
218. toms726 a library which computes recursion coefficients for orthogonal polynomials, and the abscissa and weights for related quadrature rules;
     this library is commonly called ORTHPOL;
     this is ACM TOMS algorithm 726.
219. toms738 a library which evaluates Niederreiter's quasirandom sequence;
     this is ACM TOMS algorithm 738;
220. toms772 a library which carries out computational geometry on the unit sphere;
     this library is commonly called STRIPACK;
     this is ACM TOMS algorithm 772.
221. ttyplt a library which create simple typewriter graphics images of mathematical data.
222. uncmin a library which seeks the unconstrained minimization of a scalar function of several variables.
223. uniform a library which generates uniform random numbers.
224. xerror a library which constitutes the SLATEC error message package.
225. y12m a library which implements a direct solver for sparse linear systems.
226. ziggurat, a library which samples the uniform, normal or exponential distributions, using the ziggurat method.

You can go up one level to the main web page.