我考虑过这个线性系统
Bx =（A + sI）x = b，s> 0 其中s = 10 ^ 4且A是nxn矩阵，n = 986703。
我在matlab中使用ichol命令从As计算Cholesky因子L.在此之后，我实现了CG算法来解决Bx = b的问题。 b是从B * x计算的向量，其中x = 1（n，1）。所以我知道问题的“分析”解决方案x。
当我执行CG算法时，我得到一个x_tilde并计算x_tilde和x之间的最大分量相对误差。
理论上是否可以使用双精度数将此分量相对误差降低到2 ^ -53以下？
问候