- allocatable_array_test
- alpert_rule, a C++ code which sets up an Alpert quadrature rule for functions which are regular, log(x) singular, or 1/sqrt(x) singular.
- alpert_rule_test
- anagram, a C++ code which determines anagrams of a string, by James Cherry;
- anagram_test
- analemma, a C++ code 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.
- analemma_test
- annulus_monte_carlo, a C++ code which uses the Monte Carlo method to estimate the integral of a function over the interior of a circular annulus in 2D.
- annulus_monte_carlo_test
- annulus_rule, a C++ code which computes a quadrature rule for estimating integrals of a function over the interior of a circular annulus in 2D.
- annulus_rule_test
- args, a C++ code which reports the command line arguments with which it was invoked;
- args_test
- asa005, a C++ code which evaluates the lower tail of the noncentral Student's T distribution, by BE Cooper. This is a version of Applied Statistics Algorithm 5;
- asa005_test
- asa006, a C++ code which computes the Cholesky factor of a symmetric positive definite (SPD) matrix, by Michael Healy. This is a version of Applied Statistics Algorithm 6;
- asa006_test
- asa007, a C++ code which computes the inverse of a symmetric positive definite (SPD) matrix, by Michael Healy. This is a version of Applied Statistics Algorithm 7.
- asa007_test
- asa032, a C++ code which evaluates the incomplete Gamma function, by G Bhattacharjee. This is a version of Applied Statistics Algorithm 32;
- asa032_test
- asa047, a C++ code which minimizes a scalar function of several variables using the Nelder-Mead algorithm, by R ONeill. This is a version of Applied Statistics Algorithm 47;
- asa047_test
- asa053, a C++ code which produces sample matrices from the Wishart distribution, by William Smith and Ronald Hocking. This is a version of Applied Statistics Algorithm 53.
- asa053_test
- asa058, a C++ code which carries out K-means clustering of data, by David Sparks. This is a version of Applied Statistics Algorithm 58;
- asa058_test
- asa063, a C++ code which evaluates the incomplete Beta function, by KL Majumder and G Bhattacharjee. This is a version of Applied Statistics Algorithm 63;
- asa063_test
- asa066, a C++ code which computes the cumulative density function (CDF) of the standard normal distribution, by David Hill. This is a version of Applied Statistics Algorithm 66;
- asa066_test
- asa076, a C++ code which evaluates Owen's T function, by Young and Minder. This is a version of Applied Statistics Algorithm 76;
- asa076_test
- asa082, a C++ code which computes the determinant of an orthogonal matrix; this is Applied Statistics Algorithm 82, by J C Gower.
- asa082_test
- asa091, a C++ code which computes the percentage points of the Chi-square distribution, by Best and Roberts. This is a version of Applied Statistics Algorithm 91;
- asa091_test
- asa103, a C++ code which evaluates the digamma or psi function, by Jose Bernardo. This is a version of Applied Statistics Algorithm 103;
- asa103_test
- asa109, a C++ code which inverts the incomplete Beta function, by Cran, Martin and Thomas. This is a version of Applied Statistics Algorithm 109;
- asa109_test
- asa111, a C++ code which evaluates the percentage points of the normal distribution, by Beasley and Springer. This is a version of Applied Statistics Algorithm 111;
- asa111_test
- asa113, a C++ code which clusters data using a transfer and swap algorithm, by Banfield and Bassill. This is a version of Applied Statistics Algorithm 113;
- asa113_test
- asa121, a C++ code which evaluates the trigamma function, by BE Schneider. This is a version of Applied Statistics Algorithm 121;
- asa121_test
- asa136, a C++ code which implements a clustering algorithm, by Hartigan and Wong. This is a version of Applied Statistics Algorithm 136;
- asa136_test
- asa144, a C++ code which randomly generates tables with given row and column sums, by James Boyett. This is a version of Applied Statistics Algorithm 144;
- asa144_test
- asa147, a C++ code which evaluates the incomplete Gamma function, by Chi Leung Lau. This is a version of Applied Statistics Algorithm 147;
- asa147_test
- asa152, a C++ code which evaluates the probability density function (PDF) and cumulative density function (CDF) associated with the hypergeometric distribution, by Richard Lund. This is a version of Applied Statistics Algorithm 152;
- asa152_test
- asa159, a C++ code which computes a random table with given row and column sums, by Michael Patefield. This is a version of Applied Statistics Algorithm 159;
- asa159_test
- asa172, a C++ code 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 version of Applied Statistics Algorithm 172;
- asa172_test
- asa183, a C++ code which implements a random number generator (RNG), by Wichman and Hill. This is a version of Applied Statistics Algorithm 183;
- asa183_test
- asa226, a C++ code which evaluates the noncentral incomplete Beta function, by Russell Lenth. this is a version of Applied Statistics Algorithm 226;
- asa226_test
- asa239, a C++ code which evaluates the incomplete Gamma function, by Shea. this is a version of Applied Statistics Algorithm 239;
- asa239_test
- asa241, a C++ code which computes the inverse of the normal cumulative density function (CDF), by Michael Wichura; this is a version of Applied Statistics Algorithm 241;
- asa241_test
- asa243, a C++ code which evaluates the lower tail of the noncentral Student's T distribution, by Russell Lenth; this is a version of Applied Statistics Algorithm 243;
- asa243_test
- asa245, a C++ code which evaluates the logarithm of the Gamma function, by Allan Mcleod; this is a version of Applied Statistics Algorithm 245;
- asa245_test
- asa266, a C++ code which evaluates various properties of the Dirichlet probability density function (PDF); this is a version of Applied Statistics Algorithm 266;
- asa266_test
- asa299, a C++ code which computes the lattice points (integer coordinates) in an M-dimensional simplex, by Chasalow and Brand; this is a version of Applied Statistics Algorithm 299;
- asa299_test
- asa310, a C++ code which evaluates the noncentral incomplete Beta function, by Chattamvelli and Shanmugam; this is a version of Applied Statistics Algorithm 310;
- asa310_test
- asa314, a C++ code which computes the inverse of a matrix whose elements are subject to modulo arithmetic, by Roger Payne. This is a version of Applied Statistics Algorithm 314.
- asa314_test
- atbash, a C++ code which applies the Atbash substitution cipher to a string of text.
- atbash_test
- backtrack_binary_rc, a C++ code which carries out a backtrack search for a set of binary decisions, using reverse communication.
- backtrack_binary_rc_test
- ball_grid, a C++ code which computes grid points inside a 3D ball.
- ball_grid_test
- ball_integrals, a C++ code which returns the exact value of the integral of any monomial over the interior of the unit ball in 3D.
- ball_integrals_test
- ball_monte_carlo, a C++ code which applies a Monte Carlo method to estimate integrals of a function over the interior of the unit ball in 3D;
- ball_monte_carlo_test
- bank, a C++ code which computes the check digit associated with a US Bank Routing Number check digit, or reports whether a 9-digit code is actually valid.
- bank_test
- barycentric_interp_1d, a C++ code 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.
- barycentric_interp_1d_test
- bellman_ford, a C++ code which implements the Bellman-Ford algorithm for finding the shortest distance from a given node to all other nodes in a directed graph whose edges have been assigned real-valued lengths.
- bellman_ford_test
- bernstein_polynomial, a C++ code which evaluates the Bernstein polynomials;
- bernstein_polynomial_test
- besselj, a C++ code which evaluates Bessel J functions of noninteger order;
- besselj_test
- beta_nc, a C++ code which evaluates the cumulative density function (CDF) of the noncentral Beta distribution;
- beta_nc_test
- bins, a C++ code which handles geometric data partially sorted into bins;
- bins_test
- bisection_integer, a C++ code which seeks an integer solution to the equation F(X)=0, using bisection within a user-supplied change of sign interval [A,B].
- bisection_integer_test
- bisection_rc, a C++ code 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 (RC).
- bisection_rc_test
- black_scholes, a C++ code which implements some simple approaches to the Black-Scholes option valuation theory, by Desmond Higham.
- black_scholes_test
- blas0, a C++ code which contains auxilliary functions for the Basic Linear Algebra Subprograms (BLAS).
- blas0_test
- blas1_d, a C++ code 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_d_test
- blas1_z, a C++ code which constitutes the Level 1 Basic Linear Algebra Subprograms (BLAS), using double precision complex arithmetic, by Charles Lawson, Richard Hanson, David Kincaid, Fred Krogh.
- blas1_z_test
- blas2_d, a C++ code which constitutes the Level 2 Basic Linear Algebra Subprograms (BLAS), for matrix-vector operations, using double precision real arithmetic.
- blas2_d_test
- blas3_d, a C++ code which constitutes the Level 3 Basic Linear Algebra Subprograms (BLAS) for matrix-matrix operations, using double precision arithmetic;
- blas3_d_test
- blend, a C++ code which uses transfinite interpolation to fill in data within a square or a cube based on boundary values;
- blend_test
- boost_test
- box_behnken, a C++ code which computes a Box-Behnken design, that is, a set of arguments to sample the behavior of a function of multiple parameters;
- box_behnken_test
- brownian_motion_simulation, a C++ code which simulates Brownian motion in an M-dimensional region, creating graphics files for processing by gnuplot().
- brownian_motion_simulation_test
- burgers_solution, a C++ code which evaluates exact solutions of the time-dependent 1D viscous Burgers equation.
- burgers_solution_test
- bvec, a C++ code which demonstrates how signed integers can be stored as binary vectors, and arithmetic can be performed on them.
- bvec_test
- c++_calls_c_test, a C++ code which which calls a C function.
- c++_calls_f77_test, a C++ code which calls a FORTRAN77 subroutine.
- c++_calls_f90_test, a C++ code which calls a FORTRAN90 subroutine.
- c8_complex_lib, a C++ code which implements double precision complex (C8) numbers;
- c8_complex_lib_test
- c8lib, a C++ code which implements certain elementary functions for double precision complex (C8) variables;
- c8lib_test
- caesar, a C++ code which applies a Caesar Shift Cipher to a string of text.
- caesar_test
- calendar_rd, a C++ code which computes the representation of a given date in a number of calendrical systems, by Edward Reingold and Nachum Dershowitz.
- calendar_rd_test
- calpak, a C++ code which makes various calendar calculations;
- calpak_test
- catalog, a C++ code which reads a C, C++, FORTRAN77, FORTRAN90, MATLAB, Python, or R file, and prints every line that begins with a special index tag that identifies each individualfunction. If the code has been marked up expecting this convention, it is a handy way of making a table of contents of a program file.
- catalog_test
- cauchy_principal_value, a C++ code which uses Gauss-Legendre quadrature to estimate the Cauchy Principal Value (CPV) of certain singular integrals.
- cauchy_principal_value_test
- cblas_test
- ccn_rule, a C++ code which defines a nested Clenshaw Curtis quadrature rule.
- ccn_rule_test
- ccs_io, a C++ code which reads and writes sparse linear systems stored in the Compressed Column Storage (CCS) format.
- ccs_io_test
- ccs_to_st, a C++ code which converts a sparse matrix from Compressed Column Storage (CCS) to sparse triplet (ST) format.
- ccs_to_st_test
- cdflib, a C++ code 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.
- cdflib_test
- cell, a C++ code which defines a cell array, a generalization of an array which compactly stores and retrieves vector or matrix data of varying size, such as the rows of a triangular matrix.
- cell_test
- cellular_automaton, a C++ code which demonstrates the 1D Cellular Automaton known as rule #30.
- cellular_automaton_test
- cg, a C++ code which implements a simple version of the conjugate gradient (CG) method for solving a system of linear equations of the form A*x=b, suitable for situations in which the matrix A is symmetric positive definite (SPD)
- cg_test
- cg_rc, a C++ code which implements the conjugate gradient (CG) method for solving a symmetric positive definite (SPD) sparse linear system A*x=b, using reverse communication (RC).
- cg_rc_test
- change_making, a C++ code which considers the change making problem, in which a given sum is to be formed using coins of various denominations.
- change_making_test
- chebyshev, a C++ code which computes the Chebyshev interpolant or approximant to a given function over an interval.
- chebyshev_test
- chebyshev_interp_1d, a C++ code which determines the combination of Chebyshev polynomials which interpolates a set of data, so that p(x(i)) = y(i).
- chebyshev_interp_1d_test
- chebyshev_polynomial, a C++ code which considers the Chebyshev polynomials T(i,x), U(i,x), V(i,x) and W(i,x). Functions are provided to evaluate the polynomials, determine their zeros, produce their polynomial coefficients, produce related quadrature rules, project other functions onto these polynomial bases, and integrate double and triple products of the polynomials.
- chebyshev_polynomial_test
- chebyshev_series, a C++ code which evaluates a Chebyshev series approximating a function f(x), while efficiently computing one, two or three derivatives of the series, which approximate f'(x), f''(x), and f'''(x), by Manfred Zimmer.
- chebyshev_series_test
- chebyshev1_exactness, a C++ code which tests the polynomial exactness of Gauss-Chebyshev type 1 quadrature rules.
- chebyshev1_exactness_test
- chebyshev1_rule, a C++ code which computes and writes out a Gauss-Chebyshev type 1 quadrature rule of given order.
- chebyshev1_rule_test
- chebyshev2_exactness, a C++ code which tests the polynomial exactness of Gauss-Chebyshev type 2 quadrature rules.
- chebyshev2_exactness_test
- chebyshev2_rule, a C++ code which computes and writes out a Gauss-Chebyshev type 2 quadrature rule of given order.
- chebyshev2_rule_test
- chrpak, a C++ code which manipulates characters and strings;
- chrpak_test
- circle_arc_grid, a C++ code which computes grid points along a circular arc.
- circle_arc_grid_test
- circle_integrals, a C++ code which returns the exact value of the integral of any monomial over the circumference of the unit circle in 2D.
- circle_integrals_test
- circle_monte_carlo, a C++ code which applies a Monte Carlo method to estimate the integral of a function along the circumference of the unit circle in 2D;
- circle_monte_carlo_test
- circle_rule, a C++ code which computes quadrature rules for the unit circle in 2D, that is, the circumference of the circle of radius 1 and center (0,0).
- circle_rule_test
- circle_segment, a C++ code which carries out computations associated with a circle segment, including height, angle, area, centroid, sampling, and quadrature, creating some graphics images using gnuplot().
- circle_segment_test
- cities, a C++ code which works with city-to-city distance matrices;
- cities_test
- clapack_test
- clausen, a C++ code which evaluates a Chebyshev interpolant to the Clausen function Cl2(x).
- clausen_test
- clenshaw_curtis_grid, a C++ code which sets up a Clenshaw Curtis quadrature rule in multiple dimensions;
- clenshaw_curtis_grid_test
- clenshaw_curtis_rule, a C++ code which implements a Clenshaw Curtis quadrature rule;
- clenshaw_curtis_rule_test
- cnf_io, a C++ code which reads or writes Conjunctive Normal Form (CNF) data files, used to store boolean formulas.
- cnf_io_test
- collatz_recursive, a C++ code which demonstrates recursive programming by considering the simple Collatz 3n+1 problem.
- collatz_recursive_test
- colored_noise, a C++ code which generates samples of noise obeying a 1/f^alpha power law.
- colored_noise_test
- combo, a C++ code 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.
- combo_test
- communicator_mpi, a C++ code which creates new communicators involving a subset of initial set of MPI processes in the default communicator MPI_COMM_WORLD.
- communicator_mpi_test
- compass_search, a C++ code which seeks the minimizer of a scalar function of several variables using compass search, a direct search algorithm that does not use derivatives.
- compass_search_test
- complex_numbers_test, a C++ code which demonstrates some simple features involved in the use of complex numbers in C programming.
- condition, a C++ code which implements methods of computing or estimating the condition number of a matrix.
- condition_test
- continued_fraction, a C++ code which implements some simple algorithms for dealing with simple and generalized continued fractions.
- continued_fraction_test
- cordic, a C++ code which computes a few special functions using the CORDIC algorithm.
- cordic_test
- correlation, a C++ code which contains examples of statistical correlation functions.
- correlation_test
- cosine_transform, a C++ code which demonstrates some simple properties of the discrete cosine transform (DCT).
- cosine_transform_test
- cpp_arrays_test, a C++ code which uses vectors, matrices and tensors.
- cpp_files_test, a C++ code which uses files.
- cpp_intrinsics_test, a C++ code which illustrates the use of intrinsic functions supplied by C++.
- cpp_random_test, a C++ code which uses the C++ random number generator (RNG) functions.
- cpp_test, a C++ code which illustrates some features of the C++ language;
- crc, a C++ code which computes the cyclic redundancy checksum (CRC) as used by the PNG data format.
- crc_test
- cube_arbq_rule, a C++ code which returns quadrature rules, with exactness up to total degree 15, over the interior of the symmetric cube in 3D, by Hong Xiao and Zydrunas Gimbutas.
- cube_arbq_rule_test
- cube_exactness, a C++ code which investigates the polynomial exactness of quadrature rules over the interior of a cube in 3D.
- cube_exactness_test
- cube_felippa_rule, a C++ code which returns Felippa's quadratures rules for approximating integrals over the interior of a cube in 3D.
- cube_felippa_rule_test
- cube_grid, a C++ code which computes a grid of points over the interior of a cube in 3D.
- cube_grid_test
- cube_integrals, a C++ code which returns the exact value of the integral of any monomial over the interior of the unit cube in 3D.
- cube_integrals_test
- cube_monte_carlo, a C++ code which applies a Monte Carlo method to estimate the integral of a function over the interior of the unit cube in 3D.
- cube_monte_carlo_test
- cuda_loop, a C++ code which shows how, in a CUDA program, the user's choice of block and thread factors determines the allocation of tasks to processors.
- cuda_loop_test
- cvt, a C++ code which computes an N point Centroidal Voronoi Tessellation (CVT) in M dimensions;
- cvt_test
- cvt_box, a C++ code which carries out a Centroidal Voronoi Tessellation (CVT) calculation constrained to a box region, with points projected to the boundary;
- cvt_box_test
- cycle_brent, a C++ code 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_brent_test
- cycle_floyd, a C++ code 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.
- cycle_floyd_test
- cyclic_reduction, a C++ code which solves a tridiagonal linear system using cyclic reduction;
- cyclic_reduction_test
- delaunay_tree_2d, a C++ code which computes the Delaunay triangulation of a 2D dataset, by Olivier Devillers.
- delaunay_tree_2d_test
- differ, a C++ code which determines the finite difference method (FDM) coefficients necessary to compute an approximation of given accuracy to a derivative of a given order.
- differ_test
- dijkstra, a C++ code which implements a simple version of Dijkstra's minimum distance algorithm for graphs.
- dijkstra_test
- dijkstra_openmp, a C++ code which uses OpenMP to parallelize a simple example of Dijkstra's minimum distance algorithm for graphs.
- dijkstra_openmp_test
- discrete_pdf_sample_2d, a C++ code which demonstrates how to construct a Probability Density Function (PDF) from a table of sample data, and then to use that PDF to create new samples.
- discrete_pdf_sample_2d_test
- disk_grid, a C++ code which computes grid points within the interior of a disk of user specified radius and center in 2D, using gnuplot() to create an image of the grid.
- disk_grid_test
- disk_monte_carlo, a C++ code which applies a Monte Carlo method to estimate integrals of a function over the interior of the general disk in 2D;
- disk_monte_carlo_test
- disk_rule, a C++ code which computes quadrature rules for the general disk in 2D.
- disk_rule_test
- disk01_integrals, a C++ code which returns the exact value of the integral of any monomial over the interior of the unit disk in 2D.
- disk01_integrals_test
- disk01_monte_carlo, a C++ code which applies a Monte Carlo method to estimate integrals of a function over the interior of the unit disk in 2D;
- disk01_monte_carlo_test
- disk01_quarter_monte_carlo, a C++ code which applies a Monte Carlo method to estimate integrals of a function over the interior of the unit quarter disk in 2D;
- disk01_quarter_monte_carlo_test
- disk01_rule, a C++ code which computes quadrature rules for the unit disk in 2D, that is, the interior of the circle of radius 1 and center (0,0).
- disk01_rule_test
- dislin_test, a C++ code which calls dislin(), which is a package of scientific graphics software.
- divdif, a C++ code which uses divided differences to compute the polynomial interpolant to a given set of data.
- divdif_test
- doomsday, a C++ code 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.
- doomsday_test
- dream, a C++ code which implements the DREAM algorithm for accelerating Markov Chain Monte Carlo (MCMC) convergence using differential evolution, using five user functions to define the problem, by Guannan Zhang.
- dream_test
- duel_simulation, a C++ code which simulates N repetitions of a duel between two players, each of whom has a known firing accuracy.
- duel_simulation_test
- edge, a C++ code which defines some test functions in 1D, 2D and 3D for the detection of edges or discontinuities, and using gnuplot() to display the curves or surfaces.
- edge_test
- ellipse_grid, a C++ code which computes grid points over the interior of an ellipse in 2D.
- ellipse_grid_test
- ellipse_monte_carlo, a C++ code which uses the Monte Carlo method to estimate the value of integrals over the interior of an ellipse in 2D.
- ellipse_monte_carlo_test
- ellipsoid_grid, a C++ code which computes a grid of points over the interior of an ellipsoid in 3D.
- ellipsoid_grid_test
- ellipsoid_monte_carlo, a C++ code which uses the Monte Carlo method to estimate the value of integrals over the interior of an ellipsoid in M dimensions.
- ellipsoid_monte_carlo_test
- elliptic_integral, a C++ code which evaluates complete elliptic integrals of first, second and third kind, using Carlson's elliptic integral functions.
- elliptic_integral_test
- euler, a C++ code which solves one or more ordinary differential equations (ODE) using the forward Euler method.
- euler_test
- exactness, a C++ code which investigates the exactness of quadrature rules that estimate the integral of a function with a density, such as 1, exp(-x) or exp(-x^2), over an interval such as [-1,+1], [0,+oo) or (-oo,+oo).
- exactness_test
- fair_dice_simulation, a C++ code which simulates N tosses of 2 dice, creating graphics files for processing by gnuplot().
- fair_dice_simulation_test
- fastgl, a C++ code which carries out the rapid computation of the Kth value and weight of an N point Gauss-Legendre quadrature rule, by Ignace Bogaert.
- fastgl_test
- faure, a C++ code which computes elements of the Faure M-dimensional Quasi Mont Carlo (QMC) sequence;
- faure_test
- fd_predator_prey, a C++ code which solves a pair of predator prey ordinary differential equations (ODE) using the finite difference method (FDM).
- fd_predator_prey_test
- fd1d_advection_diffusion_steady, a C++ code which applies the finite difference method (FDM) to solve the steady advection diffusion equation v*ux-k*uxx=0 in one spatial dimension, with constant velocity v and diffusivity k, writing graphics files for processing by gnuplot().
- fd1d_advection_diffusion_steady_test
- fd1d_advection_ftcs, a C++ code 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 forward time centered space (FTCS) method, writing graphics files for processing by gnuplot().
- fd1d_advection_ftcs_test
- fd1d_advection_lax, a C++ code 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_test
- fd1d_advection_lax_wendroff, a C++ code 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_advection_lax_wendroff_test
- fd1d_burgers_lax, a C++ code 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_lax_test
- fd1d_burgers_leap, a C++ code 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_burgers_leap_test
- fd1d_bvp, a C++ code which applies the finite difference method (FDM) to a two point boundary value problem (BVP) in one spatial dimension.
- fd1d_bvp_test
- fd1d_heat_explicit, a C++ code which uses the finite difference method (FDM), with explicit time stepping, to solve the time dependent 1D heat equation;
- fd1d_heat_explicit_test
- fd1d_heat_implicit, a C++ code which implements a finite difference method (FDM), implicit in time, of the time dependent 1D heat equation;
- fd1d_heat_implicit_test
- fd1d_heat_steady, a C++ code which implements a finite difference method (FDM) for the steady (time independent) 1D heat equation;
- fd1d_heat_steady_test
- fd1d_wave, a C++ code which applies the finite difference method (FDM) to solve the time-dependent wave equation utt = c * uxx in one spatial dimension.
- fd1d_wave_test
- fd2d_heat_steady, a C++ code which implements a finite difference method (FDM) for the steady (time independent) 2D heat equation;
- fd2d_heat_steady_test
- fem_basis, a C++ code which defines and evaluates 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_basis_test
- fem_io, a C++ code which reads or writes FEM files, which describe the nodes, elements, and function values associated with a finite element method (FEM) model;
- fem_io_test
- fem_to_gmsh, a C++ code which reads a pair of FEM files defining node coordinates and elements, of a 1D, 2D or 3D mesh, namely a file of node coordinates and a file of elements defined by node indices, and creates a gmsh() mesh file.
- fem_to_gmsh_test
- fem_to_triangle, a C++ code which reads FEM files defining a 2D mesh of triangles, namely a file of node coordinates and a file of elements defined by node indices, and creates a corresponding pair of node and element files for use by Jonathan Shewchuk's triangle() program.
- fem_to_triangle_test
- fem_to_xml, a C++ code which reads a pair of FEM files defining node coordinates and elements, of a 1D, 2D or 3D mesh, namely a file of node coordinates and a file of elements defined by node indices, and creates a corresponding XML file for input to DOLFIN or FENICS.
- fem_to_xml_test
- fem1d, a C++ code which applies the finite element method (FEM), using piecewise linear basis functions, to a linear two point boundary value problem (BVP) in 1D;
- fem1d_test
- fem1d_adaptive, a C++ code 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 (BVP) in 1D;
- fem1d_adaptive_test
- fem1d_bvp_linear, a C++ code which applies the finite element method (FEM), with piecewise linear elements, to a two point boundary value problem (BVP) in one spatial dimension, and compares the computed and exact solutions with the L2 and seminorm errors.
- fem1d_bvp_linear_test
- fem1d_bvp_quadratic, a C++ code which applies the finite element method (FEM), with piecewise quadratic elements, to a two point boundary value problem (BVP) in one spatial dimension, and compares the computed and exact solutions with the L2 and seminorm errors.
- fem1d_bvp_quadratic_test
- fem1d_heat_steady, a C++ code which implements the finite element method (FEM) for the steady (time independent) 1D heat equation;
- fem1d_heat_steady_test
- fem1d_lagrange, a C++ code which sets up the matrices and vectors associated with the finite element method (FEM) solution of a boundary value problem (BVP) -u''+u=f(x), using Lagrange basis polynomials.
- fem1d_lagrange_test
- fem1d_nonlinear, a C++ code which applies the finite element method (FEM), with piecewise linear basis functions, to a nonlinear two point boundary value problem (BVP) in 1D;
- fem1d_nonlinear_test
- fem1d_pack, a C++ code which contains utilities for 1D finite element method (FEM) calculations.
- fem1d_pack_test
- fem1d_pmethod, a C++ code which applies the p-method version of the finite element method (FEM) to a linear two point boundary value problem (BVP) in 1D;
- fem1d_pmethod_test
- fem1d_project, a C++ code 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_project_test
- fem1d_sample, a C++ code which samples a scalar or vector finite element function of one variable, defined by FEM files, returning interpolated values at the sample points.
- fem1d_sample_test
- fem2d_bvp_linear, a C++ code which applies the finite element method (FEM), with piecewise linear elements, to a 2D boundary value problem (BVP) in a rectangle, and compares the computed and exact solutions with the L2 and seminorm errors.
- fem2d_bvp_linear_test
- fem2d_bvp_quadratic, a C++ code which applies the finite element method (FEM), with piecewise quadratic elements, to a 2D boundary value problem (BVP) in a rectangle, and compares the computed and exact solutions with the L2 and seminorm errors.
- fem2d_bvp_quadratic_test
- fem2d_bvp_serene, a C++ code which applies the finite element method (FEM), with serendipity elements, to a 2D boundary value problem (BVP) in a rectangle, and compares the computed and exact solutions with the L2 and seminorm errors.
- fem2d_bvp_serene_test
- fem2d_heat, a C++ code which applies the finite element method (FEM) to solve the time-dependent heat equation in an arbitrary triangulated region in 2D;
- fem2d_heat_test
- fem2d_heat_rectangle, a C++ code which applies the finite element method (FEM) to solve the time dependent heat equation on a square in 2D;
- fem2d_heat_rectangle_test
- fem2d_pack, a C++ code which implements various 2D finite element method (FEM) calculations;
- fem2d_pack_test
- fem2d_poisson_rectangle, a C++ code which solves the 2D Poisson equation on a rectangle, using the finite element method (FEM), and piecewise quadratic triangular elements.
- fem2d_poisson_rectangle_test
- fem2d_poisson_rectangle_linear, a C++ code which solves the 2D Poisson equation on a rectangle, using the finite element method (FEM), and piecewise linear triangular elements.
- fem2d_poisson_rectangle_linear_test
- fem2d_project, a C++ code which projects a function f(x,y), given as data, into a given finite element space of piecewise linear triangular elements.
- fem2d_project_test
- fem2d_project_function, a C++ code which projects a function f(x,y), given as a formula, into a given finite element space of piecewise linear triangular elements.
- fem2d_project_function_test
- fem2d_sample, a C++ code which samples a finite element method (FEM) function, defined by FEM files, (three text files describing the nodes, triangles, and coefficients); at arbitrary points.
- fem2d_sample_test
- fem3d_pack, a C++ code which contains utilities for 3D finite element method (FEM) calculations.
- fem3d_pack_test
- fem3d_project, a C++ code 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_project_test
- fem3d_sample, a C++ code which evaluates a finite element function defined on a 3D tetrahedral mesh.
- fem3d_sample_test
- feynman_kac_1d, a C++ code 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_1d_test
- feynman_kac_2d, a C++ code 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_2d_test
- feynman_kac_3d, a C++ code 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.
- feynman_kac_3d_test
- fft_openmp, a C++ code which demonstrates the computation of a Fast Fourier Transform (FFT) in parallel, using OpenMP.
- fft_openmp_test
- fft_serial, a C++ code which computes a Fast Fourier Transform (FFT), and is intended as a starting point for implementing an OpenMP parallel version.
- fft_serial_test
- file_name_sequence, a C++ code which demonstrates ways to generate a sequence of filenames, which can be useful when generating a sequence of still snapshots to be animated later.
- file_name_sequence_test
- files_multiple, a C++ code which demonstrates how a program can open multiple output files at one time, and write data to any one specific file it chooses.
- files_multiple_test
- filon, a C++ code which approximates integrals in which the integrand includes an oscillatory factor of sin(k*x) or cos(k*x).
- filon_test
- filum, a C++ code which performs various operations with files;
- filum_test
- fire_simulation, a C++ code which simulates a forest fire over a rectangular array of trees, starting at a single random location. It is intended as a starting point for the development of a parallel version.
- fire_simulation_test
- floyd, a C++ code which implements Floyd's algorithm for finding the shortest distance between pairs of nodes on a directed graph.
- floyd_test
- fn, a C++ code 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.
- fn_test
- freefem_msh_io, a C++ code which reads and writes files used by the FreeFem++ finite element program to store mesh information.
- freefem_msh_io_test
- fsolve, a C++ code which solves systems of nonlinear equations, inspired by the fsolve() function in MATLAB, and based on the minpack() minimization package.
- fsolve_test
- g++_intrinsics_test, a C++ code which tests or demonstrates some of the intrinsic functions provided by the Gnu G++ compiler for C++.
- g++_quadmath_test, a C++ code which illustrates the use of quadruple precision real arithmetic, as provided on some systems by the Gnu g++ compiler for C++.
- gegenbauer_cc, a C++ code which computes the Gegenbauer weighted integral of a function f(x) using a Clenshaw-Curtis approach.
- gegenbauer_cc_test
- gegenbauer_exactness, a C++ code which tests the polynomial exactness of Gauss-Gegenbauer quadrature rules.
- gegenbauer_exactness_test
- gegenbauer_polynomial, a C++ code which evaluates the Gegenbauer polynomial and associated functions.
- gegenbauer_polynomial_test
- gegenbauer_rule, a C++ code which computes and writes out a Gauss-Gegenbauer quadrature rule of given order.
- gegenbauer_rule_test
- gen_hermite_exactness, a C++ code which tests the polynomial exactness of generalized Gauss-Hermite quadrature rules.
- gen_hermite_exactness_test
- gen_hermite_rule, a C++ code which computes and writes out a generalized Gauss-Hermite quadrature rule of given order and parameter value.
- gen_hermite_rule_test
- gen_laguerre_exactness, a C++ code which tests the polynomial exactness of generalized Gauss-Laguerre quadrature rules.
- gen_laguerre_exactness_test
- gen_laguerre_rule, a C++ code which computes and writes out a generalized Gauss-Laguerre quadrature rule of given order and parameter value.
- gen_laguerre_rule_test
- geometry, a C++ code which performs geometric calculations in 2, 3 and M dimensional space, including the computation of angles, areas, containment, distances, intersections, lengths, and volumes.
- geometry_test
- glomin, a C++ code which finds a global minimum of a scalar function of a scalar variable, without the use of derivative information, by Richard Brent.
- glomin_test
- gmsh_io, a C++ code which reads or writes some of the files created by gmsh() for 1D, 2D or 3D meshes used by the finite element method (FEM).
- gmsh_io_test
- gmsh_to_fem, a C++ code which reads a mesh data file created by gmsh() and writes a pair of node and element files that correspond to the FEM format.
- gmsh_to_fem_test
- gnuplot_test, a C++ code which uses gnuplot() to create various kinds of plots.
- gprof_test
- graphics_dislin_test, a C++ code which creates various kinds of plots using the dislin() library.
- graphics_gnuplot_test, a C++ code which creates various kinds of plots using gnuplot().
- grf_io, a C++ code which reads or writes a GRF file;
- grf_io_test
- gsl_test
- haar, a C++ code which computes the Haar transform of data.
- haar_test
- halton, a C++ code which computes elements of a Halton Quasi Monte Carlo (QMC) sequence, using a simple interface.
- halton_test
- hammersley, a C++ code which computes elements of a Hammersley Quasi Monte Carlo (QMC) sequence, using a simple interface.
- hammersley_test
- hankel_cholesky, a C++ code which computes the upper Cholesky factor R of a symmetric positive definite (SPD) matrix so that H = R' * R.
- hankel_cholesky_test
- hankel_spd, a C++ code which computes a lower triangular matrix L which is the Cholesky factor of a symmetric positive definite (SPD) Hankel matrix H, that is, H = L * L'.
- hankel_spd_test
- hdf5_test
- heat_mpi, a C++ code which uses the Message Passing Interface (MPI) to solve the 1D time dependent heat equation.
- heat_mpi_test
- heated_plate, a C++ code 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_test
- heated_plate_openmp, a C++ code which solves the steady (time independent) heat equation in a 2D rectangular region, using OpenMP to run in parallel.
- heated_plate_openmp_test
- hello, a C++ code which prints out "Hello, world!".
- hello_test
- hello_mpi, a C++ code which prints out "Hello, world!" using the Message Passing Interface (MPI) for parallel programming.
- hello_mpi_test
- hello_openmp, a C++ code which prints out "Hello, world!" within the OpenMP parallel programming environment.
- hello_openmp_test
- hermite_cubic, a C++ code which computes the value, derivatives or integral of a Hermite cubic polynomial, or manipulates an interpolating function made up of piecewise Hermite cubic polynomials.
- hermite_cubic_test
- hermite_exactness, a C++ code which tests the polynomial exactness of Gauss-Hermite quadrature rules for approximating the integral of a function with density exp(-x^2) over the interval (-oo,+oo).
- hermite_exactness_test
- hermite_integrands, a C++ code which defines test integrands for Hermite integrals with density exp(-x^2) over the interval (-oo,+oo).
- hermite_integrands_test
- hermite_interpolant, a C++ code which computes the Hermite interpolant, a polynomial that matches function values and derivatives.
- hermite_interpolant_test
- hermite_polynomial, a C++ code which evaluates the physicist's Hermite polynomial, the probabilist's Hermite polynomial, the Hermite function, and related functions.
- hermite_polynomial_test
- hermite_product_polynomial, a C++ code which defines Hermite product polynomials, creating a multivariate polynomial as the product of univariate Hermite polynomials.
- hermite_product_polynomial_test
- hermite_rule, a C++ code which computes and writes out a Gauss-Hermite quadrature rule for approximating the integral of a function with density exp(-x^2) over the interval (-oo,+oo).
- hermite_rule_test
- hexdump, a C++ code which produces a hexadecimal dump of a file;
- hexdump_test
- high_card_simulation, a C++ code which simulates a situation in which you see the cards in a deck one by one, and must select the one you think is the highest and stop; the program uses gnuplot() to create a graph of the results.
- high_card_simulation_test
- hilbert_curve, a C++ code which computes the sequence of discrete Hilbert curves whose limit is a space-filling curve.
- hilbert_curve_test
- histogram_data_2d_sample, a C++ code 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.
- histogram_data_2d_sample_test
- htmlindex, a C++ code which creates a skeleton HTML page for a marked-up C, C++ or FORTRAN file;
- htmlindex_test
- humps, a C++ code which evaluates the humps function, its first and second derivatives and its antiderivative.
- humps_test
- hyperball_integrals, a C++ code which returns the exact value of the integral of any monomial over the interior of the unit hyperball in M dimensions.
- hyperball_integrals_test
- hyperball_monte_carlo, a C++ code which applies a Monte Carlo method to estimate the integral of a function over the interior of the unit ball in M dimensions;
- hyperball_monte_carlo_test
- hyperball_volume_monte_carlo, a C++ code which applies a Monte Carlo procedure to estimate the volume of the unit ball in M dimensions;
- hyperball_volume_monte_carlo_test
- hypercube_exactness, a C++ code which tests the polynomial exactness of a quadrature rule over the interior of the unit hypercube in M dimensions.
- hypercube_exactness_test
- hypercube_grid, a C++ code which computes a grid of points over the interior of a hypercube in M dimensions.
- hypercube_grid_test
- hypercube_integrals, a C++ code which returns the exact value of the integral of any monomial over the interior of the unit hypercube in M dimensions.
- hypercube_integrals_test
- hypercube_monte_carlo, a C++ code which applies a Monte Carlo method to estimate the integral of a function over the interior of the unit hypercube in M dimensions.
- hypercube_monte_carlo_test
- hypersphere_integrals, a C++ code which returns the exact value of the integral of any monomial over the surface of the unit hypersphere in M dimensions.
- hypersphere_integrals_test
- hypersphere_monte_carlo, a C++ code which applies a Monte Carlo method to estimate the integral of a function on the surface of the unit sphere in M dimensions;
- hypersphere_monte_carlo_test
- hypersphere_properties, a C++ code which carries out various operations for an M dimensional hypersphere, including converting between Cartesian and spherical coordinates, stereographic projection, sampling the surface of the sphere, and computing the surface area and volume.
- hypersphere_properties_test
- i2_binary_to_ascii, a C++ code which converts a file of short int (I2) data from binary to ASCII format.
- i2_binary_to_ascii_test
- i4lib, a C++ code which contains many utility routines, using single precision integer (I4) arithmetic.
- i4lib_test
- i8lib, a C++ code which contains many utility routines, using double precision integer (I8) arithmetic.
- i8lib_test
- image_components, a C++ code which seeks the connected nonzero or nonblack components of an image or integer array.
- image_components_test
- image_denoise, a C++ code which applies simple filtering operations to a noisy image.
- image_denoise_test
- image_edge, a C++ code which demonstrates a simple procedure for edge detection in images.
- image_edge_test
- index, a C++ code 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.
- index_test
- interp, a C++ code which parameterizes and interpolates data;
- interp_test
- isbn, a C++ code which determines the check digit for an International Standard Book Number (ISBN) or report whether a given ISBN is valid.
- isbn_test
- ising_2d_simulation, a C++ code 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, using gnuplot() to display the initial and final configurations.
- ising_2d_simulation_test
- jacobi, a C++ code which implements the Jacobi iteration for linear systems.
- jacobi_test
- jacobi_eigenvalue, a C++ code which implements the Jacobi iteration for the iterative determination of the eigenvalues and eigenvectors of a real symmetric matrix.
- jacobi_eigenvalue_test
- jacobi_exactness, a C++ code which tests the polynomial exactness of rules for integrals with a Jacobi weight function.
- jacobi_exactness_test
- jacobi_openmp, a C++ code which illustrates the use of the OpenMP application program interface to parallelize a Jacobi iteration solving A*x=b.
- jacobi_openmp_test
- jacobi_polynomial, a C++ code which evaluates the Jacobi polynomial and associated functions.
- jacobi_polynomial_test
- jacobi_rule, a C++ code which computes and writes out a Gauss-Jacobi quadrature rule of given order.
- jacobi_rule_test
- kmeans, a C++ code which handles the K-Means problem, which organizes a set of N points in M dimensions into K clusters;
- kmeans_test
- knapsack_01, a C++ code which uses brute force to solve small versions of the 0/1 knapsack problem;
- knapsack_01_test
- kronrod, a C++ code which computes a Gauss and Gauss-Kronrod pair of quadrature rules of arbitrary order, by Robert Piessens, Maria Branders.
- kronrod_test
- lagrange_approx_1d, a C++ code 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_approx_1d_test
- lagrange_interp_1d, a C++ code 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_1d_test
- lagrange_interp_2d, a C++ code 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_2d_test
- lagrange_interp_nd, a C++ code 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).
- lagrange_interp_nd_test
- lagrange_nd, a C++ code which implements several multivariate Lagrange interpolation schemes developed by Tomas Sauer.
- lagrange_nd_test
- laguerre_exactness, a C++ code which tests the polynomial exactness of Gauss-Laguerre quadrature rules for integrating functions with density exp(-x) over [0,+oo).
- laguerre_exactness_test
- laguerre_integrands, a C++ code which implements test problems for for Gauss-Laguerre type quadrature rules for approximating the integral of a function with density exp(-x) over the interval [0,+oo).
- laguerre_integrands_test
- laguerre_polynomial, a C++ code which evaluates the Laguerre polynomial, the generalized Laguerre polynomials, and the Laguerre function.
- laguerre_polynomial_test
- laguerre_rule, a C++ code which computes and writes out a Gauss-Laguerre quadrature rule for approximating the integral of a function with density exp(-x) over the interval [0,+oo).
- laguerre_rule_test
- laplacian, a C++ code which carries out computations related to the discrete Laplacian operator, including full or sparse evaluation, evaluation for unequally spaced data sampling points, application to a set of data samples, solution of associated linear systems, eigenvalues and eigenvectors, and extension to 2D and 3D geometry.
- laplacian_test
- latin_center, a C++ code which computes Latin Center Squares of N points in M dimensions;
- latin_center_test
- latin_cover, a C++ code which produces N Latin squares which cover an NxN square, or NxN Latin cubes which cover an NxNxN cube.
- latin_cover_test
- latin_edge, a C++ code which computes Latin Edge Squares of N points in M dimensions;
- latin_edge_test
- latin_random, a C++ code which computes Latin Random Squares of N points in M dimensions;
- latin_random_test
- latinize, a C++ code which adjusts N points in M dimensions to form a Latin Hypercube;
- latinize_test
- lattice_rule, a C++ code which returns lattice rules for M-dimensional integration;
- lattice_rule_test
- lcvt, a C++ code which creates a Latinized Centroidal Voronoi Tessellation (CVT);
- lcvt_test
- lebesgue, a C++ code which is given a set of nodes in 1D, and plots the Lebesgue function, and estimates the Lebesgue constant, which measures the maximum magnitude of the potential error of Lagrange polynomial interpolation.
- lebesgue_test
- legendre_exactness, a C++ code which tests the monomial exactness of quadrature rules for the Legendre problem of integrating a function with density 1 over the interval [-1,+1].
- legendre_exactness_test
- legendre_polynomial, a C++ code which evaluates the Legendre polynomial and associated functions.
- legendre_polynomial_test
- legendre_product_polynomial, a C++ code which defines Legendre product polynomials, creating a multivariate polynomial as the product of univariate Legendre polynomials.
- legendre_product_polynomial_test
- legendre_rule, a C++ code which writes out a Gauss-Legendre quadrature rule of given order.
- legendre_rule_test
- legendre_rule_fast, a C++ code which uses a fast (order N) algorithm to compute a Gauss-Legendre quadrature rule of given order.
- legendre_rule_fast_test
- legendre_shifted_polynomial, a C++ code which evaluates the shifted Legendre polynomial, with the domain [0,1].
- legendre_shifted_polynomial_test
- levenshtein, a C++ code which returns the Levenshtein distance between two strings.
- levenshtein_test
- life_serial, a C++ code which computes a few steps of the evolution of John Conway's Game of Life, intended as a starting point for implementing a parallel version.
- life_serial_test
- line_cvt_lloyd, a C++ code which applies Lloyd's iteration repeatedly to a set of N points, to compute a Centroidal Voronoi Tessellation (CVT) over the interior of a line segment in 1D.
- line_cvt_lloyd_test
- line_fekete_rule, a C++ code which estimates the location of N Fekete points, for polynomial interpolation or quadrature, over the interior of a line segment in 1D.
- line_fekete_rule_test
- line_felippa_rule, a C++ code which returns Felippa's quadratures rules for approximating integrals over the interior of a line segment in 1D.
- line_felippa_rule_test
- line_grid, a C++ code which computes a grid of points over the interior of a line segment in 1D.
- line_grid_test
- line_integrals, a C++ code which returns the exact value of the integral of any monomial over the length of the unit line in 1D.
- line_integrals_test
- line_monte_carlo, a C++ code which applies a Monte Carlo method to estimate the integral of a function over the length of the unit line in 1D.
- line_monte_carlo_test
- line_ncc_rule, a C++ code which computes a Newton Cotes (Closed) quadrature rule, using equally spaced points, over the interior of a line segment in 1D.
- line_ncc_rule_test
- line_nco_rule, a C++ code which computes a Newton Cotes Open (NCO) quadrature rule, using equally spaced points, over the interior of a line segment in 1D.
- line_nco_rule_test
- linpack_bench, a C++ code which is the LINPACK benchmark program;
- linpack_bench_test
- linpack_d, a C++ code which factors matrices in a variety of formats, and solving the associated linear systems, using double precision real arithmetic; it includes a function for the singular value decomposition (SVD) of a rectangular matrix, by Jack Dongarra, Cleve Moler, Jim Bunch, Pete Stewart.
- linpack_d_test
- linpack_z, a C++ code which factors matrices in a variety of formats, and solving the associated linear systems, using double precision complex arithmetic; it includes a function for the singular value decomposition (SVD) of a rectangular matrix, by Jack Dongarra, Cleve Moler, Jim Bunch, Pete Stewart.
- linpack_z_test
- llsq, a C++ code which solves the simple linear least squares (LLS) problem of finding the formula of a straight line y=a*x+b which minimizes the root-mean-square (RMS) error to a set of N data points.
- llsq_test
- lobatto_polynomial, a C++ code which evaluates Lobatto polynomials similar to Legendre polynomials except that they are 0 at both endpoints.
- lobatto_polynomial_test
- local_min, a C++ code which finds a local minimum of a scalar function of a scalar variable, without the use of derivative information, by Richard Brent.
- local_min_test
- local_min_rc, a C++ code which finds a local minimum of a scalar function of a scalar variable, without the use of derivative information, using reverse communication (RC), by Richard Brent.
- local_min_rc_test
- log_normal, a C++ code which returns quantities related to the log normal Probability Distribution Function (PDF).
- log_normal_test
- log_normal_truncated_ab, a C++ code which returns quantities related to the log normal Probability Distribution Function (PDF) truncated to the interval [A,B].
- log_normal_truncated_ab_test
- lorenz_ode, a C++ code which solves the Lorenz system of ordinary differential equations (ODE), creating output files that can be displayed by gnuplot().
- lorenz_ode_test
- luhn, a C++ code which computes the Luhn check digit for a string, or verify a string, as used for error detection in credit card numbers.
- luhn_test
- machar, a C++ code which dynamically computes the values of various machine constants, by William Cody;
- machar_test
- machine, a C++ code which returns tabulated values of the constants associated with computer arithmetic.
- machine_test
- makefiles_test, a C++ code which illustrates the use of makefiles for maintaining a software project;
- mandelbrot, a C++ code which generates an ASCII Portable Pixel Map (PPM) image of the Mandelbrot set;
- mandelbrot_test
- mandelbrot_openmp, a C++ code which generates an ASCII Portable Pixel Map (PPM) image of the Mandelbrot fractal set, using OpenMP for parallel execution;
- mandelbrot_openmp_test
- mandelbrot_ppm, a C++ code which generates a binary Portable Pixel Map (PPM) image of the Mandelbrot set, by Eric Weeks.
- mandelbrot_ppm_test
- mario, a C++ code which creates a sort of "needlepoint" image of Mario, as an array of colored squares.
- mario_test
- matrix_exponential, a C++ code which demonstrates some simple approaches to the problem of computing the exponential of a matrix.
- matrix_exponential_test
- md, a C++ code which carries out a molecular dynamics simulation, and is intended as a starting point for implementing an OpenMP parallel version.
- md_test
- md_openmp, a C++ code which carries out a molecular dynamics simulation in parallel using OpenMP.
- md_openmp_test
- memory, a C++ code which declares and uses a sequence of larger and larger arrays, to see what the memory limits are on a given computer.
- memory_test
- mesh_bandwidth, a C++ code which returns the geometric bandwidth associated with a mesh of elements of any order and in a space of M dimensions.
- mesh_bandwidth_test
- metis_test
- mgmres, a C++ code which applies the restarted Generalized Minimum Residual (GMRES) algorithm to solve a sparse linear system, using Compressed Row Storage (CRS) or sparse triplet (ST) format, by Lili Ju.
- mgmres_test
- mgs, a C++ code 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.
- mgs_test
- midpoint_explicit, a C++ code which solves one or more ordinary differential equations (ODE) using the (explicit) midpoint method, also called the modified Euler method.
- midpoint_explicit_test
- midpoint_fixed, a C++ code which solves one or more ordinary differential equations (ODE) using the (implicit) midpoint method, using a fixed point iteration to solve the associated nonlinear equation.
- midpoint_fixed_test
- minesweeper, a C++ code which is a text-based implementation of the Minesweeper game, by Detelina Stoyanova.
- minpack, a C++ code which solves a system of nonlinear equations, by Jorge More, Danny Sorenson, Burton Garbow, Kenneth Hillstrom.
- minpack_test
- mixed_test, a C++ code which illustrates mixed language programming in which the main program is written in C++;
- monomial, a C++ code which enumerates, lists, ranks, unranks and randomizes multivariate monomials in a space of M dimensions, with total degree less than N, equal to N, or lying within a given range.
- monomial_test
- monomial_value, a C++ code which evaluates a monomial in M dimensions.
- monomial_value_test
- mpi_test, a C++ code which illustrates the use of the Message Passing Interface (MPI) for parallel programming;
- multigrid_poisson_1d, a C++ code which applies the multigrid method to a discretized version of the 1D Poisson equation.
- multigrid_poisson_1d_test
- multitask_mpi, a C++ code which demonstrates how to multitask, that is, to execute several unrelated and distinct tasks simultaneously, using MPI for parallel execution.
- multitask_mpi_test
- multitask_openmp, a C++ code which demonstrates how to multitask, that is, to execute several unrelated and distinct tasks simultaneously, using OpenMP for parallel execution.
- multitask_openmp_test
- mxm_openmp, a C++ code which computes a dense matrix product C=A*B, using OpenMP for parallel execution.
- mxm_openmp_test
- mxm_serial, a C++ code which sets up a matrix multiplication problem A=B*C, intended as a starting point for implementing a parallel version.
- mxm_serial_test
- naca, a C++ code which takes the parameters of certain NACA airfoils and returns the coordinates of a sequence of points that outline the wing shape. The points can be plotted with gnuplot().
- naca_test
- nas, a C++ code which runs the NASA kernel benchmark.
- nas_test
- navier_stokes_2d_exact, a C++ code which evaluates an exact solution to the incompressible time-dependent Navier-Stokes equations over an arbitrary domain in 2D.
- navier_stokes_2d_exact_test
- navier_stokes_3d_exact, a C++ code which evaluates an exact solution to the incompressible time-dependent Navier-Stokes equations over an arbitrary domain in 3D.
- navier_stokes_3d_exact_test
- nearest_interp_1d, a C++ code which interpolates a set of data using a piecewise constant interpolant defined by the nearest neighbor criterion, creating graphics files for processing by gnuplot().
- nearest_interp_1d_test
- neighbors_to_metis_graph, a C++ code 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();
- neighbors_to_metis_graph_test
- netcdf_test
- newton_interp_1d, a C++ code which finds a polynomial interpolant to data using Newton divided differences.
- newton_interp_1d_test
- niederreiter, a C++ code which computes elements of the Niederreiter Quasi Mont Carlo (QMC) M-dimensional sequence [arbitrary base];
- niederreiter_test
- niederreiter2, a C++ code which computes elements of the Niederreiter M-dimensional Quasi Mont Carlo (QMC) sequence using base 2;
- niederreiter2_test
- nint_exactness_mixed, a C++ code which measures the polynomial exactness of an M-dimensional quadrature rule based on a mixture of 1D quadrature rule factors.
- nint_exactness_mixed_test
- nintlib, a C++ code which implements approximate integration (quadrature) in M dimensions;
- nintlib_test
- normal, a C++ code which implements a random number generator (RNG) for normally distributed values.
- normal_test
- obj_io, a C++ code which reads and writes the data in an OBJ 3D graphics file;
- obj_io_test
- ode, a C++ code which solves a system of ordinary differential equations (ODE), by Shampine and Gordon;
- ode_test
- openmp_test, a C++ code which uses the OpenMP interface for parallel programming shared memory system;
- optimize_test, a C++ code which looks at the use of automatic optimization by the compiler to make the executable code run faster.
- ornstein_uhlenbeck, a C++ code 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().
- ornstein_uhlenbeck_test
- owen, a C++ code which evaluates Owen's T function;
- owen_test
- padua, a C++ code which returns the coordinates of the 2D Padua points, as well as interpolation weights or quadrature weights, and images of the points in gnuplot() graphics files.
- padua_test
- partition_problem, a C++ code which seeks solutions of the partition problem, splitting a set of integers into two subsets with equal sum.
- partition_problem_test
- patterson_rule, a C++ code which returns the points and weights of a 1D Gauss-Patterson quadrature rule of order 1, 3, 7, 15, 31, 63, 127, 255 or 511.
- patterson_rule_test
- pbma_io, a C++ code which reads or writes an ASCII Portable Bit Map (PBM) graphics file;
- pbma_io_test
- pbmb_io, a C++ code which reads or writes a binary Portable Bit Map (PBM) graphics file;
- pbmb_io_test
- pce_burgers, a C++ code 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_burgers_test
- pce_ode_hermite, a C++ code 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.
- pce_ode_hermite_test
- pdflib, a C++ code which evaluates Probability Density Functions (PDF) and produces random samples from them, including beta, binomial, chi, exponential, gamma, inverse chi, inverse gamma, multinomial, normal, scaled inverse chi, and uniform.
- pdflib_test
- persistence, a C++ code which demonstrates how to initialize, get, or set data stored within a function, using static/persistent/saved/local memory.
- persistence_test
- pgma_io, a C++ code which reads or writes an ASCII Portable Gray Map (PGM) 2D graphics file;
- pgma_io_test
- pgmb_io, a C++ code which reads or writes a binary Portable Gray Map (PGM) 2D graphics file;
- pgmb_io_test
- pink_noise, a C++ code which computes a pink noise signal obeying a 1/f power law.
- pink_noise_test
- poisson_openmp, a C++ code 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_openmp_test
- poisson_serial, a C++ code 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_serial_test
- poisson_simulation, a C++ code which simulates a Poisson process in which events randomly occur with an average waiting time of Lambda, creating graphics files for processing by gnuplot().
- poisson_simulation_test
- polpak, a C++ code which evaluates a variety of mathematical functions, including Chebyshev, Gegenbauer, Hermite, Jacobi, Laguerre, Legendre polynomials, and the Collatz sequence.
- polpak_test
- polygon_grid, a C++ code which generates a grid of points over the interior of a polygon in 2D.
- polygon_grid_test
- polygon_integrals, a C++ code which returns the exact value of the integral of any monomial over the interior of a polygon in 2D.
- polygon_integrals_test
- polygon_monte_carlo, a C++ code which applies a Monte Carlo method to estimate the integral of a function over the interior of a polygon in 2D.
- polygon_monte_carlo_test
- polygon_properties, a C++ code which computes properties of an arbitrary polygon in the plane, defined by a sequence of vertices, including interior angles, area, centroid, containment of a point, convexity, diameter, distance to a point, inradius, lattice area, nearest point in set, outradius, uniform sampling.
- polygon_properties_test
- polygon_triangulate, a C++ code which triangulates a possibly nonconvex polygon in 2D, and which can use gnuplot() to display the external edges and internal diagonals of the triangulation.
- polygon_triangulate_test
- polynomial, a C++ code which adds, multiplies, differentiates, evaluates and prints multivariate polynomials in a space of M dimensions.
- polynomial_test
- polyominoes, a C++ code which provides some utilities for manipulating polyominoes.
- polyominoes_test
- power_method, a C++ code which carries out the power method for finding a dominant eigenvalue and its eigenvector.
- power_method_test
- power_rule, a C++ code which constructs a power rule, that is, a product quadrature rule from identical 1D factor rules.
- power_rule_test
- ppma_io, a C++ code which reads or writes an ASCII Portable Pixel Map (PPM) graphics file;
- ppma_io_test
- ppmb_io, a C++ code which reads or writes a binary Portable Pixel Map (PPM) 2D pixel map graphics file;
- ppmb_io_test
- praxis, a C++ code which minimizes a scalar function of several variables, without requiring derivative information, by Richard Brent.
- praxis_test
- predator_prey_ode, a C++ code which solves a time-dependent predator-prey system of ordinary differential equations (ODE).
- predator_prey_ode_test
- prime, a C++ code which counts the number of primes between 1 and N, and is intended as a starting point for a parallel version.
- prime_test
- prime_mpi, a C++ code which counts the number of primes between 1 and N, using MPI for parallel execution.
- prime_mpi_test
- prime_openmp, a C++ code which counts the number of primes between 1 and N, using OpenMP for parallel execution.
- prime_openmp_test
- prob, a C++ code which evaluates, samples, inverts, and characterizes a number of Probability Density Functions (PDF) and Cumulative Density Functions (CDF), including anglit, arcsin, benford, birthday, bernoulli, beta_binomial, beta, binomial, bradford, burr, cardiod, cauchy, chi, chi squared, circular, cosine, deranged, dipole, dirichlet mixture, discrete, empirical, english sentence and word length, error, exponential, extreme values, f, fisk, folded normal, frechet, gamma, generalized logistic, geometric, gompertz, gumbel, half normal, hypergeometric, inverse gaussian, laplace, levy, logistic, log normal, log series, log uniform, lorentz, maxwell, multinomial, nakagami, negative binomial, normal, pareto, planck, poisson, power, quasigeometric, rayleigh, reciprocal, runs, sech, semicircular, student t, triangle, uniform, von mises, weibull, zipf.
- prob_test
- product_rule, a C++ code which creates an M-dimensional quadrature rule from a product of distinct 1d quadrature rules;
- product_rule_test
- pwl_approx_1d, a C++ code which approximates a set of data using a piecewise linear function.
- pwl_approx_1d_test
- pwl_interp_1d, a C++ code which interpolates a set of data using a piecewise linear function, creating graphics files for processing by gnuplot().
- pwl_interp_1d_test
- pwl_interp_2d, a C++ code which evaluates a piecewise linear interpolant to data defined on a regular 2D grid.
- pwl_interp_2d_test
- pwl_interp_2d_scattered, a C++ code which evaluates a piecewise linear interpolant to data which is available at an irregularly arranged set of points.
- pwl_interp_2d_scattered_test
- pwl_product_integral, a C++ code which calculates the exact value of the integral of the product of two piecewise linear functions f(x) and g(x).
- pwl_product_integral_test
- pyramid_exactness, a C++ code which investigates the polynomial exactness of a quadrature rule over the interior of the unit pyramid in 3D.
- pyramid_exactness_test
- pyramid_felippa_rule, a C++ code which returns Felippa's quadratures rules for approximating integrals over the interior of a pyramid in 3D.
- pyramid_felippa_rule_test
- pyramid_grid, a C++ code which computes a grid of points over the interior of the unit pyramid in 3D;
- pyramid_grid_test
- pyramid_integrals, a C++ code which returns the exact value of the integral of any monomial over the interior of the unit pyramid in 3D.
- pyramid_integrals_test
- pyramid_monte_carlo, a C++ code which applies a Monte Carlo method to estimate integrals of a function over the interior of the unit pyramid in 3D;
- pyramid_monte_carlo_test
- pyramid_rule, a C++ code which computes a quadrature rule for a pyramid.
- pyramid_rule_test
- qr_solve, a C++ code which computes the linear least squares (LLS) solution of a system A*x=b.
- qr_solve_test
- quad_mesh, a C++ code which handles meshes of quadrilaterals over a 2D region;
- quad_mesh_test
- quad_mesh_rcm, a C++ code which computes the Reverse Cuthill-McKee (RCM) reordering for nodes in a mesh of 4-node quadrilaterals.
- quad_mesh_rcm_test
- quad_mpi, a C++ code which applies a quadrature rule to estimate an integral, working in parallel by using the Message Passing Interface (MPI);
- quad_mpi_test
- quad_openmp, a C++ code which applies a quadrature rule to estimate an integral over a 1D interval, using OpenMP for parallel execution;
- quad_openmp_test
- quad_serial, a C++ code which applies a quadrature rule to estimate an integral, intended as a starting point for parallelization exercises.
- quad_serial_test
- quad2d_openmp, a C++ code which applies a product quadrature rule to estimate an integral over a 2D rectangle, using OpenMP for parallel execution;
- quad2d_openmp_test
- quad2d_serial, a C++ code which applies a product quadrature rule to estimate an integral over a 2D rectangle, intended as a starting point for parallelization exercises.
- quad2d_serial_test
- quadmom, a C++ code which computes a Gaussian quadrature rule for a weight function rho(x) based on the Golub-Welsch procedure that only requires knowledge of the moments of rho(x).
- quadmom_test
- quadrature_golub_welsch, a C++ code which computes the points and weights of a Gaussian quadrature rule using the Golub-Welsch procedure.
- quadrature_golub_welsch_test
- quadrature_least_squares, a C++ code which computes weights for sub-interpolatory quadrature rules, that is, it estimates integrals by integrating a polynomial that approximates the function data in a least squares sense.
- quadrature_least_squares_test
- quadrature_weights_vandermonde, a C++ code which computes the weights of a quadrature rule using the Vandermonde matrix, assuming that the points have been specified.
- quadrature_weights_vandermonde_test
- quadrature_weights_vandermonde_2d, a C++ code which computes the weights of a 2D quadrature rule using the Vandermonde matrix, assuming that the points have been specified.
- quadrature_weights_vandermonde_2d_test
- quadrule, a C++ code which defines quadrature rules for approximating integrals;
- quadrule_test
- quality, a C++ code which computes quality of distribution measures for pointsets in M dimensions;
- quality_test
- quaternions, a C++ code which carries out some simple arithmetic operations for quaternions.
- quaternions_test
- r83, a C++ code which contains linear algebra routines for r83 matrices (real, 64 bit, tridiagonal 3xN format).
- r83_test
- r8col, a C++ code which contains utility routines for R8COL's, that is, double precision real MxN arrays, considered as N column vectors, each of length M. The data may be thought of as a matrix of multiple columns, and many operations will be carried out columnwise.
- r8col_test
- r8ge, a C++ code which contains linear algebra routines for R8GE matrices (real, 64 bit, General format).
- r8ge_test
- r8lib, a C++ code which contains many utility routines, using double precision real (R8) arithmetic.
- r8lib_test
- r8lt, a C++ code which contains linear algebra routines for R8LT matrices (real, 64 bit, Lower Triangular).
- r8lt_test
- r8poly, a C++ code which contains a number of utilities for polynomials with R8 coefficients, that is, using double precision or 64 bit real arithmetic.
- r8poly_test
- r8row, a C++ code which contains utility routines for R8ROW's, that is, double precision real MxN arrays, considered as M row vectors, each of length N. The data may be thought of as a matrix of multiple rows, and many operations will be carried out rowwise.
- r8row_test
- r8st, a C++ code which contains linear algebra routines for R8ST matrices (real, 64 bit, Sparse Triplet (ST) (I,J,AIJ) format).
- r8st_test
- r8ut, a C++ code which contains linear algebra routines for R8UT matrices (real, 64 bit, Upper Triangular).
- r8ut_test
- rand48_test, a C++ code which demonstrates the use of the RAND48 family of random number generators (RNG) available in the C/C++ standard library.
- randlc, a C++ code which implements a random number generator (RNG) used by the NAS Benchmark programs.
- randlc_test
- random_data, a C++ code which uses a random number generator (RNG) to sample points distribute according to various probability density functions (PDF), for various spatial dimensions, and geometries, including the M-dimensional cube, ellipsoid, simplex and sphere.
- random_data_test
- random_mpi, a C++ code which demonstrates one way to generate the same sequence of random numbers for both sequential execution and parallel execution under MPI.
- random_mpi_test
- random_openmp, a C++ code which illustrates how a parallel program using OpenMP can generate multiple distinct streams of random numbers.
- random_openmp_test
- random_sorted, a C++ code which generates vectors of random values which are already sorted.
- random_sorted_test
- ranlib, a C++ code which produces random samples from Probability Density Functions (PDF), 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.
- ranlib_test
- rbf_interp_1d, a C++ code which defines and evaluates radial basis function (RBF) interpolants to 1D data.
- rbf_interp_1d_test
- rbf_interp_2d, a C++ code which defines and evaluates radial basis function (RBF) interpolants to 2D data.
- rbf_interp_2d_test
- rbf_interp_nd, a C++ code which defines and evaluates radial basis function (RBF) interpolants to M-dimensional data.
- rbf_interp_nd_test
- rcm, a C++ code which applies the Reverse Cuthill McKee (RCM) algorithm for reordering the nodes of a graph, and reducing the bandwidth of a corresponding sparse matrix;
- rcm_test
- reactor simulation, a C++ code 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."
- reactor simulation_test
- ring_mpi, a C++ code which uses the MPI parallel programming environment, and measures the time necessary to copy a set of data around a ring of processes.
- ring_mpi_test
- rk4, a C++ code which applies the fourth order Runge-Kutta (RK) algorithm to estimate the solution of an ordinary differential equation (ODE) at the next time step.
- rk4_test
- rkf45, a C++ code which applies the Runge-Kutta-Fehlberg (RKF) algorithm to solve a system of ordinary differential equations (ODE);
- rkf45_test
- rnglib, a C++ code which implements a random number generator (RNG) with splitting facilities, allowing multiple independent streams to be computed, by L'Ecuyer and Cote.
- rnglib_test
- root_rc, a C++ code which seeks a solution of a scalar nonlinear equation f(x)=0, using reverse communication (RC), by Gaston Gonnet.
- root_rc_test
- roots_rc, a C++ code which seeks solutions of a system of nonlinear equations, using reverse communication (RC), by Gaston Gonnet.
- roots_rc_test
- rot13, a C++ code which makes a copy of a file which has been encoded using the ROT13 coding and a ROT5 coding for digits.
- rot13_test
- row_echelon_integer, a C++ code which carries out the exact computation of the integer row echelon form (IREF) and integer reduced row echelon form (IRREF) of an integer matrix.
- row_echelon_integer_test
- satisfy, a C++ code which demonstrates, for a particular circuit, an exhaustive search for solutions of the circuit satisfy problem.
- satisfy_test
- satisfy_mpi, a C++ code 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_mpi_test
- satisfy_openmp, a C++ code which demonstrates, for a particular circuit, an exhaustive search for solutions of the circuit satisfy problem, using OpenMP for parallel execution.
- satisfy_openmp_test
- schedule_openmp, a C++ code which demonstrates the default, static, and dynamic methods of scheduling loop iterations in OpenMP to avoid work imbalance.
- schedule_openmp_test
- sde, a C++ code which illustrates properties of stochastic ordinary differential equations (SODE), and common algorithms for their analysis, including the Euler method, the Euler-Maruyama method, and the Milstein method, making graphics files for processing and display by gnuplot(), by Desmond Higham;
- sde_test
- search_mpi, a C++ code which searches integers between A and B for a solution J such that F(J)=C, carried out in parallel using MPI.
- search_mpi_test
- search_serial, a C++ code 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.
- search_serial_test
- set_theory, a C++ code which demonstrates various set theoretic operations using several models of a set.
- set_theory_test
- sftpack, a C++ code which implements the slow Fourier transform, intended as a teaching tool and comparison with the Fast Fourier Transform (FFT).
- sftpack_test
- sgefa_openmp, a C++ code which compares a standard linear algebra solver against a revised version which can be run in parallel with OpenMP.
- sgefa_openmp_test
- shallow_water_1d, a C++ code which simulates the evolution of a 1D fluid governed by the time-dependent shallow water equations.
- shallow_water_1d_test
- shepard_interp_1d, a C++ code which defines and evaluates Shepard interpolants to 1D data, based on inverse distance weighting.
- shepard_interp_1d_test
- shepard_interp_2d, a C++ code which defines and evaluates Shepard interpolants to 2D data, based on inverse distance weighting.
- shepard_interp_2d_test
- shepard_interp_nd, a C++ code which defines and evaluates Shepard interpolants to M-dimensional data, based on inverse distance weighting.
- shepard_interp_nd_test
- simple_ga, a C++ code which is a simple example of a genetic algorithm, by Dennis Cormier and Sita Raghavan.
- simple_ga_test
- simplex_coordinates, a C++ code which computes the Cartesian coordinates of the vertices of a regular simplex in M dimensions.
- simplex_coordinates_test
- simplex_gm_rule, a C++ code which defines Grundmann-Moeller quadrature rules over the interior of a triangle in 2D, a tetrahedron in 3D, or over the interior of the simplex in M dimensions.
- simplex_gm_rule_test
- simplex_grid, a C++ code which generates a regular grid of points over the interior of an arbitrary simplex in M dimensions.
- simplex_grid_test
- simplex_integrals, a C++ code which returns the exact value of the integral of any monomial over the interior of the unit simplex in M dimensions.
- simplex_integrals_test
- simplex_monte_carlo, a C++ code which uses the Monte Carlo method to estimate an integral over the interior of the unit simplex in M dimensions.
- simplex_monte_carlo_test
- sine_transform, a C++ code which demonstrates some simple properties of the discrete sine transform.
- sine_transform_test
- snakes_and_ladders, a C++ code which simulates a one-player version of the game of Snakes and Ladders.
- snakes_and_ladders_test
- sobol, a C++ code which computes elements of the Sobol M-dimensional Quasi Mont Carlo (QMC) sequence, extended to a maximal spatial dimension of 1111, by Bennett Fox.
- sobol_test
- solve, a C++ code which implements a linear solver which makes it easy to create doubly-dimensioned arrays and solve associated linear systems.
- solve_test
- sor, a C++ code which implements a simple version of the Successive Over-Relaxation (SOR) method for the iterative solution of a linear system of equations.
- sor_test
- sort_rc, a C++ code which sorts a list of any kind of objects, using reverse communication (RC).
- sort_rc_test
- sparse_count, a C++ code which analyzes sparse grids in which a single family of 1D quadrature rules is used for all spatial dimensions, with a variety of growth rules.
- sparse_count_test
- sparse_display, a C++ code which reads information defining a matrix of numbers and displays the sparsity pattern or location of the nonzero elements using gnuplot(). This operation is already available in the built-in MATLAB spy() command.
- sparse_display_test
- sparse_grid_cc, a C++ code which creates sparse grids based on Clenshaw-Curtis rules.
- sparse_grid_cc_test
- sparse_grid_gl, a C++ code which creates sparse grids based on Gauss-Legendre rules.
- sparse_grid_gl_test
- sparse_grid_hermite, a C++ code which creates sparse grids based on Gauss-Hermite rules.
- sparse_grid_hermite_test
- sparse_grid_hw, a C++ code 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_hw_test
- sparse_grid_laguerre, a C++ code which creates sparse grids based on Gauss-Laguerre rules.
- sparse_grid_laguerre_test
- sparse_interp_nd, a C++ code which defines a sparse interpolant to a function f(x) of a M-dimensional argument.
- sparse_interp_nd_test
- sphere_exactness, a C++ code which tests the polynomial exactness of a quadrature rule over the surface of the unit sphere in 3D;
- sphere_exactness_test
- sphere_fibonacci_grid, a C++ code which uses a Fibonacci spiral to create a grid of points on the surface of the unit sphere in 3D.
- sphere_fibonacci_grid_test
- sphere_grid, a C++ code 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 surface of the unit sphere in 3D.
- sphere_grid_test
- sphere_integrals, a C++ code which returns the exact value of the integral of any monomial over the surface of the unit sphere in 3D.
- sphere_integrals_test
- sphere_lebedev_rule, a C++ code which computes Lebedev quadrature rules for the surface of the unit sphere in 3D;
- sphere_lebedev_rule_test
- sphere_llq_grid, a C++ code which uses longitudes and latitudes to create grids of points, lines, and quadrilaterals on the surface of the unit sphere in 3D.
- sphere_llq_grid_test
- sphere_llt_grid, a C++ code which uses longitudes and latitudes to create grids of points, lines, and triangles on the surface of the unit sphere in 3D.
- sphere_llt_grid_test
- sphere_monte_carlo, a C++ code which applies a Monte Carlo method to estimate the integral of a function over the surface of the unit sphere in 3D;
- sphere_monte_carlo_test
- sphere_quad, a C++ code which uses triangulation to approximate an integral over the surface of the unit sphere in 3D;
- sphere_quad_test
- sphere_stereograph, a C++ code which computes the stereographic mapping between points on the plane Z=1 and the surface of the unit sphere in 3D.
- sphere_stereograph_test
- spiral_data, a C++ code which computes a velocity vector field that satisfies the continuity equation, writing the data to a file that can be plotted by gnuplot().
- spiral_data_test
- spline, a C++ code which interpolates or approximates data by piecewise polynomial functions;
- spline_test
- square_arbq_rule, a C++ code which returns quadrature rules, with exactness up to total degree 20, over the interior of the symmetric square in 2D, by Hong Xiao and Zydrunas Gimbutas.
- square_arbq_rule_test
- square_exactness, a C++ code which investigates the polynomial exactness of quadrature rules for f(x,y) over the interior of a square in 2D.
- square_exactness_test
- square_felippa_rule, a C++ code which returns Felippa's quadratures rules for approximating integrals over the interior of a square in 2D.
- square_felippa_rule_test
- square_grid, a C++ code which computes a grid of points over the interior of a square in 2D.
- square_grid_test
- square_integrals, a C++ code which returns the exact value of the integral of any monomial over the interior of the unit square in 2D.
- square_integrals_test
- square_minimal_rule, a C++ code which returns "almost minimal" quadrature rules, with exactness up to total degree 55, over the interior of the symmetric square in 2D, by Mattia Festa and Alvise Sommariva.
- square_minimal_rule_test
- square_monte_carlo, a C++ code which applies a Monte Carlo method to estimate the integral of a function over the interior of the unit square in 2D.
- square_monte_carlo_test
- square_symq_rule, a C++ code which returns symmetric quadrature rules, with exactness up to total degree 20, over the interior of the symmetric square in 2D, by Hong Xiao and Zydrunas Gimbutas.
- square_symq_rule_test
- st_io, a C++ code which reads and writes sparse linear systems stored in the Sparse Triplet (ST) format.
- st_io_test
- st_to_ccs, a C++ code which converts sparse matrix data from Sparse Triplet (ST) format to Compressed Column Storage (CCS) format;
- st_to_ccs_test
- star_discrepancy, a C++ code which computes bounds on the star discrepancy of a set of N points in M dimensions, stored as a table file, by Eric Thiemard.
- star_discrepancy_test
- stiff_ode, a C++ code which considers an ordinary differential equation (ODE) which is an example of a stiff ODE.
- stiff_ode_test
- stla_io, a C++ code which reads and writes ASCII stereolithography (STL) 3d graphics files;
- stla_io_test
- stochastic_diffusion, functions which implement several versions of a stochastic diffusivity coefficient.
- stochastic_diffusion_test
- stochastic_heat2d, a C++ code which implements a finite difference method (FDM) for the steady (time independent) 2D heat equation, with a stochastic heat diffusivity coefficient.
- stochastic_heat2d_test
- stochastic_rk, a C++ code which applies a Runge-Kutta (RK) scheme to a stochastic ordinary differential equation (SODE).
- stochastic_rk_test
- stokes_2d_exact, a C++ code which evaluates exact solutions to the incompressible steady Stokes equations over the unit square in 2D.
- stokes_2d_exact_test
- string_pde, a C++ code which sets up and solves the partial differential equations (PDE) describing a vibrating string, creating files that can be displayed by gnuplot().
- string_pde_test
- stroud, a C++ code 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.
- stroud_test
- subpak, a C++ code which contains many utility routines;
- subpak_test
- subset, a C++ code which enumerates, generates, randomizes, ranks and unranks combinatorial objects including combinations, compositions, Gray codes, index sets, partitions, permutations, polynomials, subsets, and Young tables. Backtracking routines are included to solve some combinatorial problems.
- subset_test
- subset_sum, a C++ code which seeks solutions of the subset sum problem.
- subset_sum_test
- subset_sum_serial, a C++ code which seeks solutions of the subset sum problem, intended as a starting point for a parallel computing approach.
- subset_sum_serial_test
- svd_basis, a C++ code which applies the singular value decomposition (SVD) to a set of vectors, to extract the dominant modes;
- svd_basis_test
- svd_snowfall, a C++ code which reads a file containing historical snowfall data and analyzes the data with the Singular Value Decomposition (SVD), displaying the results using gnuplot().
- svd_snowfall_test
- svd_test
- svd_truncated, a C++ code 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.
- svd_truncated_test
- table_io, a C++ code which reads and writes a simple dataset stored as a table file;
- table_io_test
- task_division, a C++ code 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.
- task_division_test
- test_approx, a C++ code which implements test problems for approximation, provided as a set of (x,y) data.
- test_approx_test
- test_eigen, a C++ code which implements test matrices for eigenvalue analysis.
- test_eigen_test
- test_int, a C++ code which implements test problems for approximate integration in one dimension.
- test_int_test
- test_int_2d, a C++ code which implements test problems for approximate integration in two dimensions.
- test_int_2d_test
- test_interp, a C++ code which defines test problems for interpolation, provided as a set of (x,y) data.
- test_interp_test
- test_interp_1d, a C++ code which defines test problems for interpolation of data y(x), which depends on a 1D argument.
- test_interp_1d_test
- test_interp_2d, a C++ code which defines test problems for interpolation of data z(x,y), depending on a 2D argument.
- test_interp_2d_test
- test_interp_nd, a C++ code which defines test problems for interpolation of data z(x), depending on an M-dimensional argument.
- test_interp_nd_test
- test_lls, a C++ code which implements linear least squares (LLS) test problems of the form A*x=b.
- test_lls_test
- test_mat, a C++ code which defines test matrices for which the condition number, determinant, eigenvalues, eigenvectors, inverse, null vectors, P*L*U factorization or linear system solution are known. Examples include the Fibonacci, Hilbert, Redheffer, Vandermonde, Wathen and Wilkinson matrices.
- test_mat_test
- test_matrix_exponential, a C++ code which defines a set of test cases for computing the matrix exponential.
- test_matrix_exponential_test
- test_min, a C++ code which implements test problems for minimization of a scalar function of a scalar variable.
- test_min_test
- test_nint, a C++ code which implements test problems for approximate integration (quadrature) in M dimensions.
- test_nint_test
- test_opt_con, a C++ code 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_opt_con_test
- test_optimization, a C++ code which implements test problems for optimization of a scalar function of several variables, as described by Molga and Smutnicki.
- test_optimization_test
- test_partial_digest, a C++ code which generates cases of the partial digest problem.
- test_partial_digest_test
- test_values, a C++ code which returns a small set of values for a variety of mathematical functions, suitable for a rough test of accuracy for algorithms;
- test_values_test
- test_zero, a C++ code which implements test problems for the solution of a single nonlinear equation in one variable.
- test_zero_test
- tester, a BASH script which runs the test programs.
- testpack, a C++ code which demonstrates the use of a set of test integrand functions chosen by Alan Genz.
- testpack_test
- tet_mesh, a C++ code which works with tetrahedral meshes in 3D;
- tet_mesh_test
- tet_mesh_boundary, a C++ code 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_boundary_test
- tet_mesh_l2q, a C++ code which reads information about a 4-node (linear) tet mesh and creates data defining a corresponding 10-node (quadratic) tet mesh;
- tet_mesh_l2q_test
- tet_mesh_q2l, a C++ code which reads information about a 10-node (quadratic) tet mesh and creates data defining a corresponding 4-node (linear) tet mesh;
- tet_mesh_q2l_test
- tet_mesh_quad, a C++ code which estimates the integral of a function over a region defined by a tetrahedral mesh.
- tet_mesh_quad_test
- tet_mesh_quality, a C++ code which reads a node coordinate file and a tetrahedron file and returns quality measures for the tetrahedronal mesh;
- tet_mesh_quality_test
- tet_mesh_rcm, a C++ code 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_rcm_test,
- tet_mesh_refine, a C++ code which refines a tet mesh;
- tet_mesh_refine_test
- tet_mesh_tet_neighbors, a C++ code which reads information about a tet mesh and writes out a list of the tetrahedrons that are adjacent to a given tetrahedron;
- tet_mesh_tet_neighbors_test
- tet_mesh_to_gmsh, a C++ code which converts a tet mesh to gmsh() format;
- tet_mesh_to_gmsh_test
- tet_mesh_to_xml, a C++ code which converts a tet mesh to XML format;
- tet_mesh_to_xml_test
- tet_mesh_volumes, a C++ code which computes the volume of each tetrahedron in a tet mesh;
- tet_mesh_volumes_test
- tetrahedron_arbq_rule, a C++ code which returns quadrature rules, with exactness up to total degree 15, over the interior of a tetrahedron in 3D, by Hong Xiao and Zydrunas Gimbutas.
- tetrahedron_arbq_rule_test
- tetrahedron_exactness, a C++ code which tests the polynomial exactness of a quadrature rule over the interior of a tetrahedron in 3D.
- tetrahedron_exactness_test
- tetrahedron_felippa_rule, a C++ code which returns Felippa's quadratures rules for approximating integrals over the interior of a tetrahedron in 3D.
- tetrahedron_felippa_rule_test
- tetrahedron_grid, a C++ code which computes a grid of points over the interior of a tetrahedron in 3D.
- tetrahedron_grid_test
- tetrahedron_integrals, a C++ code which returns the exact value of the integral of any monomial over the interior of the unit tetrahedron in 3D.
- tetrahedron_integrals_test
- tetrahedron_keast_rule, a C++ code which defines ten quadrature rules, of degree of exactness 0 through 8, over the interior of the tetrahedron in 3D.
- tetrahedron_keast_rule_test
- tetrahedron_monte_carlo, a C++ code which uses the Monte Carlo method to estimate an integral over the interior of a general tetrahedron in 3D.
- tetrahedron_monte_carlo_test
- tetrahedron_ncc_rule, a C++ code which defines Newton-Cotes closed (NCC) quadrature rules over the interior of a tetrahedron in 3D.
- tetrahedron_ncc_rule_test
- tetrahedron_nco_rule, a C++ code which defines Newton-Cotes open (NCO) quadrature rules over the interior of a tetrahedron in 3D.
- tetrahedron_nco_rule_test
- tetrahedron_properties, a C++ code which computes properties including the centroid, circumsphere, dihedral angles, edge lengths, face angles, face areas, insphere, quality, solid angles, and volume of a tetrahedron in 3D.
- tetrahedron_properties_test
- tetrahedron01_monte_carlo, a C++ code which uses the Monte Carlo method to estimate an integral over the interior of the unit tetrahedron in 3D.
- tetrahedron01_monte_carlo_test
- three_body_ode, a C++ code 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.
- three_body_ode_test
- tiler_3d, a C++ code which blends a pattern repeatedly into a region in 3D;
- tiler_3d_test
- timer_test, a C++ code which measures the CPU time used by a portion of a user's program;
- timestamp, a C++ code which prints the current YMDHMS date as a timestamp, and other timekeeping operations;
- timestamp_test
- toeplitz_cholesky, a C++ code which computes the Cholesky factorization of a symmetric positive definite (SPD) Toeplitz matrix.
- toeplitz_cholesky_test
- toms097, a C++ code which computes the distance between all pairs of nodes in a directed graph with weighted edges, using Floyd's algorithm. This is a version of ACM TOMS algorithm 97.
- toms097_test
- toms112, a C++ code which determines whether a point is contained in a polygon, by Moshe Shimrat. This is a version of ACM TOMS algorithm 112.
- toms112_test
- toms178, a C++ code which seeks the minimizer of a scalar function of several variables using the Hooke-Jeeves method, by Arthur Kaupe. This is a version of ACM TOMS algorithm 178.
- toms178_test
- toms179, a C++ code which calculates the incomplete Beta ratio, by Oliver Ludwig. This is a version of ACM TOMS algorithm 179.
- toms179_test
- toms243, a C++ code which evaluates the logarithm of a complex value, by David Collens. This is a version of ACM TOMS algorithm 243.
- toms243_test
- toms291, a C++ code which approximates the logarithm of the Gamma function; this is a version of ACM TOMS algorithm 291.
- toms291_test
- toms443, a C++ code which evaluates Lambert's W function. This is a version of ACM TOMS algorithm 443.
- toms443_test
- toms446, a C++ code which manipulates Chebyshev series for interpolation and approximation. This is a version of ACM TOMS algorithm 446, by Roger Broucke.
- toms446_test
- toms462, a C++ code 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 version of ACM TOMS algorithm 462.
- toms462_test
- toms515, a C++ code which selects subsets of size K from a set of size N, This is a version of ACM TOMS Algorithm 515, by Bill Buckles, Matthew Lybanon.
- toms515_test
- toms577, a C++ code which evaluates Carlson's elliptic integral functions RC, RD, RF and RJ. This is a version of ACM TOMS algorithm 577.
- toms577_test
- toms655, a C++ code which computes the weights for interpolatory quadrature rule; this library is commonly called IQPACK, by Sylvan Elhay and Jaroslav Kautsky. This is a version of ACM TOMS algorithm 655.
- toms655_test
- toms743, a C++ code which evaluates Lambert's W function. This is a version of ACM TOMS algorithm 743, by Barry, Barry and Culligan-Hensley.
- toms743_test
- toms886, a C++ code which defines the Padua points for interpolation in a 2D region, including the rectangle, triangle, and ellipse, by Marco Caliari, Stefano de Marchi, Marco Vianello. This is a version of ACM TOMS algorithm 886.
- toms886_test
- toms917, a C++ code which evaluates the Wright Omega function, by Piers Lawrence, Robert Corless, David Jeffrey. This is a version of ACM TOMS algorithm 917.
- toms917_test
- treepack, a C++ code which carries out computations on trees, a simple kind of graph that is minimally connected.
- treepack_test
- tri_surface_io, a C++ code which reads and writes the 3D graphics information in a tri_surface file;
- tri_surface_io_test
- triangle_analyze, a C++ code which reads a triangle defined in a file, and uses triangle_properties() to compute angles, area, centroid, circumcircle, edge lengths, incircle, orientation, orthocenter, and quality.
- triangle_analyze_test
- triangle_dunavant_rule, a C++ code which defines Dunavant quadrature rules over the interior of a triangle in 2D.
- triangle_dunavant_rule_test
- triangle_exactness, a C++ code which tests the polynomial exactness of a quadrature rule over the interior of a triangle in 2D.
- triangle_exactness_test
- triangle_fekete_rule, a C++ code which sets up one of seven Fekete rules for interpolation or quadrature over the interior of a triangle in 2D.
- triangle_fekete_rule_test
- triangle_felippa_rule, a C++ code which returns Felippa's quadratures rules for approximating integrals over the interior of a triangle in 2D.
- triangle_felippa_rule_test
- triangle_grid, a C++ code which computes a grid of points over the interior of a triangle in 2D.
- triangle_grid_test
- triangle_histogram, a C++ code which computes histograms of data over the interior of the unit triangle in 2D.
- triangle_histogram_test
- triangle_integrals, a C++ code which returns the exact value of the integral of any polynomial over the interior of an arbitrary triangle in 2D.
- triangle_integrals_test
- triangle_interpolate, a C++ code which shows how vertex data can be interpolated at any point in the interior of a triangle.
- triangle_interpolate_test
- triangle_io, a C++ code which reads or writes files created by Jonathan Shewchuk's triangle() program.
- triangle_io_test
- triangle_lyness_rule, a C++ code which returns Lyness-Jespersen quadrature rules over the interior of a triangle in 2D;
- triangle_lyness_rule_test
- triangle_monte_carlo, a C++ code which uses the Monte Carlo method to estimate an integral over the interior of a general triangle in 2D.
- triangle_monte_carlo_test
- triangle_ncc_rule, a C++ code which defines Newton-Cotes closed (NCC) quadrature rules over the interior of a triangle in 2D.
- triangle_ncc_rule_test
- triangle_nco_rule, a C++ code which defines Newton-Cotes open (NCO) quadrature rules over the interior of a triangle in 2D.
- triangle_nco_rule_test
- triangle_properties, a C++ code which computes properties, including angles, area, centroid, circumcircle, edge lengths, incircle, orientation, orthocenter, and quality, of a triangle in 2D.
- triangle_properties_test
- triangle_svg, a C++ code which uses Scalable Vector Graphics (SVG) to plot a triangle and any number of points, to illustrate quadrature rules and sampling techniques.
- triangle_svg_test
- triangle_symq_rule, a C++ code which returns efficient symmetric quadrature rules, with exactness up to total degree 50, over the interior of a triangle in 2D, by Hong Xiao and Zydrunas Gimbutas.
- triangle_symq_rule_test
- triangle_to_fem, a C++ code 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_fem_test
- triangle_to_xml, a C++ code which reads the NODE and ELE files created by triangle() to describe a triangular mesh in 2D, and writes out a corresponding XML mesh file for use by DOLFIN or FENICS.
- triangle_to_xml_test
- triangle_twb_rule, a C++ code which generates the points and weights of quadrature rules over the interior of a triangle in 2D, determined by Taylor, Wingate, and Bos.
- triangle_twb_rule_test
- triangle_wandzura_rule, a C++ code which returns quadrature rules of exactness 5, 10, 15, 20, 25 and 30 over the interior of the triangle in 2D.
- triangle_wandzura_rule_test
- triangle01_integrals, a C++ code which returns the integral of any monomial over the interior of the unit triangle in 2D.
- triangle01_integrals_test
- triangle01_monte_carlo, a C++ code which uses the Monte Carlo method to estimate an integral over the interior of the unit triangle in 2D.
- triangle01_monte_carlo_test
- triangulation, a C++ code which computes the triangulation of a set of points in the plane, and performs various operations using a triangulation, including searching a Delaunay triangulation to find which triangle contains a given point.
- triangulation_test
- triangulation_boundary_nodes, a C++ code which reads data defining a triangulation, determines which nodes lie on the boundary, and writes their coordinates to a file;
- triangulation_boundary_nodes_test
- triangulation_corner, a C++ code which modifies triangulations in which one or more triangles have more than one boundary edge.
- triangulation_corner_test
- triangulation_delaunay_discrepancy, a C++ code which measures the amount by which a triangulation fails the local Delaunay test;
- triangulation_delaunay_discrepancy_test
- triangulation_histogram, a C++ code which computes histograms of data over a triangulation.
- triangulation_histogram_test
- triangulation_l2q, a C++ code which reads information about a 3-node (linear) triangulation and creates data defining a corresponding 6-node (quadratic) triangulation;
- triangulation_l2q_test
- triangulation_mask, a C++ code which reads a triangulation and calls a user-supplied routine to consider each triangle for deletion;
- triangulation_mask_test
- triangulation_node_to_element, a C++ code which reads files describing a set of nodes, their triangulation, and the value of one or more quantities at each node, and outputs a file that averages the quantities for each element. This operation in effect creates an "order1" finite element model of the data.
- triangulation_node_to_element_test
- triangulation_orient, a C++ code which ensures that the triangles in an order 3 or order 6 triangulation have positive orientation;
- triangulation_orient_test
- triangulation_plot, a C++ code which plots the nodes and elements of a triangulation as a PostScript (PS) file;
- triangulation_plot_test
- triangulation_q2l, a C++ code which reads information about a 6-node triangulation and creates data defining a corresponding 3-node triangulation;
- triangulation_q2l_test
- triangulation_quad, a C++ code 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_quad_test
- triangulation_quality, a C++ code which reads information about a triangulation and computes various quality measures;
- triangulation_quality_test
- triangulation_refine, a C++ code which refines a triangulation;
- triangulation_refine_test
- triangulation_svg, a C++ code which creates an SVG image of a triangulation, which can be displayed by a web browser.
- triangulation_svg_test
- triangulation_t3_to_t4, a C++ code which reads information about a 3-node triangulation and creates data defining a corresponding 4-node triangulation (vertices + centroid);
- triangulation_t3_to_t4_test
- triangulation_triangle_neighbors, a C++ code which reads data defining a triangulation, finds the three neighbor triangles for each triangle, and writes them to a file;
- triangulation_triangle_neighbors_test
- truncated_normal, a C++ code which works with the truncated normal distribution over [A,B], or [A,+oo) or (-oo,B], returning the probability density function (PDF), the cumulative density function (CDF), the inverse CDF, the mean, the variance, and sample values.
- truncated_normal_test
- truncated_normal_rule, a C++ code which computes a quadrature rule for a normal probability density function (PDF), sometimes called a Gaussian distribution, that has been truncated to [A,+oo), (-oo,B] or [A,B].
- truncated_normal_rule_test
- tsp_brute, a C++ code which reads a file of city-to-city distances and solves a (small!) traveling salesperson problem, using brute force.
- tsp_brute_test
- ubvec, a C++ code which demonstrates how nonnegative integers can be stored as unsigned binary vectors, and arithmetic can be performed on them.
- ubvec_test
- umfpack_test
- unicycle, a C++ code which considers permutations containing a single cycle, sometimes called cyclic permutations.
- unicycle_test
- uniform, a C++ code which implements uniform random number generators (RNG) for a variety of arithmetic types.
- uniform_test
- upc, a C++ code which computes the check digit of a uniform product code (UPC), or reports whether a full 12-digit UPC is actually valid.
- upc_test
- valgrind_test
- van_der_corput, a C++ code which evaluates the van der Corput 1-dimensional Quasi Mont Carlo (QMC) sequence, using a simple interface.
- van_der_corput_test
- vandermonde, a C++ code which implements the Bjork-Pereyra algorithm for accurate solution of linear systems involving the Vandermonde matrix.
- vandermonde_test
- vandermonde_approx_1d, a C++ code 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_1d_test
- vandermonde_approx_2d, a C++ code 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_approx_2d_test
- vandermonde_interp_1d, a C++ code 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, creating graphics files for processing by gnuplot().
- vandermonde_interp_1d_test
- vandermonde_interp_2d, a C++ code 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.
- vandermonde_interp_2d_test
- walker_sample, a C++ code which efficiently samples a discrete probability vector using Walker sampling.
- walker_sample_test
- walsh, a C++ code which implements versions of the Walsh and Haar transforms.
- walsh_test
- wathen, a C++ code which compares storage schemes (full, banded, sparse triplet ) and solution strategies: Linpack full, Linpack banded, conjugate gradient (CG), for linear systems involving the Wathen matrix, which can arise when solving a problem using the finite element method (FEM).
- wathen_test
- wave_mpi, a C++ code which uses finite differences and the Message Passing Interface (MPI) to estimate a solution to the wave equation.
- wave_mpi_test
- wavelet, a C++ code which does some simple wavelet calculations;
- wavelet_test
- wedge_exactness, a C++ code which investigates the polynomial exactness of a quadrature rule over the interior of the unit wedge in 3D.
- wedge_exactness_test
- wedge_felippa_rule, a C++ code which returns quadrature rules for approximating integrals over the interior of the unit wedge in 3D.
- wedge_felippa_rule_test
- wedge_grid, a C++ code which computes a grid of points over the interior of the unit wedge in 3D.
- wedge_grid_test
- wedge_integrals, a C++ code which returns the exact value of the integral of any monomial over the interior of the unit wedge in 3D.
- wedge_integrals_test
- wedge_monte_carlo, a C++ code which uses the Monte Carlo method to estimate an integral over the interior of the unit wedge in 3D.
- wedge_monte_carlo_test
- weekday, a C++ code which determines the day of the week corresponding to a given date, such as 14 October 1066, Julian calendar, ... which was a Saturday.
- weekday_test
- wishart, a C++ code which produces sample matrices from the Wishart or Bartlett distributions, useful for sampling random covariance matrices.
- wishart_test
- wtime, a C++ code which returns a reading of the wall clock time in seconds.
- wtime_test
- zero, a C++ code which seeks a solution of a scalar nonlinear equation f(x)=0, by Richard Brent.
- zero_test
- zero_rc, a C++ code which seeks a solution of a scalar nonlinear equation f(x)=0, using reverse communication (RC), by Richard Brent.
- zero_rc_test
- ziggurat, a C++ code which samples the uniform, normal or exponential distributions, using the ziggurat method, by Marsaglia and Tsang.
- ziggurat_test
- ziggurat_openmp, a C++ code which demonstrates how ziggurat() can be used as a random number generator (RNG) in an OpenMP parallel program.
- ziggurat_openmp_test