代码之家 › 专栏 › 技术社区 › idladle

如何在Fortran中计算数组元素的乘积之和?

fortran

idladle · 技术社区 · 7 年前

鉴于 A(I,J,K) 具有 I =1,2,3和 J =1,2,3,我想取下面的和

B(K) = c(1)*c(1)*A(1,1,K) + c(1)*c(2)*A(1,2,K) + c(1)*c(3)*A(1,3,K) + &
       c(2)*c(1)*A(2,1,K) + c(2)*c(2)*A(2,2,K) + c(2)*c(3)*A(2,3,K) + &
       c(3)*c(1)*A(3,1,K) + c(3)*c(2)*A(3,2,K) + c(3)*c(3)*A(3,3,K)

对于较大的 我 和 J 具有 c(I) 和 c(J) 为每个 我 或 J 。

如何更有效地编写此代码?我认为 DO WHILE 循环也许是答案,但我是个初学者,不知道怎么做。有人能给我一个提示吗?

我的尝试:

DO K = 1,100
   J = 1.d0
   DO WHILE (J .LE. 3)
      I = 1.d0
      DO WHILE (I .LE. 3)
         A(I,J,K) = c(I)*c(J)*A(I,J,K) ! + ???
         I = I + 1.d0
      END DO
   END DO
END DO

1 回复 | 直到 7 年前

L. Scott Johnson user3836982 7 年前

只需对J使用do循环,我就像你对K所做的那样。将总和累加到B(K)中,从0开始。

DO K = 1,100
B(K) = 0
   DO J = 1,3
      DO I = 1,3
         B(K) = B(K) + c(I)*c(J)*A(I,J,K) 
      END DO
   END DO
END DO

推荐文章

Fib · 查找R包中使用的Fortran文件

3 年前

ethanmorton · Fortran/MPI编译错误隐式NONE语句不能跟在属性声明后面

3 年前

Holmz · 使用ifort 2020和OneAPI更改字符串行为?

3 年前

imronuke · nvfortran编译器有用标志的建议

3 年前

chaemin · fortran中的某些计算机没有将数组的初始值指定为零。为什么?

3 年前

R. N · gfortran指针过程类型问题

7 年前

Gavin Chandler · 错误LNK1104:无法打开文件“DFOR”。lib的Visual Studio

7 年前

ptb · 将Torvalds的“好品味”应用于Fortran链表

7 年前

veryreverie · Fortran 2008中带有非多态过程的重载延迟过程

7 年前

Fabricio Fernández · 我怎样才能阅读。Fortran程序中的原始文件?

7 年前