21 int main(
int argc,
char* argv[])
25 Scalar*
a = internal::aligned_new<Scalar>(size2);
26 Scalar*
b = internal::aligned_new<Scalar>(size2+4)+1;
27 Scalar*
c = internal::aligned_new<Scalar>(size2);
37 for (
int k=0; k<10; ++k)
45 for (
int innersize =
size; innersize>2 ; --innersize)
47 if (size2%innersize==0)
49 int outersize = size2/innersize;
54 for (
int k=0; k<3; ++k)
60 std::cout << innersize <<
" x " << outersize <<
" " <<
timer.
value() <<
"s " << (double(size2*
REPEAT)/
timer.
value())/(1024.*1024.*1024.) <<
" GFlops\n";
68 for (
int k=0; k<3; ++k)
81 for (
int k=0; k<
REPEAT; ++k)
87 for (
int k=0; k<
REPEAT; ++k)
96 for (
int k=0; k<
REPEAT; ++k)
97 for (
int i=0;
i<
size;
i+=PacketSize*8)