cg_rc


cg_rc, a C code which implements the conjugate gradient method for solving a positive definite sparse linear system A*x=b, using reverse communication (RC).

Licensing:

The computer code and data files described and made available on this web page are distributed under the MIT license

Languages:

cg_rc is available in a C version and a C++ version and a FORTRAN90 version and a MATLAB version and a Python version.

Related Data and Programs:

backtrack_binary_rc, a C code which carries out a backtrack search for a set of binary decisions, using reverse communication (RC).

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).

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 positive definite (only real, positive eigenvalues) and symmetric.

cg_rc_test

CSPARSE, a C code which carries out the direct solution of sparse linear systems, by Timothy Davis.

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.

MGMRES, a C code which applies the restarted GMRES algorithm to solve a sparse linear system.

MULTIGRID_POISSON_1D, a C code which applies the multigrid method to a discretized version of the 1D Poisson equation.

ROOT_RC, a C code which seeks a solution of a scalar nonlinear equation f(x) = 0, or a system of nonlinear equations, using reverse communication (RC), by Gaston Gonnet.

ROOTS_RC, a C code which seeks a solution of a system of nonlinear equations f(x) = 0, using reverse communication (RC), by Gaston Gonnet.

SORT_RC, a C code which can sort a list of any kind of objects, using reverse communication (RC).

ZERO_RC, a C code which seeks solutions of a scalar nonlinear equation f(x) = 0, using reverse communication (RC).

Reference:

  1. Richard Barrett, Michael Berry, Tony Chan, James Demmel, June Donato, Jack Dongarra, Victor Eijkhout, Roidan Pozo, Charles Romine, Henk van der Vorst,
    Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods,
    SIAM, 1994,
    ISBN: 0898714710,
    LC: QA297.8.T45.
  2. Jonathan Shewchuk,
    An introduction to the conjugate gradient method without the agonizing pain,
    Edition 1.25, August 1994.

Source Code:


Last revised on 13 June 2019.