Lines Matching refs:subdiag
24 ComputationInfo computeFromTridiagonal_impl(DiagType& diag, SubDiagType& subdiag, const Index maxIt…
238 …omputeFromTridiagonal(const RealVectorType& diag, const SubDiagonalType& subdiag , int options=Com…
393 static void tridiagonal_qr_step(RealScalar* diag, RealScalar* subdiag, Index start, Index end, Scal…
451 ::computeFromTridiagonal(const RealVectorType& diag, const SubDiagonalType& subdiag , int options)
457 m_subdiag = subdiag;
482 ComputationInfo computeFromTridiagonal_impl(DiagType& diag, SubDiagType& subdiag, const Index maxIt…
501 …if (internal::isMuchSmallerThan(abs(subdiag[i]),(abs(diag[i])+abs(diag[i+1])),precision) || abs(su…
502 subdiag[i] = 0;
505 while (end>0 && subdiag[end-1]==RealScalar(0))
517 while (start>0 && subdiag[start-1]!=0)
520 …r_step<MatrixType::Flags&RowMajorBit ? RowMajor : ColMajor>(diag.data(), subdiag.data(), start, en…
808 static void tridiagonal_qr_step(RealScalar* diag, RealScalar* subdiag, Index start, Index end, Scal…
812 RealScalar e = subdiag[end-1];
823 RealScalar e2 = numext::abs2(subdiag[end-1]);
830 RealScalar z = subdiag[start];
837 RealScalar sdk = rot.s() * diag[k] + rot.c() * subdiag[k];
838 RealScalar dkp1 = rot.s() * subdiag[k] + rot.c() * diag[k+1];
840 …diag[k] = rot.c() * (rot.c() * diag[k] - rot.s() * subdiag[k]) - rot.s() * (rot.c() * subdiag[k] -…
842 subdiag[k] = rot.c() * sdk - rot.s() * dkp1;
846 subdiag[k - 1] = rot.c() * subdiag[k-1] - rot.s() * z;
848 x = subdiag[k];
852 z = -rot.s() * subdiag[k+1];
853 subdiag[k + 1] = rot.c() * subdiag[k+1];