老版本fortran做循环manova代码

时间:2017-10-27 23:29:55

标签: c++ fortran

我试图复制1971年出版的多变量统计代码。未指定Fortran版本。我将int移植到C ++并且不理解以下代码段; 来自Cooley和Lohnes - 多变量数据分析

DO 12 NS = 1, NG
READ(5, TIT)        (V(J),  J = 1, M)
DO 12   J = 1, M
U(J) =  U(J) + V(J)
T(J) =  T(J) + V(J)
DO 12   K = 1, M
A(J, K) = A(J,K) + V(J) * V(K)

12 C(J,K)= C(J,K)+ V(J)* V(K)

具体来说,我不明白V(J)如何在第二行填充。非常感谢Fortran大师的任何帮助。 BTW TIT是一个包含项目描述的数组,因此我不明白为什么要在这里阅读。

1 个答案:

答案 0 :(得分:0)

V(J), J=1,M基本上是一个隐式循环,所以它大致等同于

for (i=1; i <= m; i++)
    infile >> v[i];