对称正定矩阵分解，计算结果错误，实在找不到原因了

pasuka

matlab给出的结果，lz答案哪里错了？！
15.279446979521216                   0                   0                   0
   7.450681962022638   4.805271927871565                   0                   0
  16.758610461700346   0.534147362099914   0.579449729096019                   0
   9.494433940864122   5.112903632921073   5.217080776328907   6.142467643894371

估计是双精度和单精度的浮点数误差
btw，宋某某那本书翻过几页，真是一本烂书，估计是混职称用的

pasuka

likm1110 发表于 2015-3-27 09:40
但是老师让我把那本书所有程序都打一遍我也没办法。还得边打便给他改错我也是醉了。。。。。另外就是我手 ...

cholesky.zip (1.11 KB, 下载次数: 8)

2015-3-27 12:08 上传

点击文件名下载附件

单精度和双精度的问题，给你改了代码，自己慢慢琢磨吧
cholesky_test.f90 (1.27 KB, 下载次数: 7)

2015-3-27 12:06 上传

点击文件名下载附件