|
我想对WIJGUU这个矩阵求逆,前面我已经定义了COMPLEX*16 WIJGUU(NK/2,NK/2),并且open了文件,据我了解,前面这个子程序是LU分解,后面这个子程序是求逆,两个子程序都是联合起来使用的,我在运行中时,运算老是直接到100处,从而导致我后面算的数据write不到我open的文件中(100是我do循环的终点,我后续的运算必须在循环内,不能把100移到write上面),这两句是我根据其他的程序仿写的,不知道为什么它能把数据写到文件中我这个写不到文件中?这个IF(INFO1.NE.0)和IF(INFO2.NE.0)中的INFO1和INFO2是什么意思呢?
|
|