Package no.uib.cipr.matrix
Class LowerSPDPackMatrix
java.lang.Object
no.uib.cipr.matrix.AbstractMatrix
no.uib.cipr.matrix.LowerSymmPackMatrix
no.uib.cipr.matrix.LowerSPDPackMatrix
- All Implemented Interfaces:
Iterable<MatrixEntry>,Matrix
Lower symmetrical positive definite packed matrix. Same layout as
LowerSymmPackMatrix. This
class does not enforce the SPD property, but serves as a tag so that more
efficient algorithms can be used in the solvers.-
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
ConstructorsConstructorDescriptionLowerSPDPackMatrix(int n) Constructor for LowerSPDPackMatrixConstructor for LowerSPDPackMatrixLowerSPDPackMatrix(Matrix A, boolean deep) Constructor for LowerSPDPackMatrix -
Method Summary
Modifier and TypeMethodDescriptioncopy()Creates a deep copy of the matrixdouble[]getData()Returns the matrix contents.y = alpha*A*x + yA = alpha*x*yT + A.A = alpha*x*yT + alpha*y*xT + A.A=B.X = A\B.x = A\b.transMultAdd(double alpha, Vector x, Vector y) y = alpha*AT*x + yTransposes the matrix in-place.transSolve(Matrix B, Matrix X) X = AT\B.transSolve(Vector b, Vector x) x = AT\b.zero()Zeros all the entries in the matrix, while preserving any underlying structure.Methods inherited from class no.uib.cipr.matrix.LowerSymmPackMatrix
add, get, setMethods 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, iterator, max, max, mult, mult, mult, mult, multAdd, multAdd, multAdd, norm, norm1, normF, normInf, numColumns, numRows, rank1, rank1, rank1, rank1, rank1, rank2, rank2, rank2, scale, set, toString, transABmult, transABmult, transABmultAdd, transABmultAdd, transAmult, transAmult, transAmultAdd, transAmultAdd, transBmult, transBmult, transBmultAdd, transBmultAdd, transMult, transMult, transMultAdd, transpose, transRank1, transRank1, transRank2, transRank2Methods inherited from class java.lang.Object
clone, equals, finalize, getClass, hashCode, notify, notifyAll, wait, wait, waitMethods inherited from interface java.lang.Iterable
forEach, spliterator
-
Constructor Details
-
LowerSPDPackMatrix
public LowerSPDPackMatrix(int n) Constructor for LowerSPDPackMatrix- Parameters:
n- Size of the matrix. Since the matrix must be square, this equals both the number of rows and columns
-
LowerSPDPackMatrix
Constructor for LowerSPDPackMatrix- Parameters:
A- Matrix to copy contents from. Only the entries of the relevant part are copied
-
LowerSPDPackMatrix
Constructor for LowerSPDPackMatrix- 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 Details
-
copy
Description copied from interface:MatrixCreates a deep copy of the matrix- Specified by:
copyin interfaceMatrix- Overrides:
copyin classLowerSymmPackMatrix- Returns:
- A
-
solve
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 -
multAdd
Description copied from interface:Matrixy = alpha*A*x + y- Specified by:
multAddin interfaceMatrix- Overrides:
multAddin classAbstractMatrix- Parameters:
x- Vector of sizeA.numColumns()y- Vector of sizeA.numRows()- Returns:
- y
-
transMultAdd
Description copied from interface:Matrixy = alpha*AT*x + y- Specified by:
transMultAddin interfaceMatrix- Overrides:
transMultAddin classAbstractMatrix- Parameters:
x- Vector of sizeA.numRows()y- Vector of sizeA.numColumns()- Returns:
- y
-
rank1
Description copied from interface:MatrixA = alpha*x*yT + A. The matrix must be square, and the vectors of the same length- Specified by:
rank1in interfaceMatrix- Overrides:
rank1in classAbstractMatrix- Returns:
- A
-
rank2
Description copied from interface:MatrixA = alpha*x*yT + alpha*y*xT + A. The matrix must be square, and the vectors of the same length- Specified by:
rank2in interfaceMatrix- Overrides:
rank2in classAbstractMatrix- Returns:
- A
-
solve
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
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
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
-
transpose
Description copied from interface:MatrixTransposes the matrix in-place. In most cases, the matrix must be square for this to work.- Specified by:
transposein interfaceMatrix- Overrides:
transposein classAbstractMatrix- Returns:
- This matrix
-
getData
public double[] getData()Returns the matrix contents. Ordering depends on the underlying storage assumptions -
set
Description copied from interface:MatrixA=B. The matrices must be of the same size- Specified by:
setin interfaceMatrix- Overrides:
setin classAbstractMatrix- Returns:
- A
-
zero
Description copied from interface:MatrixZeros all the entries in the matrix, while preserving any underlying structure. Useful for general, unstructured matrices.- Specified by:
zeroin interfaceMatrix- Overrides:
zeroin classAbstractMatrix- Returns:
- A
-