GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-05-01 16:24:28 +08:00
Code Issues Packages Projects Releases Wiki Activity

Product kernel: skip loop on columns if there is no remaining rows

Browse Source
This commit is contained in:
Gael Guennebaud 2014-04-25 16:54:30 +02:00
parent 6f64b0b487
commit f9d2f3903e
1 changed files with 1 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
Eigen/src/Core/products/GeneralBlockPanelKernel.h
View File

@ -1255,6 +1255,7 @@ void gebp_kernel<LhsScalar,RhsScalar,Index,mr,nr,ConjugateLhs,ConjugateRhs>
} }
} }
//---------- Process remaining rows, 1 at once ---------- //---------- Process remaining rows, 1 at once ----------
if(peeled_mc1<rows)
{ {
// loop on each panel of the rhs // loop on each panel of the rhs
for(Index j2=0; j2<packet_cols4; j2+=nr) for(Index j2=0; j2<packet_cols4; j2+=nr)
@ -1288,13 +1289,11 @@ void gebp_kernel<LhsScalar,RhsScalar,Index,mr,nr,ConjugateLhs,ConjugateRhs>
straits.loadLhsUnaligned(blB+0*SwappedTraits::LhsProgress, A0); straits.loadLhsUnaligned(blB+0*SwappedTraits::LhsProgress, A0);
straits.loadLhsUnaligned(blB+1*SwappedTraits::LhsProgress, A1); straits.loadLhsUnaligned(blB+1*SwappedTraits::LhsProgress, A1);
straits.loadRhsQuad(blA+0*spk, B_0); straits.loadRhsQuad(blA+0*spk, B_0);
straits.loadRhsQuad(blA+1*spk, B_1); straits.loadRhsQuad(blA+1*spk, B_1);
straits.madd(A0,B_0,C0,B_0); straits.madd(A0,B_0,C0,B_0);
straits.madd(A1,B_1,C1,B_1); straits.madd(A1,B_1,C1,B_1);
straits.loadLhsUnaligned(blB+2*SwappedTraits::LhsProgress, A0); straits.loadLhsUnaligned(blB+2*SwappedTraits::LhsProgress, A0);
straits.loadLhsUnaligned(blB+3*SwappedTraits::LhsProgress, A1); straits.loadLhsUnaligned(blB+3*SwappedTraits::LhsProgress, A1);
straits.loadRhsQuad(blA+2*spk, B_0); straits.loadRhsQuad(blA+2*spk, B_0);