dgemm-naive.c

void dgemm( int m, int n, float *A, float *C )
{
  for( int i = 0; i < m; i++ )
    for( int k = 0; k < n; k++ ) 
      for( int j = 0; j < m; j++ ) 
	C[i+j*m] += A[i+k*m] * A[j+k*m];
}