Lines Matching refs:RealScalar
94 typedef typename NumTraits<Scalar>::Real RealScalar; typedef
103 typedef typename internal::plain_col_type<MatrixType, RealScalar>::type RealVectorType;
391 template<int StorageOrder,typename RealScalar, typename Scalar, typename Index>
393 static void tridiagonal_qr_step(RealScalar* diag, RealScalar* subdiag, Index start, Index end, Scal…
433 RealScalar scale = mat.cwiseAbs().maxCoeff();
434 if(scale==RealScalar(0)) scale = RealScalar(1);
494 typedef typename DiagType::RealScalar RealScalar;
495 const RealScalar considerAsZero = (std::numeric_limits<RealScalar>::min)();
496 const RealScalar precision = RealScalar(2)*NumTraits<RealScalar>::epsilon();
505 while (end>0 && subdiag[end-1]==RealScalar(0))
806 template<int StorageOrder,typename RealScalar, typename Scalar, typename Index>
808 static void tridiagonal_qr_step(RealScalar* diag, RealScalar* subdiag, Index start, Index end, Scal…
811 RealScalar td = (diag[end-1] - diag[end])*RealScalar(0.5);
812 RealScalar e = subdiag[end-1];
818 RealScalar mu = diag[end];
819 if(td==RealScalar(0))
823 RealScalar e2 = numext::abs2(subdiag[end-1]);
824 RealScalar h = numext::hypot(td,e);
825 …if(e2==RealScalar(0)) mu -= (e / (td + (td>RealScalar(0) ? RealScalar(1) : RealScalar(-1)))) * (e …
826 else mu -= e2 / (td + (td>RealScalar(0) ? h : -h));
829 RealScalar x = diag[start] - mu;
830 RealScalar z = subdiag[start];
833 JacobiRotation<RealScalar> rot;
837 RealScalar sdk = rot.s() * diag[k] + rot.c() * subdiag[k];
838 RealScalar dkp1 = rot.s() * subdiag[k] + rot.c() * diag[k+1];