/**
* <p>
* Puts all the real eigenvectors into the columns of a matrix. If an eigenvalue is imaginary
* then the corresponding eigenvector will have zeros in its column.
* </p>
*
* @param eig An eigenvalue decomposition which has already decomposed a matrix.
* @return An m by m matrix containing eigenvectors in its columns.
*/
public static DMatrixRMaj createMatrixV(EigenDecomposition_F64 <DMatrixRMaj> eig)
{
int N = eig.getNumberOfEigenvalues();
DMatrixRMaj V = new DMatrixRMaj(N, N);
for (int i = 0; i < N; i++)
{
Complex_F64 c = eig.getEigenvalue(i);
if (c.isReal())
{
DMatrixRMaj v = eig.getEigenVector(i);
if (v != null)
{
for (int j = 0; j < N; j++)
{
V.set(j, i, v.get(j, 0));
}
}
}
}
return(V);
}
/**
* <p>
* A diagonal matrix where real diagonal element contains a real eigenvalue. If an eigenvalue
* is imaginary then zero is stored in its place.
* </p>
*
* @param eig An eigenvalue decomposition which has already decomposed a matrix.
* @return A diagonal matrix containing the eigenvalues.
*/
public static DMatrixRMaj createMatrixD(EigenDecomposition_F64 <DMatrixRMaj> eig)
{
int N = eig.getNumberOfEigenvalues();
DMatrixRMaj D = new DMatrixRMaj(N, N);
for (int i = 0; i < N; i++)
{
Complex_F64 c = eig.getEigenvalue(i);
if (c.isReal())
{
D.set(i, i, c.real);
}
}
return(D);
}
public int getNumberOfEigenvalues()
{
return(symmetric ? symmetricAlg.getNumberOfEigenvalues() :
generalAlg.getNumberOfEigenvalues());
}
