Gromacs  2016.4
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
List of all members | Public Attributes
gmx_sparsematrix Struct Reference

#include <gromacs/linearalgebra/sparsematrix.h>


Sparse matrix storage format.

This structure specifies a storage format for a sparse matrix. The memory requirements are only proportional to the number of nonzero elements, and it provides a reasonably fast way to perform matrix-vector multiplications.

The data format is very similar to a neighborlist. It is optimized for fast access, but it is difficult to add entries. If you are constructing a matrix you should either do it in exactly the order specified here, or use some other more flexible intermediate structure.

The index array is of size nrow+1. All non-zero matrix elements on row i are stored in positions index[i] through index[i+1]-1 in the arrays column and value. The column array contains the column index for each entry, in ascending order, and the corresponding position in the value array contains the floating point matrix element.

index[nrow] should be equal to the total number of elements stored.

Thus, to find the value of matrix element [5,4] you should loop over positions index[5] to index[6]-1 in column until you either find the value 4, or a higher value (meaning the element was zero).

It is fairly easy to construct the matrix on-the-fly if you can do it row-by-row.

IMPORTANT: If compressed_symmetric is set to TRUE, you should only store EITHER the upper OR lower triangle (and the diagonal), and the other half is assumed to be symmetric. Otherwise, if compressed_symmetric==FALSE, no symmetry is implied and all elements should be stored.

The symmetry compression saves us a factor 2 both in storage and matrix multiplication CPU-time, which can be very useful for huge eigenproblems.

If you are unsure, just set compressed_symmetric to FALSE and list all elements. If you enable it but still list all elements (both upper and lower triangle) you will be sorry...

Internally, the sparse data is stored as a separate list for each row, where the list element is a structure with a column and (floating-point) data value. This makes it possible, although not completely transparent, to update values in random access order. The drawback is that the structure will allocate nrow memory regions. The matrix data could be stored in a single contiguous array with indices for each row, but then we could only insert elements at the end without copying the entire matrix.

After you have

In other words: Not perfect, but it works.

Public Attributes

gmx_bool compressed_symmetric
 Store half elements and assume symmetry. More...
int nrow
 Number of rows in matrix.
int * ndata
 Number of entries on each row (list)
int * nalloc
 Allocated entry list length for each row.
gmx_sparsematrix_entry_t ** data
 data[i] is a list with entries on row i

Member Data Documentation

gmx_bool gmx_sparsematrix::compressed_symmetric

Store half elements and assume symmetry.

The documentation for this struct was generated from the following file: