Package no.uib.cipr.matrix
Class UpperSymmDenseMatrix
java.lang.Object
no.uib.cipr.matrix.AbstractMatrix
no.uib.cipr.matrix.UpperSymmDenseMatrix
- All Implemented Interfaces:
Iterable<MatrixEntry>,Matrix
- Direct Known Subclasses:
UpperSPDDenseMatrix
Upper symmetrix dense matrix. It has the same storage layout as the
DenseMatrix, but only refers to
elements above 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
ConstructorsConstructorDescriptionUpperSymmDenseMatrix(int n) Constructor for UpperSymmDenseMatrixConstructor for UpperSymmDenseMatrixUpperSymmDenseMatrix(Matrix A, boolean deep) Constructor for UpperSymmDenseMatrix -
Method Summary
Modifier and TypeMethodDescriptionvoidadd(int row, int column, double value) A(row,column) += valuecopy()Creates a deep copy of the matrixdoubleget(int row, int column) ReturnsA(row,column)double[]getData()Returns the matrix contents.C = alpha*A*B + Cy = alpha*A*x + yA = alpha*C*CT + A.A = alpha*x*yT + A.A = alpha*B*CT + alpha*C*BT + A.A = alpha*x*yT + alpha*y*xT + A.voidset(int row, int column, double value) A(row,column) = valueA=B.X = A\B.x = A\b.transAmultAdd(double alpha, Matrix B, Matrix C) C = alpha*AT*B + CtransMultAdd(double alpha, Vector x, Vector y) y = alpha*AT*x + yTransposes the matrix in-place.transRank1(double alpha, Matrix C) A = alpha*CT*C + AThe matrices must be square and of the same sizetransRank2(double alpha, Matrix B, Matrix C) A = alpha*BT*C + alpha*CT*B + A.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.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, 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
-
UpperSymmDenseMatrix
public UpperSymmDenseMatrix(int n) Constructor for UpperSymmDenseMatrix- Parameters:
n- Size of the matrix. Since the matrix must be square, this equals both the number of rows and columns
-
UpperSymmDenseMatrix
Constructor for UpperSymmDenseMatrix- Parameters:
A- Matrix to copy. It must be a square matrix, and only the upper triangular part is copied
-
UpperSymmDenseMatrix
Constructor for UpperSymmDenseMatrix- Parameters:
A- Matrix to copy. It must be a square matrix, and only the upper triangular part is copieddeep- If false, a shallow copy is made. In that case,Amust be a dense matrix
-
-
Method Details
-
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
Description copied from interface:MatrixCreates a deep copy of the matrix- Specified by:
copyin interfaceMatrix- Overrides:
copyin classAbstractMatrix- Returns:
- A
-
multAdd
Description copied from interface:MatrixC = alpha*A*B + C- Specified by:
multAddin interfaceMatrix- Overrides:
multAddin classAbstractMatrix- Parameters:
B- Matrix such thatB.numRows() == A.numColumns()andB.numColumns() == C.numColumns()C- Matrix such thatC.numRows() == A.numRows()andB.numColumns() == C.numColumns()- Returns:
- C
-
transAmultAdd
Description copied from interface:MatrixC = alpha*AT*B + C- Specified by:
transAmultAddin interfaceMatrix- Overrides:
transAmultAddin classAbstractMatrix- Parameters:
B- Matrix such thatB.numRows() == A.numRows()andB.numColumns() == C.numColumns()C- Matrix such thatC.numRows() == A.numColumns()andB.numColumns() == C.numColumns()- Returns:
- C
-
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
-
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*C*CT + A. The matrices must be square and of the same size- Specified by:
rank1in interfaceMatrix- Overrides:
rank1in classAbstractMatrix- Returns:
- A
-
transRank1
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
Description copied from interface:MatrixA = alpha*B*CT + alpha*C*BT + A. This matrix must be square- Specified by:
rank2in interfaceMatrix- Overrides:
rank2in classAbstractMatrix- Parameters:
B- 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
Description copied from interface:MatrixA = alpha*BT*C + alpha*CT*B + A. This matrix must be square- Specified by:
transRank2in interfaceMatrix- Overrides:
transRank2in classAbstractMatrix- Parameters:
B- 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
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
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
-