Fortran计算奇异矩阵问题

xinrenxuef · 发表于 2020-6-16 18:08:39

编写的程序中，有一段代码是计算矩阵的逆矩阵的，但是计算到求逆矩阵的时候就会提示这个矩阵的逆矩阵是奇异矩阵， Singular matrix from CMatrixInv with D =
(0.000000000000000E+000,0.000000000000000E+000)，但是我在代码的下面输入write（*，*）命令让其显示该矩阵时，就会直接显示矩阵，也不会在提示奇异矩阵，请教各位大佬这是怎么回事？

li913 · 发表于 2020-6-16 19:56:23

代码有bug，只能自己调试

necrohan · 发表于 2020-6-17 08:14:36

近似奇异矩阵求逆有专门的算法，另找吧。输出和奇异矩阵没任何关系。

xinrenxuef · 发表于 2020-6-17 10:57:34

necrohan 发表于 2020-6-17 08:14
近似奇异矩阵求逆有专门的算法，另找吧。输出和奇异矩阵没任何关系。

刚开始用Fortran也不太懂，这是个1024次循环变量的，变量在其他数值时都没事，但是循环到4，5，11时就出现了这个问题，让变量在4时，单独计算也会出现这样的显示，但是这时候在这个代码前后在输入write命令，让它显示矩阵，这个时候矩阵就会正常显示，也不会显示出错误，所以一直没搞懂怎么回事，另外虽然都这样显示，但是每次都能计算出结果。有没有什么能查找一下问题，请大佬赐教

necrohan · 发表于 2020-6-18 08:03:31

这不是fortran的问题，是算法的问题

[数值问题] Fortran计算奇异矩阵问题

规矩勋章