Lines Matching refs:m_matT
113 EigenSolver() : m_eivec(), m_eivalues(), m_isInitialized(false), m_realSchur(), m_matT(), m_tmp() {}
127 m_matT(size, size),
152 m_matT(matrix.rows(), matrix.cols()),
306 MatrixType m_matT;
376 m_matT = m_realSchur.matrixT();
385 if (i == matrix.cols() - 1 || m_matT.coeff(i+1, i) == Scalar(0))
387 m_eivalues.coeffRef(i) = m_matT.coeff(i, i);
392 Scalar p = Scalar(0.5) * (m_matT.coeff(i, i) - m_matT.coeff(i+1, i+1));
393 Scalar z = sqrt(abs(p * p + m_matT.coeff(i+1, i) * m_matT.coeff(i, i+1)));
394 m_eivalues.coeffRef(i) = ComplexScalar(m_matT.coeff(i+1, i+1) + p, z);
395 m_eivalues.coeffRef(i+1) = ComplexScalar(m_matT.coeff(i+1, i+1) + p, -z);
443 norm += m_matT.row(j).segment((std::max)(j-1,Index(0)), size-(std::max)(j-1,Index(0))).cwiseAbs().sum();
463 m_matT.coeffRef(n,n) = 1.0;
466 Scalar w = m_matT.coeff(i,i) - p;
467 Scalar r = m_matT.row(i).segment(l,n-l+1).dot(m_matT.col(n).segment(l, n-l+1));
480 m_matT.coeffRef(i,n) = -r / w;
482 m_matT.coeffRef(i,n) = -r / (eps * norm);
486 Scalar x = m_matT.coeff(i,i+1);
487 Scalar y = m_matT.coeff(i+1,i);
490 m_matT.coeffRef(i,n) = t;
492 m_matT.coeffRef(i+1,n) = (-r - w * t) / x;
494 m_matT.coeffRef(i+1,n) = (-lastr - y * t) / lastw;
498 Scalar t = abs(m_matT.coeff(i,n));
500 m_matT.col(n).tail(size-i) /= t;
510 if (abs(m_matT.coeff(n,n-1)) > abs(m_matT.coeff(n-1,n)))
512 m_matT.coeffRef(n-1,n-1) = q / m_matT.coeff(n,n-1);
513 m_matT.coeffRef(n-1,n) = -(m_matT.coeff(n,n) - p) / m_matT.coeff(n,n-1);
517 std::complex<Scalar> cc = cdiv<Scalar>(0.0,-m_matT.coeff(n-1,n),m_matT.coeff(n-1,n-1)-p,q);
518 m_matT.coeffRef(n-1,n-1) = numext::real(cc);
519 m_matT.coeffRef(n-1,n) = numext::imag(cc);
521 m_matT.coeffRef(n,n-1) = 0.0;
522 m_matT.coeffRef(n,n) = 1.0;
525 Scalar ra = m_matT.row(i).segment(l, n-l+1).dot(m_matT.col(n-1).segment(l, n-l+1));
526 Scalar sa = m_matT.row(i).segment(l, n-l+1).dot(m_matT.col(n).segment(l, n-l+1));
527 Scalar w = m_matT.coeff(i,i) - p;
541 m_matT.coeffRef(i,n-1) = numext::real(cc);
542 m_matT.coeffRef(i,n) = numext::imag(cc);
547 Scalar x = m_matT.coeff(i,i+1);
548 Scalar y = m_matT.coeff(i+1,i);
555 m_matT.coeffRef(i,n-1) = numext::real(cc);
556 m_matT.coeffRef(i,n) = numext::imag(cc);
559 m_matT.coeffRef(i+1,n-1) = (-ra - w * m_matT.coeff(i,n-1) + q * m_matT.coeff(i,n)) / x;
560 m_matT.coeffRef(i+1,n) = (-sa - w * m_matT.coeff(i,n) - q * m_matT.coeff(i,n-1)) / x;
564 cc = cdiv(-lastra-y*m_matT.coeff(i,n-1),-lastsa-y*m_matT.coeff(i,n),lastw,q);
565 m_matT.coeffRef(i+1,n-1) = numext::real(cc);
566 m_matT.coeffRef(i+1,n) = numext::imag(cc);
572 Scalar t = (max)(abs(m_matT.coeff(i,n-1)),abs(m_matT.coeff(i,n)));
574 m_matT.block(i, n-1, size-i, 2) /= t;
591 m_tmp.noalias() = m_eivec.leftCols(j+1) * m_matT.col(j).segment(0, j+1);
