Package no.uib.cipr.matrix
Class BandMatrix
- java.lang.Object
-
- no.uib.cipr.matrix.AbstractMatrix
-
- no.uib.cipr.matrix.BandMatrix
-
- All Implemented Interfaces:
java.lang.Iterable<MatrixEntry>,Matrix
public class BandMatrix extends AbstractMatrix
Banded matrix. The banded matrix is a useful sparse structure for many kinds of direct computations, however it should only be used if the band is sufficiently narrow as wide bands actually wastes both memory and compute time. The matrixa11 a12 a21 a22 a23 a31 a32 a33 a34 a42 a43 a44 a45 a53 a54 a55 has two lower diagonals and one upper diagonal. It will be stored in the array
a11 a21 a31 a21 a22 a32 a42 a23 a33 a43 a53 a34 a44 a54 a45 a55 Empty cells are allocated, but never referenced.
-
-
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 BandMatrix(int n, int kl, int ku)Constructor for BandMatrixBandMatrix(Matrix A, int kl, int ku)Constructor for BandMatrixBandMatrix(Matrix A, int kl, int ku, boolean deep)Constructor for BandMatrix
-
Method Summary
All Methods Instance Methods Concrete Methods Modifier and Type Method Description voidadd(int row, int column, double value)A(row,column) += valueBandMatrixcopy()Creates a deep copy of the matrixdoubleget(int row, int column)ReturnsA(row,column)double[]getData()Returns the matrix contentsjava.util.Iterator<MatrixEntry>iterator()VectormultAdd(double alpha, Vector x, Vector y)y = alpha*A*x + yintnumSubDiagonals()Returns the number of lower diagonalsintnumSuperDiagonals()Returns the number of upper diagonalsvoidset(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.VectortransMultAdd(double alpha, Vector x, Vector y)y = alpha*AT*x + yMatrixtranspose()Transposes the matrix in-place.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, max, max, mult, mult, mult, mult, multAdd, multAdd, multAdd, norm, norm1, normF, normInf, numColumns, numRows, rank1, rank1, rank1, rank1, rank1, rank1, rank2, 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, transRank2, transSolve, transSolve
-
-
-
-
Constructor Detail
-
BandMatrix
public BandMatrix(int n, int kl, int ku)Constructor for BandMatrix- Parameters:
n- Size of the matrix. Since the matrix must be square, this equals both the number of rows and columnskl- Number of bands above the main diagonal (superdiagonals)ku- Number of bands below the main diagonal (subdiagonals)
-
BandMatrix
public BandMatrix(Matrix A, int kl, int ku)
Constructor for BandMatrix- Parameters:
A- Matrix to copy contents from. Only the parts ofAthat lie within the allocated band are copied over, the rest is ignoredkl- Number of bands above the main diagonal (superdiagonals)ku- Number of bands below the main diagonal (subdiagonals)
-
BandMatrix
public BandMatrix(Matrix A, int kl, int ku, boolean deep)
Constructor for BandMatrix- Parameters:
A- Matrix to copy contents from. Only the parts ofAthat lie within the allocated band are copied over, the rest is ignoredkl- Number of bands above the main diagonal (superdiagonals)ku- Number of bands below the main diagonal (subdiagonals)deep- True for a deep copy. For shallow copies,Amust be a banded matrix
-
-
Method Detail
-
copy
public BandMatrix copy()
Description copied from interface:MatrixCreates a deep copy of the matrix- Specified by:
copyin interfaceMatrix- Overrides:
copyin classAbstractMatrix- Returns:
- A
-
zero
public Matrix zero()
Description copied from interface:MatrixZeros all the entries in the matrix, while preserving any underlying structure. Useful for general, unstructured matrices.
-
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
-
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
-
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
-
numSubDiagonals
public int numSubDiagonals()
Returns the number of lower diagonals
-
numSuperDiagonals
public int numSuperDiagonals()
Returns the number of upper diagonals
-
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
-
iterator
public java.util.Iterator<MatrixEntry> iterator()
- Specified by:
iteratorin interfacejava.lang.Iterable<MatrixEntry>- Overrides:
iteratorin classAbstractMatrix
-
-