C++ for Scientists - Technische Universität Dresden

182 CHAPTER 5. META-PROGRAMMING
};
}
C(i + Index0, k + Max1)= tmp.value;
next().update(tmp.sub, C, i, k);
template
struct mult block
{
template
void operator()(Tmp& tmp, const Matrix& A, const Matrix& B, unsigned i, unsigned j, unsigned k)
{
std::cout ≪ ”tmp.” ≪ tmp.bs ≪ ”+= A[” ≪ i + Max0 ≪ ”][” ≪ j ≪ ”] ∗ B[” ≪ j ≪ ”][” ≪
k + Max1 ≪ ”]\n”;
tmp.value+= A(i + Max0, j) ∗ B(j, k + Max1);
}
};
template
void update(const Tmp& tmp, Matrix& C, unsigned i, unsigned k)
{
std::cout ≪ ”C[” ≪ i + Max0 ≪ ”][” ≪ k + Max1 ≪ ”]= tmp.” ≪ tmp.bs ≪ ”\n”;
C(i + Max0, k + Max1)= tmp.value;
}
In order to verify that all operations are performed, we log them completely but look here only
at tmp.4 and tmp.3:
tmp.4+= A[1][0] * B[0][0]
tmp.3+= A[1][0] * B[0][1]
tmp.4+= A[1][1] * B[1][0]
tmp.3+= A[1][1] * B[1][1]
tmp.4+= A[1][2] * B[2][0]
tmp.3+= A[1][2] * B[2][1]
tmp.4+= A[1][3] * B[3][0]
tmp.3+= A[1][3] * B[3][1]
C[1][0]= tmp.4
C[1][1]= tmp.3
tmp.4+= A[3][0] * B[0][0]
tmp.3+= A[3][0] * B[0][1]
tmp.4+= A[3][1] * B[1][0]
tmp.3+= A[3][1] * B[1][1]
tmp.4+= A[3][2] * B[2][0]
tmp.3+= A[3][2] * B[2][1]
tmp.4+= A[3][3] * B[3][0]
tmp.3+= A[3][3] * B[3][1]
C[3][0]= tmp.4
C[3][1]= tmp.3