Fortran ASSOCIATE 语法允许索引访问吗？

Question

有没有一种很好的方法来编写 Fortran ASSOCIATE 语句，将其变成

FORALL (i = 2:n-2)
    v(:,i) = v(:,i) + MATMUL(A, &
             c(2)*u(:,i-2) + c(1)*u(:,i-1) + c(0)*u(:,i) + c(1)*u(:,i+1) + c(2)*u(:,i+2))
END FORALL

如下所示的内容

ASSOCIATE ( U => ..., V => ...)
    FORALL (i = 2:n-2)
        V(i) = V(i) + MATMUL(A, &
               c(2)*U(i-2) + c(1)*U(i-1) + c(0)*U(i) + c(1)*U(i+1) + c(2)*U(i+2))
    END FORALL
END ASSOCIATE

我盯着 Adams 等人的 The Fortran 2003 Handbook 第 8.2 节，但我不知道如何编写关联名称 =>选择器 构造以允许对associate-name 进行索引访问。

显然，我想要的几行内容是多余的。我有一堆我想浓缩的东西。

Answer 1

除非我误读了事情，否则我认为这是不可能的。 规范说（第 8.1.4.3 节）：

在 SELECT TYPE 或 ASSOCIATE 构造中，每个关联实体与其所属的实体具有相同的等级
关联的选择器。

据我所知，您需要一个 1 级关联实体 (V)，并且需要一个 2 级关联实体
关联的选择器（用于保存v）。

Answer 2

如果目标是使代码更短/更好，我认为实现这一目标的最佳方法不需要 ASSOCIATE 构造：

forall(i=2:n-2)
  v(:,i) = v(:,i) + MATMUL(A,MATMUL(u(:,i-2:i+2),c([2,1,0,1,2]))
end forall

只要 5 大小的系数数组不改变，就可以预初始化为 a

real :: c5(5)
c5 = c([2,1,0,1,2])

然后在单行中运行，

forall(i=2:n-2) v(:,i) = v(:,i) + MATMUL(A,MATMUL(u(:,i-2:i+2),c5))

注意它必须是 lbound(u,2)<=0。