我将Eigen-c ++用于线性弹性的有限元应用。我在MATLAB中开发了一个原型代码并对其进行了测试。现在我试图用Eigen-c ++将它翻译成c ++。
我使用稀疏矩阵乘法，不幸的是我有一些问题;可能精确。我在计算中使用double个数字，我观察到MATLAB的最终矩阵与Eigen-c ++的最终矩阵之间存在很大差异。
即，MATLAB中1.6397029e+13的值为1.63974e+13，带有eigen-c ++。这是数百万的差异，并产生错误的矩阵解决方案。
我该如何解决这个问题？