maspack.matrix

Class SparseCRSMatrix

java.lang.Object

maspack.matrix.MatrixBase

maspack.matrix.SparseMatrixBase

maspack.matrix.SparseCRSMatrix

All Implemented Interfaces:

java.lang.Cloneable, LinearTransformNd, Matrix, SparseMatrix, Clonable

public class SparseCRSMatrix
extends SparseMatrixBase
implements Clonable

Implements a square sparse matrix composed of MatrixBlocks.

Nested Class Summary

Nested classes/interfaces inherited from interface maspack.matrix.Matrix

Matrix.Partition, Matrix.WriteFormat

Field Summary

Fields inherited from interface maspack.matrix.Matrix

INDEFINITE, POSITIVE_DEFINITE, SPD, SYMMETRIC

Constructor Summary

Constructors

Constructor and Description
SparseCRSMatrix() Creates an empty SparseCRSMatrix with 0 x 0 size.
SparseCRSMatrix(double[] vals, int[] cols, int[] rowOffs) Creates a sparse matrix from CRS data.
SparseCRSMatrix(int numRows, int numCols) Creates an empty sparse matrix with a specified size.
SparseCRSMatrix(Matrix M) Creates a sparse matrix from an existing general matrix.

Method Summary

Modifier and Type	Method and Description
void	add(int i, int j, double v)
void	add(SparseCRSMatrix M) Adds this matrix to M and places the result in this matrix.
void	add(SparseCRSMatrix M1, SparseCRSMatrix M2) Adds matrix M1 to M2 and places the result in this matrix.
void	addRow(DynamicIntArray colIdxs, DynamicDoubleArray values)
void	addRow(int[] colIdxs, double[] values, int numVals)
void	addScaledRow(VectorNd v, double s, int i) Adds s times the i-th row to the vector v.
void	checkConsistency() Check that the internal structures of this matrix are consistent.
SparseCRSMatrix	clone()
int	colSize() Number of columns in the matrix associated with this transformation.
void	combine(double a, SparseCRSMatrix M1, double b, SparseCRSMatrix M2) Forms the combination a M1 + b M2 and places the result in this matrix.
double	dotRow(int i, VectorNd v) Returns the dot product of the i-th row with a vector v.
double	dotRowSelf(int i) Returns the dot product of the i-th row with itself.
double	get(int i, int j) Gets a single element of this matrix.
void	get(MatrixNd M) Copies the contents of this matrix into a dense MatrixNd.
int	getCCSIndices(int[] rowIdxs, int[] colOffs, Matrix.Partition part, int numRows, int numCols) Gets the compressed column storage (CCS) indices for a principal sub-matrix of this matrix delimited by the first numRows rows and the first numCols columns.
int	getCCSValues(double[] vals, Matrix.Partition part, int numRows, int numCols) Gets the compressed column storage (CCS) values for a principal sub-matrix of this matrix delimited by the first numRows rows and the first numCols columns.
int	getCRSIndices(int[] colIdxs, int[] rowOffs, Matrix.Partition part, int numRows, int numCols) Gets the compressed row storage (CRS) indices for a principal sub-matrix of this matrix delimited by the first numRows rows and the first numCols columns.
int	getCRSValues(double[] vals, Matrix.Partition part, int numRows, int numCols) Gets the compressed row storage (CRS) values for a principal sub-matrix of this matrix delimited by the first numRows rows and the first numCols columns.
boolean	isFixedSize() Returns true if this matrix is of fixed size.
double	maxNorm() Returns the max-norm of this matrix.
void	mul(MatrixNd MR, MatrixNd M1) Computes the product
void	mul(SparseCRSMatrix M) Multiplies this matrix by M and places the result in this matrix.
void	mul(SparseCRSMatrix M1, SparseCRSMatrix M2) Multiplies matrix M1 by M2 and places the result in this matrix.
void	mulLeft(MatrixNd MR, MatrixNd M1) Computes the product
void	mulTranspose(SparseCRSMatrix M) Multiplies this matrix by the transpose of M and places the result in this matrix.
void	mulTransposeBoth(SparseCRSMatrix M1, SparseCRSMatrix M2) Multiplies the transpose of matrix M1 by the transpose of M2 and places the result in this matrix.
void	mulTransposeLeft(SparseCRSMatrix M1, SparseCRSMatrix M2) Multiplies the transpose of matrix M1 by M2 and places the result in this matrix.
void	mulTransposeRight(MatrixNd MR, MatrixNd M1) Computes the product
void	mulTransposeRight(SparseCRSMatrix M1, SparseCRSMatrix M2) Multiplies matrix M1 by the transpose of M2 and places the result in this matrix.
void	negate() Negates this matrix in place.
void	negate(SparseCRSMatrix M) Sets this matrix to the negative of another matrix M.
int	numNonZeroVals() Returns the number of non-zero values in this matrix object.
int	numNonZeroVals(Matrix.Partition part, int numRows, int numCols) Returns the number of non-zero values for a specified partition of a principal sub-matrix of this matrix delimited by the first numRows rows and the first numCols columns.
void	removeRow(int idx)
void	removeRows(int[] idxs)
int	rowSize() Number of rows in the matrix associated with this transformation.
void	scale(double s) Scales the elements of this matrix by s.
void	scale(double s, SparseCRSMatrix M) Scales the elements of matrix M1 by s and places the results in this matrix.
void	scaledAdd(double s, SparseCRSMatrix M1, SparseCRSMatrix M2) Computes s M1 + M2 and places the result in this matrix.
void	set(double[] values, int[] indices, int numVals) Sets the elements of this matrix from an array of doubles and a list of indices.
void	set(int i, int j, double val) Sets a single element of this matrix.
void	set(Matrix M) Sets the values of this matrix from an existing general matrix M.
void	set(SparseCRSMatrix M) Sets the values of this matrix from an existing SparseCRSMatrix M.
void	setCRSValues(double[] vals, int[] colIdxs, int[] rowOffs, int numVals, int numRows, Matrix.Partition part) Sets the contents of this matrix given a set of values in compressed row storage (CRS).
void	setDiagonal(double[] vals) Sets this matrix to a diagonal matrix whose diagonal elements are specified by an array.
void	setIdentity() Sets this matrix to the identity matrix.
void	setRandom() Sets this matrix to have a random sparsity structure and values, with between one and ten nonzero elements per row.
void	setRandomValues() Sets the nonzero elements of this matrix to random values in the range -0.5 (inclusive) to 0.5 (exclusive).
void	setRandomValues(double lower, double upper, java.util.Random generator) Sets the nonzero elements of this matrix to random values in a specified range.
void	setSize(int numRows, int numCols) Sets the size of this matrix.
void	setZero() Sets the elements of this matrix to zero.
void	sub(SparseCRSMatrix M) Subtracts M from this matrix.
void	sub(SparseCRSMatrix M1, SparseCRSMatrix M2) Subtracts M1 from M2 and places the result in this matrix.
java.lang.String	toString() Returns a String representation of this sparse matrix.
java.lang.String	toString(NumberFormat fmt) Returns a String representation of this sparse matrix.
java.lang.String	toString(java.lang.String fmtStr) Returns a String representation of this sparse matrix.
void	transpose() Sets this matrix its transpose, changing its size and sparsity structure if necessary.
void	transpose(SparseCRSMatrix M) Sets this matrix to the transpose of M, changing the size and sparsity structure of this matrix if necessary.
void	zeroColumn(int j)
void	zeroRow(int i)

Methods inherited from class maspack.matrix.SparseMatrixBase

setCCSValues, setRandom, setRandom, writeToFileCRS

Methods inherited from class maspack.matrix.MatrixBase

containsNaN, determinant, epsilonEquals, equals, frobeniusNorm, frobeniusNormSquared, get, get, getCCSIndices, getCCSIndices, getCCSValues, getCCSValues, getColumn, getColumn, getColumn, getCRSIndices, getCRSIndices, getCRSValues, getCRSValues, getDefaultFormat, getRow, getRow, getRow, getSize, getSubMatrix, hasNaN, idString, infinityNorm, isSymmetric, isWritable, mul, mul, mul, mulAdd, mulAdd, mulAdd, mulTranspose, mulTranspose, mulTranspose, mulTransposeAdd, mulTransposeAdd, mulTransposeAdd, oneNorm, scan, scan, setCRSValues, setDefaultFormat, trace, write, write, write, write, write, write, write, writeToFile

Methods inherited from class java.lang.Object

equals, getClass, hashCode, notify, notifyAll, wait, wait, wait

Methods inherited from interface maspack.matrix.Matrix

determinant, epsilonEquals, equals, frobeniusNorm, frobeniusNormSquared, get, getCCSIndices, getCCSIndices, getCCSValues, getCCSValues, getColumn, getColumn, getColumn, getCRSIndices, getCRSIndices, getCRSValues, getCRSValues, getRow, getRow, getRow, getSize, getSubMatrix, infinityNorm, isSymmetric, mul, mul, mul, mulAdd, mulAdd, mulAdd, mulTranspose, mulTranspose, mulTranspose, mulTransposeAdd, mulTransposeAdd, mulTransposeAdd, oneNorm, scan, trace, write, write, write

Constructor Detail

SparseCRSMatrix

public SparseCRSMatrix()

Creates an empty SparseCRSMatrix with 0 x 0 size.

SparseCRSMatrix

public SparseCRSMatrix(int numRows,
                       int numCols)

Creates an empty sparse matrix with a specified size.

Parameters:

numRows - number of rows

numCols - number of columns

SparseCRSMatrix

public SparseCRSMatrix(double[] vals,
                       int[] cols,
                       int[] rowOffs)

Creates a sparse matrix from CRS data. The number of columns is determined from the column data.

Parameters:

vals - if non-null, gives all the nonzero values, in row order

cols - column index for each value, in row order.

rowOffs - gives the starting offset of each row into the cols and vals arrays

Throws:

java.lang.IllegalArgumentException - if the specified storgae format is inconsistent

SparseCRSMatrix

public SparseCRSMatrix(Matrix M)

Creates a sparse matrix from an existing general matrix.

Parameters:

M - matrix to copy

Method Detail

setSize

public void setSize(int numRows,
                    int numCols)

Description copied from class: MatrixBase

Sets the size of this matrix. This operation is only supported if isFixedSize returns false.

Specified by:

setSize in interface Matrix

Overrides:

setSize in class MatrixBase

Parameters:

numRows - new row size

numCols - new column size

See Also:

Matrix.isFixedSize()

rowSize

public int rowSize()

Number of rows in the matrix associated with this transformation. This will equal the size of the output vector in mul.

Specified by:

rowSize in interface LinearTransformNd

Specified by:

rowSize in interface Matrix

Specified by:

rowSize in class MatrixBase

Returns:

number of rows

colSize

public int colSize()

Number of columns in the matrix associated with this transformation. This will equal the size of the input vector in mul.

Specified by:

colSize in interface LinearTransformNd

Specified by:

colSize in interface Matrix

Specified by:

colSize in class MatrixBase

Returns:

number of columns

isFixedSize

public boolean isFixedSize()

Returns true if this matrix is of fixed size. If this matrix is not of fixed size, then it can be resized dynamically, either explicitly using setSize, or implicitly when used as a result for various matrix operations.

Specified by:

isFixedSize in interface Matrix

Overrides:

isFixedSize in class MatrixBase

Returns:

true if this matrix is of fixed size

See Also:

Matrix.setSize(int, int)

get

public double get(int i,
                  int j)

Gets a single element of this matrix.

Specified by:

get in interface Matrix

Specified by:

get in class MatrixBase

Parameters:

i - element row index

j - element column index

Returns:

element value

get

public void get(MatrixNd M)

Copies the contents of this matrix into a dense MatrixNd.

Parameters:

M - matrix to receive the contents of this matrix

set

public void set(int i,
                int j,
                double val)

Sets a single element of this matrix. The row and column indices must lie within the current matrix dimenions. If the specified location does not currently contain an explicit element, a new element is added. Otherwise, the previous element is replaced.

Parameters:

i - element row index

j - element column index

val - new element value

set

public void set(double[] values,
                int[] indices,
                int numVals)

Sets the elements of this matrix from an array of doubles and a list of indices. Each value should be associated with two (i,j) index values, which, for the k-th value, are given by indices[k*2] and indices[k*2+1], respectively. All non-specified elements are set to zero.

Specified by:

set in interface Matrix

Parameters:

values - nonzero values for specific matrix locations

indices - i and j indices for each explicit value

numVals - number of explicit values

set

public void set(Matrix M)

Sets the values of this matrix from an existing general matrix M.

Specified by:

set in interface Matrix

Overrides:

set in class SparseMatrixBase

Parameters:

M - matrix used to set values

set

public void set(SparseCRSMatrix M)

Sets the values of this matrix from an existing SparseCRSMatrix M.

Parameters:

M - matrix used to set values

setIdentity

public void setIdentity()

Sets this matrix to the identity matrix. If this matrix matrix is not square, then element (i,j) is set to 1 if i and j are equal, and 0 otherwise.

This method does not trim matrix capacity.

setDiagonal

public void setDiagonal(double[] vals)

Sets this matrix to a diagonal matrix whose diagonal elements are specified by an array.

This method does not trim matrix capacity.

Parameters:

vals - diagonal elements for this matrix

Throws:

ImproperSizeException - if the length of diag does not equal the minimum matrix dimension

setZero

public void setZero()

Sets the elements of this matrix to zero. This method does not trim matrix capacity.

setRandom

public void setRandom()

Sets this matrix to have a random sparsity structure and values, with between one and ten nonzero elements per row.

setRandomValues

public void setRandomValues()

Sets the nonzero elements of this matrix to random values in the range -0.5 (inclusive) to 0.5 (exclusive).

setRandomValues

public void setRandomValues(double lower,
                            double upper,
                            java.util.Random generator)

Sets the nonzero elements of this matrix to random values in a specified range.

Parameters:

lower - lower random value (inclusive)

upper - upper random value (exclusive)

generator - random number generator

mul

public void mul(SparseCRSMatrix M)

Multiplies this matrix by M and places the result in this matrix. This matrix is resized if necessary.

Parameters:

M - right-hand matrix

Throws:

ImproperSizeException - if this matrix and M do not conform

mul

public void mul(SparseCRSMatrix M1,
                SparseCRSMatrix M2)
         throws ImproperSizeException

Multiplies matrix M1 by M2 and places the result in this matrix. This matrix is resized if necessary.

Parameters:

M1 - left-hand matrix

M2 - right-hand matrix

Throws:

ImproperSizeException - if this matrix and M1 and M2 do not conform.

mulTranspose

public void mulTranspose(SparseCRSMatrix M)

Multiplies this matrix by the transpose of M and places the result in this matrix. This matrix is resized if necessary.

Parameters:

M - right-hand matrix

Throws:

ImproperSizeException - if this matrix and the transpose of M do not conform

mulTransposeRight

public void mulTransposeRight(SparseCRSMatrix M1,
                              SparseCRSMatrix M2)
                       throws ImproperSizeException

Multiplies matrix M1 by the transpose of M2 and places the result in this matrix. This matrix is resized if necessary.

Parameters:

M1 - left-hand matrix

M2 - right-hand matrix

Throws:

ImproperSizeException - if M1 and the transpose M2 do not comform

mulTransposeLeft

public void mulTransposeLeft(SparseCRSMatrix M1,
                             SparseCRSMatrix M2)
                      throws ImproperSizeException

Multiplies the transpose of matrix M1 by M2 and places the result in this matrix. This matrix is resized if necessary.

Parameters:

M1 - left-hand matrix

M2 - right-hand matrix

Throws:

ImproperSizeException - if M1 and the transpose M2 do not comform

mulTransposeBoth

public void mulTransposeBoth(SparseCRSMatrix M1,
                             SparseCRSMatrix M2)
                      throws ImproperSizeException

Multiplies the transpose of matrix M1 by the transpose of M2 and places the result in this matrix. This matrix is resized if necessary.

Parameters:

M1 - left-hand matrix

M2 - right-hand matrix

Throws:

ImproperSizeException - if M1 and the transpose M2 do not comform

dotRowSelf

public double dotRowSelf(int i)

Returns the dot product of the i-th row with itself.

Parameters:

i - row index

Returns:

dot product of row i with itself

dotRow

public double dotRow(int i,
                     VectorNd v)

Returns the dot product of the i-th row with a vector v.

Parameters:

i - row index

v - vector to take dot product with

Returns:

dot product of row i with v

addScaledRow

public void addScaledRow(VectorNd v,
                         double s,
                         int i)

Adds s times the i-th row to the vector v.

Parameters:

v - vector to add scaled row to

s - scaling factor

i - row index

add

public void add(SparseCRSMatrix M)

Adds this matrix to M and places the result in this matrix.

Parameters:

M - right-hand matrix

Throws:

ImproperSizeException - if this matrix and M have different sizes.

add

public void add(SparseCRSMatrix M1,
                SparseCRSMatrix M2)

Adds matrix M1 to M2 and places the result in this matrix. This matrix is resized if necessary.

Parameters:

M1 - left-hand matrix

M2 - right-hand matrix

Throws:

ImproperSizeException - if matrices M1 and M2 have different sizes.

scaledAdd

public void scaledAdd(double s,
                      SparseCRSMatrix M1,
                      SparseCRSMatrix M2)

Computes s M1 + M2 and places the result in this matrix. This matrix is resized if necessary.

Parameters:

s - scaling factor

M1 - matrix to be scaled

M2 - matrix to be added

Throws:

ImproperSizeException - if matrices M1 and M2 have different sizes

sub

public void sub(SparseCRSMatrix M)

Subtracts M from this matrix.

Parameters:

M - matrix to subtract

Throws:

ImproperSizeException - if this matrix and M have different sizes.

sub

public void sub(SparseCRSMatrix M1,
                SparseCRSMatrix M2)

Subtracts M1 from M2 and places the result in this matrix. This matrix is resized if necessary.

Parameters:

M1 - left-hand matrix

M2 - right-hand matrix

Throws:

ImproperSizeException - if matrices M1 and M2 have different sizes.

combine

public void combine(double a,
                    SparseCRSMatrix M1,
                    double b,
                    SparseCRSMatrix M2)

Forms the combination a M1 + b M2 and places the result in this matrix. This matrix is resized if necessary.

Parameters:

a - scaling factor for M1

M1 - left-hand matrix

b - scaling factor for M2

M2 - right-hand matrix

Throws:

ImproperSizeException - if matrices M1 and M2 have different sizes.

negate

public void negate()

Negates this matrix in place.

negate

public void negate(SparseCRSMatrix M)

Sets this matrix to the negative of another matrix M. This matrix is resized if necessary.

Parameters:

M - matrix to be

scale

public void scale(double s)

Scales the elements of this matrix by s.

Specified by:

scale in interface SparseMatrix

Parameters:

s - scaling factor

scale

public void scale(double s,
                  SparseCRSMatrix M)

Scales the elements of matrix M1 by s and places the results in this matrix. This matrix is resized if necessary.

Parameters:

s - scaling factor

M - matrix to be scaled

transpose

public void transpose()

Sets this matrix its transpose, changing its size and sparsity structure if necessary.

transpose

public void transpose(SparseCRSMatrix M)

Sets this matrix to the transpose of M, changing the size and sparsity structure of this matrix if necessary.

Parameters:

M - matrix to take the transpose of

mul

public void mul(MatrixNd MR,
                MatrixNd M1)

Computes the product

 MR = MS M1
 

where MS is this matrix, and stores the result in MR, which is sized appropriately. If MR == M1 and the required size of MR does not equal that of M1, then an internal copy of M1 is made for computational purposes and M1 is resized as required.

mulLeft

public void mulLeft(MatrixNd MR,
                    MatrixNd M1)

Computes the product

 MR = M1 MS
 

where MS is this matrix, and stores the result in MR, which is sized appropriately. If MR == M1 and the required size of MR does not equal that of M1, then an internal copy of M1 is made for computational purposes and M1 is resized as required.

mulTransposeRight

public void mulTransposeRight(MatrixNd MR,
                              MatrixNd M1)

Computes the product

 MR = MS M1^T
 

here MS is this matrix, and stores the result in MR, which will be sized appropriately. If MR == M1, then an internal copy will be made of M1 for computational purposes and M1 will be resized as needed to store the result.

toString

public java.lang.String toString()

Returns a String representation of this sparse matrix. Each explicit element is output in row-major order as a tuple ( i j value), one tuple per line. Elements which are not printed are zero.

Overrides:

toString in class MatrixBase

Returns:

String representation of this matrix

See Also:

MatrixBase.toString(String)

toString

public java.lang.String toString(java.lang.String fmtStr)

Returns a String representation of this sparse matrix. The output format is the same as that produced by toString, except that elements values are themselves formatted using a C printf style format string. For a description of the format string syntax, see NumberFormat.

Specified by:

toString in interface Matrix

Overrides:

toString in class MatrixBase

Parameters:

fmtStr - numeric format specification

Returns:

String representation of this matrix

See Also:

NumberFormat

toString

public java.lang.String toString(NumberFormat fmt)

Returns a String representation of this sparse matrix. The output format is the same as that produced by toString, except that elements values are themselves formatted using a C printf style as decribed by the parameter NumberFormat. When called numerous times, this routine can be more efficient than toString(String), because the NumberFormat does not need to be recreated each time from a specification string.

Specified by:

toString in interface Matrix

Overrides:

toString in class MatrixBase

Parameters:

fmt - numeric format

Returns:

String representation of this matrix

numNonZeroVals

public int numNonZeroVals()

Description copied from class: MatrixBase

Returns the number of non-zero values in this matrix object. For dense matrices, this will simply be the product of the number of rows times the number of columns.

Specified by:

numNonZeroVals in interface Matrix

Overrides:

numNonZeroVals in class MatrixBase

Returns:

number of non-zero values

numNonZeroVals

public int numNonZeroVals(Matrix.Partition part,
                          int numRows,
                          int numCols)

Description copied from class: MatrixBase

Returns the number of non-zero values for a specified partition of a principal sub-matrix of this matrix delimited by the first numRows rows and the first numCols columns. If the matrix is dense and the partition is Full, then this will simply be the product of the number of rows times the number of columns.

Specified by:

numNonZeroVals in interface Matrix

Overrides:

numNonZeroVals in class MatrixBase

Parameters:

part - matrix parition to be examined

numRows - number of rows delimiting the sub-matrix

numCols - number of columns delimiting the sub-matrix

Returns:

number of non-zero values

setCRSValues

public void setCRSValues(double[] vals,
                         int[] colIdxs,
                         int[] rowOffs,
                         int numVals,
                         int numRows,
                         Matrix.Partition part)

Sets the contents of this matrix given a set of values in compressed row storage (CRS). The argument part specifies what portion of the matrix is supplied. If it is Full, then the entire matrix is supplied. If it is Matrix.Partition.UpperTriangular, then only the upper triangular part is assumed to be supplied, and the lower triangular part is set from its transpose (matrix bounds permitting). For Partition.UpperTriangular, supplied entries which are not upper triangular will either be ignored or generate an exception.

All specified data must fit within the current matrix bounds; the matrix is not resized.

CRS data takes the form of three arrays:

vals

An array giving the values of all the non-zero elements in the sparse matrix, in row-major order.

colIdxs

An array giving the column indices of all the non-zero elements in the sparse matrix, in row-major order. The indices are 1-based, so that 1 denotes the first column.

rowOffs

An array of size nrows+1, where nrows is the number of matrix rows, giving the offsets into vals and colIdxs corresponding to the first non-zero element in each row. All values are 1-based, so that first offset value is 1, the second offset value is n+1, where n is the number of non-zero elements in the first row, etc. The final value is set to nvals+1, where nvals is the number of non-zero elements in the matrix.

CRS index data is 1-based because this is the most standard form used by solvers and matrix formats.

Specified by:

setCRSValues in interface Matrix

Overrides:

setCRSValues in class SparseMatrixBase

Parameters:

vals - non-zero element values. This array must have a length equal at least to nvals.

colIdxs - column indices for each non-zero element. This array must have a length equal at least to nvals.

rowOffs - row start offsets into vals and colIdxs. This array must have a length equal at least to nrows+1, where nrows is the number of rows in this matrix.

numVals - number of non-zero values

numRows - number of specified rows. Indicates the maximum value of rowOffs that will be used; will not resize matrix.

part - must be either Full or UpperTriangular.

getCRSIndices

public int getCRSIndices(int[] colIdxs,
                         int[] rowOffs,
                         Matrix.Partition part,
                         int numRows,
                         int numCols)

Gets the compressed row storage (CRS) indices for a principal sub-matrix of this matrix delimited by the first numRows rows and the first numCols columns. Indices are 1-based and supplied in row-major order. For a detailed description of the CRS format, see setCRSValues(). Some matrix types may place restrictions on the sub-matrix; for instance, a block-structured matrix may require that the sub-matrix be block aligned.

Specified by:

getCRSIndices in interface Matrix

Overrides:

getCRSIndices in class MatrixBase

Parameters:

colIdxs - returns the column indices of each non-zero element, in row-major order. This array must have a length equal at least to the number of non-zero elements.

rowOffs - returns the row start offsets into colIdxs, followed by nvals+1, where nvals is the number of non-zero elements in the sub-matrix. This array must have a length equal at least to nrows+1, where nrows is the number of rows in the sub-matrix.

part - specifies what portion of the sub-matrix to store; must be either Full or UpperTriangular

numRows - number of rows delimiting the sub-matrix

numCols - number of columns delimiting the sub-matrix

Returns:

number of non-zero elements

getCRSValues

public int getCRSValues(double[] vals,
                        Matrix.Partition part,
                        int numRows,
                        int numCols)

Gets the compressed row storage (CRS) values for a principal sub-matrix of this matrix delimited by the first numRows rows and the first numCols columns. Values are supplied in row-major order. For a detailed decsription of the CRS format, see setCRSValues(). Some matrix types may place restrictions on the sub-matrix; for instance, a block-structured matrix may require that the sub-matrix be block aligned.

Specified by:

getCRSValues in interface Matrix

Overrides:

getCRSValues in class MatrixBase

Parameters:

vals - returns the value of each non-zero element. This array must have a length equal at least to the number of non-zero elements in the sub-matrix.

part - specifies what portion of the sub-matrix to store; must be either Full or UpperTriangular

numRows - number of rows delimiting the sub-matrix

numCols - number of columns delimiting the sub-matrix

Returns:

number of non-zero elements

getCCSIndices

public int getCCSIndices(int[] rowIdxs,
                         int[] colOffs,
                         Matrix.Partition part,
                         int numRows,
                         int numCols)

Gets the compressed column storage (CCS) indices for a principal sub-matrix of this matrix delimited by the first numRows rows and the first numCols columns. Indices are 1-based and supplied in column-major order. For a detailed decsription of the CCS format, see setCCSValues(). Some matrix types may place restrictions on the sub-matrix; for instance, a block-structured matrix may require that the sub-matrix be block aligned.

Specified by:

getCCSIndices in interface Matrix

Overrides:

getCCSIndices in class MatrixBase

Parameters:

rowIdxs - returns the row indices of each non-zero element, in column-major order. This array must have a length equal at least to the number of non-zero elements.

colOffs - returns the column start offsets into rowIdxs, followed by nvals+1, where nvals is the number of non-zero elements in the sub-matrix. This array must have a length equal at least to ncol+1, where ncols is the number of columns in the sub-matrix.

part - specifies what portion of the sub-matrix to store; must be either Full or LowerTriangular

numRows - number of rows delimiting the sub-matrix

numCols - number of columns delimiting the sub-matrix

Returns:

number of non-zero elements

getCCSValues

public int getCCSValues(double[] vals,
                        Matrix.Partition part,
                        int numRows,
                        int numCols)

Gets the compressed column storage (CCS) values for a principal sub-matrix of this matrix delimited by the first numRows rows and the first numCols columns. Values are supplied in column-major order. For a detailed decsription of the CCS format, see setCCSValues(). Some matrix types may place restrictions on the sub-matrix; for instance, a block-structured matrix may require that the sub-matrix be block aligned.

Specified by:

getCCSValues in interface Matrix

Overrides:

getCCSValues in class MatrixBase

Parameters:

vals - returns the value of each non-zero element. This array must have a length equal at least to the number of non-zero elements in the sub-matrix.

part - specifies what portion of the sub-matrix to store; must be either Full or LowerTriangular

numRows - number of rows delimiting the sub-matrix

numCols - number of columns delimiting the sub-matrix

Returns:

number of non-zero elements

checkConsistency

public void checkConsistency()

Check that the internal structures of this matrix are consistent. Used for testing.

Specified by:

checkConsistency in interface Matrix

clone

public SparseCRSMatrix clone()

Specified by:

clone in interface Clonable

Overrides:

clone in class java.lang.Object

maxNorm

public double maxNorm()

Description copied from interface: Matrix

Returns the max-norm of this matrix. This is equal to the maximum of the absolute values of all entries (i.e. element-wise, with p=q=inf)

Specified by:

maxNorm in interface Matrix

Overrides:

maxNorm in class MatrixBase

Returns:

max norm of this matrix

zeroRow

public void zeroRow(int i)

zeroColumn

public void zeroColumn(int j)

add

public void add(int i,
                int j,
                double v)

removeRow

public void removeRow(int idx)

removeRows

public void removeRows(int[] idxs)

addRow

public void addRow(int[] colIdxs,
                   double[] values,
                   int numVals)

addRow

public void addRow(DynamicIntArray colIdxs,
                   DynamicDoubleArray values)