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()
152 m_matT(matrix.rows(), matrix.cols()),
293 MatrixType m_matT; variable
359 m_matT = m_realSchur.matrixT(); in compute()
368 if (i == matrix.cols() - 1 || m_matT.coeff(i+1, i) == Scalar(0)) in compute()
370 m_eivalues.coeffRef(i) = m_matT.coeff(i, i); in compute()
375 Scalar p = Scalar(0.5) * (m_matT.coeff(i, i) - m_matT.coeff(i+1, i+1)); in compute()
376 … Scalar z = internal::sqrt(internal::abs(p * p + m_matT.coeff(i+1, i) * m_matT.coeff(i, i+1))); in compute()
377 m_eivalues.coeffRef(i) = ComplexScalar(m_matT.coeff(i+1, i+1) + p, z); in compute()
378 m_eivalues.coeffRef(i+1) = ComplexScalar(m_matT.coeff(i+1, i+1) + p, -z); in compute()
424 …norm += m_matT.row(j).segment((std::max)(j-1,Index(0)), size-(std::max)(j-1,Index(0))).cwiseAbs().… in doComputeEigenvectors()
444 m_matT.coeffRef(n,n) = 1.0; in doComputeEigenvectors()
447 Scalar w = m_matT.coeff(i,i) - p; in doComputeEigenvectors()
448 Scalar r = m_matT.row(i).segment(l,n-l+1).dot(m_matT.col(n).segment(l, n-l+1)); in doComputeEigenvectors()
461 m_matT.coeffRef(i,n) = -r / w; in doComputeEigenvectors()
463 m_matT.coeffRef(i,n) = -r / (eps * norm); in doComputeEigenvectors()
467 Scalar x = m_matT.coeff(i,i+1); in doComputeEigenvectors()
468 Scalar y = m_matT.coeff(i+1,i); in doComputeEigenvectors()
471 m_matT.coeffRef(i,n) = t; in doComputeEigenvectors()
473 m_matT.coeffRef(i+1,n) = (-r - w * t) / x; in doComputeEigenvectors()
475 m_matT.coeffRef(i+1,n) = (-lastr - y * t) / lastw; in doComputeEigenvectors()
479 Scalar t = internal::abs(m_matT.coeff(i,n)); in doComputeEigenvectors()
481 m_matT.col(n).tail(size-i) /= t; in doComputeEigenvectors()
491 if (internal::abs(m_matT.coeff(n,n-1)) > internal::abs(m_matT.coeff(n-1,n))) in doComputeEigenvectors()
493 m_matT.coeffRef(n-1,n-1) = q / m_matT.coeff(n,n-1); in doComputeEigenvectors()
494 m_matT.coeffRef(n-1,n) = -(m_matT.coeff(n,n) - p) / m_matT.coeff(n,n-1); in doComputeEigenvectors()
498 std::complex<Scalar> cc = cdiv<Scalar>(0.0,-m_matT.coeff(n-1,n),m_matT.coeff(n-1,n-1)-p,q); in doComputeEigenvectors()
499 m_matT.coeffRef(n-1,n-1) = internal::real(cc); in doComputeEigenvectors()
500 m_matT.coeffRef(n-1,n) = internal::imag(cc); in doComputeEigenvectors()
502 m_matT.coeffRef(n,n-1) = 0.0; in doComputeEigenvectors()
503 m_matT.coeffRef(n,n) = 1.0; in doComputeEigenvectors()
506 Scalar ra = m_matT.row(i).segment(l, n-l+1).dot(m_matT.col(n-1).segment(l, n-l+1)); in doComputeEigenvectors()
507 Scalar sa = m_matT.row(i).segment(l, n-l+1).dot(m_matT.col(n).segment(l, n-l+1)); in doComputeEigenvectors()
508 Scalar w = m_matT.coeff(i,i) - p; in doComputeEigenvectors()
522 m_matT.coeffRef(i,n-1) = internal::real(cc); in doComputeEigenvectors()
523 m_matT.coeffRef(i,n) = internal::imag(cc); in doComputeEigenvectors()
528 Scalar x = m_matT.coeff(i,i+1); in doComputeEigenvectors()
529 Scalar y = m_matT.coeff(i+1,i); in doComputeEigenvectors()
536 m_matT.coeffRef(i,n-1) = internal::real(cc); in doComputeEigenvectors()
537 m_matT.coeffRef(i,n) = internal::imag(cc); in doComputeEigenvectors()
540 … m_matT.coeffRef(i+1,n-1) = (-ra - w * m_matT.coeff(i,n-1) + q * m_matT.coeff(i,n)) / x; in doComputeEigenvectors()
541 m_matT.coeffRef(i+1,n) = (-sa - w * m_matT.coeff(i,n) - q * m_matT.coeff(i,n-1)) / x; in doComputeEigenvectors()
545 cc = cdiv(-lastra-y*m_matT.coeff(i,n-1),-lastsa-y*m_matT.coeff(i,n),lastw,q); in doComputeEigenvectors()
546 m_matT.coeffRef(i+1,n-1) = internal::real(cc); in doComputeEigenvectors()
547 m_matT.coeffRef(i+1,n) = internal::imag(cc); in doComputeEigenvectors()
553 Scalar t = (max)(internal::abs(m_matT.coeff(i,n-1)),internal::abs(m_matT.coeff(i,n))); in doComputeEigenvectors()
555 m_matT.block(i, n-1, size-i, 2) /= t; in doComputeEigenvectors()
572 m_tmp.noalias() = m_eivec.leftCols(j+1) * m_matT.col(j).segment(0, j+1); in doComputeEigenvectors()