如何正确比较 matlab/octave 和 C++ 之间的结果

Question

我正在将一些在 Matlab/Octave 中可用的代码写入 C++ 代码。我只有八度，所以从现在开始我只会说八度。 我想正确比较八度音程代码和 C++ 代码之间的结果。我正在编写的算法将 2D 矩阵作为输入，并输出另一个 2D 矩阵。

为了比较结果，我使用函数save A.mat A 从八度音阶写入输入矩阵A，并带有默认选项。这将创建一个 ascii 文件 A.mat，其开头类似于

# Created by Octave 3.8.1, Tue May 27 12:12:53 2014 CEST <remi@desktop>
# name: values
# type: matrix
# rows: 25
# columns: 5 
43.0656 6.752420000000001 68.39323 35.75617 98.85446
...

我使用 octave 运行算法，并同样保存输出矩阵 B。

在我的 C++ 代码中，我使用以下代码加载矩阵 A 和 B：

// I opened the file A.mat with std::ifstream infile(filename);
// and read the first lines starting by # and loaded the matrix dimensions
std::string buffer;
double* matBuffer = new double[rows*cols];
while (std::getline(infile, buffer)) {
  std::istringstream iss(buffer);
  while (iss >> *matBuffer) {
    matBuffer++;
  }
}

然后我使用从 A.mat 读取的值运行 C++ 代码，并通过计算 B 读取与 B 计算的系数上的 mean squared error(MSE) 将结果与从 B.mat 读取的值进行比较。

但是，通过这样的设计，我可以期望 MSE 在 C++ 和 octave 代码之间为 0 吗？当然，我用同一台机器在 octave 和 C++ 上进行计算。但是由于在文件中写入/读取矩阵而导致的精度损失呢？另外，我假设倍频程矩阵的系数默认以双精度存储，这是正确的吗？

Answer 1

我可以期望 C++ 和 octave 代码之间的 MSE 为 0 吗？

我不这么认为，因为转换的层次很多，精度损失是很难避免的。

另外，我假设倍频程矩阵的系数默认存储为双精度，这样对吗？

Octave 使用 double 精度来表示值的内部表示，但在以 ASCII 格式存储值时，精度可能会有所损失。

我建议您尝试使用二进制格式来存储值，这将排除任何精度问题。您可以使用 HDF5 格式使用

save -hdf5 A.mat A

然后您可以使用HDF5 API 读取您的 CPP 应用程序中的值。