}
}
if (m != j) {
if (its == maxIter)
throw new InvalidMatrixException(
new MaxIterationsExceededException(maxIter));
its++;
double q = (realEigenvalues[j + 1] - realEigenvalues[j]) / (2 * e[j]);
double t = FastMath.sqrt(1 + q * q);
if (q < 0.0) {
q = realEigenvalues[m] - realEigenvalues[j] + e[j] / (q - t);