Armadillo 保存双精度 ascii

Question

我想用 armadillo 库以 ascii 格式保存一些复杂的双精度矩阵，但不知何故只保存了 6 位有效数字...我知道这种截断（或舍入，尚未检查）可以通过使用二进制格式来克服，但我想使用 ascii 以便我可以手动验证和修改一些值。在下面提供的 mwe 中，我希望矩阵 A 和 C 对于大约 16 位有效数字是相同的。是否有一个秘密开关告诉犰狳以更高的精度以ascii格式保存复杂矩阵？谢谢

#include <iostream>
#include <armadillo>

using namespace std;
using namespace arma;

int main()
{
cout.precision(20);
cx_mat A = randu<cx_mat>(2,2);
A(0,0) = 1.2345678901234567890;
A.raw_print(cout);
A.save("A.arma_ascii", arma_ascii);
cout << endl;
cx_mat C;
C.load("A.arma_ascii", arma_ascii);
C.raw_print(cout);
return 0;
}

EDIT 该问题只出现在复数矩阵中。对于真实的犰狳，将精度设置为 diskio_meat.hpp。似乎解决此限制的唯一方法是修补犰狳或单独保存复数和虚数部分。

编辑 2 问题在犰狳库 (git commit) 中得到解决，仅在此处发布我的问题 2 小时后。干得好！

Answer 1

我认为问题在于 cx_mat 的基础数据类型是 complex<double>，因此您不能期望超过 15-17 位有效小数位的精度，请参阅 https://en.wikipedia.org/wiki/Double-precision_floating-point_format。在 Armadillo 代码中，他们使用 std::ofstream 作为输出，默认显示 6 位数字，就像您在 A.arma_ascii 文件中一样。

Armadillo 保存双精度 ascii

Armadillo save double precision ascii

c++

armadillo