Package no.uib.cipr.matrix
Class DenseMatrix
- java.lang.Object
-
- no.uib.cipr.matrix.AbstractMatrix
-
- no.uib.cipr.matrix.DenseMatrix
-
- All Implemented Interfaces:
java.lang.Iterable<MatrixEntry>,Matrix
public class DenseMatrix extends AbstractMatrix
Dense matrix. It is a good all-round matrix structure, with fast access and efficient algebraic operations. The matrixa11 a12 a13 a14 a21 a22 a23 a24 a31 a32 a33 a34 a41 a42 a43 a44 is stored column major in a single array, as follows:
a11 a21 a31 a41 a12 a22 a32 a42 a13 a23 a33 a43 a14 a24 a34 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 DenseMatrix(double[][] values)Constructor for DenseMatrix.DenseMatrix(int numRows, int numColumns)Constructor for DenseMatrixDenseMatrix(MatrixVectorReader r)Constructor for DenseMatrixDenseMatrix(Matrix A)Constructor for DenseMatrixDenseMatrix(Matrix A, boolean deep)Constructor for DenseMatrixDenseMatrix(Vector x)Constructor for DenseMatrix.DenseMatrix(Vector[] x)Constructor for DenseMatrix.DenseMatrix(Vector x, boolean deep)Constructor for DenseMatrix.
-
Method Summary
All Methods Instance Methods Concrete Methods Modifier and Type Method Description voidadd(int row, int column, double value)A(row,column) += valueDenseMatrixcopy()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, Vector x, Vector y)A = alpha*x*yT + 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.MatrixtransABmultAdd(double alpha, Matrix B, Matrix C)C = alpha*AT*BT + CMatrixtransAmultAdd(double alpha, Matrix B, Matrix C)C = alpha*AT*B + CMatrixtransBmultAdd(double alpha, Matrix B, Matrix C)C = alpha*A*BT + CVectortransMultAdd(double alpha, Vector x, Vector y)y = alpha*AT*x + yMatrixtransSolve(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, rank1, rank2, rank2, rank2, rank2, scale, set, toString, transABmult, transABmult, transABmultAdd, transAmult, transAmult, transAmultAdd, transBmult, transBmult, transBmultAdd, transMult, transMult, transMultAdd, transpose, transpose, transRank1, transRank1, transRank2, transRank2
-
-
-
-
Constructor Detail
-
DenseMatrix
public DenseMatrix(MatrixVectorReader r) throws java.io.IOException
Constructor for DenseMatrix- Parameters:
r- Reader to get the matrix from- Throws:
java.io.IOException
-
DenseMatrix
public DenseMatrix(int numRows, int numColumns)Constructor for DenseMatrix- Parameters:
numRows- Number of rowsnumColumns- Number of columns
-
DenseMatrix
public DenseMatrix(Matrix A)
Constructor for DenseMatrix- Parameters:
A- Matrix to copy. A deep copy is made
-
DenseMatrix
public DenseMatrix(Matrix A, boolean deep)
Constructor for DenseMatrix- Parameters:
A- Matrix to copy contents fromdeep- If true,Ais copied, else a shallow copy is made and the matrices share underlying storage. For this,Amust be a dense matrix
-
DenseMatrix
public DenseMatrix(Vector x, boolean deep)
Constructor for DenseMatrix. Builds the matrix from a vector- Parameters:
x- Vector to copy from. This will form this matrix' single columndeep- If true, x is copied, if false, the internal storage of this matrix is the same as that of the vector. In that case,xmust be aDenseVector
-
DenseMatrix
public DenseMatrix(Vector x)
Constructor for DenseMatrix. Builds the matrix from a vector- Parameters:
x- The vector which forms this matrix' single column. It is copied, not referenced
-
DenseMatrix
public DenseMatrix(Vector[] x)
Constructor for DenseMatrix. Builds the matrix from vectors. Each vector will correspond to a column of the matrix- Parameters:
x- Vectors which forms the columns of this matrix. Every vector must have the same size
-
DenseMatrix
public DenseMatrix(double[][] values)
Constructor for DenseMatrix. Copies from the passed array- Parameters:
values- Arrays to copy from. Every sub-array must have the same size
-
-
Method Detail
-
copy
public DenseMatrix 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
-
transBmultAdd
public Matrix transBmultAdd(double alpha, Matrix B, Matrix C)
Description copied from interface:MatrixC = alpha*A*BT + C- Specified by:
transBmultAddin interfaceMatrix- Overrides:
transBmultAddin classAbstractMatrixB- Matrix such thatB.numRows() == A.numRows()andB.numColumns() == C.numColumns()C- Matrix such thatC.numRows() == A.numColumns()andB.numColumns() == C.numColumns()- Returns:
- C
-
transABmultAdd
public Matrix transABmultAdd(double alpha, Matrix B, Matrix C)
Description copied from interface:MatrixC = alpha*AT*BT + C- Specified by:
transABmultAddin interfaceMatrix- Overrides:
transABmultAddin classAbstractMatrixB- Matrix such thatB.numColumns() == A.numRows()andB.numRows() == C.numColumns()C- Matrix such thatC.numRows() == A.numColumns()andB.numRows() == 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
-
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
-
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
-
getData
public double[] getData()
Returns the matrix contents. Ordering depends on the underlying storage assumptions
-
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
-
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
-
-