LOOP快还是直接运算快

Kong · 发表于 2020-7-17 11:22:53

我有一个数组 fields_f(i, j, k) 和 f_f_k(i, j, k) ，想请教一下这样的运算

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

fields_f = fields_f + f_f_k * dt

直接利用fortran的特性这样算快呢，还是通过以下的方法利用我计算机的48个线程各自循环算快呢

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

!$omp parallel do default(shared) private(i,j,k)
LOOP
fields_f(i, j, k) = fields_f(i, j, k) + f_f_k(i, j, k) * dt
ENDLOOP
!$omp end parallel do

Kong · 发表于 2020-7-17 11:36:09

好像应该发在算法区的

liudy02 · 发表于 2020-7-17 23:43:42

把计算量搞大点，自己加点计时代码测试下不就什么都清楚了……

[数值问题] LOOP快还是直接运算快