代码

const int n = 1024;
void matrix_multiply(int _, const double* A, const double* B, double* C) {
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			double& x = C[(i << 10) | j];
			x = 0;
			for (int k = 0; k < n; k++) {
				x += A[(i << 10) | k] * B[(k << 10) | j];
			}
		}
	}
}

评测结果