//Sort Eigenvalues, Eigenvectors
public static void SortEigen(ref Vector eigenvalues, ref MMatrix eigenvectors)
{
if (eigenvalues.Elements.Length != eigenvectors.col)
throw new MMatrixException("Eigenvectors do not match eigenvalues.");
double temp;
for (int i = 0; i < eigenvalues.Elements.Length; ++i)
for (int j = i + 1; j < eigenvalues.Elements.Length; ++j)
if (eigenvalues[i] < eigenvalues[j])
{
temp = eigenvalues[i];
eigenvalues[i] = eigenvalues[j];
eigenvalues[j] = temp;
eigenvectors.InterchangeColumn(i, j);
}
}
