/external/eigen/Eigen/src/Core/ |
D | StableNorm.h | 57 typedef typename Derived::RealScalar RealScalar; in blueNorm_impl() typedef 63 static RealScalar b1, b2, s1m, s2m, rbig, relerr; in blueNorm_impl() 67 RealScalar eps; in blueNorm_impl() 76 … ibeta = std::numeric_limits<RealScalar>::radix; // base for floating-point numbers in blueNorm_impl() 77 …it = std::numeric_limits<RealScalar>::digits; // number of base-beta digits in m… in blueNorm_impl() 78 iemin = std::numeric_limits<RealScalar>::min_exponent; // minimum exponent in blueNorm_impl() 79 iemax = std::numeric_limits<RealScalar>::max_exponent; // maximum exponent in blueNorm_impl() 80 rbig = (std::numeric_limits<RealScalar>::max)(); // largest floating-point number in blueNorm_impl() 83 b1 = RealScalar(pow(RealScalar(ibeta),RealScalar(iexp))); // lower boundary of midrange in blueNorm_impl() 85 b2 = RealScalar(pow(RealScalar(ibeta),RealScalar(iexp))); // upper boundary of midrange in blueNorm_impl() [all …]
|
D | ConditionEstimator.h | 21 return (v_abs.array() == static_cast<typename Vector::RealScalar>(0)) in run() 30 return (v.array() < static_cast<typename Vector::RealScalar>(0)) 56 typename Decomposition::RealScalar rcond_invmatrix_L1_norm_estimate(const Decomposition& dec) 60 typedef typename Decomposition::RealScalar RealScalar; 62 typedef typename internal::plain_col_type<MatrixType, RealScalar>::type RealVector; 84 RealScalar lower_bound = v.template lpNorm<1>(); 91 RealScalar old_lower_bound = lower_bound; 133 Scalar alternating_sign(RealScalar(1)); 136 …v[i] = alternating_sign * static_cast<RealScalar>(RealScalar(1) + (RealScalar(i) / (RealScalar(n -… 140 const RealScalar alternate_lower_bound = (2 * v.template lpNorm<1>()) / (3 * RealScalar(n)); [all …]
|
D | Dot.h | 125 RealScalar z = n.squaredNorm(); 127 if(z>RealScalar(0)) 144 RealScalar z = squaredNorm(); 146 if(z>RealScalar(0)) 168 RealScalar w = n.cwiseAbs().maxCoeff(); 169 RealScalar z = (n/w).squaredNorm(); 170 if(z>RealScalar(0)) 190 RealScalar w = cwiseAbs().maxCoeff(); 191 RealScalar z = (derived()/w).squaredNorm(); 192 if(z>RealScalar(0)) [all …]
|
D | Fuzzy.h | 23 static bool run(const Derived& x, const OtherDerived& y, const typename Derived::RealScalar& prec) in run() 35 static bool run(const Derived& x, const OtherDerived& y, const typename Derived::RealScalar&) 45 static bool run(const Derived& x, const OtherDerived& y, const typename Derived::RealScalar& prec) 55 static bool run(const Derived& x, const OtherDerived&, const typename Derived::RealScalar&) 65 … bool run(const Derived& x, const typename Derived::RealScalar& y, const typename Derived::RealSca… 75 …c bool run(const Derived& x, const typename Derived::RealScalar&, const typename Derived::RealScal… 105 const RealScalar& prec 127 const RealScalar& prec 147 const RealScalar& prec
|
/external/eigen/unsupported/Eigen/CXX11/src/Tensor/ |
D | TensorFFT.h | 64 typedef typename NumTraits<typename XprTraits::Scalar>::Real RealScalar; 65 typedef typename std::complex<RealScalar> ComplexScalar; 67 …typedef typename conditional<FFTResultType == RealPart || FFTResultType == ImagPart, RealScalar, C… 92 typedef typename Eigen::NumTraits<Scalar>::Real RealScalar; 93 typedef typename std::complex<RealScalar> ComplexScalar; 94 …l::conditional<FFTResultType == RealPart || FFTResultType == ImagPart, RealScalar, ComplexScalar>:… 124 typedef typename Eigen::NumTraits<Scalar>::Real RealScalar; 125 typedef typename std::complex<RealScalar> ComplexScalar; 129 …l::conditional<FFTResultType == RealPart || FFTResultType == ImagPart, RealScalar, ComplexScalar>:… 212 buf[i] = MakeComplex<internal::is_same<InputScalar, RealScalar>::value>()(m_impl.coeff(i)); [all …]
|
/external/eigen/unsupported/Eigen/src/MatrixFunctions/ |
D | MatrixPower.h | 42 typedef typename MatrixType::RealScalar RealScalar; typedef 51 … MatrixPowerParenthesesReturnValue(MatrixPower<MatrixType>& pow, RealScalar p) : m_pow(pow), m_p(p) in MatrixPowerParenthesesReturnValue() 68 const RealScalar m_p; 95 typedef typename MatrixType::RealScalar RealScalar; typedef 96 typedef std::complex<RealScalar> ComplexScalar; 101 RealScalar m_p; 104 void compute2x2(ResultType& res, RealScalar p) const; 109 static ComplexScalar computeSuperDiag(const ComplexScalar&, const ComplexScalar&, RealScalar p); 110 static RealScalar computeSuperDiag(RealScalar, RealScalar, RealScalar p); 124 MatrixPowerAtomic(const MatrixType& T, RealScalar p); [all …]
|
D | MatrixLogarithm.h | 27 typedef typename NumTraits<Scalar>::Real RealScalar; typedef 28 static const int value = std::numeric_limits<RealScalar>::digits<= 24? 5: // single precision 29 std::numeric_limits<RealScalar>::digits<= 53? 7: // double precision 30 std::numeric_limits<RealScalar>::digits<= 64? 8: // extended precision 31 std::numeric_limits<RealScalar>::digits<=106? 10: // double-double 40 typedef typename MatrixType::RealScalar RealScalar; in matrix_log_compute_2x2() typedef 58 else if ((abs(A(0,0)) < RealScalar(0.5)*abs(A(1,1))) || (abs(A(0,0)) > 2*abs(A(1,1)))) in matrix_log_compute_2x2() 65 …ingNumber = static_cast<int>(ceil((imag(logA11 - logA00) - RealScalar(EIGEN_PI)) / RealScalar(2*EI… in matrix_log_compute_2x2() 134 typedef typename NumTraits<typename MatrixType::Scalar>::Real RealScalar; in matrix_log_compute_pade() typedef 139 const RealScalar nodes[][maxPadeDegree] = { in matrix_log_compute_pade() [all …]
|
/external/eigen/Eigen/src/misc/ |
D | RealSvd2x2.h | 18 template<typename MatrixType, typename RealScalar, typename Index> 20 JacobiRotation<RealScalar> *j_left, in real_2x2_jacobi_svd() 21 JacobiRotation<RealScalar> *j_right) in real_2x2_jacobi_svd() 25 Matrix<RealScalar,2,2> m; in real_2x2_jacobi_svd() 28 JacobiRotation<RealScalar> rot1; in real_2x2_jacobi_svd() 29 RealScalar t = m.coeff(0,0) + m.coeff(1,1); in real_2x2_jacobi_svd() 30 RealScalar d = m.coeff(1,0) - m.coeff(0,1); in real_2x2_jacobi_svd() 32 if(abs(d) < (std::numeric_limits<RealScalar>::min)()) in real_2x2_jacobi_svd() 34 rot1.s() = RealScalar(0); in real_2x2_jacobi_svd() 35 rot1.c() = RealScalar(1); in real_2x2_jacobi_svd() [all …]
|
/external/eigen/test/ |
D | stable_norm.cpp | 26 typedef typename NumTraits<Scalar>::Real RealScalar; in stable_norm() typedef 34 ibeta = std::numeric_limits<RealScalar>::radix; // base for floating-point numbers in stable_norm() 35 … it = std::numeric_limits<RealScalar>::digits; // number of base-beta digits in mantissa in stable_norm() 36 iemin = std::numeric_limits<RealScalar>::min_exponent; // minimum exponent in stable_norm() 37 iemax = std::numeric_limits<RealScalar>::max_exponent; // maximum exponent in stable_norm() 42 Scalar inf = std::numeric_limits<RealScalar>::infinity(); in stable_norm() 43 if(NumTraits<Scalar>::IsComplex && (numext::isnan)(inf*RealScalar(1)) ) in stable_norm() 48 … mess up complex*real product, " << inf << " * " << 1.0 << " = " << inf*RealScalar(1) << std::endl; in stable_norm() 59 while(numext::abs2(factor)<RealScalar(1e-4)) in stable_norm() 61 Scalar big = factor * ((std::numeric_limits<RealScalar>::max)() * RealScalar(1e-4)); in stable_norm() [all …]
|
D | svd_fill.h | 25 typedef typename MatrixType::RealScalar RealScalar; typedef 28 RealScalar s = std::numeric_limits<RealScalar>::max_exponent10/4; 29 s = internal::random<RealScalar>(1,s); 30 Matrix<RealScalar,Dynamic,1> d = Matrix<RealScalar,Dynamic,1>::Random(diagSize); 32 d(k) = d(k)*pow(RealScalar(10),internal::random<RealScalar>(-s,s)); 68 samples << 0, four_denorms<RealScalar>(), 69 …RealScalar(1)/NumTraits<RealScalar>::highest(), RealScalar(1)/NumTraits<RealScalar>::highest(), (s…
|
/external/eigen/blas/ |
D | level2_cplx_impl.h | 19 …GEN_BLAS_FUNC(hemv)(const char *uplo, const int *n, const RealScalar *palpha, const RealScalar *pa… in EIGEN_BLAS_FUNC() 20 … const RealScalar *px, const int *incx, const RealScalar *pbeta, RealScalar *py, const int *incy) in EIGEN_BLAS_FUNC() 105 int EIGEN_BLAS_FUNC(hpr)(char *uplo, int *n, RealScalar *palpha, RealScalar *px, int *incx, RealSca… in EIGEN_BLAS_FUNC() 107 typedef void (*functype)(int, Scalar*, const Scalar*, RealScalar); in EIGEN_BLAS_FUNC() 117 RealScalar alpha = *palpha; in EIGEN_BLAS_FUNC() 149 …AS_FUNC(hpr2)(char *uplo, int *n, RealScalar *palpha, RealScalar *px, int *incx, RealScalar *py, i… in EIGEN_BLAS_FUNC() 197 int EIGEN_BLAS_FUNC(her)(char *uplo, int *n, RealScalar *palpha, RealScalar *px, int *incx, RealSca… in EIGEN_BLAS_FUNC() 209 RealScalar alpha = *reinterpret_cast<RealScalar*>(palpha); in EIGEN_BLAS_FUNC() 219 if(alpha==RealScalar(0)) in EIGEN_BLAS_FUNC() 244 …AS_FUNC(her2)(char *uplo, int *n, RealScalar *palpha, RealScalar *px, int *incx, RealScalar *py, i… in EIGEN_BLAS_FUNC() [all …]
|
D | level1_cplx_impl.h | 13 typedef RealScalar result_type; 15 inline RealScalar operator() (const Scalar& a) const { return numext::norm1(a); } in EIGEN_EMPTY_STRUCT_CTOR() 28 RealScalar EIGEN_CAT(EIGEN_CAT(REAL_SCALAR_SUFFIX,SCALAR_SUFFIX),asum_)(int *n, RealScalar *px, int… 40 int EIGEN_BLAS_FUNC(dotcw)(int *n, RealScalar *px, int *incx, RealScalar *py, int *incy, RealScalar… 63 int EIGEN_BLAS_FUNC(dotuw)(int *n, RealScalar *px, int *incx, RealScalar *py, int *incy, RealScalar… 84 RealScalar EIGEN_CAT(EIGEN_CAT(REAL_SCALAR_SUFFIX,SCALAR_SUFFIX),nrm2_)(int *n, RealScalar *px, int… 97 …X,REAL_SCALAR_SUFFIX),rot_)(int *n, RealScalar *px, int *incx, RealScalar *py, int *incy, RealScal… 103 RealScalar c = *pc; 104 RealScalar s = *ps; 120 …T(EIGEN_CAT(SCALAR_SUFFIX,REAL_SCALAR_SUFFIX),scal_)(int *n, RealScalar *palpha, RealScalar *px, i… [all …]
|
D | level1_impl.h | 12 …EIGEN_BLAS_FUNC(axpy)(const int *n, const RealScalar *palpha, const RealScalar *px, const int *inc… in EIGEN_BLAS_FUNC() 29 int EIGEN_BLAS_FUNC(copy)(int *n, RealScalar *px, int *incx, RealScalar *py, int *incy) in EIGEN_BLAS_FUNC() 54 int EIGEN_CAT(EIGEN_CAT(i,SCALAR_SUFFIX),amax_)(int *n, RealScalar *px, int *incx) in EIGEN_CAT() 65 int EIGEN_CAT(EIGEN_CAT(i,SCALAR_SUFFIX),amin_)(int *n, RealScalar *px, int *incx) in EIGEN_CAT() 76 int EIGEN_BLAS_FUNC(rotg)(RealScalar *pa, RealScalar *pb, RealScalar *pc, RealScalar *ps) in EIGEN_BLAS_FUNC() 83 RealScalar* c = pc; in EIGEN_BLAS_FUNC() 106 if (ab > aa && *c!=RealScalar(0)) in EIGEN_BLAS_FUNC() 113 RealScalar norm,scale; in EIGEN_BLAS_FUNC() 114 if(abs(a)==RealScalar(0)) in EIGEN_BLAS_FUNC() 116 *c = RealScalar(0); in EIGEN_BLAS_FUNC() [all …]
|
/external/eigen/unsupported/Eigen/src/LevenbergMarquardt/ |
D | LevenbergMarquardt.h | 117 typedef typename JacobianType::RealScalar RealScalar; typedef 151 m_ftol = sqrt(NumTraits<RealScalar>::epsilon()); in resetParameters() 152 m_xtol = sqrt(NumTraits<RealScalar>::epsilon()); in resetParameters() 158 void setXtol(RealScalar xtol) { m_xtol = xtol; } in setXtol() 161 void setFtol(RealScalar ftol) { m_ftol = ftol; } in setFtol() 164 void setGtol(RealScalar gtol) { m_gtol = gtol; } in setGtol() 167 void setFactor(RealScalar factor) { m_factor = factor; } in setFactor() 170 void setEpsilon (RealScalar epsfcn) { m_epsfcn = epsfcn; } in setEpsilon() 179 RealScalar xtol() const {return m_xtol; } in xtol() 182 RealScalar ftol() const {return m_ftol; } in ftol() [all …]
|
D | LMonestep.h | 25 RealScalar temp, temp1,temp2; in minimizeOneStep() 26 RealScalar ratio; in minimizeOneStep() 27 RealScalar pnorm, xnorm, fnorm1, actred, dirder, prered; in minimizeOneStep() 131 temp = RealScalar(.5); in minimizeOneStep() 133 temp = RealScalar(.5) * dirder / (dirder + RealScalar(.5) * actred); in minimizeOneStep() 134 if (RealScalar(.1) * fnorm1 >= m_fnorm || temp < RealScalar(.1)) in minimizeOneStep() 137 m_delta = temp * (std::min)(m_delta, pnorm / RealScalar(.1)); in minimizeOneStep() 139 } else if (!(m_par != 0. && ratio < RealScalar(.75))) { in minimizeOneStep() 140 m_delta = pnorm / RealScalar(.5); in minimizeOneStep() 141 m_par = RealScalar(.5) * m_par; in minimizeOneStep() [all …]
|
/external/eigen/unsupported/Eigen/src/IterativeSolvers/ |
D | MINRES.h | 33 typename Dest::RealScalar& tol_error) in minres() 36 typedef typename Dest::RealScalar RealScalar; in minres() typedef 41 const RealScalar rhsNorm2(rhs.squaredNorm()); in minres() 53 …const RealScalar threshold2(tol_error*tol_error*rhsNorm2); // convergence threshold (compared to r… in minres() 59 RealScalar residualNorm2(v_new.squaredNorm()); in minres() 63 RealScalar beta_new2(v_new.dot(w_new)); in minres() 65 RealScalar beta_new(sqrt(beta_new2)); in minres() 66 const RealScalar beta_one(beta_new); in minres() 70 RealScalar c(1.0); // the cosine of the Givens rotation in minres() 71 RealScalar c_old(1.0); in minres() [all …]
|
/external/eigen/unsupported/Eigen/src/Polynomials/ |
D | PolynomialSolver.h | 35 typedef typename NumTraits<Scalar>::Real RealScalar; typedef 36 typedef std::complex<RealScalar> RootType; 70 const RealScalar& absImaginaryThreshold = NumTraits<Scalar>::dummy_precision() ) const 86 RealScalar norm2 = numext::abs2( m_roots[0] ); in selectComplexRoot_withRespectToNorm() 89 const RealScalar currNorm2 = numext::abs2( m_roots[i] ); in selectComplexRoot_withRespectToNorm() 117 inline const RealScalar& selectRealRoot_withRespectToAbsRealPart( 120 const RealScalar& absImaginaryThreshold = NumTraits<Scalar>::dummy_precision() ) const 125 RealScalar abs2(0); 139 const RealScalar currAbs2 = m_roots[i].real() * m_roots[i].real(); 158 inline const RealScalar& selectRealRoot_withRespectToRealPart( [all …]
|
/external/eigen/Eigen/src/SparseCore/ |
D | SparseSparseProductWithPruning.h | 20 …(const Lhs& lhs, const Rhs& rhs, ResultType& res, const typename ResultType::RealScalar& tolerance) in sparse_sparse_product_with_pruning_impl() 89 typedef typename ResultType::RealScalar RealScalar; 91 static void run(const Lhs& lhs, const Rhs& rhs, ResultType& res, const RealScalar& tolerance) 102 typedef typename ResultType::RealScalar RealScalar; 103 static void run(const Lhs& lhs, const Rhs& rhs, ResultType& res, const RealScalar& tolerance) 116 typedef typename ResultType::RealScalar RealScalar; 117 static void run(const Lhs& lhs, const Rhs& rhs, ResultType& res, const RealScalar& tolerance) 129 typedef typename ResultType::RealScalar RealScalar; 130 static void run(const Lhs& lhs, const Rhs& rhs, ResultType& res, const RealScalar& tolerance) 149 typedef typename ResultType::RealScalar RealScalar; [all …]
|
/external/eigen/Eigen/src/QR/ |
D | ColPivHouseholderQR.h | 60 typedef typename MatrixType::RealScalar RealScalar; 67 typedef typename internal::plain_row_type<MatrixType, RealScalar>::type RealRowVectorType; 233 typename MatrixType::RealScalar absDeterminant() const; 247 typename MatrixType::RealScalar logAbsDeterminant() const; 259 RealScalar premultiplied_threshold = abs(m_maxpivot) * threshold(); 353 ColPivHouseholderQR& setThreshold(const RealScalar& threshold) 378 RealScalar threshold() const 384 … : NumTraits<Scalar>::epsilon() * RealScalar(m_qr.diagonalSize()); 403 RealScalar maxPivot() const { return m_maxpivot; } 442 RealScalar m_prescribedThreshold, m_maxpivot; [all …]
|
/external/eigen/unsupported/test/ |
D | matrix_function.cpp | 20 return ((a-b).array().abs() < test_precision<typename Type1::RealScalar>()).all(); in test_isApprox_abs() 30 typedef typename MatrixType::RealScalar RealScalar; in randomMatrixWithRealEivals() typedef 33 diag(i, i) = Scalar(RealScalar(internal::random<int>(0,2))) in randomMatrixWithRealEivals() 34 + internal::random<Scalar>() * Scalar(RealScalar(0.01)); in randomMatrixWithRealEivals() 84 typedef typename MatrixType::RealScalar RealScalar; in run() typedef 88 diag(i, i) = Scalar(RealScalar(internal::random<Index>(-1, 1))) * imagUnit in run() 89 + internal::random<Scalar>() * Scalar(RealScalar(0.01)); in run() 102 typedef typename NumTraits<Scalar>::Real RealScalar; in testMatrixExponential() typedef 103 typedef std::complex<RealScalar> ComplexScalar; in testMatrixExponential() 112 typedef typename NumTraits<Scalar>::Real RealScalar; in testMatrixLogarithm() typedef [all …]
|
/external/eigen/Eigen/src/Jacobi/ |
D | Jacobi.h | 37 typedef typename NumTraits<Scalar>::Real RealScalar; typedef 66 bool makeJacobi(const RealScalar& x, const Scalar& y, const RealScalar& z); 83 bool JacobiRotation<Scalar>::makeJacobi(const RealScalar& x, const Scalar& y, const RealScalar& z) in makeJacobi() 87 typedef typename NumTraits<Scalar>::Real RealScalar; in makeJacobi() typedef 88 RealScalar deno = RealScalar(2)*abs(y); in makeJacobi() 89 if(deno < (std::numeric_limits<RealScalar>::min)()) in makeJacobi() 97 RealScalar tau = (x-z)/deno; in makeJacobi() 98 RealScalar w = sqrt(numext::abs2(tau) + RealScalar(1)); in makeJacobi() 99 RealScalar t; in makeJacobi() 100 if(tau>RealScalar(0)) in makeJacobi() [all …]
|
/external/eigen/Eigen/src/Eigenvalues/ |
D | SelfAdjointEigenSolver.h | 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)) [all …]
|
/external/eigen/Eigen/src/SVD/ |
D | BDCSVD.h | 79 typedef typename NumTraits<typename MatrixType::Scalar>::Real RealScalar; 80 typedef typename NumTraits<RealScalar>::Literal Literal; 96 typedef Matrix<RealScalar, Dynamic, Dynamic, ColMajor> MatrixXr; 97 typedef Matrix<RealScalar, Dynamic, 1> VectorType; 98 typedef Array<RealScalar, Dynamic, 1> ArrayXr; 186 …atic RealScalar secularEq(RealScalar x, const ArrayRef& col0, const ArrayRef& diag, const IndicesR… 246 const RealScalar considerZero = (std::numeric_limits<RealScalar>::min)(); 262 RealScalar scale = matrix.cwiseAbs().maxCoeff(); 283 RealScalar a = abs(m_computed.coeff(i, i)); 399 const RealScalar considerZero = (std::numeric_limits<RealScalar>::min)(); [all …]
|
/external/eigen/Eigen/src/IterativeLinearSolvers/ |
D | LeastSquareConjugateGradient.h | 30 typename Dest::RealScalar& tol_error) in least_square_conjugate_gradient() 34 typedef typename Dest::RealScalar RealScalar; in least_square_conjugate_gradient() typedef 38 RealScalar tol = tol_error; in least_square_conjugate_gradient() 46 RealScalar rhsNorm2 = (mat.adjoint()*rhs).squaredNorm(); in least_square_conjugate_gradient() 54 RealScalar threshold = tol*tol*rhsNorm2; in least_square_conjugate_gradient() 55 RealScalar residualNorm2 = normal_residual.squaredNorm(); in least_square_conjugate_gradient() 67 …RealScalar absNew = numext::real(normal_residual.dot(p)); // the square of the absolute value of … in least_square_conjugate_gradient() 84 RealScalar absOld = absNew; in least_square_conjugate_gradient() 86 …RealScalar beta = absNew / absOld; // calculate the Gram-Schmidt value used to create… in least_square_conjugate_gradient() 160 typedef typename MatrixType::RealScalar RealScalar;
|
/external/eigen/Eigen/src/Cholesky/ |
D | LDLT.h | 62 typedef typename NumTraits<typename MatrixType::Scalar>::Real RealScalar; typedef 217 RealScalar rcond() const in rcond() 224 LDLT& rankUpdate(const MatrixBase<Derived>& w, const RealScalar& alpha=1); 279 RealScalar m_l1_norm; 298 typedef typename MatrixType::RealScalar RealScalar; 308 if (numext::real(mat.coeff(0,0)) > static_cast<RealScalar>(0) ) sign = PositiveSemiDef; 309 else if (numext::real(mat.coeff(0,0)) < static_cast<RealScalar>(0)) sign = NegativeSemiDef; 361 RealScalar realAkk = numext::real(mat.coeffRef(k,k)); 362 bool pivot_is_valid = (abs(realAkk) > RealScalar(0)); 384 if (realAkk < static_cast<RealScalar>(0)) sign = Indefinite; [all …]
|