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

条件Fortran循环的优化

blas lapack fortran optimization

sponce · 技术社区 · 7 年前

我想优化以下代码块的速度:

DO i=1, dim1
  DO j=1, dim2
    DO k=1, dim3
      IF (A(k,j,i)>0) &
        B(k,j,i) = exp(C(k))/A(k,j,i)
    ENDDO
  ENDDO
ENDDO

非常重要的是, A INTEGER B 和 C 是 COMPLEX

有两个问题: 1) 如何用BLAS/LAPACK调用来替换这个?问题是条件。 2) exp的计算很慢。如何加快速度?

2 回复 | 直到 7 年前

RussF 7 年前

我做了几个测试 idim[1-3] 各种排列的 [40,40,1000] 并发现使用一个临时数组作为指数并保持原始循环的顺序比提供的另一个答案快2倍或更多。你的理解可能因编译器等的不同而不同。

d=exp(c)
DO i=1, dim1
  DO j=1, dim2
    DO k=1, dim3
      IF (A(k,j,i)>0) &
        B(k,j,i) = d(k)/A(k,j,i)
    ENDDO
  ENDDO
ENDDO

Yves Daoust 7 年前

DO k=1, dim3
  expCk= exp(C(k))
  DO i=1, dim1
    DO j=1, dim2
      IF (A(k,j,i)>0) &
        B(k,j,i) = expCk/A(k,j,i)
    ENDDO
  ENDDO
ENDDO

推荐文章

Mr. W · 如何生成一个重复的整数,如二进制中的1001001001,时间复杂度为O(n)?

11 月前

giantjenga · 优化整数向量到二进制向量的转换

11 月前

Costinean Sebastian · 如何优化使用for循环处理大型数组的MATLAB脚本的执行时间?"

11 月前

Daniel Lobo · 使用约束进行优化

11 月前

Imyaf · 在什么条件下,在组装时对“if”和“if-else”进行比较?

1 年前

Sergio · python中大量数字的乘法

1 年前

Kiran Raj R · 为什么尽管使用了useMemo和useCallback,我的React组件却意外地重新渲染?

1 年前

Sergey Dev · 临时表与表变量

1 年前

JoReyner · 单字符元素的空间差异:列表<字符>与列表<字符串>

1 年前

John · 减少C中的内存消耗++

1 年前