umfpack_test

umfpack_testK, FORTRAN90 programs which call umfpack to solve sparse linear systems.

The UMFPACK library, written by Timothy Davis, is part of the SuiteSparse package, which is available from http://www.cise.ufl.edu/research/sparse/ .

Calling UMFPACK requires linking additional libraries from the SuiteSparse collection. The executable shell scripts given here suggest how this might be done.

UMFPACK requires that the sparse matrix be stored in a common format known as "compressed column storage". UMFPACK uses the C indexing convention, in which the first element of a vector has index 0.

As an example, the matrix whose full representation is:

       11 12  0  0  0
       21  0 23  0 25
        0 32 33 34  0
        0  0 43  0  0
        0 52 53  0 55
      

         #  R   CC   X
        --  -   --  --
         0  0    0  11
         1  1       21

         2  0    2  12
         3  2       32
         4  4       52

         5  1    5  23
         6  2       33
         7  3       43
         8  4       53

         9  2    9  34

        10  1   10  25
        11  4       55

        12  -   12
      

Calling a C library from FORTRAN90 is often difficult, and the details can vary depending on the computer operating system and compilers used. Moreover, calling UMFPACK requires linking additional libraries from the SuiteSparse collection. The executable shell scripts given here suggest how this might be done.

Licensing:

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

Languages:

umfpack_test is available in a C version and a C++ version and a FORTRAN90 version and a MATLAB version.

Related Data and Programs:

CC, a data directory which contains examples of the compressed column (CC) sparse matrix format;

CG_RC, a FORTRAN90 library which implements the conjugate gradient method for solving a positive definite sparse linear system A*x=b, using reverse communication.

HB_IO, a FORTRAN90 library which reads and writes sparse linear systems stored in the Harwell Boeing (HB) format for sparse matrices.

MGMRES, a FORTRAN90 library which applies the restarted Generalized Minimum Residual (GMRES) algorithm to solve a sparse linear system, by Lili Ju.

ST_IO, a FORTRAN90 library which reads and writes sparse linear systems stored in the Sparse Triplet (ST) format.

SUPERLU, FORTRAN90 programs which illustrate the usage of the SUPERLU library, which applies a fast direct solution method to solve sparse linear systems, by James Demmel, John Gilbert, and Xiaoye Li.

Reference:

1. Timothy Davis,
   Algorithm 832: UMFPACK, an unsymmetric-pattern multifrontal method,
   ACM Transactions on Mathematical Software,
   Volume 30, Number 2, June 2004, pages 196-199.
2. Timothy Davis,
   Direct Methods for Sparse Linear Systems,
   SIAM, 2006,
   ISBN: 0898716136,
   LC: QA188.D386.

Source Code:

• umf4_f77wrapper.c, a C library of interfaces that can be used to allow a FORTRAN program to invoke UMFPACK functions. This library must be compiled and linked along with the user calling program.

Examples and Tests:

The SIMPLE example demonstrates the use of UMFPACK on a 5x5 "sparse" matrix.

• umfpack_simple.f90, the source code.
• umfpack_simple_output.txt, the output file.

The WATHEN example demonstrates the use of UMFPACK on a sparse matrix that is an instance of the Wathen finite element mass matrix. The matrix is created using the sparse triplet (ST) format, and must be converted to the compressed column (CC) format before calling UMFPACK. This demonstrates why the 5x5 SIMPLE example is not representative of the typical use of UMFPACK.

• umfpack_wathen.f90, the source code.
• umfpack_wathen_output.txt, the output file.

The WEST example demonstrates the use of UMFPACK on a sparse matrix in the compressed column (CC) format, stored in three files. The matrix is read from the files, and then UMFPACK is called to handle an associated linear system.

• umfpack_west.f90, the source code.
• umfpack_west_output.txt, the output file.
• west_acc.txt, a file containing the nonzero matrix values.
• west_ccc.txt, a file containing the compressed column indices.
• west_icc.txt, a file containing the row indices of the nonzero matrix values.

You can go up one level to the FORTRAN90 source codes.