C++ Source Codes

anagram, a program which determines anagrams of a string, by James Cherry;
analemma, a program which evaluates the equation of time, a formula for the difference between the uniform 24 hour day and the actual position of the sun, creating data files that can be plotted with gnuplot, based on a C program by Brian Tung.
ann, a library which computes Approximate Nearest Neighbors, by David Mount and Sunil Arya;
ann_test, a program which uses the ann library to approximate the nearest neighbors of a set of points stored in a file;
ann_to_fig, a program which plots the search tree used by the ann library, while seeking the approximate nearest neighbors of a set of points;
args, a program which reports the command line arguments with which it was invoked;
asa005, a library which evaluates the lower tail of the noncentral Student's T distribution, by BE Cooper. This is a C++ version of Applied Statistics Algorithm 5;
asa006, a library which computes the Cholesky factor of a positive definite symmetric matrix, by Michael Healy. This is a C++ version of Applied Statistics Algorithm 6;
asa007, a library which computes the inverse of a positive definite symmetric matrix, by Michael Healy. This is a C++ version of Applied Statistics Algorithm 7.
asa032, a library which evaluates the incomplete Gamma function, by G Bhattacharjee. This is a C++ version of Applied Statistics Algorithm 32;
asa047, a library which minimizes a scalar function of several variables using the Nelder-Mead algorithm, by R ONeill. This is a C++ version of Applied Statistics Algorithm 47;
asa058, a library which carries out K-means clustering of data, by David Sparks. This is a C++ version of Applied Statistics Algorithm 58;
asa063, a library which evaluates the incomplete Beta function, by KL Majumder and G Bhattacharjee. This is a C++ version of Applied Statistics Algorithm 63;
asa066, a library which computes the cumulative density function (CDF) of the standard normal distribution, by David Hill. This is a C++ version of Applied Statistics Algorithm 66;
asa076, a library which evaluates Owen's T function, by Young and Minder. This is a C++ version of Applied Statistics Algorithm 76;
asa091, a library which computes the percentage points of the Chi-square distribution, by Best and Roberts. This is a C++ version of Applied Statistics Algorithm 91;
asa103, a library which evaluates the digamma or psi function, by Jose Bernardo. This is a C++ version of Applied Statistics Algorithm 103;
asa109, a library which inverts the incomplete Beta function, by Cran, Martin and Thomas. This is a C++ version of Applied Statistics Algorithm 109;
asa111, a library which evaluates the percentage points of the normal distribution, by Beasley and Springer. This is a C++ version of Applied Statistics Algorithm 111;
asa113, a library which clusters data using a transfer and swap algorithm, by Banfield and Bassill. This is a C++ version of Applied Statistics Algorithm 113;
asa121, a library which evaluates the trigamma function, by BE Schneider. This is a C++ version of Applied Statistics Algorithm 121;
asa136, a library which implements a clustering algorithm, by Hartigan and Wong. This is a C++ version of Applied Statistics Algorithm 136;
asa144, a library which randomly generates tables with given row and column sums, by James Boyett. This is a C++ version of Applied Statistics Algorithm 144;
asa147, a library which evaluates the incomplete Gamma function, by Chi Leung Lau. This is a C++ version of Applied Statistics Algorithm 147;
asa152, a library which evaluates the probability density function (PDF) and cumulative density function (CDF) associated with the hypergeometric distribution, by Richard Lund. This is a C++ version of Applied Statistics Algorithm 152;
asa159, a library which computes a random table with given row and column sums, by Michael Patefield. This is a C++ version of Applied Statistics Algorithm 159;
asa172, a library which generates all M-dimensional indices in a given range, simulating the behavior of an arbitrary number of nested loops, by OFlaherty and MacKenzie. This is a C++ version of Applied Statistics Algorithm 172;
asa183, a library which implements a random number generator (RNG), by Wichman and Hill. This is a C++ version of Applied Statistics Algorithm 183;
asa226, a library which evaluates the noncentral incomplete Beta function, by Russell Lenth. this is a C++ version of Applied Statistics Algorithm 226;
asa239, a library which evaluates the incomplete Gamma function, by Shea. this is a C++ version of Applied Statistics Algorithm 239;
asa241, a library which computes the inverse of the normal cumulative density function (CDF), by Michael Wichura; this is a C++ version of Applied Statistics Algorithm 241;
asa243, a library which evaluates the lower tail of the noncentral Student's T distribution, by Russell Lenth; this is a C++ version of Applied Statistics Algorithm 243;
asa245, a library which evaluates the logarithm of the Gamma function, by Allan Mcleod; this is a C++ version of Applied Statistics Algorithm 245;
asa299, a library which computes the lattice points (integer coordinates) in an M-dimensional simplex, by Chasalow and Brand; this is a C++ version of Applied Statistics Algorithm 299;
asa310, a library which evaluates the noncentral incomplete Beta function, by Chattamvelli and Shanmugam; this is a C++ version of Applied Statistics Algorithm 310;
ball_grid, a library which computes grid points inside a 3D ball.
ball_volume_monte_carlo, a program which applies a Monte Carlo procedure to estimate the volume of the M-dimensional unit ball;
barycentric_interp_1d, a library which defines and evaluates the barycentric Lagrange polynomial p(x) which interpolates a set of data, so that p(x(i)) = y(i). The barycentric approach means that very high degree polynomials can safely be used.
bernstein_polynomial, a library which evaluates the Bernstein polynomials;
beta_nc, a library which evaluates the cumulative density function (CDF) of the noncentral Beta distribution;
bins, a library which handles geometric data partially sorted into bins;
bisection_integer, a library which seeks an integer solution to the equation F(X)=0, using bisection within a user-supplied change of sign interval [A,B].
bisection_rc, a library which seeks a solution to the equation F(X)=0 using bisection within a user-supplied change of sign interval [A,B]. The procedure is written using reverse communication.
bits_to_ui, a program which reads a text file of bits and writes a binary file in which groups of 32 bits have been interpreted as an unsigned integer;
black_scholes, a library which implements some simple approaches to the Black-Scholes option valuation theory, by Desmond Higham.
blas1, a library which contains the Basic Linear Algebra Subprograms (BLAS) for level 1 (vector-vector operations), for single and double precision, and for real and complex arithmetic, by Charles Lawson, Richard Hanson, David Kincaid, Fred Krogh.
blas1_c, a library which constitutes the Level 1 Basic Linear Algebra Subprograms (BLAS), using single precision complex arithmetic, by Charles Lawson, Richard Hanson, David Kincaid, Fred Krogh.
blas1_d, a library which constitutes the Level 1 Basic Linear Algebra Subprograms (BLAS), using double precision real arithmetic, by Charles Lawson, Richard Hanson, David Kincaid, Fred Krogh.
blas1_s, a library which constitutes the Level 1 Basic Linear Algebra Subprograms (BLAS), using single precision real arithmetic, by Charles Lawson, Richard Hanson, David Kincaid, Fred Krogh.
blas1_z, a library which constitutes the Level 1 Basic Linear Algebra Subprograms (BLAS), using double precision complex arithmetic, by Charles Lawson, Richard Hanson, David Kincaid, Fred Krogh.
blend, a library which uses transfinite interpolation to fill in data within a square or a cube based on boundary values;
bmp_io, a library which reads and writes graphics information stored in a Microsoft BMP file;
bmp_to_ppma, a program which converts a Microsoft BMP file to an ASCII Portable Pixel Map (PPM) file;
bmp_to_ppmb, a program which converts a Microsoft BMP file to a binary Portable Pixel Map (PPM) file;
boost, programs which illustrate the use of some of the Boost libraries;
box_behnken, a library which computes a Box-Behnken design, that is, a set of arguments to sample the behavior of a function of multiple parameters;
brent, a library which contains routines for finding zeroes or minima of a scalar function of a scalar variable, without the use of derivative information, including a reverse communication option, by Richard Brent. The C++ version includes some customizations by John Denker.
brent_old, a library which contains routines for finding zeroes or minima of a scalar function of a scalar variable, without the use of derivative information, including a reverse communication option, by Richard Brent.
brownian_motion_simulation, a library which simulates Brownian motion in an M-dimensional region, creating graphics files for processing by gnuplot.
burgers_solution, a library which evaluates an exact solution of the time-dependent 1D viscous Burgers equation.
c++_calls_f77, programs which illustrate a C++ program calling a FORTRAN77 subroutine.
c++_calls_f90, programs which illustrate a C++ program calling a FORTRAN90 subroutine.
c4lib, a library which implements a basic level of single precision complex (C4) arithmetic.
c8_complex_lib, a C++ class for double precision complex (C8) numbers;
c8lib, a library which implements certain elementary functions for double precision complex (C8) variables;
calendar_rd, a program which computes the representation of a given date in a number of calendrical systems, by Edward Reingold and Nachum Dershowitz.
calpak, a library which makes various calendar calculations;
catalog, a program which extracts routine definition lines from a marked-up C, C++, FORTRAN77 or FORTRAN90 file;
caustic_opengl, a program which computes a caustic curve and displays it using OpenGL.
ccn_rule, a program which defines a nested Clenshaw Curtis quadrature rule.
ccvt_box, a program which carries out a Centroidal Voronoi Tessellation (CVT) calculation constrained to a box region, with points projected to the boundary;
cdflib, a library which evaluates the cumulative density function (CDF), inverse CDF, and certain other inverse functions, for distributions including beta, binomial, chi-square, noncentral chi-square, F, noncentral F, gamma, negative binomial, normal, Poisson, and students T, by Barry Brown, James Lovato, Kathy Russell.
cell, a library which defines a cell array, a generalization of an array which can compactly store and retrieve vector or matrix data of varying size, such as the rows of a triangular matrix.
cellular_automaton, a program which demonstrates the 1D Cellular Automaton known as rule #30.
cg_rc, a library which implements the conjugate gradient method for solving a positive definite sparse linear system A*x=b, using reverse communication.
chebyshev, a library which computes the Chebyshev interpolant/approximant to a given function over an interval.
chebyshev_interp_1d, a library which determines the combination of Chebyshev polynomials which interpolates a set of data, so that p(x(i)) = y(i).
chebyshev_polynomial, a library which evaluates the Chebyshev polynomial and associated functions.
chebyshev1_rule, a program which computes and writes out a Gauss-Chebyshev type 1 quadrature rule of given order.
chebyshev2_rule, a program which computes and writes out a Gauss-Chebyshev type 2 quadrature rule of given order.
chrpak, a library which can manipulate characters and strings;
circle_arc_grid, a program which computes grid points along a circular arc.
circle_grid, a library which computes grid points inside a circle.
cities, a library which works with city-to-city distance matrices;
clenshaw_curtis_rule, a program which implements a Clenshaw Curtis quadrature rule;
cnf_io, a library which reads or writes Conjunctive Normal Form (CNF) data files, used to store boolean formulas.
collatz_recursive, a library which demonstrates recursive programming by considering the simple Collatz 3n+1 problem.
colored_noise, a library which generates samples of noise obeying a 1/f^alpha power law.
combination_lock, a program which simulates the process of determining the secret combination of a lock.
combo, a library which includes routines for ranking, unranking, enumerating and randomly selecting balanced sequences, cycles, graphs, Gray codes, subsets, partitions, permutations, restricted growth functions, Pruefer codes and trees.
communicator_mpi, a program which creates new communicators involving a subset of initial set of MPI processes in the default communicator MPI_COMM_WORLD.
compass_search, a library which seeks the minimizer of a scalar function of several variables using compass search, a direct search algorithm that does not use derivatives.
complex_numbers, a program which demonstrates some simple features involved in the use of complex numbers in C programming.
condition, a library which implements methods of computing or estimating the condition number of a matrix.
cordic, a library which computes a few special functions using the CORDIC algorithm.
correlation, a library which contains examples of statistical correlation functions.
cpp, programs which illustrate some features of the C++ language;
cpp_arrays, programs which illustrate the use of vectors, matrices and tensors.
cpp_files, programs which illustrate the use of files.
cpp_intrinsics, a program which illustrates the use of intrinsic functions supplied by C++.
cpp_random, programs which illustrate the use of the C++ random number generator (RNG) functions.
cr2crlf, a program which converts carriage returns to carriage returns + linefeeds in a file;
cr2lf, a program which converts carriage returns to linefeeds in a file;
crc, a library which computes the cyclic redundancy checksum (CRC) as used by the PNG data format.
crrm, a program which removes carriage returns from a file;
cvt, a library which computes an N point CVT Centroidal Voronoi Tessellation (CVT) in M dimensions;
cvt_dataset, a program which creates a Centroidal Voronoi Tessellation (CVT) dataset;
cycle_brent, a library which carries out an iterated function evaluation to determine the nearest element of a cycle, and the cycle's length, using Brent's method.
cycle_floyd, a library which carries out an iterated function evaluation, and seeks to determine the nearest element of a cycle, and the cycle's length, using Floyd's method.
cyclic_reduction, a library which solves a tridiagonal linear system using cyclic reduction;
deblank, a program which removes all blank lines from a file;
decomment, a program which removes from a file every line that begins with "#" (a common convention for indicating comments);
delaunay_tree_2d, a program which computes the Delaunay triangulation of a 2D dataset, by Olivier Devillers.
detroff, a program which removes from a file every occurrence of the troff overprint pair character+backspace, which can make a MAN page printable;
diaphony, a program which reads a file of N points in M dimensions and computes its diaphony, a measure of point dispersion.
dijkstra, a program which runs a simple example of Dijkstra's minimum distance algorithm for graphs.
dijkstra_openmp, a program which uses OpenMP to parallelize a simple example of Dijkstra's minimum distance algorithm for graphs.
discrete_pdf_sample_2d, a program which demonstrates how to construct a Probability Density Function (PDF) from sample data over a 2D domain, and then to use that PDF to create new samples.
dislin, programs which illusrate the use of DISLIN, a scientific plotting package;
divdif, a library which uses divided differences to compute the polynomial interpolant to a given set of data.
doomsday, a library which is given the year, month and day of a date, and uses John Conway's doomsday algorithm to determine the corresponding day of the week.
dream, a library which implements the DREAM algorithm for accelerating Markov Chain Monte Carlo (MCMC) convergence using differential evolution, by Guannan Zhang.
duel_simulation, a program which simulates N repetitions of a duel between two players, each of whom has a known firing accuracy.
dunavant, a library which defines Dunavant quadrature rules in a triangle;
eispack, a library which carries out eigenvalue computations. It includes a function to compute the singular value decomposition (SVD) of a rectangular matrix. This library was superseded by LAPACK;
ellipse_grid, a library which computes grid points inside an ellipse.
ellipsoid_grid, a library which computes grid points inside a 3D ellipsoid.
emps, a program which expands a compressed MPS linear programming file, by David Gay.
faure, a library which computes elements of the Faure M-dimensional quasirandom sequence;
faure_dataset, a program which creates a Faure quasirandom dataset;
fd_predator_prey, a program which solves a pair of predator prey ordinary differential equations (ODE's) using the finite difference method (FDM).
fd1d_advection_ftcs, a program which applies the finite difference method (FDM) to solve the time-dependent advection equation ut = - c * ux in one spatial dimension, with a constant velocity, using the FTCS method, forward time difference, centered space difference, writing graphics files for processing by gnuplot.
fd1d_advection_lax, a program which applies the finite difference method (FDM) to solve the time-dependent advection equation ut = - c * ux in one spatial dimension, with a constant velocity, using the Lax method to approximate the time derivative, writing graphics files for processing by gnuplot.
fd1d_advection_lax_wendroff, a program which applies the finite difference method (FDM) to solve the time-dependent advection equation ut = - c * ux in one spatial dimension, with a constant velocity, using the Lax-Wendroff method to approximate the time derivative, writing graphics files for processing by gnuplot.
fd1d_burgers_lax, a program which applies the finite difference method (FDM) and the Lax Wendroff method to solve the non-viscous time-dependent Burgers equation in one spatial dimension.
fd1d_burgers_leap, a program which applies the finite difference method (FDM) and the leapfrog approach to solve the non-viscous time-dependent Burgers equation in one spatial dimension.
fd1d_bvp, a program which applies the finite difference method (FDM) to a two point boundary value problem in one spatial dimension.
fd1d_heat_explicit, a library which uses the finite difference method (FDM), with explicit time stepping, to solve the time dependent 1D heat equation;
fd1d_heat_implicit, a program which implements a finite difference method (FDM), implicit in time, of the time dependent 1D heat equation;
fd1d_heat_steady, a program which implements a finite difference method (FDM) for the steady (time independent) 1D heat equation;
fd1d_wave, a program which applies the finite difference method (FDM) to solve the time-dependent wave equation utt = c * uxx in one spatial dimension.
fair_dice_simulation, a program which simulates N tosses of 2 dice, making a histogram of the results.
fekete, a library which sets up one of seven Fekete rules for interpolation or quadrature in a triangle;
felippa, a library which defines quadrature rules for lines, triangles, quadrilaterals, pyramids, wedges, tetrahedrons and hexahedrons.
fem_basis, a library which can define and evaluate finite element method (FEM) basis functions for any degree in an M-dimensional simplex (1D interval, 2D triangle, 3D tetrahedron, and higher dimensional generalizations.)
fem_io, a library which reads or writes FEM files (three text files describing a finite element model);
fem_to_mesh, a program which reads a pair of FEM files defining node coordinates and elements, and creates a corresponding MESH file.
fem1d, a program which applies the finite element method (FEM), using piecewise linear basis functions, to a linear two point boundary value problem in 1D;
fem1d_adaptive, a program which uses an adaptive mesh when applying the finite element method (FEM), with piecewise linear basis functions, to a linear two point boundary value problem in 1D;
fem1d_bvp_linear, a program which applies the finite element method (FEM), with piecewise linear elements, to a two point boundary value problem in one spatial dimension, and compares the computed and exact solutions with the L2 and seminorm errors.
fem1d_heat_steady, a program which implements the finite element method (FEM) for the steady (time independent) 1D heat equation;
fem1d_nonlinear, a program which applies the finite element method (FEM), with piecewise linear basis functions, to a nonlinear two point boundary value problem in 1D;
fem1d_pack, a library which contains utilities for 1D finite element method (FEM) calculations.
fem1d_pmethod, a program which applies the p-method version of the finite element method (FEM) to a linear two point boundary value problem in 1D;
fem1d_project, a program which projects data into a finite element space, including the least squares approximation of data, or the projection of a finite element solution from one mesh to another.
fem1d_sample, a program which samples a scalar or vector finite element function of one variable, defined by FEM files, returning interpolated values at the sample points.
fem2d_heat, a program which applies the finite element method (FEM) to solve the time-dependent heat equation in an arbitrary triangulated region in 2D;
fem2d_heat_rectangle, a program which applies the finite element method (FEM) to solve the time dependent heat equation on a square in 2D;
fem2d_heat_square, a library which defines the geometry of a square region, as well as boundary and initial conditions for a given heat problem, and is called by fem2d_heat as part of a solution procedure.
fem2d_navier_stokes, a program which applies the finite element method (FEM) to solve the steady incompressible Navier Stokes equations in an arbitrary triangulated region in 2D;
fem2d_navier_stokes_cavity, a library which contains the user-supplied routines necessary to run fem2d_navier_stokes on the cavity problem.
fem2d_navier_stokes_channel, a library which contains the user-supplied routines necessary to run fem2d_navier_stokes on the channel problem.
fem2d_navier_stokes_inout, a library which contains the user-supplied routines necessary to run fem2d_navier_stokes on the inout problem.
fem2d_pack, a library which implements various 2D finite element method (FEM) calculations;
fem2d_poisson, a program which applies the finite element method (FEM) to solve Poisson's equation in an arbitrary triangulated region in 2D;
fem2d_poisson_cg, a program which solves Poisson's equation on a triangulated region, using the finite element method (FEM), sparse storage, and a conjugate gradient solver.
fem2d_poisson_cg_baffle, a library which defines the geometry of a channel with 13 hexagonal baffles, as well as boundary conditions for a given Poisson problem, and is called by fem2d_poisson_cg as part of a solution procedure.
fem2d_poisson_cg_ell, a library which defines the geometry of an L-shaped region, as well as boundary conditions for a given Poisson problem, and is called by fem2d_poisson_cg as part of a solution procedure.
fem2d_poisson_cg_lake, a library which defines the geometry of a lake-shaped region, as well as boundary conditions for a given Poisson problem, and is called by fem2d_poisson_cg as part of a solution procedure.
fem2d_poisson_ell, a library which defines the geometry of an L-shaped region, as well as boundary conditions for a given Poisson problem, and is called by fem2d_poisson as part of a solution procedure.
fem2d_poisson_lake, a library which defines the geometry of a lake-shaped region, as well as boundary conditions for a given Poisson problem, and is called by fem2d_poisson as part of a solution procedure.
fem2d_poisson_rectangle, a program which solves the 2D Poisson equation on a rectangle, using the finite element method (FEM), and piecewise quadratic triangular elements.
fem2d_poisson_rectangle_linear, a program which solves the 2D Poisson equation on a rectangle, using the finite element method (FEM), and piecewise linear triangular elements.
fem2d_poisson_sparse, a program which uses the finite element method to solve Poisson's equation on an arbitrary triangulated region in 2D; (This is a version of fem2d_poisson which replaces the banded storage and direct solver by a sparse storage format and an iterative solver.)
fem2d_poisson_sparse_baffle, a library which defines the geometry of a rectangular channel with 13 baffles, as well as boundary conditions for a given Poisson problem, and is called by fem2d_poisson_sparse as part of a solution procedure.
fem2d_poisson_sparse_ell, a library which defines the geometry of an L-shaped region, as well as boundary conditions for a given Poisson problem, and is called by fem2d_poisson_sparse as part of a solution procedure.
fem2d_poisson_sparse_lake, a library which defines the geometry of a lake-shaped region, as well as boundary conditions for a given Poisson problem, and is called by fem2d_poisson_sparse as part of a solution procedure.
fem2d_project, a program which projects a function f(x,y), given as data, into a given finite element space of piecewise linear triangular elements.
fem2d_project_function, a program which projects a function f(x,y), given as a formula, into a given finite element space of piecewise linear triangular elements.
fem2d_sample, a program which samples a finite element function, defined by FEM files, (three text files describing the nodes, triangles, and coefficients); at arbitrary points.
fem2d_stokes, a program which applies the finite element method (FEM) to solve the steady incompressible Stokes's equations in an arbitrary triangulated region in 2D;
fem2d_stokes_cavity, a library which contains the user-supplied routines necessary to run fem2d_stokes on the cavity problem.
fem2d_stokes_channel, a library which contains the user-supplied routines necessary to run fem2d_stokes on the channel problem.
fem2d_stokes_inout, a library which contains the user-supplied routines necessary to run fem2d_stokes on the inout problem.
fem3d_pack, a library which contains utilities for 3D finite element method (FEM) calculations.
fem3d_project, a program which projects a function f(x,y,z), given as a data, into a given space of piecewise linear tetrahedrons associated with the finite element method (FEM).
fem3d_sample, a program library which evaluates a finite element function defined on a 3D tetrahedral mesh.
fern_opengl, a program which uses OpenGL to display the Barnsley fractal fern.
feynman_kac_1d, a program which demonstrates the use of the Feynman-Kac algorithm to solve Poisson's equation in a 1D interval by averaging stochastic paths to the boundary.
feynman_kac_2d, a program which demonstrates the use of the Feynman-Kac algorithm to solve Poisson's equation in a 2D ellipse by averaging stochastic paths to the boundary.
feynman_kac_3d, a program which demonstrates the use of the Feynman-Kac algorithm to solve Poisson's equation in a 3D ellipsoid by averaging stochastic paths to the boundary.
fft_openmp, a program which demonstrates the computation of a Fast Fourier Transform (FFT) in parallel, using OpenMP.
fft_serial, a program which computes a Fast Fourier Transform (FFT), and is intended as a starting point for implementing an OpenMP parallel version.
file_name_sequence, a program which demonstrates ways to generate a sequence of filenames, which can be useful when generating a sequence of still snapshots to be animated later.
files_multiple, a program which demonstrates how a program can open multiple output files at one time, and write data to any one specific file it chooses.
filum, a library which performs various operations with files;
floyd, a library which implements Floyd's algorithm for finding the shortest distance between pairs of nodes on a directed graph.
fn, a library which evaluates elementary and special functions using Chebyshev polynomials; functions include Airy, Bessel I, J, K and Y, beta, confluent hypergeometric, error, gamma, log gamma, Pochhammer, Spence; integrals include hyperbolic cosine, cosine, Dawson, exponential, logarithmic, hyperbolic sine, sine; by Wayne Fullerton.
fsu_dakota, a library which implements a variety of methods for the generation of quasirandom numbers, submitted for inclusion in the Sandia Labs DAKOTA package;
g++_intrinsics, a program which tests or demonstrates some of the intrinsic functions provided by the Gnu G++ compiler for C++.
g++_quadmath, a program which illustrates the use of quadruple precision real arithmetic, as provided on some systems by the Gnu g++ compiler for C++.
gegenbauer_rule, a program which computes and writes out a Gauss-Gegenbauer quadrature rule of given order.
gen_hermite_rule, a program which computes and writes out a generalized Gauss-Hermite quadrature rule of given order and parameter value.
gen_laguerre_rule, a program which computes and writes out a generalized Gauss-Laguerre quadrature rule of given order and parameter value.
geometry, a library which performs geometric calculations in 2, 3 and M dimensional space, including the computation of angles, areas, containment, distances, intersections, lengths, and volumes.
geompack, a library which computes the Delaunay triangulation of 2D data, by Barry Joe;
gm_rule, a library which defines Grundmann-Moeller quadratures rules for an M-dimensional simplex.
gnuplot, programs which illustrate how a program can write data and command files so that gnuplot can create plots of the program results.
gnuplot_i_examples, programs which demonstrate the use of the gnuplot_i library for interactive runtime gnuplot graphics.
gprof, programs which illustrate the use of the gprof program performance monitor;
graphics_examples, programs which illustrate how various kinds of data can be displayed and analyzed graphically.
grf_display_opengl, a program which reads a GRF file defining a mathematical graph and displays it in the OpenGL graphics window.
grf_io, a library which reads or writes a GRF file;
grid, a library which generates N random points on an M-dimensional uniform grid;
grid_dataset, a program which creates a file of grid values;
grid_to_bmp, a program which reads a text file of data on a rectangular grid and creates a BMP file containing a color image of the data.
gsl, programs which illustrate the use of the Gnu Scientific Library (GSL);
haar, a library which computes the Haar transform of data.
halton, a library which computes elements of the Halton M-dimensional quasirandom sequence;
halton_dataset, a program which creates a Halton quasirandom sequence;
hammersley, a library which computes elements of the Hammersley M-dimensional quasirandom sequence;
hammersley_dataset, a program which creates a Hammersley quasirandom sequence;
hb_io, a library which reads and writes matrices in the Harwell Boeing (HB) sparse matrix format;
heat_mpi, a program which demonstrates the use of MPI, by solving the 1D time dependent heat equation.
heated_plate, a program which solves the steady state heat equation in a 2D rectangular region, and is intended as a starting point for implementing an OpenMP parallel version.
heated_plate_openmp, a program which solves the steady (time independent) heat equation in a 2D rectangular region, using OpenMP to run in parallel.
hello, a program which prints out "Hello, world!".
hello_mpi, a program which prints out "Hello, world!" using MPI for parallel programming.
hello_openmp, a program which prints out "Hello, world!" within the OpenMP parallel programming environment.
hermite, a library which computes the Hermite interpolant, a polynomial that matches function values and derivatives.
hermite_cubic, a library which can compute the value, derivatives or integral of a Hermite cubic polynomial, or manipulate an interpolating function made up of piecewise Hermite cubic polynomials.
hermite_polynomial, a library which evaluates the physicist's Hermite polynomial, the probabilist's Hermite polynomial, the Hermite function, and related functions.
hermite_rule, a program which computes and writes out a Gauss-Hermite quadrature rule of given order.
hex_grid, a library which creates a hexagonal grid of points in the unit square or an arbitrary rectangle;
hex_grid_dataset, a program which creates a hexagonal grid of points, using the hex_grid library;
hexdump, a program which produces a hexadecimal dump of a file;
htmlindex, a program which creates a skeleton HTML page for a marked-up C, C++ or FORTRAN file;
i2_binary_to_ascii, a program which converts a file of short int (I2) data from binary to ASCII format.
i4lib, a library which contains many utility routines, using single precision integer (I4) arithmetic.
i8lib, a library which contains many utility routines, using double precision integer (I8) arithmetic.
ice_io, a library which reads or writes ICE grid files stored in netcdf format.
ice_to_mesh, a program which reads a netcdf file containing an ICE dataset and rewrites the information as a MESH file.
ihs, a library which computes elements of the Improved Hypercube Sampling (IHS) M-dimensional quasirandom sequence;
ihs_dataset, a program which creates an Improved Hypercube Sampling (IHS) dataset;
image_components, a library which seeks the connected nonzero or nonblack components of an image or integer array.
image_denoise, a library which applies simple filtering operations to a noisy image.
image_edge, a library which demonstrates a simple procedure for edge detection in images.
index, a library which converts a M-dimensional vector index to a one-dimensional vector index; it can handle zero and one based indexing schemes, as well as column major and row major conventions.
int_exactness, a program which tests the polynomial exactness of a 1D quadrature rule for a finite interval;
int_exactness_chebyshev1, a program which tests the polynomial exactness of Gauss-Chebyshev type 1 quadrature rules.
int_exactness_chebyshev2, a program which tests the polynomial exactness of Gauss-Chebyshev type 2 quadrature rules.
int_exactness_gegenbauer, a program which tests the polynomial exactness of Gauss-Gegenbauer quadrature rules.
int_exactness_gen_hermite, a program which tests the polynomial exactness of generalized Gauss-Hermite quadrature rules.
int_exactness_gen_laguerre, a program which tests the polynomial exactness of generalized Gauss-Laguerre quadrature rules.
int_exactness_hermite, a program which tests the polynomial exactness of Gauss-Hermite quadrature rules.
int_exactness_jacobi, a program which tests the polynomial exactness of Gauss-Jacobi quadrature rules.
int_exactness_laguerre, a program which tests the polynomial exactness of Gauss-Laguerre quadrature rules.
int_exactness_legendre, a program which tests the polynomial exactness of Gauss-Legendre quadrature rules.
ising_2d_simulation, a program which simulates the evolution of a 2D array of positive and negative charges, each of which is likely to flip to be in agreement with neighbors.
ivcon, a program which converts between various 3D graphic file formats;
jacobi, a program which implements the Jacobi iteration for linear systems.
jacobi_polynomial, a library which evaluates the Jacobi polynomial and associated functions.
jacobi_rule, a program which computes and writes out a Gauss-Jacobi quadrature rule of given order.
keast, a library which defines Keast quadrature rules in a tetrahedron;
kmeans, a library which handles the K-Means problem, which organizes a set of N points in M dimensions into K clusters;
kronrod, a library which can compute a Gauss and Gauss-Kronrod pair of quadrature rules of arbitrary order, by Robert Piessens, Maria Branders.
lagrange_approx_1d, a library which defines and evaluates the Lagrange polynomial p(x) of degree m which approximates a set of nd data points (x(i),y(i)).
lagrange_interp_1d, a library which defines and evaluates the Lagrange polynomial p(x) which interpolates a set of data depending on a 1D argument, so that p(x(i)) = y(i).
lagrange_interp_2d, a library which defines and evaluates the Lagrange polynomial p(x,y) which interpolates a set of data depending on a 2D argument that was evaluated on a product grid, so that p(x(i),y(j)) = z(i,j).
lagrange_interp_nd, a library which defines and evaluates the Lagrange polynomial p(x) which interpolates a set of data depending on a M-dimensional argument that was evaluated on a product grid, so that p(x(i)) = z(i).
laguerre_polynomial, a library which evaluates the Laguerre polynomial, the generalized Laguerre polynomials, and the Laguerre function.
laguerre_rule, a program which computes and writes out a Gauss-Laguerre quadrature rule of given order.
latin_center, a library which computes Latin Center Squares of N points in M dimensions;
latin_center_dataset, a program which creates a Latin Center Hypercube dataset;
latin_cover, a library which produces N Latin squares which cover an NxN square, or NxN Latin cubes which cover an NxNxN cube.
latin_edge, a library which computes Latin Edge Squares of N points in M dimensions;
latin_edge_dataset, a program which creates a Latin Edge Hypercube dataset;
latin_random, a library which computes Latin Random Squares of N points in M dimensions;
latin_random_dataset, a program which creates a Latin Random Hypercube dataset;
latinize, a library which adjusts N points in M dimensions to form a Latin Hypercube;
lattice_rule, a library which returns lattice rules for M-dimensional integration;
lcvt, a library which creates a Latinized Centroidal Voronoi Tessellation (CVT);
lcvt_dataset, a program which creates a Latinized Centroidal Voronoi Tessellation (CVT);
legendre_polynomial, a library which evaluates the Legendre polynomial and associated functions.
legendre_rule, a program which writes out a Gauss-Legendre quadrature rule of given order.
legendre_rule_fast, a program which uses a fast (order N) algorithm to compute a Gauss-Legendre quadrature rule of given order.
lf2cr, a program which converts linefeeds to carriage returns in a file;
lf2crlf, a program which converts linefeeds to carriage return + linefeeds in a file;
lfrm, a program which removes all linefeed characters from a file;
life_opengl, a program which uses OpenGL to display the evolution of John Conway's Game of Life.
lights_out_opengl, a program which sets up a Lights Out game and allows the user to solve it, using the OpenGL graphics window.
linpack, a library which solves linear systems for a variety of matrix storage schemes, real or complex arithmetic, and single or double precision, by Jack Dongarra, Cleve Moler, Jim Bunch, Pete Stewart.
linpack_bench, a program which is a version of the LINPACK benchmark program;
linpack_c, a library which factors matrices in a variety of formats, and solving the associated linear systems, using single precision complex arithmetic; it includes a function to compute the singular value decomposition (SVD) of a rectangular matrix. This library is by Jack Dongarra, Cleve Moler, Jim Bunch, Pete Stewart.
linpack_d, a library which factors matrices in a variety of formats, and solving the associated linear systems, using double precision real arithmetic; it includes a function to compute the singular value decomposition (SVD) of a rectangular matrix. This library is by Jack Dongarra, Cleve Moler, Jim Bunch, Pete Stewart.
linpack_s, a library which factors matrices in a variety of formats, and solving the associated linear systems, using single precision real arithmetic; it includes a function to compute the singular value decomposition (SVD) of a rectangular matrix. This library is by Jack Dongarra, Cleve Moler, Jim Bunch, Pete Stewart.
linpack_z, a library which factors matrices in a variety of formats, and solving the associated linear systems, using double precision complex arithmetic; it includes a function to compute the singular value decomposition (SVD) of a rectangular matrix. This library is by Jack Dongarra, Cleve Moler, Jim Bunch, Pete Stewart.
linplus, a library which carries out operations such as matrix-vector products, matrix factorization, direct and iterative linear solves for matrices in a variety of formats, including banded, border-banded, circulant, lower triangular, pentadiagonal, sparse, symmetric, toeplitz, tridiagonal, upper triangular and vandermonde formats.
lissajous_opengl, a program which computes a Lissajous curve and displays it using OpenGL.
llsq, a library which solves the simple linear least squares problem of finding the formula of a straight line y=a*x+b which minimizes the root-mean-square error to a set of N data points.
lyness_rule, a library which returns Lyness-Jespersen quadrature rules for the triangle.
machar, a library which dynamically computes the values of various machine constants, by William Cody;
machine, a library which returns tabulated values of the constants associated with computer arithmetic
makefiles, programs which illustrate the use of makefiles for maintaining a software project;
mandelbrot, a program which generates an ASCII Portable Pixel Map (PPM) image of the Mandelbrot set;
mandelbrot_openmp, a program which generates an ASCII Portable Pixel Map (PPM) image of the Mandelbrot fractal set, using OpenMP for parallel execution;
mandelbrot_ppm, a program which generates a binary Portable Pixel Map (PPM) image of the Mandelbrot set, by Eric Weeks.
matrix_exponential, a library which demonstrates some simple approaches to the problem of computing the exponential of a matrix.
md, a program which carries out a molecular dynamics simulation, and is intended as a starting point for implementing an OpenMP parallel version.
md_openmp, a program which carries out a molecular dynamics simulation in parallel using OpenMP.
memory_test, a program which tests the memory available on the computer by declaring and using big vectors.
mesh_bandwidth, a program which returns the geometric bandwidth associated with a mesh of elements of any order and in a space of M dimensions.
mesh_display_opengl, a program which reads data defining a polygonal mesh and displays it using OpenGL.
mesh_io, a library which can read or write MESH files, which can be used to define the geometry of a finite element mesh in 2D or 3D, using triangles, quadrilaterals, tetrahedrons or hexahedrons.
mesh_to_ice, a program which reads the ICE dataset information from a MESH files and rewrites it to a netcdf file.
mgmres, a library which applies the restarted GMRES algorithm to a sparse linear system, by Lili Ju;
mgs, a library which is an example of legacy code, which the student is encouraged to study, understand, and document. This example is due to Diane O'Leary.
minesweeper, a program which is a text-based implementation of the Minesweeper game, by Detelina Stoyanova.
minpack, a library which solves linear and nonlinear least squares problems.
mixed, programs which illustrate mixed language programming in which the main program is written in C++;
monte_carlo_rule, a program which generates a dataset of N random M-dimensional points, regards it as a quadrature rule for the unit hypercube, and writes out three files of information.
mpas_draw, a program which reads an MPAS netcdf file and allows the interactive user to select grid and solution information to be displayed using OpenGL, by Doug Jacobsen, John Burkardt, Geoff Womeldorff.
mpas_grid_display_opengl, programs which can read an MPAS netcdf grid file and display the primary polygonal mesh or the dual triangular mesh, using OpenGL.
mpi, programs which illustrate the use of the MPI Message Passing Interface library, for parallel programming;
mpi_stubs, a library which may be used if a parallel program is to be compiled, loaded and run in a serial machine;
multigrid_poisson_1d, a library which applies the multigrid method to a discretized version of the 1D Poisson equation.
multitask_mpi, a program which demonstrates how to multitask, that is, to execute several unrelated and distinct tasks simultaneously, using MPI for parallel execution.
multitask_openmp, a program which demonstrates how to multitask, that is, to execute several unrelated and distinct tasks simultaneously, using OpenMP for parallel execution.
mxm, a program which sets up a matrix multiplication problem A=B*C of arbitrary size, and compares the time required for IJK, IKJ, JIK, JKI, KIJ and KJI orderings of the loops.
mxm_openmp, a program which computes a dense matrix product C=A*B, using OpenMP for parallel execution.
mxm_serial, a program which sets up a matrix multiplication problem A=B*C, intended as a starting point for implementing a parallel version.
mxv, a program which compares the performance of (for i, for j) loops and (for j, for i) loops for computing the product of an MxN matrix A and an N vector X.
mxv_openmp, a program which compares the performance of the matrix multiplication problem y=A*x, with and without parallelization by OpenMP.
my_calendar, a program which checks a calendar file and prints any line containing the current date;
nas, a program which runs the NASA kernel benchmark.
nast2d, a program which applies the 2D finite volume method to solve transient Navier Stokes flow, by Griebel, Dornseifer and Neunhoffer;
ncc_tetrahedron, a library which defines Newton-Cotes closed quadrature rules in a tetrahedron;
ncc_triangle, a library which defines Newton-Cotes closed quadrature rules in a triangle;
nco_tetrahedron, a library which defines Newton-Cotes open quadrature rules in a tetrahedron;
nco_triangle, a library which defines Newton-Cotes open quadrature rules in a triangle;
nearest_interp_1d, a library which interpolates a set of data using a piecewise constant interpolant defined by the nearest neighbor criterion.
neighbors_to_metis_graph, a program which reads a file describing the neighbor structure of a triangulation or tet mesh, and writes a metis graph file suitable for input to one of the family of programs based on metis;
netcdf, programs which read and write netcdf files.
netcdf_mpas, a library which reads MPAS grid data stored as a netcdf file.
niederreiter, a library which computes elements of the Niederreiter quasirandom M-dimensional sequence [arbitrary base];
niederreiter2, a library which computes elements of the Niederreiter M-dimensional quasirandom sequence using base 2;
niederreiter2_dataset, a program which creates a Niederreiter quasirandom dataset with base 2;
nint_exactness, a program which tests the polynomial exactness of an M-dimensional quadrature rule defined over a finite rectangular product region.
nint_exactness_mixed, a program which measures the polynomial exactness of an M-dimensional quadrature rule based on a mixture of 1D quadrature rule factors.
nintlib, a library which implements approximate integration (quadrature) in M dimensions;
normal, a library which implements a random number generator (RNG) for normally distributed values.
normal_dataset, a program which creates a multivariate normal random dataset;
obj_io, a library which reads and writes the data in an OBJ 3D graphics file;
ode, a library which solves a system of ordinary differential equations, by Shampine and Gordon;
off2plc, a program which reads an OFF file describing a surface in 3D, and creates an equivalent PLC file, by Alex Rand.
opengl, programs which illustrate the use of the OpenGL graphics library;
openmp, programs which illustrate the use of the OpenMP interface for parallel programming shared memory system;
openmp_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.
optimize, a program which looks at the use of automatic optimization by the compiler to make the executable code run faster.
ornstein_uhlenbeck, a library which approximates solutions of the Ornstein-Uhlenbeck stochastic ordinary differential equation (SODE) using the Euler method, the Euler-Maruyama method, and the Milstein method, and creating graphics files for processing by gnuplot.
owens, a library which evaluates Owen's T function;
pages, a program which counts the number of pages (sets of 60 lines) in a file;
partition_problem, a library which seeks solutions of the partition problem, splitting a set of integers into two subsets with equal sum.
patterson_rule, a program which looks up and writes out a Gauss-Patterson quadrature rule of given order.
pbma_io, a library which reads or writes an ASCII Portable Bit Map (PBM) graphics file;
pbma_to_pbmb, a program which converts an ASCII Portable Bit Map (PBM) file to a binary PBM file;
pbmb_io, a library which reads or writes a binary Portable Bit Map (PBM) graphics file;
pbmb_to_pbma, a program which converts a binary Portable Bit Map (PBM) file to an ASCII PBM file;
pce_burgers, a program which defines and solves the time-dependent viscous Burgers equation, with uncertain viscosity, using a polynomial chaos expansion in terms of Hermite polynomials, by Gianluca Iaccarino.
pce_ode_hermite, a program which sets up a simple scalar ordinary differential equation (ODE) for exponential decay with an uncertain decay rate, using a polynomial chaos expansion in terms of Hermite polynomials.
pcl_read, a program which reads a PCL file and extracts the numeric data;
pgma_io, a library which reads or writes an ASCII Portable Gray Map (PGM) 2D graphics file;
pgma_to_pgmb, a program which converts an ASCII Portable Gray Map (PGM) file to a binary PGM file;
pgmb_io, a library which reads or writes a binary Portable Gray Map (PGM) 2D graphics file;
pgmb_to_pgma, a program which converts a binary Portable Gray Map (PGM) file to an ASCII PGM file;
piecewise_linear_product_integral, a library which calculates the exact value of the integral of the product of two piecewise linear functions f(x) and g(x).
pink_noise, a library which computes a pink noise signal obeying a 1/f power law.
png_crc, a program which computes the Cyclic Redundancy Checksum (CRC) for a PNG graphics file (or any other file, for that matter);
point_merge, a library which considers N points in M dimensional space, and counts or indexes the unique or tolerably unique items.
poisson_openmp, a program which computes an approximate solution to the Poisson equation in a rectangle, using the Jacobi iteration to solve the linear system, and OpenMP to carry out the Jacobi iteration in parallel.
poisson_serial, a program which computes an approximate solution to the Poisson equation in a rectangle, and is intended as the starting point for the creation of a parallel version.
poisson_simulation, a library which simulates a Poisson process in which events randomly occur with an average waiting time of Lambda.
polpak, a library which evaluates a variety of mathematical functions, including Chebyshev, Gegenbauer, Hermite, Jacobi, Laguerre, Legendre polynomials, and the Collatz sequence.
polygon_moments, a library which computes arbitrary moments of a polygon.
polygonal_surface_display_opengl, a program which displays a surface in 3D described as a set of polygons using OpenGL;
power_method, a library which carries out the power method for finding a dominant eigenvalue and its eigenvector.
power_rule, a program which constructs a power rule, that is, a product quadrature rule from identical 1D factor rules.
ppma_io, a library which reads or writes an ASCII Portable Pixel Map (PPM) graphics file;
ppma_to_bmp, a program which converts an ASCII Portable Pixel Map (PPM) graphics file to a Microsoft BMP file;
ppma_to_ppmb, a program which converts an ASCII Portable Pixel Map (PPM) graphics file to binary PPM format;
ppmb_io, a library which reads or writes a binary Portable Pixel Map (PPM) 2D pixel map graphics file;
ppmb_to_bmp, a program which converts a binary Portable Pixel Map (PPM) file to a Microsoft BMP file;
ppmb_to_ppma, a program which converts a binary Portable Pixel Map (PPM) file to ASCII PPM format;
prime_mpi, a program which counts the number of primes between 1 and N, using MPI for parallel execution.
prime_openmp, a program which counts the number of primes between 1 and N, using OpenMP for parallel execution.
prime_serial, a program which counts the number of primes between 1 and N, and is intended as a starting point for a parallel version.
prob, a library which evaluates the probability density function (PDF), cumulative density function (CDF), mean, variance and samples of various probability distributions;
product_rule, a program which creates an M-dimensional quadrature rule from a product of distinct 1d quadrature rules;
pthreads, programs which illustrate the use of the POSIX thread library to carry out parallel program execution.
pwl_approx_1d, a library which approximates a set of data using a piecewise linear function.
pwl_interp_1d, a library which interpolates a set of data using a piecewise linear function.
pwl_interp_2d, a library which evaluates a piecewise linear interpolant to data defined on a regular 2D grid.
pwl_interp_2d_scattered, a library which evaluates a piecewise linear interpolant to data which is available at an irregularly arranged set of points.
pyramid_exactness, a program which investigates the polynomial exactness of a quadrature rule for the pyramid.
pyramid_rule, a program which computes a quadrature rule for a pyramid.
qr_solve, a library which computes the least squares solution of a linear system A*x=b.
quad_mesh, a library which handles meshes of quadrilaterals over a 2D region;
quad_mesh_rcm, a program which computes the Reverse Cuthill-McKee (RCM) reordering for nodes in a mesh of 4-node quadrilaterals.
quad_mpi, a program which applies a quadrature rule to estimate an integral, working in parallel by using MPI;
quad_openmp, a program which applies a quadrature rule to estimate an integral over a 1D interval, using OpenMP for parallel execution;
quad_serial, a program which applies a quadrature rule to estimate an integral, intended as a starting point for parallelization exercises.
quad2d_openmp, a program which applies a product quadrature rule to estimate an integral over a 2D rectangle, using OpenMP for parallel execution;
quad2d_serial, a program which applies a product quadrature rule to estimate an integral over a 2D rectangle, intended as a starting point for parallelization exercises.
quadrature_test, a program which reads the definition of an M-dimensional quadrature rule from three files, applies the rule to a number of test integrals, and prints the results.
quadrature_test_2d, a program which reads the definition of a 2D quadrature rule from three files, applies the rule to a number of test integrals, and prints the results.
quadrule, a library which defines quadrature rules for approximating integrals;
quality, a library which computes quality of distribution measures for pointsets in M dimensions;
r4lib, a library which contains many utility routines, using single precision real (R4) arithmetic.
r8lib, a library which contains many utility routines, using double precision real (R8) arithmetic.
rand48, a program which demonstrates the use of the RAND48 family of random number generators (RNG's) available in the C/C++ standard library.
randlc, a library which implements a random number generator (RNG) used by the NAS Benchmark programs.
random_data, a library which generates sample points for various probability distributions, spatial dimensions, and geometries;
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.
random_openmp, a program which illustrates how a parallel program using OpenMP can generate multiple distinct streams of random numbers.
ranlib, a library which produces random samples from Probability Density Functions (PDF's), including Beta, Chi-square Exponential, F, Gamma, Multivariate normal, Noncentral chi-square, Noncentral F, Univariate normal, random permutations, Real uniform, Binomial, Negative Binomial, Multinomial, Poisson and Integer uniform, by Barry Brown and James Lovato.
rbf_interp_1d, a library which defines and evaluates radial basis function (RBF) interpolants to 1D data.
rbf_interp_2d, a library which defines and evaluates radial basis function (RBF) interpolants to 2D data.
rbf_interp_nd, a library which defines and evaluates radial basis function (RBF) interpolants to M-dimensional data.
rcm, a library which applies the Reverse Cuthill McKee (RCM) algorithm for reordering the nodes of a graph, and reducing the bandwidth of a corresponding sparse matrix;
reactor simulation, a program which is a simple Monte Carlo simulation of the shielding effect of a slab of a certain thickness in front of a neutron source. This program was provided as an example with the book "Numerical Methods and Software."
recomment, a program which converts C style comments to C++ style;
reword, a program which makes a copy of a file in which every line has the same number of words;
ring_mpi, a program which uses the MPI parallel programming environment, and measures the time necessary to copy a set of data around a ring of processes.
rk4, a library which applies the fourth order Runge-Kutta (RK) algorithm to estimate the solution of an ordinary differential equation (ODE) at the next time step.
rkf45, a library which applies the Runge-Kutta-Fehlberg (RKF) algorithm to solve a system of ordinary differential equations (ODEs);
rnglib, a library which implements a random number generator (RNG) with splitting facilities, allowing multiple independent streams to be computed, by L'Ecuyer and Cote.
rot13, a program which makes a copy of a file which has been encoded using the ROT13 coding and a ROT5 coding for digits.
rotating_cube_display_opengl, a program which displays a rotating color cube in 3D, using OpenGL;
rsites, a program which generates random pointsets with integer coordinates for testing computational geometry algorithms;
sandia_cubature, a library which implements quadrature rules for certain M-dimensional regions and weight functions.
sandia_rules, a library which can produce a standard 1D Gauss quadrature rule of Chebyshev, Gegenbauer, generalized Hermite, generalized Laguerre, Hermite, Jacobi, Laguerre, and Legendre types.
sandia_rules2, a library which contains a very small selection of functions which serve as an interface between sandia_sgmg or sandia_sgmga and sandia_rules.
sandia_sgmg, a library which creates sparse grids based a mixture of 1D factor rules, and which experiments with the idea of using linear growth rate for the 1D quadrature rules. This is a version of sparse_grid_mixed_growth that uses a different procedure for supplying the parameters needed to evaluate certain quadrature rules.
sandia_sgmga, a library which creates sparse grids based on a mixture of 1D quadrature rules, allowing anisotropic weights for each dimension. This is a version of sgmga that uses a different procedure for supplying the parameters needed to evaluate certain quadrature rules.
sandia_sgmgg, a library which explores a generalized construction method for sparse grids.
sandia_sparse, a library which can produce an M-dimensional sparse grid, based on a variety of 1D quadrature rules; only isotropic grids are generated, that is, the same rule is used in each dimension, and the same maximum order is used in each dimension.
satisfy, a program which demonstrates, for a particular circuit, an exhaustive search for solutions of the circuit satisfy problem.
satisfy_mpi, a program which demonstrates, for a particular circuit, an exhaustive search for solutions of the circuit satisfy problem, using MPI to perform the calculation in parallel.
satisfy_openmp, a program which demonstrates, for a particular circuit, an exhaustive search for solutions of the circuit satisfy problem, using OpenMP for parallel execution.
schedule_openmp, a program which demonstrates the default, static, and dynamic methods of scheduling loop iterations in OpenMP to avoid work imbalance.
scvt_mpi, a program which determines a Centroidal Voronoi Tessellation (CVT) on a spherical surface, using triangle for triangulation, and the boost mpi and serialization libraries, by Doug Jacobsen;
sde, a library which illustrates properties of stochastic ordinary differential equations (SODE's), and common algorithms for their analysis, by Desmond Higham;
search_mpi, a program which searches integers between A and B for a solution J such that F(J)=C, carried out in parallel using MPI.
search_serial, a program which searches the integers from A to B for a value J such that F(J) = C. this version of the program is intended as a starting point for a parallel approach.
set_theory, a library which demonstrates various set theoretic operations using several models of a set.
sftpack, a library which implements the slow Fourier transform, intended as a teaching tool and comparison with the Fast Fourier Transform (FFT).
sgefa_openmp, a program which compares a standard linear algebra solver against a revised version which can be run in parallel with OpenMP.
sgmg, a library which creates sparse grids based a mixture of 1D factor rules, and which experiments with the idea of using linear growth rate for the 1D quadrature rules.
sgmga, a library which creates sparse grids based on a mixture of 1D quadrature rule families, allowing user specified growth rates for each family, allowing anisotropic weights for each dimension.
shallow_water_1d, a program which simulates the evolution of a 1D fluid governed by the time-dependent shallow water equations.
shepard_interp_1d, a library which defines and evaluates Shepard interpolants to 1D data, based on inverse distance weighting.
shepard_interp_2d, a library which defines and evaluates Shepard interpolants to 2D data, based on inverse distance weighting.
shepard_interp_nd, a library which defines and evaluates Shepard interpolants to M-dimensional data, based on inverse distance weighting.
simple_ga, a program which is a simple example of a genetic algorithm, by Dennis Cormier and Sita Raghavan.
simplex_coordinates, a library which computes the Cartesian coordinates of the vertices of a regular simplex in M dimensions.
sine_transform, a library which demonstrates some simple properties of the discrete sine transform.
sobol, a library which computes elements of the Sobol M-dimensional quasirandom sequence, extended to a maximal spatial dimension of 1111, by Bennett Fox.
sobol_dataset, a program which creates a Sobol quasirandom dataset;
sor, a library which implements a simple version of the Successive Over-Relaxation (SOR) method for the iterative solution of a linear system of equations.
sparse_grid_cc, a library which creates sparse grids based on Clenshaw-Curtis rules.
sparse_grid_cc_dataset, a program which creates a sparse grid dataset based on Clenshaw-Curtis rules.
sparse_grid_gl, a library which creates sparse grids based on Gauss-Legendre rules.
sparse_grid_gl_dataset, a program which creates a sparse grid dataset based on Gauss-Legendre rules.
sparse_grid_hermite, a library which creates sparse grids based on Gauss-Hermite rules.
sparse_grid_hermite_dataset, a program which creates a sparse grid dataset based on Gauss-Hermite rules.
sparse_grid_hw, a library which creates sparse grids based on Gauss-Legendre, Gauss-Hermite, Gauss-Patterson, or a nested variation of Gauss-Hermite rules, by Florian Heiss and Viktor Winschel.
sparse_grid_laguerre, a library which creates sparse grids based on Gauss-Laguerre rules.
sparse_grid_laguerre_dataset, a program which creates a sparse grid dataset based on Gauss-Laguerrre rules.
sparse_grid_mixed, a library which creates sparse grids based on a mixture of 1D factor rules.
sparse_grid_mixed_dataset, a program which creates a sparse grid dataset based on a mixture of 1D rules.
sparse_grid_open, a library which creates sparse grids based on open rules (Fejer 2, Gauss-Patterson, Newton-Cotes-Open).
sparse_grid_open_dataset, a program which creates a sparse grid dataset based on open rules (Fejer 2, Gauss-Patterson, Newton-Cotes-Open).
sparse_grid_pure, a library which for the analysis and creation of sparse grids in which a single family of 1D quadrature rules is used for all spatial dimensions.
sparse_interp_nd, a library which can be used to define a sparse interpolant to a function f(x) of a M-dimensional argument.
sphere_exactness, a program which tests the polynomial exactness of a quadrature rule for the unit sphere;
sphere_grid, a library which provides a number of ways of generating grids of points, or of points and lines, or of points and lines and faces, over the unit sphere.
sphere_lebedev_rule, a library which computes Lebedev quadrature rules for the unit sphere;
sphere_monte_carlo, a library which applies a Monte Carlo method to estimate the integral of a function over the surface of the unit sphere in 3D;
sphere_quad, a library which approximates an integral over the surface of the unit sphere by applying a triangulation to the surface;
sphere_stereograph, a library which computes the stereographic mapping between points on the unit sphere and points on the plane Z = 1; a generalized mapping is also available.
sphere_triangle_quad, a library which estimates the integral of a function over a spherical triangle.
sphere_voronoi_display_opengl, a program which displays a sphere and randomly selected generator points, and then gradually colors in points in the sphere that are closest to each generator.
sphere_xyz_display_opengl, a program which reads XYZ information defining points in 3D that lie (presumably) on the unit sphere, and displays an image of the unit sphere and the points, using OpenGL;
spline, a library which interpolates or approximates data by piecewise polynomial functions;
spring_ode, a program which shows how line printer graphics can be used to make a crude illustration of a solution of the ordinary differential equation (ODE) that describes the motion of a weight attached to a spring.
spring_ode2, a program which shows how gnuplot graphics can be used to illustrate a solution of the ordinary differential equation (ODE) that describes the motion of a weight attached to a spring.
star_discrepancy, a program which computes bounds on the star discrepancy of a set of N points in M dimensions, stored as a table file, by Eric Thiemard.
stla_io, a library which reads and writes ASCII stereolithography (STL) 3d graphics files;
stochastic_rk, a library which applies a Runge-Kutta (RK) scheme to a stochastic ordinary differential equation (SODE).
string_simulation, a program which simulates the behavior of a vibrating string, creating files that can be displayed by gnuplot.
stroud, a library which defines quadrature rules for a variety of M-dimensional regions, including the interior of the square, cube and hypercube, the pyramid, cone and ellipse, the hexagon, the M-dimensional octahedron, the circle, sphere and hypersphere, the triangle, tetrahedron and simplex, and the surface of the circle, sphere and hypersphere.
subpak, a library which contains many utility routines;
subset, a library which ranks, unranks and generates random subset, combinations, permutations and other combinatorial objects;
subset_sum, a program which seeks solutions of the subset sum problem.
sum_million, a program which sums the integers from 1 to 1,000,000, as a demonstration of how to rate a computer's speed;
svd_basis, a program which applies the singular value decomposition (SVD) to a set of vectors, to extract the dominant modes;
svd_demo, a program which demonstrates the calculation of the singular value decomposition (SVD) and some of its properties;
svd_snowfall, a program which reads a file containing historical snowfall data and analyzes the data with the Singular Value Decomposition (SVD), displaying the results using gnuplot.
svd_truncated, a program which demonstrates the computation of the reduced or truncated Singular Value Decomposition (SVD) that is useful for cases when one dimension of the matrix is much smaller than the other.
table_border, a program which reads a table file of 2D data, and adds data representing a boundary layer of zero values;
table_delaunay, a program which computes the Delaunay triangulation of a set of points stored as a table file;
table_io, a library which reads and writes a simple dataset stored as a table file;
table_latinize, a program which reads a dataset of N points in M dimensions from a table file, adjusts data to form a Latin Hypercube;
table_quality, a program which reads a dataset of N points in the M-dimensional unit hypercube from a table file, and analyzes the quality of dispersion;
table_unborder, a program which reads a table file of 2D data and deletes the data representing the boundary layer;
table_voronoi, a program which reads a set of 2D points from a table file, and prints information about the Voronoi diagram;
tanh_quad, a library which sets up the tanh quadrature rule;
tanh_sinh_rule, a program which computes and writes out a tanh-sinh quadrature rule of given order.
task_division, a library which implements a simple procedure for smoothly dividing T tasks among P processors; such a method can be useful in MPI and other parallel environments, particularly when T is not an exact multiple of P, and when the processors can be indexed starting from 0 or from 1.
test_approx, a library which implements test problems for approximation, provided as a set of (x,y) data.
test_eigen, a library which implements test matrices for eigenvalue analysis.
test_int, a library which implements test problems for approximate integration (quadrature) in one dimension.
test_int_2d, a library which implements test problems for approximate integration (quadrature) in two dimensions.
test_int_hermite, a library which defines test integrands for Hermite integrals with interval (-oo,+oo) and weight exp(-x^2).
test_int_laguerre, a library which implements test problems for approximate integration over a semi-infinite interval, often carried out using Gauss-Laguerre quadrature.
test_interp, a library which defines test problems for interpolation, provided as a set of (x,y) data.
test_interp_1d, a library which defines test problems for interpolation of data y(x), which depends on a 1D argument.
test_interp_2d, a library which defines test problems for interpolation of data z(x,y), depending on a 2D argument.
test_interp_nd, a library which defines test problems for interpolation of data z(x), depending on an M-dimensional argument.
test_ls, a library which implements linear least squares test problems of the form A*x=b.
test_mat, a library which implements test matrices.
test_matrix_exponential, a library which defines a set of test cases for computing the matrix exponential.
test_min, a library which implements test problems for minimization of a scalar function of a scalar variable.
test_nearest, a library which tests the time complexity of various procedures for solving the nearest neighbor problem.
test_nint, a library which implements test problems for approximate integration (quadrature) in M dimensions.
test_opt_con, a library which defines test problems for the minimization of a scalar function of several variables, with the search constrained to lie within a specified hyper-rectangle.
test_optimization, a library which implements test problems for optimization of a scalar function of several variables, as described by Molga and Smutnicki.
test_tri_int, a library which implements test functions for integration over a triangle in 2 dimensions.
test_values, a library which returns a small set of values for a variety of mathematical functions, suitable for a rough test of accuracy for algorithms;
test_zero, a library which implements test problems for the solution of a single nonlinear equation in one variable.
testpack, a program which demonstrates the use of a set of test integrand functions chosen by Alan Genz.
tet_mesh, a library which works with tetrahedral meshes in 3D;
tet_mesh_boundary, a program which reads information defining a tetrahedral mesh of points in 3D, and determines the triangular faces that form the boundary of the mesh; it writes files containing the nodes and elements defining this tri_surface data.
tet_mesh_display_opengl, a program which reads two files describing a tetrahedral mesh in 3D and displays it using OpenGL.
tet_mesh_l2q, a program which reads information about a 4-node (linear) tet mesh and creates data defining a corresponding 10-node (quadratic) tet mesh;
tet_mesh_q2l, a program which reads information about a 10-node (quadratic) tet mesh and creates data defining a corresponding 4-node (linear) tet mesh;
tet_mesh_quad, a program which estimates the integral of a function over a region defined by a tetrahedral mesh.
tet_mesh_quality, a program which reads a node coordinate file and a tetrahedron file and returns quality measures for the tetrahedronal mesh;
tet_mesh_rcm, a program which computes the Reverse Cuthill McKee (RCM) reordering for the nodes of a tet_mesh that uses 4-node or 10-node tetrahedrons;
tet_mesh_refine, a program which refines a tet mesh;
tet_mesh_tet_neighbors, a program which reads information about a tet mesh and writes out a list of the tetrahedrons that are adjacent to a given tetrahedron;
tet_mesh_volumes, a program which computes the volume of each tetrahedron in a tet mesh;
tetrahedron_exactness, a program which tests the polynomial exactness of a quadrature rule for the tetrahedron;
tetrahedron_grid, a library which computes a grid of points inside a 3D tetrahedron.
tetrahedron_monte_carlo, a program which uses the Monte Carlo method to estimate an integral over a tetrahedron.
tetrahedron_properties, a program which prints out the properties of a tetrahedron whose vertex coordinates are read from a file;
tga_io, a library which reads a TARGA (TGA) graphics file;
three_body_simulation, a program which simulates the behavior of three planets, constrained to lie in a plane, and moving under the influence of gravity, by Walter Gander and Jiri Hrebicek.
tiler_3d, a library which blends a pattern repeatedly into a region in 3D;
timer, programs which measure the CPU time used by a portion of a user's program;
timestamp, a library which prints the current YMDHMS date as a timestamp, and other timekeeping operations;
toeplitz_cholesky, a library which computes the Cholesky factorization of a nonnegative definite symmetric Toeplitz matrix.
toms178, a library which seeks the minimizer of a scalar function of several variables using the Hooke-Jeeves method, by Arthur Kaupe; this is a C++ version of ACM TOMS algorithm 178.
toms179, a library which calculates the incomplete Beta ratio, by Oliver Ludwig; this is a C++ version of ACM TOMS algorithm 179.
toms291, a library which approximates the logarithm of the Gamma function; this is a C++ version of ACM TOMS algorithm 291.
toms446, a library which manipulates Chebyshev series for interpolation and approximation; this is a C++ version of ACM TOMS algorithm 446, by Roger Broucke.
toms462, a library which evaluates the upper right tail of the bivariate normal distribution; that is, the probability that normal variables X and Y with correlation R will satisfy H <= X and K <= Y; this is a C++ version of ACM TOMS algorithm 462.
toms655, a library which computes the weights for interpolatory quadrature rule; this library is commonly called IQPACK; this is a C++ version of ACM TOMS algorithm 655.
tri_surface_display_opengl, a program which displays a 3D triangulated surface or TRI_SURFACE using OpenGL;
tri_surface_io, a library which reads and writes the 3D graphics information in a tri_surface file;
triangle_exactness, a program which tests the polynomial exactness of a quadrature rule for triangles;
triangle_grid, a library which computes a triangular grid of points.
triangle_histogram, a program which computes histograms of data on the unit triangle.
triangle_io, a library which can read or write some of the files created by Jonathan Shewchuk's triangle program.
triangle_monte_carlo, a program which uses the Monte Carlo method to estimate an integral over a triangle.
triangle_properties, a program which computes properties of a triangle whose vertex coordinates are read from a file.
triangle_to_fem, a program which reads the NODE and ELE files created by the triangle program to describe a triangular mesh, and writes a corresponding pair of node and element files in the 2D FEM format.
triangle_to_mesh, a program which reads the NODE and ELE files created by the triangle program to describe a triangular mesh, and writes a corresponding MESH file for input to medit.
triangulation, a library which computes the triangulation of a set of points in the plane, to search it, plot it, and manipulate it in other ways;
triangulation_boundary_nodes, a program which reads data defining a triangulation, determines which nodes lie on the boundary, and writes their coordinates to a file;
triangulation_corner, a program which modifies triangulations in which one or more triangles have more than one boundary edge.
triangulation_delaunay_discrepancy, a program which measures the amount by which a triangulation fails the local Delaunay test;
triangulation_display_opengl, a program which reads two files describing a triangulation and displays it using OpenGL.
triangulation_histogram, a program which computes histograms of data over a triangulation.
triangulation_l2q, a program which reads information about a 3-node (linear) triangulation and creates data defining a corresponding 6-node (quadratic) triangulation;
triangulation_mask, a program which reads a triangulation and calls a user-supplied routine to consider each triangle for deletion;
triangulation_orient, a program which ensures that the triangles in an order 3 or order 6 triangulation have positive orientation;
triangulation_plot, a program which plots the nodes and elements of a triangulation as a PostScript (PS) file;
triangulation_q2l, a program which reads information about a 6-node triangulation and creates data defining a corresponding 3-node triangulation;
triangulation_quad, a program which reads information about a triangulation and the value of a function at the nodes and estimates the integral of the function over the triangulated region.
triangulation_quality, a program which reads information about a triangulation and computes various quality measures;
triangulation_rcm, a program which computes the Reverse Cuthill McKee (RCM) reordering for the nodes of a triangulation that uses 3-node or 6-node triangles;
triangulation_refine, a program which refines a triangulation;
triangulation_t3_to_t4, a program which reads information about a 3-node triangulation and creates data defining a corresponding 4-node triangulation (vertices + centroid);
triangulation_triangle_neighbors, a program which reads data defining a triangulation, finds the three neighbor triangles for each triangle, and writes them to a file;
uncontrol, a program which removes control characters from a file;
unicycle, a library which considers permutations containing a single cycle, sometimes called cyclic permutations.
uniform, a library which implements uniform random number generators (RNG's) for a variety of arithmetic types.
uniform_dataset, a program which creates a file of uniform pseudorandom values;
valgrind, programs which illustrate the use of the VALGRIND suite of programs which includes a memory leak detector.
van_der_corput, a library which evaluates the van der Corput 1-dimensional quasirandom sequence;
van_der_corput_dataset, a program which creates a file of quasirandom van der Corput values;
vandermonde_approx_1d, a library which finds a polynomial approximant to data y(x) of a 1D argument by setting up and solving an overdetermined linear system for the polynomial coefficients involving the Vandermonde matrix.
vandermonde_approx_2d, a library which finds a polynomial approximant p(x,y) to data z(x,y) of a 2D argument by setting up and solving an overdetermined linear system for the polynomial coefficients involving the Vandermonde matrix.
vandermonde_interp_1d, a library which finds a polynomial interpolant to data y(x) of a 1D argument by setting up and solving a linear system for the polynomial coefficients involving the Vandermonde matrix.
vandermonde_interp_2d, a library which finds a polynomial interpolant to data z(x,y) of a 2D argument by setting up and solving a linear system for the polynomial coefficients involving the Vandermonde matrix.
vector_read, a library which reads, one line at a time, numeric data, storing it in a vector returned to the user;
walsh, a library which implements versions of the Walsh and Haar transforms.
wandzura, a library which returns a Wandzura rule for quadrature in a triangle;
wavelet, a library which does some simple wavelet calculations;
weekday, a library which determines the day of the week corresponding to a given date, such as 14 October 1066, Julian calendar, ...which was a Saturday.
width, a program which returns the length and location of the longest line in a text file;
wrap, a program which wraps a file so that no line is longer than 80 characters;
wrap2, a program which wraps a text file so that no line is longer than a given maximum; the wrapping may be done slightly earlier than at the maximum, in order to avoid breaking words;
wtime, a library which returns a reading of the wall clock time in seconds.
x_window, programs which illustrate the use of the X Window graphical terminal system;
xy_display_opengl, a program which reads XY information defining points in 2D and displays an image using OpenGL;
xy_io, a library which reads and writes XY, XYL and XYF files defining points, lines, and faces in 2D.
xyf_display_opengl, a program which reads XYF information defining points and faces in 2D and displays an image using OpenGL;
xyl_display_opengl, a program which reads XYL information defining points and lines in 2D and displays an image using OpenGL;
xyz_display_opengl, a program which reads XYZ information defining points in 3D and displays an image using OpenGL;
xyz_io, a library which reads and writes XYZ, XYZL and XYZF files defining points, lines and faces in 3D.
xyzf_display_opengl, a program which reads XYZF information defining points and faces in 3D and displays an image using OpenGL;
xyzl_display_opengl, a program which reads XYZL information defining points and lines in 3D and displays an image using OpenGL;
zero_rc, a library which seeks solutions of a scalar nonlinear equation f(x)=0, or a system of nonlinear equations, using reverse communication.
ziggurat, a library which samples the uniform, normal or exponential distributions, using the ziggurat method, by Marsaglia and Tsang.
ziggurat_openmp, a program which demonstrates how the ziggurat library can be used as a random number generator (RNG) in an OpenMP parallel program.

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

Last revised on 27 March 2013.