Lines Matching refs:m_matT
84 : m_matT(size, size), in m_matT() function
105 : m_matT(matrix.rows(),matrix.cols()),
147 return m_matT; in matrixT()
227 MatrixType m_matT;
261 m_matT.setZero(matrix.rows(),matrix.cols()); in compute()
276 m_matT *= scale; in compute()
286 m_matT = matrixH; in computeFromHessenberg()
293 m_workspaceVector.resize(m_matT.cols()); in computeFromHessenberg()
300 Index iu = m_matT.cols() - 1; in computeFromHessenberg()
315 m_matT.coeffRef(iu,iu) = m_matT.coeff(iu,iu) + exshift; in computeFromHessenberg()
317 m_matT.coeffRef(iu, iu-1) = Scalar(0); in computeFromHessenberg()
355 const Index size = m_matT.cols(); in computeNormOfT()
361 norm += m_matT.col(j).segment(0, (std::min)(size,j+2)).cwiseAbs().sum(); in computeNormOfT()
373 Scalar s = abs(m_matT.coeff(res-1,res-1)) + abs(m_matT.coeff(res,res)); in findSmallSubdiagEntry()
374 if (abs(m_matT.coeff(res,res-1)) <= NumTraits<Scalar>::epsilon() * s) in findSmallSubdiagEntry()
387 const Index size = m_matT.cols(); in splitOffTwoRows()
391 Scalar p = Scalar(0.5) * (m_matT.coeff(iu-1,iu-1) - m_matT.coeff(iu,iu)); in splitOffTwoRows()
392 …Scalar q = p * p + m_matT.coeff(iu,iu-1) * m_matT.coeff(iu-1,iu); // q = tr^2 / 4 - det = discr/4 in splitOffTwoRows()
393 m_matT.coeffRef(iu,iu) += exshift; in splitOffTwoRows()
394 m_matT.coeffRef(iu-1,iu-1) += exshift; in splitOffTwoRows()
401 rot.makeGivens(p + z, m_matT.coeff(iu, iu-1)); in splitOffTwoRows()
403 rot.makeGivens(p - z, m_matT.coeff(iu, iu-1)); in splitOffTwoRows()
405 m_matT.rightCols(size-iu+1).applyOnTheLeft(iu-1, iu, rot.adjoint()); in splitOffTwoRows()
406 m_matT.topRows(iu+1).applyOnTheRight(iu-1, iu, rot); in splitOffTwoRows()
407 m_matT.coeffRef(iu, iu-1) = Scalar(0); in splitOffTwoRows()
413 m_matT.coeffRef(iu-1, iu-2) = Scalar(0); in splitOffTwoRows()
422 shiftInfo.coeffRef(0) = m_matT.coeff(iu,iu); in computeShift()
423 shiftInfo.coeffRef(1) = m_matT.coeff(iu-1,iu-1); in computeShift()
424 shiftInfo.coeffRef(2) = m_matT.coeff(iu,iu-1) * m_matT.coeff(iu-1,iu); in computeShift()
431 m_matT.coeffRef(i,i) -= shiftInfo.coeff(0); in computeShift()
432 Scalar s = abs(m_matT.coeff(iu,iu-1)) + abs(m_matT.coeff(iu-1,iu-2)); in computeShift()
452 m_matT.coeffRef(i,i) -= s; in computeShift()
467 const Scalar Tmm = m_matT.coeff(im,im); in initFrancisQRStep()
470 v.coeffRef(0) = (r * s - shiftInfo.coeff(2)) / m_matT.coeff(im+1,im) + m_matT.coeff(im,im+1); in initFrancisQRStep()
471 v.coeffRef(1) = m_matT.coeff(im+1,im+1) - Tmm - r - s; in initFrancisQRStep()
472 v.coeffRef(2) = m_matT.coeff(im+2,im+1); in initFrancisQRStep()
476 const Scalar lhs = m_matT.coeff(im,im-1) * (abs(v.coeff(1)) + abs(v.coeff(2))); in initFrancisQRStep()
477 …const Scalar rhs = v.coeff(0) * (abs(m_matT.coeff(im-1,im-1)) + abs(Tmm) + abs(m_matT.coeff(im+1,i… in initFrancisQRStep()
490 const Index size = m_matT.cols(); in performFrancisQRStep()
500 v = m_matT.template block<3,1>(k,k-1); in performFrancisQRStep()
509 m_matT.coeffRef(k,k-1) = -m_matT.coeff(k,k-1); in performFrancisQRStep()
511 m_matT.coeffRef(k,k-1) = beta; in performFrancisQRStep()
514 m_matT.block(k, k, 3, size-k).applyHouseholderOnTheLeft(ess, tau, workspace); in performFrancisQRStep()
515 m_matT.block(0, k, (std::min)(iu,k+3) + 1, 3).applyHouseholderOnTheRight(ess, tau, workspace); in performFrancisQRStep()
521 Matrix<Scalar, 2, 1> v = m_matT.template block<2,1>(iu-1, iu-2); in performFrancisQRStep()
528 m_matT.coeffRef(iu-1, iu-2) = beta; in performFrancisQRStep()
529 m_matT.block(iu-1, iu-1, 2, size-iu+1).applyHouseholderOnTheLeft(ess, tau, workspace); in performFrancisQRStep()
530 m_matT.block(0, iu-1, iu+1, 2).applyHouseholderOnTheRight(ess, tau, workspace); in performFrancisQRStep()
538 m_matT.coeffRef(i,i-2) = Scalar(0); in performFrancisQRStep()
540 m_matT.coeffRef(i,i-3) = Scalar(0); in performFrancisQRStep()