代码

#define idx(i, j) ((i) * n + (j))
void sgemm(const double *A, const double *B, double *C, int m, int n, int k)
{
    for (int i = 0; i < m; i++)
    {
        for (int j = 0; j < n; j++)
        {
            C[i * n + j] = 0;
            for (int l = 0; l < k; l++)
            {
                C[i * n + j] += A[i * k + l] * B[l * n + j];
            }
        }
    }
}
void matrix_multiply(int n, const double *A, const double *B, double *C) {
	sgemm(A,B,C,1024,1024,1024);
}

评测结果