Class matrix

The matrix class includes methods to create and manipulate dense and (in the future) sparse matrices in ways required for BEM and FEM analysis^4.1 with the following methods:

• matCreate is the basic constructor, it allocates memory for the matrix structure and its data.
• matDuplicate is an alternate constructor which creates a new matrix object identical to an existing one.
• matDestroy is the basic destructor, which frees all memory associated with a matrix object.

Please see the file matrix.h (installed in prefix/Julian) for error returns.

The only current implementation of the compute-intensive matrix methods consists of wrappers to the native system Basic Linear Algebra Subroutines, or BLAS. Alternative dgemm-based implementations of the dgetrf and dtrsm routines and also dgemm itself (for LU-decomposition, back substitution and matrix multiplication respectively) are provided, which make callbacks to external routines to indicate progress toward completion of their operations. These routines are currently optimized for single processor operation, but at some point may have their backend replaced by subroutines which store and operate on matrix components in parallel.