Package no.uib.cipr.matrix
Class LowerSymmDenseMatrix
- java.lang.Object
-
- no.uib.cipr.matrix.AbstractMatrix
-
- no.uib.cipr.matrix.LowerSymmDenseMatrix
-
- All Implemented Interfaces:
java.lang.Iterable<MatrixEntry>,Matrix
- Direct Known Subclasses:
LowerSPDDenseMatrix
public class LowerSymmDenseMatrix extends AbstractMatrix
Lower symmetric dense matrix. It has the same storage layout as theDenseMatrix, but only refers to elements below or on the main diagonal. The remaining elements are never accessed nor changed, and is known only by symmetry.
-
-
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 LowerSymmDenseMatrix(int n)Constructor for LowerSymmDenseMatrixLowerSymmDenseMatrix(Matrix A)Constructor for LowerSymmDenseMatrixLowerSymmDenseMatrix(Matrix A, boolean deep)Constructor for LowerSymmDenseMatrix
-
Method Summary
All Methods Instance Methods Concrete Methods Modifier and Type Method Description voidadd(int row, int column, double value)A(row,column) += valueLowerSymmDenseMatrixcopy()Creates a deep copy of the matrixdoubleget(int row, int column)ReturnsA(row,column)double[]getData()Returns the matrix contents.MatrixmultAdd(double alpha, Matrix B, Matrix C)C = alpha*A*B + CVectormultAdd(double alpha, Vector x, Vector y)y = alpha*A*x + yMatrixrank1(double alpha, Matrix C)A = alpha*C*CT + A.Matrixrank1(double alpha, Vector x, Vector y)A = alpha*x*yT + A.Matrixrank2(double alpha, Matrix B, Matrix C)A = alpha*B*CT + alpha*C*BT + A.Matrixrank2(double alpha, Vector x, Vector y)A = alpha*x*yT + alpha*y*xT + A.voidset(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.MatrixtransAmultAdd(double alpha, Matrix B, Matrix C)C = alpha*AT*B + CVectortransMultAdd(double alpha, Vector x, Vector y)y = alpha*AT*x + yMatrixtranspose()Transposes the matrix in-place.MatrixtransRank1(double alpha, Matrix C)A = alpha*CT*C + AThe matrices must be square and of the same sizeMatrixtransRank2(double alpha, Matrix B, Matrix C)A = alpha*BT*C + alpha*CT*B + A.MatrixtransSolve(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, iterator, max, max, mult, mult, mult, mult, multAdd, multAdd, norm, norm1, normF, normInf, numColumns, numRows, rank1, rank1, rank1, rank1, rank2, rank2, scale, set, toString, transABmult, transABmult, transABmultAdd, transABmultAdd, transAmult, transAmult, transAmultAdd, transBmult, transBmult, transBmultAdd, transBmultAdd, transMult, transMult, transMultAdd, transpose, transRank1, transRank2
-
-
-
-
Constructor Detail
-
LowerSymmDenseMatrix
public LowerSymmDenseMatrix(int n)
Constructor for LowerSymmDenseMatrix- Parameters:
n- Size of the matrix. Since the matrix must be square, this equals both the number of rows and columns
-
LowerSymmDenseMatrix
public LowerSymmDenseMatrix(Matrix A)
Constructor for LowerSymmDenseMatrix- Parameters:
A- Matrix to copy. It must be a square matrix, and only the lower triangular part is copied
-
LowerSymmDenseMatrix
public LowerSymmDenseMatrix(Matrix A, boolean deep)
Constructor for LowerSymmDenseMatrix- Parameters:
A- Matrix to copy. It must be a square matrix, and only the lower triangular part is copieddeep- If false, a shallow copy is made. In that case,Amust be a dense matrix
-
-
Method Detail
-
add
public void add(int row, int column, double value)Description copied from interface:MatrixA(row,column) += value
-
get
public double get(int row, int column)Description copied from interface:MatrixReturnsA(row,column)
-
set
public void set(int row, int column, double value)Description copied from interface:MatrixA(row,column) = value
-
copy
public LowerSymmDenseMatrix copy()
Description copied from interface:MatrixCreates a deep copy of the matrix- Specified by:
copyin interfaceMatrix- Overrides:
copyin classAbstractMatrix- Returns:
- A
-
multAdd
public Matrix multAdd(double alpha, Matrix B, Matrix C)
Description copied from interface:MatrixC = alpha*A*B + C- Specified by:
multAddin interfaceMatrix- Overrides:
multAddin classAbstractMatrixB- Matrix such thatB.numRows() == A.numColumns()andB.numColumns() == C.numColumns()C- Matrix such thatC.numRows() == A.numRows()andB.numColumns() == C.numColumns()- Returns:
- C
-
transAmultAdd
public Matrix transAmultAdd(double alpha, Matrix B, Matrix C)
Description copied from interface:MatrixC = alpha*AT*B + C- Specified by:
transAmultAddin interfaceMatrix- Overrides:
transAmultAddin classAbstractMatrixB- Matrix such thatB.numRows() == A.numRows()andB.numColumns() == C.numColumns()C- Matrix such thatC.numRows() == A.numColumns()andB.numColumns() == C.numColumns()- Returns:
- C
-
rank1
public Matrix rank1(double alpha, Vector x, Vector y)
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
public Matrix rank2(double alpha, Vector x, Vector y)
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
-
multAdd
public Vector multAdd(double alpha, Vector x, Vector y)
Description copied from interface:Matrixy = alpha*A*x + y- Specified by:
multAddin interfaceMatrix- Overrides:
multAddin classAbstractMatrixx- Vector of sizeA.numColumns()y- Vector of sizeA.numRows()- Returns:
- y
-
transMultAdd
public Vector transMultAdd(double alpha, Vector x, Vector y)
Description copied from interface:Matrixy = alpha*AT*x + y- Specified by:
transMultAddin interfaceMatrix- Overrides:
transMultAddin classAbstractMatrixx- Vector of sizeA.numRows()y- Vector of sizeA.numColumns()- Returns:
- y
-
rank1
public Matrix rank1(double alpha, Matrix C)
Description copied from interface:MatrixA = alpha*C*CT + A. The matrices must be square and of the same size- Specified by:
rank1in interfaceMatrix- Overrides:
rank1in classAbstractMatrix- Returns:
- A
-
transRank1
public Matrix transRank1(double alpha, Matrix C)
Description copied from interface:MatrixA = alpha*CT*C + AThe matrices must be square and of the same size- Specified by:
transRank1in interfaceMatrix- Overrides:
transRank1in classAbstractMatrix- Returns:
- A
-
rank2
public Matrix rank2(double alpha, Matrix B, Matrix C)
Description copied from interface:MatrixA = alpha*B*CT + alpha*C*BT + A. This matrix must be square- Specified by:
rank2in interfaceMatrix- Overrides:
rank2in classAbstractMatrixB- Matrix with the same number of rows asAand the same number of columns asCC- Matrix with the same number of rows asAand the same number of columns asB- Returns:
- A
-
transRank2
public Matrix transRank2(double alpha, Matrix B, Matrix C)
Description copied from interface:MatrixA = alpha*BT*C + alpha*CT*B + A. This matrix must be square- Specified by:
transRank2in interfaceMatrix- Overrides:
transRank2in classAbstractMatrixB- Matrix with the same number of rows asCand the same number of columns asAC- Matrix with the same number of rows asBand the same number of columns asA- Returns:
- A
-
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
-
transpose
public Matrix 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
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
-
-