Lines Matching full:mpreal
16 #include <mpreal.h>
30 * You can find a copy of MPFR C++ that is known to be compatible in the unsupported/test/mpreal folder.
43 mpreal::set_default_prec(256);
45 typedef Matrix<mpreal,Dynamic,Dynamic> MatrixXmp;
46 typedef Matrix<mpreal,Dynamic,1> VectorXmp;
60 template<> struct NumTraits<mpfr::mpreal>
61 : GenericNumTraits<mpfr::mpreal>
73 typedef mpfr::mpreal Real;
74 typedef mpfr::mpreal NonInteger;
76 inline static Real highest (long Precision = mpfr::mpreal::get_default_prec()) { return mpfr::maxval(Precision); }
77 inline static Real lowest (long Precision = mpfr::mpreal::get_default_prec()) { return -mpfr::maxval(Precision); }
80 inline static Real Pi (long Precision = mpfr::mpreal::get_default_prec()) { return mpfr::const_pi(Precision); }
81 inline static Real Euler (long Precision = mpfr::mpreal::get_default_prec()) { return mpfr::const_euler(Precision); }
82 inline static Real Log2 (long Precision = mpfr::mpreal::get_default_prec()) { return mpfr::const_log2(Precision); }
83 inline static Real Catalan (long Precision = mpfr::mpreal::get_default_prec()) { return mpfr::const_catalan(Precision); }
85 inline static Real epsilon (long Precision = mpfr::mpreal::get_default_prec()) { return mpfr::machine_epsilon(Precision); }
90 unsigned int weak_prec = ((mpfr::mpreal::get_default_prec()-1) * 90) / 100;
97 template<> inline mpfr::mpreal random<mpfr::mpreal>()
102 template<> inline mpfr::mpreal random<mpfr::mpreal>(const mpfr::mpreal& a, const mpfr::mpreal& b)
104 return a + (b-a) * random<mpfr::mpreal>();
107 inline bool isMuchSmallerThan(const mpfr::mpreal& a, const mpfr::mpreal& b, const mpfr::mpreal& eps)
112 inline bool isApprox(const mpfr::mpreal& a, const mpfr::mpreal& b, const mpfr::mpreal& eps)
117 inline bool isApproxOrLessThan(const mpfr::mpreal& a, const mpfr::mpreal& b, const mpfr::mpreal& eps)
122 template<> inline long double cast<mpfr::mpreal,long double>(const mpfr::mpreal& x)
125 template<> inline double cast<mpfr::mpreal,double>(const mpfr::mpreal& x)
128 template<> inline long cast<mpfr::mpreal,long>(const mpfr::mpreal& x)
131 template<> inline int cast<mpfr::mpreal,int>(const mpfr::mpreal& x)
134 // Specialize GEBP kernel and traits for mpreal (no need for peeling, nor complicated stuff)
135 // This also permits to directly call mpfr's routines and avoid many temporaries produced by mpreal
137 class gebp_traits<mpfr::mpreal, mpfr::mpreal, false, false>
140 typedef mpfr::mpreal ResScalar;
151 struct gebp_kernel<mpfr::mpreal,mpfr::mpreal,Index,mr,nr,ConjugateLhs,ConjugateRhs>
153 typedef mpfr::mpreal mpreal;
156 void operator()(mpreal* res, Index resStride, const mpreal* blockA, const mpreal* blockB, Index rows, Index depth, Index cols, mpreal alpha,
157 Index strideA=-1, Index strideB=-1, Index offsetA=0, Index offsetB=0, mpreal* /*unpackedB*/ = 0)
159 mpreal acc1, acc2, tmp;
167 mpreal *C1 = res + j*resStride;
168 mpreal *C2 = res + (j+1)*resStride;
171 mpreal *B = const_cast<mpreal*>(blockB) + j*strideB + offsetB*actual_nr;
172 mpreal *A = const_cast<mpreal*>(blockA) + i*strideA + offsetA;
177 mpfr_mul(tmp.mpfr_ptr(), A[k].mpfr_ptr(), B[0].mpfr_ptr(), mpreal::get_default_rnd());
178 mpfr_add(acc1.mpfr_ptr(), acc1.mpfr_ptr(), tmp.mpfr_ptr(), mpreal::get_default_rnd());
181 mpfr_mul(tmp.mpfr_ptr(), A[k].mpfr_ptr(), B[1].mpfr_ptr(), mpreal::get_default_rnd());
182 mpfr_add(acc2.mpfr_ptr(), acc2.mpfr_ptr(), tmp.mpfr_ptr(), mpreal::get_default_rnd());
188 mpfr_mul(acc1.mpfr_ptr(), acc1.mpfr_ptr(), alpha.mpfr_ptr(), mpreal::get_default_rnd());
189 mpfr_add(C1[i].mpfr_ptr(), C1[i].mpfr_ptr(), acc1.mpfr_ptr(), mpreal::get_default_rnd());
192 mpfr_mul(acc2.mpfr_ptr(), acc2.mpfr_ptr(), alpha.mpfr_ptr(), mpreal::get_default_rnd());
193 mpfr_add(C2[i].mpfr_ptr(), C2[i].mpfr_ptr(), acc2.mpfr_ptr(), mpreal::get_default_rnd());
