GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-09-12 17:33:15 +08:00
Code Issues Packages Projects Releases Wiki Activity

Fixed MPRealSupport

Browse Source
This commit is contained in:
Christoph Hertzberg 2015-02-28 16:41:00 +01:00
parent 33f40b2883
commit 682196e9fc
2 changed files with 18 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
unsupported/Eigen/MPRealSupport
View File

@ -141,20 +141,32 @@ int main()
public: public:
typedef mpfr::mpreal ResScalar; typedef mpfr::mpreal ResScalar;
enum { enum {
Vectorizable = false,
LhsPacketSize = 1,
RhsPacketSize = 1,
ResPacketSize = 1,
NumberOfRegisters = 1,
nr = 1, nr = 1,
mr = 1, mr = 1,
LhsProgress = 1, LhsProgress = 1,
RhsProgress = 1 RhsProgress = 1
}; };
typedef ResScalar LhsPacket;
typedef ResScalar RhsPacket;
typedef ResScalar ResPacket;
}; };
template<typename Index, bool ConjugateLhs, bool ConjugateRhs>
struct gebp_kernel<mpfr::mpreal,mpfr::mpreal,Index,1,1,ConjugateLhs,ConjugateRhs>
template<typename Index, typename DataMapper, bool ConjugateLhs, bool ConjugateRhs>
struct gebp_kernel<mpfr::mpreal,mpfr::mpreal,Index,DataMapper,1,1,ConjugateLhs,ConjugateRhs>
{ {
typedef mpfr::mpreal mpreal; typedef mpfr::mpreal mpreal;
EIGEN_DONT_INLINE EIGEN_DONT_INLINE
void operator()(mpreal* res, Index resStride, const mpreal* blockA, const mpreal* blockB, Index rows, Index depth, Index cols, mpreal alpha, void operator()(const DataMapper& res, const mpreal* blockA, const mpreal* blockB,
Index rows, Index depth, Index cols, const mpreal& alpha,
Index strideA=-1, Index strideB=-1, Index offsetA=0, Index offsetB=0) Index strideA=-1, Index strideB=-1, Index offsetA=0, Index offsetB=0)
{ {
if(rows==0 || cols==0 || depth==0) if(rows==0 || cols==0 || depth==0)
@ -170,8 +182,6 @@ int main()
{ {
for(Index j=0; j<cols; ++j) for(Index j=0; j<cols; ++j)
{ {
mpreal *C1 = res + j*resStride;
const mpreal *A = blockA + i*strideA + offsetA; const mpreal *A = blockA + i*strideA + offsetA;
const mpreal *B = blockB + j*strideB + offsetB; const mpreal *B = blockB + j*strideB + offsetB;
@ -183,7 +193,7 @@ int main()
} }
mpfr_mul(acc1.mpfr_ptr(), acc1.mpfr_srcptr(), alpha.mpfr_srcptr(), mpreal::get_default_rnd()); mpfr_mul(acc1.mpfr_ptr(), acc1.mpfr_srcptr(), alpha.mpfr_srcptr(), mpreal::get_default_rnd());
mpfr_add(C1[i].mpfr_ptr(), C1[i].mpfr_srcptr(), acc1.mpfr_srcptr(), mpreal::get_default_rnd()); mpfr_add(res(i,j).mpfr_ptr(), res(i,j).mpfr_srcptr(), acc1.mpfr_srcptr(), mpreal::get_default_rnd());
} }
} }
} }

3
unsupported/test/mpreal/mpreal.h
View File

@ -57,7 +57,8 @@
#include <limits> #include <limits>
// Options // Options
#define MPREAL_HAVE_INT64_SUPPORT // Enable int64_t support if possible. Available only for MSVC 2010 & GCC. // FIXME HAVE_INT64_SUPPORT leads to clashes with long int and int64_t on some systems.
//#define MPREAL_HAVE_INT64_SUPPORT // Enable int64_t support if possible. Available only for MSVC 2010 & GCC.
#define MPREAL_HAVE_MSVC_DEBUGVIEW // Enable Debugger Visualizer for "Debug" builds in MSVC. #define MPREAL_HAVE_MSVC_DEBUGVIEW // Enable Debugger Visualizer for "Debug" builds in MSVC.
#define MPREAL_HAVE_DYNAMIC_STD_NUMERIC_LIMITS // Enable extended std::numeric_limits<mpfr::mpreal> specialization. #define MPREAL_HAVE_DYNAMIC_STD_NUMERIC_LIMITS // Enable extended std::numeric_limits<mpfr::mpreal> specialization.
// Meaning that "digits", "round_style" and similar members are defined as functions, not constants. // Meaning that "digits", "round_style" and similar members are defined as functions, not constants.