Re: SIMD benchmark

[a]

On Wednesday, 18 January 2012 at 01:50:00 UTC, Timon Gehr wrote:
Anyway, I was after a general matrix*matrix multiplication, 
where the operands can get arbitrarily large and where any 
potential use of __restrict is rendered unnecessary by array 
vector ops.

Here you go. But I agree there are use cases for restrict where 
vector operations don't help

void matmul(A,B,C)(A a, B b, C c, size_t n, size_t m, size_t l)
{
   for(size_t i = 0; i < n; i++)
   {
       c[i*l..i*l + l] = 0;
       for(size_t j = 0; j < m; j++)
           c[i*l..i*l + l] += a[i*m + j] * b[j*l..j*l + l];
   }
}