代码

#define n 1024
void matrix_multiply(int _, const double* A, const double* B, double* C) {
	double s0, s1, s2, s3;
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			s0 = s1 = s2 = s3 = 0;
			const double* a = A + (i << 10);
			const double* b = B + j;
			for (int k = 0; k < n; k += 4) {
				s0 += (*a) * (*b);
				s1 += (*a + 1) * (*b + (1 << 10));
				s2 += (*a + 2) * (*b + (2 << 10));
				s3 += (*a + 3) * (*b + (3 << 10));
				a += 4;
				b += 4 << 10;
			}
			C[(i << 10) | j] = s0 + s1 + s2 + s3;
		}
	}
}

评测结果