Package no.uib.cipr.matrix
Class SymmTridiagMatrix
- java.lang.Object
-
- no.uib.cipr.matrix.AbstractMatrix
-
- no.uib.cipr.matrix.SymmTridiagMatrix
-
- All Implemented Interfaces:
java.lang.Iterable<MatrixEntry>,Matrix
- Direct Known Subclasses:
SPDTridiagMatrix
public class SymmTridiagMatrix extends AbstractMatrix
Symmetrical tridiagonal matrix. Storage as forTridiagMatrix, but only one off-diagonal array is stored.
-
-
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 SymmTridiagMatrix(double[] diag, double[] offDiag)Constructor for SymmTridiagMatrixSymmTridiagMatrix(double[] diag, double[] offDiag, int n)Constructor for SymmTridiagMatrixSymmTridiagMatrix(int n)Constructor for SymmTridiagMatrixSymmTridiagMatrix(Matrix A)Constructor for SymmTridiagMatrixSymmTridiagMatrix(Matrix A, boolean deep)Constructor for SymmTridiagMatrix
-
Method Summary
All Methods Instance Methods Concrete Methods Modifier and Type Method Description voidadd(int row, int column, double value)A(row,column) += valueSymmTridiagMatrixcopy()Creates a deep copy of the matrixdoubleget(int row, int column)ReturnsA(row,column)double[]getDiagonal()Returns the diagonal entries.double[]getOffDiagonal()Returns the off diagonal entries.java.util.Iterator<MatrixEntry>iterator()voidset(int row, int column, double value)A(row,column) = valueMatrixsolve(Matrix B, Matrix X)X = A\B.Vectorsolve(Vector b, Vector x)x = A\b.Matrixtranspose()Transposes the matrix in-place.MatrixtransSolve(Matrix B, Matrix X)X = AT\B.VectortransSolve(Vector b, Vector x)x = AT\b.SymmTridiagMatrixzero()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, mult, multAdd, multAdd, multAdd, multAdd, norm, norm1, normF, normInf, numColumns, numRows, rank1, rank1, rank1, rank1, rank1, rank1, rank2, rank2, rank2, rank2, scale, set, set, toString, transABmult, transABmult, transABmultAdd, transABmultAdd, transAmult, transAmult, transAmultAdd, transAmultAdd, transBmult, transBmult, transBmultAdd, transBmultAdd, transMult, transMult, transMultAdd, transMultAdd, transpose, transRank1, transRank1, transRank2, transRank2
-
-
-
-
Constructor Detail
-
SymmTridiagMatrix
public SymmTridiagMatrix(double[] diag, double[] offDiag, int n)Constructor for SymmTridiagMatrix- Parameters:
diag- Main diagonaloffDiag- Offdiagonals, both upper and lowern- Size of the matrix. The main diagonal must be at least as long as n, and the off diagonal part must be at least as long as n-1
-
SymmTridiagMatrix
public SymmTridiagMatrix(double[] diag, double[] offDiag)Constructor for SymmTridiagMatrix- Parameters:
diag- Main diagonaloffDiag- Offdiagonals. Must be one shorter than diag
-
SymmTridiagMatrix
public SymmTridiagMatrix(int n)
Constructor for SymmTridiagMatrix- Parameters:
n- Size of the matrix. Since the matrix must be square, this equals both the number of rows and columns.ncannot be zero
-
SymmTridiagMatrix
public SymmTridiagMatrix(Matrix A)
Constructor for SymmTridiagMatrix- Parameters:
A- Matrix to copy contents from. Only main and the superdiagonal is copied over
-
SymmTridiagMatrix
public SymmTridiagMatrix(Matrix A, boolean deep)
Constructor for SymmTridiagMatrix- Parameters:
A- Matrix to copy contents from. Only main and the superdiagonal is copied over. It must be square and cannot have any zero dimension lengthsdeep- True for a deep copy. For shallow copiesAmust be aSymmTridiagMatrix
-
-
Method Detail
-
getDiagonal
public double[] getDiagonal()
Returns the diagonal entries. Length equaln
-
getOffDiagonal
public double[] getOffDiagonal()
Returns the off diagonal entries. Length equaln-1
-
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
-
get
public double get(int row, int column)Description copied from interface:MatrixReturnsA(row,column)- Specified by:
getin interfaceMatrix- Overrides:
getin 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
-
copy
public SymmTridiagMatrix copy()
Description copied from interface:MatrixCreates a deep copy of the matrix- Specified by:
copyin interfaceMatrix- Overrides:
copyin classAbstractMatrix- Returns:
- A
-
zero
public SymmTridiagMatrix 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
-
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
-
iterator
public java.util.Iterator<MatrixEntry> iterator()
- Specified by:
iteratorin interfacejava.lang.Iterable<MatrixEntry>- Overrides:
iteratorin classAbstractMatrix
-
-