Package no.uib.cipr.matrix
Class LowerTriangPackMatrix
- java.lang.Object
-
- no.uib.cipr.matrix.AbstractMatrix
-
- no.uib.cipr.matrix.LowerTriangPackMatrix
-
- All Implemented Interfaces:
java.lang.Iterable<MatrixEntry>,Matrix
- Direct Known Subclasses:
UnitLowerTriangPackMatrix
public class LowerTriangPackMatrix extends AbstractMatrix
Lower triangular packed matrix. In contrast withLowerTriangDenseMatrix, this matrix exploits the sparsity by only storing about half the matrix. As such, the triangular matrixa11 a21 a22 a31 a32 a33 a41 a42 a43 a44 is packed as follows:
a11 a21 a31 a41 a22 a32 a42 a33 a43 a44
-
-
Nested Class Summary
-
Nested classes/interfaces inherited from interface no.uib.cipr.matrix.Matrix
Matrix.Norm
-
-
Field Summary
-
Fields inherited from class no.uib.cipr.matrix.AbstractMatrix
numColumns, numRows
-
-
Constructor Summary
Constructors Constructor Description LowerTriangPackMatrix(int n)Constructor for LowerTriangPackMatrixLowerTriangPackMatrix(Matrix A)Constructor for LowerTriangPackMatrixLowerTriangPackMatrix(Matrix A, boolean deep)Constructor for LowerTriangPackMatrix
-
Method Summary
All Methods Instance Methods Concrete Methods Modifier and Type Method Description voidadd(int row, int column, double value)A(row,column) += valueLowerTriangPackMatrixcopy()Creates a deep copy of the matrixdoubleget(int row, int column)ReturnsA(row,column)double[]getData()Returns the matrix contents.java.util.Iterator<MatrixEntry>iterator()Vectormult(double alpha, Vector x, Vector y)y = alpha*A*xvoidset(int row, int column, double value)A(row,column) = valueMatrixset(Matrix B)A=B.Matrixsolve(Matrix B, Matrix X)X = A\B.Vectorsolve(Vector b, Vector x)x = A\b.VectortransMult(double alpha, Vector x, Vector y)y = alpha*AT*xMatrixtransSolve(Matrix B, Matrix X)X = AT\B.VectortransSolve(Vector b, Vector x)x = AT\b.Matrixzero()Zeros all the entries in the matrix, while preserving any underlying structure.-
Methods inherited from class no.uib.cipr.matrix.AbstractMatrix
add, add, check, checkMultAdd, checkMultAdd, checkRank1, checkRank1, checkRank2, checkRank2, checkSize, checkSolve, checkSolve, checkTransABmultAdd, checkTransAmultAdd, checkTransBmultAdd, checkTransMultAdd, checkTranspose, checkTranspose, checkTransRank1, checkTransRank2, isSquare, max, max, mult, mult, mult, multAdd, multAdd, multAdd, multAdd, norm, norm1, normF, normInf, numColumns, numRows, rank1, rank1, rank1, rank1, rank1, rank1, rank2, rank2, rank2, rank2, scale, set, toString, transABmult, transABmult, transABmultAdd, transABmultAdd, transAmult, transAmult, transAmultAdd, transAmultAdd, transBmult, transBmult, transBmultAdd, transBmultAdd, transMult, transMultAdd, transMultAdd, transpose, transpose, transRank1, transRank1, transRank2, transRank2
-
-
-
-
Constructor Detail
-
LowerTriangPackMatrix
public LowerTriangPackMatrix(int n)
Constructor for LowerTriangPackMatrix- Parameters:
n- Size of the matrix. Since the matrix must be square, this equals both the number of rows and columns
-
LowerTriangPackMatrix
public LowerTriangPackMatrix(Matrix A)
Constructor for LowerTriangPackMatrix- Parameters:
A- Matrix to copy contents from. Only the entries of the relevant part are copied
-
LowerTriangPackMatrix
public LowerTriangPackMatrix(Matrix A, boolean deep)
Constructor for LowerTriangPackMatrix- Parameters:
A- Matrix to copy contents from. Only the entries of the relevant part are copieddeep- True if the copy is deep, else false (giving a shallow copy). For shallow copies,Amust be a packed matrix
-
-
Method Detail
-
add
public void add(int row, int column, double value)Description copied from interface:MatrixA(row,column) += value- Specified by:
addin interfaceMatrix- Overrides:
addin classAbstractMatrix
-
set
public void set(int row, int column, double value)Description copied from interface:MatrixA(row,column) = value- Specified by:
setin interfaceMatrix- Overrides:
setin classAbstractMatrix
-
get
public double get(int row, int column)Description copied from interface:MatrixReturnsA(row,column)- Specified by:
getin interfaceMatrix- Overrides:
getin classAbstractMatrix
-
copy
public LowerTriangPackMatrix copy()
Description copied from interface:MatrixCreates a deep copy of the matrix- Specified by:
copyin interfaceMatrix- Overrides:
copyin classAbstractMatrix- Returns:
- A
-
mult
public Vector mult(double alpha, Vector x, Vector y)
Description copied from interface:Matrixy = alpha*A*x- Specified by:
multin interfaceMatrix- Overrides:
multin classAbstractMatrixx- Vector of sizeA.numColumns()y- Vector of sizeA.numRows()- Returns:
- y
-
transMult
public Vector transMult(double alpha, Vector x, Vector y)
Description copied from interface:Matrixy = alpha*AT*x- Specified by:
transMultin interfaceMatrix- Overrides:
transMultin classAbstractMatrixx- Vector of sizeA.numRows()y- Vector of sizeA.numColumns()- Returns:
- y
-
solve
public Matrix solve(Matrix B, Matrix X)
Description copied from interface:MatrixX = A\B. Not all matrices support this operation, those that do not throwUnsupportedOperationException. Note that it is often more efficient to use a matrix decomposition and its associated solver- Specified by:
solvein interfaceMatrix- Overrides:
solvein classAbstractMatrix- Parameters:
B- Matrix with the same number of rows asA, and the same number of columns asXX- Matrix with a number of rows equalA.numColumns(), and the same number of columns asB- Returns:
- X
-
solve
public Vector solve(Vector b, Vector x)
Description copied from interface:Matrixx = A\b. Not all matrices support this operation, those that do not throwUnsupportedOperationException. Note that it is often more efficient to use a matrix decomposition and its associated solver- Specified by:
solvein interfaceMatrix- Overrides:
solvein classAbstractMatrix- Parameters:
b- Vector of sizeA.numRows()x- Vector of sizeA.numColumns()- Returns:
- x
-
transSolve
public Matrix transSolve(Matrix B, Matrix X)
Description copied from interface:MatrixX = AT\B. Not all matrices support this operation, those that do not throwUnsupportedOperationException. Note that it is often more efficient to use a matrix decomposition and its associated transpose solver- Specified by:
transSolvein interfaceMatrix- Overrides:
transSolvein classAbstractMatrix- Parameters:
B- Matrix with a number of rows equalA.numColumns(), and the same number of columns asXX- Matrix with the same number of rows asA, and the same number of columns asB- Returns:
- X
-
transSolve
public Vector transSolve(Vector b, Vector x)
Description copied from interface:Matrixx = AT\b. Not all matrices support this operation, those that do not throwUnsupportedOperationException. Note that it is often more efficient to use a matrix decomposition and its associated solver- Specified by:
transSolvein interfaceMatrix- Overrides:
transSolvein classAbstractMatrix- Parameters:
b- Vector of sizeA.numColumns()x- Vector of sizeA.numRows()- Returns:
- x
-
iterator
public java.util.Iterator<MatrixEntry> iterator()
- Specified by:
iteratorin interfacejava.lang.Iterable<MatrixEntry>- Overrides:
iteratorin classAbstractMatrix
-
getData
public double[] getData()
Returns the matrix contents. Ordering depends on the underlying storage assumptions
-
set
public Matrix set(Matrix B)
Description copied from interface:MatrixA=B. The matrices must be of the same size- Specified by:
setin interfaceMatrix- Overrides:
setin classAbstractMatrix- Returns:
- A
-
-