测测你的双精度矩阵乘法-1k 排行榜

题目描述

给定两个 n * n 大小的双精度矩阵 A 和 B，要求计算 C = AB。

接口

void matrix_multiply(int n, const double *A, const double *B, double *C);

数据范围

n 等于 1024

注意事项

A 和 B 矩阵内的所有元素在 [-1, 1] 之内均匀分布。

你的答案被认为是正确的，当且仅当你给出的矩阵 C 中每个元素的值与标准答案的值相差不超过 3 * n * n * DBL_EPSILON。其中 DBL_EPSILON 大约为 2e-16，其准确值可在 C 语言头文件 float.h 中找到。

在 matrix_multiply 函数中，A, B, C 这三个变量的地址都是 4096 字节对齐的。即，表达式 (long) A % 4096 == 0 && (long) B % 4096 == 0 && (long) C % 4096 == 0 为真。