Lines Matching refs:m_matT
113 …EigenSolver() : m_eivec(), m_eivalues(), m_isInitialized(false), m_realSchur(), m_matT(), m_tmp() … in EigenSolver()
127 m_matT(size, size), in EigenSolver()
153 m_matT(matrix.rows(), matrix.cols()),
317 MatrixType m_matT; variable
395 m_matT = m_realSchur.matrixT(); in compute()
404 if (i == matrix.cols() - 1 || m_matT.coeff(i+1, i) == Scalar(0)) in compute()
406 m_eivalues.coeffRef(i) = m_matT.coeff(i, i); in compute()
418 Scalar p = Scalar(0.5) * (m_matT.coeff(i, i) - m_matT.coeff(i+1, i+1)); in compute()
423 Scalar t0 = m_matT.coeff(i+1, i); in compute()
424 Scalar t1 = m_matT.coeff(i, i+1); in compute()
432 m_eivalues.coeffRef(i) = ComplexScalar(m_matT.coeff(i+1, i+1) + p, z); in compute()
433 m_eivalues.coeffRef(i+1) = ComplexScalar(m_matT.coeff(i+1, i+1) + p, -z); in compute()
468 …norm += m_matT.row(j).segment((std::max)(j-1,Index(0)), size-(std::max)(j-1,Index(0))).cwiseAbs().… in doComputeEigenvectors()
488 m_matT.coeffRef(n,n) = Scalar(1); in doComputeEigenvectors()
491 Scalar w = m_matT.coeff(i,i) - p; in doComputeEigenvectors()
492 Scalar r = m_matT.row(i).segment(l,n-l+1).dot(m_matT.col(n).segment(l, n-l+1)); in doComputeEigenvectors()
505 m_matT.coeffRef(i,n) = -r / w; in doComputeEigenvectors()
507 m_matT.coeffRef(i,n) = -r / (eps * norm); in doComputeEigenvectors()
511 Scalar x = m_matT.coeff(i,i+1); in doComputeEigenvectors()
512 Scalar y = m_matT.coeff(i+1,i); in doComputeEigenvectors()
515 m_matT.coeffRef(i,n) = t; in doComputeEigenvectors()
517 m_matT.coeffRef(i+1,n) = (-r - w * t) / x; in doComputeEigenvectors()
519 m_matT.coeffRef(i+1,n) = (-lastr - y * t) / lastw; in doComputeEigenvectors()
523 Scalar t = abs(m_matT.coeff(i,n)); in doComputeEigenvectors()
525 m_matT.col(n).tail(size-i) /= t; in doComputeEigenvectors()
535 if (abs(m_matT.coeff(n,n-1)) > abs(m_matT.coeff(n-1,n))) in doComputeEigenvectors()
537 m_matT.coeffRef(n-1,n-1) = q / m_matT.coeff(n,n-1); in doComputeEigenvectors()
538 m_matT.coeffRef(n-1,n) = -(m_matT.coeff(n,n) - p) / m_matT.coeff(n,n-1); in doComputeEigenvectors()
542 …ComplexScalar cc = ComplexScalar(Scalar(0),-m_matT.coeff(n-1,n)) / ComplexScalar(m_matT.coeff(n-1,… in doComputeEigenvectors()
543 m_matT.coeffRef(n-1,n-1) = numext::real(cc); in doComputeEigenvectors()
544 m_matT.coeffRef(n-1,n) = numext::imag(cc); in doComputeEigenvectors()
546 m_matT.coeffRef(n,n-1) = Scalar(0); in doComputeEigenvectors()
547 m_matT.coeffRef(n,n) = Scalar(1); in doComputeEigenvectors()
550 Scalar ra = m_matT.row(i).segment(l, n-l+1).dot(m_matT.col(n-1).segment(l, n-l+1)); in doComputeEigenvectors()
551 Scalar sa = m_matT.row(i).segment(l, n-l+1).dot(m_matT.col(n).segment(l, n-l+1)); in doComputeEigenvectors()
552 Scalar w = m_matT.coeff(i,i) - p; in doComputeEigenvectors()
566 m_matT.coeffRef(i,n-1) = numext::real(cc); in doComputeEigenvectors()
567 m_matT.coeffRef(i,n) = numext::imag(cc); in doComputeEigenvectors()
572 Scalar x = m_matT.coeff(i,i+1); in doComputeEigenvectors()
573 Scalar y = m_matT.coeff(i+1,i); in doComputeEigenvectors()
580 m_matT.coeffRef(i,n-1) = numext::real(cc); in doComputeEigenvectors()
581 m_matT.coeffRef(i,n) = numext::imag(cc); in doComputeEigenvectors()
584 … m_matT.coeffRef(i+1,n-1) = (-ra - w * m_matT.coeff(i,n-1) + q * m_matT.coeff(i,n)) / x; in doComputeEigenvectors()
585 m_matT.coeffRef(i+1,n) = (-sa - w * m_matT.coeff(i,n) - q * m_matT.coeff(i,n-1)) / x; in doComputeEigenvectors()
589 …cc = ComplexScalar(-lastra-y*m_matT.coeff(i,n-1),-lastsa-y*m_matT.coeff(i,n)) / ComplexScalar(last… in doComputeEigenvectors()
590 m_matT.coeffRef(i+1,n-1) = numext::real(cc); in doComputeEigenvectors()
591 m_matT.coeffRef(i+1,n) = numext::imag(cc); in doComputeEigenvectors()
596 Scalar t = numext::maxi<Scalar>(abs(m_matT.coeff(i,n-1)),abs(m_matT.coeff(i,n))); in doComputeEigenvectors()
598 m_matT.block(i, n-1, size-i, 2) /= t; in doComputeEigenvectors()
615 m_tmp.noalias() = m_eivec.leftCols(j+1) * m_matT.col(j).segment(0, j+1); in doComputeEigenvectors()