diff options
| author |  Gael Guennebaud <g.gael@free.fr> | 2014-04-25 11:48:22 +0200 |
|---|---|---|
| committer | Gael Guennebaud <g.gael@free.fr> | 2014-04-25 11:48:22 +0200 |
| commit | c9788d55b910c4bd8fb6cdc7978695f6b5d97e4c (patch) | |
| tree | 351775be9471d3b3dab4a3dd63244e77046148cb | |
| parent | ae4d9434e23e62d9403e570c20aeb3b8b44a2dd3 (diff) | |
Disable 3pX4 kernel on Altivec: despite this platform has 32 registers, this version seems significantly slower.
| -rw-r--r-- | Eigen/src/Core/products/GeneralBlockPanelKernel.h | 26 |
1 files changed, 13 insertions, 13 deletions
diff --git a/Eigen/src/Core/products/GeneralBlockPanelKernel.h b/Eigen/src/Core/products/GeneralBlockPanelKernel.h index e76c12c39..41c46c67a 100644 --- a/Eigen/src/Core/products/GeneralBlockPanelKernel.h +++ b/Eigen/src/Core/products/GeneralBlockPanelKernel.h @@ -188,7 +188,7 @@ public: nr = 4, // register block size along the M direction (currently, this one cannot be modified) -#ifdef EIGEN_HAS_FUSED_MADD +#if defined(EIGEN_HAS_FUSED_MADD) && !defined(EIGEN_VECTORIZE_ALTIVEC) // we assume 16 registers mr = 3*LhsPacketSize, #else @@ -296,7 +296,7 @@ public: NumberOfRegisters = EIGEN_ARCH_DEFAULT_NUMBER_OF_REGISTERS, nr = 4, -#ifdef EIGEN_HAS_FUSED_MADD +#if defined(EIGEN_HAS_FUSED_MADD) && !defined(EIGEN_VECTORIZE_ALTIVEC) // we assume 16 registers mr = 3*LhsPacketSize, #else @@ -759,29 +759,29 @@ void gebp_kernel<LhsScalar,RhsScalar,Index,mr,nr,ConjugateLhs,ConjugateRhs> for(Index k=0; k<peeled_kc; k+=pk) { EIGEN_ASM_COMMENT("begin gegp micro kernel 3p x 4"); - RhsPacket B_0; + RhsPacket B_0, T0; LhsPacket A2; - + #define EIGEN_GEBGP_ONESTEP(K) \ internal::prefetch(blA+(3*K+16)*LhsProgress); \ traits.loadLhs(&blA[(0+3*K)*LhsProgress], A0); \ traits.loadLhs(&blA[(1+3*K)*LhsProgress], A1); \ traits.loadLhs(&blA[(2+3*K)*LhsProgress], A2); \ traits.loadRhs(&blB[(0+4*K)*RhsProgress], B_0); \ - traits.madd(A0, B_0, C0, B_0); \ - traits.madd(A1, B_0, C4, B_0); \ + traits.madd(A0, B_0, C0, T0); \ + traits.madd(A1, B_0, C4, T0); \ traits.madd(A2, B_0, C8, B_0); \ traits.loadRhs(&blB[1+4*K*RhsProgress], B_0); \ - traits.madd(A0, B_0, C1, B_0); \ - traits.madd(A1, B_0, C5, B_0); \ + traits.madd(A0, B_0, C1, T0); \ + traits.madd(A1, B_0, C5, T0); \ traits.madd(A2, B_0, C9, B_0); \ traits.loadRhs(&blB[2+4*K*RhsProgress], B_0); \ - traits.madd(A0, B_0, C2, B_0); \ - traits.madd(A1, B_0, C6, B_0); \ + traits.madd(A0, B_0, C2, T0); \ + traits.madd(A1, B_0, C6, T0); \ traits.madd(A2, B_0, C10, B_0); \ traits.loadRhs(&blB[3+4*K*RhsProgress], B_0); \ - traits.madd(A0, B_0, C3 , B_0); \ - traits.madd(A1, B_0, C7, B_0); \ + traits.madd(A0, B_0, C3 , T0); \ + traits.madd(A1, B_0, C7, T0); \ traits.madd(A2, B_0, C11, B_0) internal::prefetch(blB+(48+0)); @@ -802,7 +802,7 @@ void gebp_kernel<LhsScalar,RhsScalar,Index,mr,nr,ConjugateLhs,ConjugateRhs> // process remaining peeled loop for(Index k=peeled_kc; k<depth; k++) { - RhsPacket B_0; + RhsPacket B_0, T0; LhsPacket A2; EIGEN_GEBGP_ONESTEP(0); blB += 4*RhsProgress; |