JSci.maths.matrices.ComplexSquareMatrix Maven / Gradle / Ivy
Go to download
Show more of this group Show more artifacts with this name
Show all versions of jsci Show documentation
Show all versions of jsci Show documentation
JSci is a set of open source Java packages. The aim is to encapsulate scientific methods/principles in the most natural way possible. As such they should greatly aid the development of scientific based software.
It offers: abstract math interfaces, linear algebra (support for various matrix and vector types), statistics (including probability distributions), wavelets, newtonian mechanics, chart/graph components (AWT and Swing), MathML DOM implementation, ...
Note: some packages, like javax.comm, for the astro and instruments package aren't listed as dependencies (not available).
The newest version!
/* AUTO-GENERATED */
package JSci.maths.matrices;
import JSci.maths.ArrayMath;
import JSci.maths.Complex;
import JSci.maths.ComplexMapping;
import JSci.maths.LinearMath;
import JSci.maths.DimensionException;
import JSci.maths.MaximumIterationsExceededException;
import JSci.maths.vectors.AbstractComplexVector;
import JSci.maths.vectors.ComplexVector;
import JSci.maths.groups.AbelianGroup;
/**
* The ComplexSquareMatrix class provides an object for encapsulating square matrices containing complex numbers.
* @version 2.2
* @author Mark Hale
*/
public class ComplexSquareMatrix extends AbstractComplexSquareMatrix {
/**
* Arrays containing the elements of the matrix.
*/
protected final double matrixRe[][],matrixIm[][];
/**
* Constructs a matrix by wrapping two arrays.
* @param arrayRe an array of real values
* @param arrayIm an array of imaginary values
* @exception MatrixDimensionException If the array is not square.
*/
public ComplexSquareMatrix(final double arrayRe[][],final double arrayIm[][]) {
super(arrayRe.length);
if(!ArrayMath.isSquare(arrayRe))
throw new MatrixDimensionException("Array is not square.");
if(!ArrayMath.isSquare(arrayIm))
throw new MatrixDimensionException("Array is not square.");
matrixRe=arrayRe;
matrixIm=arrayIm;
}
/**
* Constructs an empty matrix.
* @param size the number of rows/columns
*/
public ComplexSquareMatrix(final int size) {
this(new double[size][size], new double[size][size]);
}
/**
* Constructs a matrix from an array.
* @param array an assigned value
* @exception MatrixDimensionException If the array is not square.
*/
public ComplexSquareMatrix(final Complex array[][]) {
this(array.length);
for(int i=0;i=0 && i=0 && j=0 && i=0 && j=0 && i=0 && j=0 && i=0 && j=0 && i=0 && j-norm.
* @author Taber Smith
*/
public double infNorm() {
double result=0.0,tmpResult;
for(int i=0;iresult)
result=tmpResult;
}
return result;
}
/**
* Returns the Frobenius or Hilbert-Schmidt (l2) norm.
* @jsci.planetmath FrobeniusMatrixNorm
* @author Taber Smith
*/
public double frobeniusNorm() {
double result=0.0;
for(int j,i=0;imax) {
max=tmp;
pivotrow=i;
}
}
// swap row j with pivotrow
if(pivotrow!=j) {
double[] tmprow = arrayURe[j];
arrayURe[j] = arrayURe[pivotrow];
arrayURe[pivotrow] = tmprow;
tmprow = arrayUIm[j];
arrayUIm[j] = arrayUIm[pivotrow];
arrayUIm[pivotrow] = tmprow;
int k=pivot[j];
pivot[j]=pivot[pivotrow];
pivot[pivotrow]=k;
// update parity
pivot[N]=-pivot[N];
}
// divide by pivot
double tmpRe=arrayURe[j][j];
double tmpIm=arrayUIm[j][j];
if(Math.abs(tmpRe)