彭国伦---《Fortran95程序设计》----------纠错

FLY · 发表于 2014-7-22 14:20:56

感觉第143页EX0714中第七行和第九行的注释错了，第七行为下三角，第九行为上三角，如果我说错了，勿喷。

FLY · 发表于 2014-7-23 10:20:10

fcode 发表于 2014-7-23 07:09
这个没有错哦。Fortran 是列优先的，所以 a(i,j) 一般称呼 I 为列，J 为行，即第 J 列，第 I 行。
I>J ...

额，我还是有一些疑问。个人理解：
a(i,j)中，i表示的行，j表示的列，fortran中列优先指的是在内存里储存数据的时候按照第一列、第二列的顺序储存，比如对于a(5,5)这个数组按照a(1,1),a(2,1).....a(5,1),a(2,2)...a(5,2)...这样的顺序储存的。所以i表示的行，j表示的列，i>j，下半部分啊。而且在矩阵论里面不也是这样来定义的么？
求解答。

FLY · 发表于 2014-7-23 10:25:35

第455页中，逆阵的求解的结果感觉有点不对，逆阵中第二行第一列-0.63应该是-0.063，最后一行最后一列-0.63也应该是-0.063。我用的MKL中sgetrf和sgetri函数算出来是这个结果

[通用算法] 彭国伦---《Fortran95程序设计》----------纠错

评分