优化代码的fortran Do循环索引问题

文章/答案/技术大牛

发布

1回答

matrix、optimization、indexing、fortran、column-major-order

首先，我的英语不好。对不起。就我所知道。Fortran地址是主要列。我的旧Fortran代码在很长一段时间内都没有优化。我尝试更改我的Fortran90代码索引以获得更快的速度。一个代码几乎是三维矩阵。(i，j，k)和几乎Do-loop大约是i和j。i和j的大小大约是2000~3000，k正好是2，它意味着x，y。我的旧<em

浏览 36提问于2019-03-06得票数 0

回答已采纳

1回答

gfortran优化在运行时导致fortran do变量循环错误

optimization、iterator、fortran、do-loops

我已经写了一个fortran例程，它使用了一些遗留的fortran 77代码来处理有限元。然而，对于特定的网格，当-O优化标志打开时，一个重要的do循环迭代器正在以某种方式被修改，尽管fortran理应禁止这样做。我使用gfortran4.5编译了这段代码，并启用了-fcheck=do运行时检查，它验证了我上面提到的内容。只有当优化打开并直接指向

浏览 4提问于2011-04-08得票数 1

2回答

为什么Fortran DO循环索引大于循环后的上限？

loops、fortran

我正在研究一个用fortran编写的空气质量模型的代码，并对fortran从do循环中传递变量的方式有一些疑问。这个非常简单的例子说明了我的意思： IMPLICIT NONE WRITE(*,*) 'L = ', L ! WRITE(*,*) 'I

浏览 1提问于2012-03-27得票数 5

回答已采纳

1回答

fortran90:寻址内存位置

c、pointers、fortran90、memory-address

我有一个fortran90代码要优化。现在，我想访问外部循环中结构的内存位置，然后访问嵌套循环中最深的结构。就像这样：do i = 1, N jj = some other integer c= a(ii, jj)%b(i)enddo第二个for

浏览 1提问于2014-10-13得票数 0

2回答

Fortran forall限制

loops、fortran、fortran95

我还发现write语句在forall块中是被禁止的，我怀疑read语句也是如此。除了有时在顺序无关紧要的时候替换do循环之外，这个构造到底是为了什么？我认为这会使编码更加清晰和优雅，特别是在操作顺序不重要时显示出来，但它似乎对在forall中可以完成的操作有很大的限制。这些限制的原因是什么，即它们保护/防止用户搞砸的是什么？现在，在我正在处理的代码</em

浏览 1提问于2011-12-22得票数 12

回答已采纳

1回答

标量到阵列的内禀分配

fortran

来自Fortran 2008规范，7.2.1.3.5 A:do i=1,3enddointeger, dimension(3) :: xC:x = 1 什么是为数组分配标量值的最佳实践(性能和

浏览 2提问于2014-09-09得票数 3

回答已采纳

1回答

fortran编译器优化包含语句中的do循环

fortran

我有一个用于HPC的大型Fortran代码，其中有几个方程是在计算网格上求解的。因此，在源代码中的许多地方，都将嵌套的do循环循环在entair计算域上循环，shuch为 do j = 1,number_of_point_in_yend doen

浏览 1提问于2022-10-11得票数 0

2回答

复制数组:循环与数组操作

arrays、fortran

我现在和Fortran一起工作很长时间，但是我有一个问题，我找不到令人满意的答案。如果我有两个数组，并希望将其中一个复制到另一个数组中：... do j=0,100enddoreal,dimension(0:100,0:100) :: array1,array2array1 = array2 计算时间也有很大的

浏览 3提问于2015-08-03得票数 3

回答已采纳

2回答

Matlab中的Fortran逻辑索引

matlab、fortran

我正在尝试编写Fortran 95代码来模仿我在MATLAB中所做的事情，但我很难访问数组索引。代码比下面显示的代码复杂得多，但这就是要点。我更喜欢避免do循环。例如--> Matlab命令。如何在fortran中存储和使用"a“数组中的那些索引，并对其他数组中的相同索引进行操作？

浏览 0提问于2017-09-30得票数 0

1回答

Fortran编译错误(预条件共轭梯度法)：我试图编译fortran代码并得到错误

loops、compiler-errors、fortran

我试着编译下面的代码，得到了下面的错误，伙计们，你们能帮我解决这个错误吗？代码： DIMENSION aa(*),ja(*),ia(*),al(*) DO 1 K=1,DO 3 i=1,n h(i)=g(i)

浏览 10提问于2021-12-09得票数 -1

1回答

在fortran中高效构造结构矩阵

performance、matrix、fortran、fortran90

离开Fortran已经几年了，现在我不得不拿起它，重新开始工作。我想用entry(i,j)构造一个f(x_i,y_j)形式的矩阵，其中f是两个变量的函数，例如f(x,y)=cos(x-y)。在Matlab或Python(Numpy)中，有一些有效的方法来处理这种特定的问题。我不知道Fortran是否有这样的优化。顺便说一句，在Fortran中，向量化操作比do/for循环</em

浏览 3提问于2015-03-17得票数 0

回答已采纳

2回答

将矩阵打印到文件中

fortran

我写了这几行WRITE(105,*) DFGRD1DOi=1,3END DO 然而，答案并不相同。它们是彼此的转置。由于我无法访问DFGRD1的实际值，因此我不知道哪一个是正确的。如果有人能帮我一下，我将不胜感激。

浏览 1提问于2019-01-30得票数 1

1回答

在非并行的依赖时间的do循环中使用OpenMP段并行

multithreading、parallel-processing、fortran、openmp

关于OpenMP，我有一个简短的问题。通常我们可以像这样做一个并行的部分(用fortran编写，有两个部分)：现在我真正想要在do循环中运行fortran代码块A和B，这个循环本身不应该并行化，因为这个do循环是一个依赖时间的循环，每一步都依赖于前一步的结果

浏览 2提问于2014-09-17得票数 1

1回答

Intel Fortran矢量化:向量循环成本高于标量

fortran、vectorization、simd、intel-fortran

我正在使用Intel Fortran 15测试和优化遗留代码，我有一个简单的循环： G1(lev)%D(ir) = 0.d0LOOP BEGIN at MLFMATransl

浏览 3提问于2015-07-22得票数 0

回答已采纳

2回答

多处理具有多个参数的for循环，这些参数取决于循环索引

python、python-multiprocessing

你好，我想多处理一个for循环，在那里我计算4个不同的独立函数。唯一的问题是这些函数有多个参数。其中一些参数是我使用循环索引选择正确元素的列表。下面是一个示例代码 #do heavy calculationsd

浏览 0提问于2017-10-12得票数 0

1回答

Fortran宽度第一次搜索运行缓慢

performance、fortran、shortest-path、breadth-first-search

背景：I编写了一个Fortran子程序，该子程序从邻接矩阵计算无加权有向网络的最短路径长度矩阵。子程序采用宽度优先搜索算法。最短路径长度矩阵的对角线元素是最短的环路长度。但是，这个子程序在MATLAB中的运行速度略快于Fortran实现中的运行速度(在MATLAB中为44秒，在相同的1000个节点网络中为46秒)。这对我来说没有意义，因为我的理解是，对于基于循环的操作，Fortran</em

浏览 2提问于2015-04-13得票数 1

回答已采纳

2回答

Cuda Fortran* 4D数组*

arrays、cuda、fortran、pgi

我的代码被全局内存中我的4D数组访问减慢了。我正在访问的4D数组是只从设备读取的，并且大小在运行时是已知的。我的第一个想法是关于指针，但不熟悉指针fortran。因为这只会使我的代码</em

浏览 1提问于2013-09-23得票数 0

1回答

如何在DO循环中优化Fortran* IF-ELSE？*

performance、if-statement、optimization、fortran、do-loops

下面的函数在我的程序中用得最多。r_n=size(R,1); do i=1,4; ar=ABS(R(i)) sqrt(-7.0_dp+12.0_dp*ar-4.0_dp*ar*ar))*0.125_dp end doR, Del：长度为4的向量所以我想提高这个函数的速度。此外，它还处于<

浏览 0提问于2019-03-13得票数 1

1回答

Fortran中循环退出后索引变量的值

loops、indexing、fortran

在下面的do-循环中，索引i从1运行到n - 1。如果我在循环退出后直接打印i的值，那么我可以确认gfortran输出了n的值(在下面的情况下= 10 )。这种行为总是按照Fortran标准来保证吗？(也就是说，对于所有Fortran编译器/选项，我们都能期望这种行为吗？)program main integer :: i, n

浏览 1提问于2020-06-06得票数 2

回答已采纳

1回答

我不知道为什么循环中会更新一个变量数组

arrays、types、fortran

看看下面的代码，integer :: i, j real, dimensionhw1 如果您查看包含L_D的行，它基本上只是在根据Mi数组值计算其值之后更新L_D数组。问题是，Mi数组在第一次第i次循环后以某种方式进行了更新，即使我没有任何更新Mi数组本身的行。例如，当j=0时，我应该得到Mi=0.5，当j =1时，我应该得到Mi=0.6，依此类推。但是在第一次第i次<e

浏览 0提问于2018-02-27得票数 0

点击加载更多