HB_READ
Simplified Harwell-Boeing Sparse Matrix File Reader

HB_READ is a FORTRAN90 library which reads files in the Harwell Boeing (HB) sparse matrix format. This is a simplified interface intended to handle only the most common format, complex unsymmetric assembled (CUA) or real unsymmetric assembled (RUA).

Licensing:

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

Languages:

HB_READ is available in a FORTRAN90 version.

Related Data and Programs:

DLAP_IO, a FORTRAN90 library which reads and writes matrices in a variety of formats, and includes a routine for reading matrices in HB format.

HB, a dataset directory which contains a collection of sparse matrices in the HB format.

HB_IO, a FORTRAN90 library which reads and writes sparse linear systems stored in the Harwell Boeing (HB) format for sparse matrices. This library is intended to handle the full range of HB formats.

HB_TO_MSM, a MATLAB program which reads an HB file and converts the matrix to MATLAB's sparse matrix format.

HB_TO_ST, a FORTRAN90 program which converts the sparse matrix information stored in a Harwell-Boeing file into a sparse triplet file.

HBSMC, a dataset directory which contains the Harwell Boeing Sparse Matrix Collection;

LINPLUS, a FORTRAN90 library which contains some routines to manipulate matrices stored in the HB sparse matrix format.

MM_IO, a FORTRAN90 library which reads and writes matrices in the Matrix Market format.

MSM_TO_HB, a MATLAB program which takes a matrix in MATLAB's sparse matrix format and writes it out in the HB format.

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

SUPERLU, FORTRAN90 programs which illustrate how to use 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. Iain Duff, Roger Grimes, John Lewis,
   User's Guide for the Harwell-Boeing Sparse Matrix Collection,
   Technical Report TR/PA/92/86, CERFACS,
   October 1992.
2. Iain Duff, Roger Grimes, John Lewis,
   The Rutherford-Boeing Sparse Matrix Collection,
   Technical Report RAL-TR-97-031,
   Rutherford Appleton Laboratory, 1997.
3. Iain Duff, Roger Grimes, John Lewis,
   Sparse Matrix Test Problems,
   ACM Transactions on Mathematical Software,
   Volume 15, pages 1-14, March 1989.
4. http://math.nist.gov/MatrixMarket/data/Harwell-Boeing/ the Harwell-Boeing area of the MatrixMarket web site.

Source Code:

• hb_read.f90, the source code.

Examples and Tests:

• hb_read_test.f90, a sample calling program.
• hb_read_test.txt, the output file.

Various Harwell-Boeing matrices used by the tests include:

• cavity_big_rua.txt, a real unsymmetric assembled 5030x5030 matrix stored in a Harwell Boeing matrix file, with 132463 nonzeros.
• cavity_small_rua.txt, a real unsymmetric assembled 92x92 matrix stored in a Harwell Boeing matrix file, with 1363 nonzeros.
• cg05_cua.txt, a complex unsymmetric assembled 25x25 matrix stored in a Harwell Boeing matrix file, with 105 nonzeros.
• cg20_cua.txt, a complex unsymmetric assembled 400x400 matrix stored in a Harwell Boeing matrix file, with 1920 nonzeros.
• cmat_cua.txt, a complex unsymmetric assembled 1280x1280 matrix stored in a Harwell Boeing matrix file, with 22778 nonzeros.
• g05_rua.txt, a real unsymmetric assembled 25x25 matrix stored in a Harwell Boeing matrix file, with 105 nonzeros.
• g10_rua.txt, a real unsymmetric assembled 100x100 matrix stored in a Harwell Boeing matrix file, with 460 nonzeros.
• g20_rua.txt, a real unsymmetric assembled 400x400 matrix stored in a Harwell Boeing matrix file, with 1920 nonzeros.

List of Routines:

• C4_HB_DATA_READ reads the data for a C4 Harwell-Boeing file.
• C4_HB_HEADER_PRINT prints the header for a C4 Harwell-Boeing file.
• C4_HB_HEADER_READ reads the header for a C4 Harwell-Boeing file.
• C4_HB_QUICK_PRINT prints a sparse matrix in C4 Harwell-Boeing format.
• C4_HB_QUICK_READ reads a sparse matrix in C4 Harwell-Boeing format.
• C4_HB_STRUCTURE_PRINT prints the structure of a C4 Harwell-Boeing matrix.
• C4_HB_VALUES_PRINT prints the values of a C4 Harwell-Boeing matrix.
• C4MAT_PRINT_SOME prints some of a C4MAT.
• C8_HB_READ reads a sparse matrix in C8 Harwell-Boeing format.
• C8_HB_DATA_READ reads the data for a C8 Harwell-Boeing file.
• C8_HB_HEADER_PRINT prints the header for a C8 Harwell-Boeing file.
• C8_HB_HEADER_READ reads the header for a C8 Harwell-Boeing file.
• C8_HB_QUICK_PRINT prints a sparse matrix in C8 Harwell-Boeing format.
• C8_HB_QUICK_READ reads a sparse matrix in C8 Harwell-Boeing format.
• C8_HB_STRUCTURE_PRINT prints the structure of a C8 Harwell-Boeing matrix.
• C8_HB_VALUES_PRINT prints the values of a C8 Harwell-Boeing matrix.
• C8MAT_PRINT prints a C8MAT.
• C8MAT_PRINT_SOME prints some of a C8MAT.
• GET_UNIT returns a free FORTRAN unit number.
• R4_HB_DATA_READ reads the data for an R4 Harwell-Boeing file.
• R4_HB_HEADER_PRINT prints the header for an R4 Harwell-Boeing file.
• R4_HB_HEADER_READ reads the header for an R4 Harwell-Boeing file.
• R4_HB_QUICK_PRINT prints a sparse matrix in R4 Harwell-Boeing format.
• R4_HB_QUICK_READ reads a sparse matrix in R4 Harwell-Boeing format.
• R4_HB_STRUCTURE_PRINT prints the structure of an R4 Harwell-Boeing matrix.
• R4_HB_VALUES_PRINT prints the values of an R4 Harwell-Boeing matrix.
• R4MAT_PRINT prints an R4MAT.
• R4MAT_PRINT_SOME prints some of an R4MAT.
• R8_HB_DATA_READ reads the data for an R8 Harwell-Boeing file.
• R8_HB_HEADER_PRINT prints the header for an R8 Harwell-Boeing file.
• R8_HB_HEADER_READ reads the header for an R8 Harwell-Boeing file.
• R8_HB_QUICK_PRINT prints a sparse matrix in R8 Harwell-Boeing format.
• R8_HB_QUICK_READ reads a sparse matrix in R8 Harwell-Boeing format.
• R8_HB_STRUCTURE_PRINT prints the structure of an R8 Harwell-Boeing matrix.
• R8_HB_VALUES_PRINT prints the values of an R8 Harwell-Boeing matrix.
• R8MAT_PRINT prints an R8MAT.
• R8MAT_PRINT_SOME prints some of an R8MAT.
• TIMESTAMP prints the current YMDHMS date as a time stamp.

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