一个数组运算的问题

andy8496 · 发表于 2017-8-4 18:27:18

遇到一个数组片段的问题。

代码1：
do i=1,size(ia_index,2)
  ra_val(ia_index(1,i)) = ra_val(ia_index(1,i)) + ra_add(i)
  ra_val(ia_index(2,i)) = ra_val(ia_index(2,i)) + ra_add(i)
enddo

代码2：
  ra_val(ia_index(1,:)) = ra_val(ia_index(1,:)) + ra_add
  ra_val(ia_index(2,:)) = ra_val(ia_index(2,:)) + ra_add

我一直以为代码1和2得到的结果是一样的，今天遇到一个问题发现并不一样，代码1的结果是正确的，当ia_index(1,:)表示的某个下标索引出现1-2次的时候，两段代码结果一致，当出现3次的时候，代码1可以累加，代码2中ra_val索引指向的成员就只有最后一个值，没有累加。不知何故！

事实上，我最初使用代码2之前，由于不放心，是做过实验的：

[Fortran] 纯文本查看 复制代码

1program main
2 implicit none
3   integer :: ia(5),ib(9)
4   ia = [1,5,5,5,2]
5   ib = 0
6   ib(ia) = ib(ia) + [1,1,1,1,1]
7   write(*,*) ib
8   stop
9 end program

结果见附件图片，无误，能够叠加我才这么写的，真是郁闷了...
还请大伙指点，先谢谢了！

fcode · 发表于 2017-8-5 04:48:18

It's hard to check why ? May be we need more code !

andy8496 · 发表于 2017-8-5 10:21:45

我编了一组数据：

[Fortran] 纯文本查看 复制代码

01program main
02implicit none
03  real :: ra_val1(11),ra_val2(11),ra_add(9)
04  integer :: ia_index(2,9)
05   
06  integer :: i
07    
08  ra_add = [0.014000,0.016500,0.016500,0.014000,0.014000,0.014750,0.014750,0.010750,0.010750]
09   
10  ia_index(:,1) = [ 1  , 2  ]
11  ia_index(:,2) = [ 1  , 9  ]
12  ia_index(:,3) = [ 1  , 11 ]
13  ia_index(:,4) = [ 2  , 3  ]
14  ia_index(:,5) = [ 3  , 4  ]
15  ia_index(:,6) = [ 4  , 8  ]
16  ia_index(:,7) = [ 4  , 10 ]
17  ia_index(:,8) = [ 5  , 6  ]
18  ia_index(:,9) = [ 6  , 7  ]
19   
20 
21   
22  ra_val1 = 0.
23  ra_val2 = 0.
24   
25  ! 写法1
26  do i=1,9
27    ra_val1(ia_index(1,i)) = ra_val1(ia_index(1,i)) + ra_add(i)
28    ra_val1(ia_index(2,i)) = ra_val1(ia_index(2,i)) + ra_add(i)
29  enddo
30   
31  ! 写法2
32  ra_val2(ia_index(1,:)) = ra_val2(ia_index(1,:)) + ra_add
33  ra_val2(ia_index(2,:)) = ra_val2(ia_index(2,:)) + ra_add
34   
35   
36  write(*,"(A6,2A12)") "No.","Val1","Val2"
37  do i=1,11
38    write(*,"(I6,2F12.6)") i,ra_val1(i),ra_val2(i)
39  enddo
40   
41  stop
42end program

andy8496 · 发表于 2017-8-5 10:23:01

本帖最后由 andy8496 于 2017-8-5 10:27 编辑

fcode 发表于 2017-8-5 04:48
It's hard to check why ? May be we need more code !

上面的代码为什么写法1和2结果不同？

理解不了，我认为结果应该是一样的，都是写法1的结果。还请老大赐教！

维尼猴 · 发表于 2017-8-5 15:52:06

本帖最后由维尼猴于 2017-8-5 15:58 编辑

哇塞，竟然能这样操作，但感觉可能和缓存的处理有关系吧

我试了一下，比如加个中间变量 temp

[Fortran] 纯文本查看 复制代码

`1`	`temp` `=` `ia_index(1,:)`

`2`	`ra_val3(temp` `)` `=` `ra_val3(temp` `)` `+` `ra_add`

这样处理之后是可以累加的，可能是你数组声明成  2,9 的问题，内存不够连续？？
来我试试  9,2 的

诶………………我试了一下，失败了，把  ia_index 交换了维数也没有得到对应效果

维尼猴 · 发表于 2017-8-5 16:00:59

我补充测试了一下，似乎下标数组不能取切片

[Fortran] 纯文本查看 复制代码

1program main
2implicit none
3  integer :: ia(5),ib(9)
4  ia = [1,5,5,5,2]
5  ib = 0
6  ib(ia) = ib(    ia  (:)   ) + [1,1,1,1,1]
7  write(*,*) ib
8  stop
9end program

只是在你原来例子上给ia加了个切片操作，就不能累加了

andy8496 · 发表于 2017-8-5 21:29:22

唉，我都忘了我哪些地方用了这种写法，但是肯定不止目前发现的这个位置。真是郁闷了！
要是能通过什么编译选项能设置解决就好了！
谢了！

维尼猴 · 发表于 2017-8-6 16:07:59

andy8496 发表于 2017-8-5 21:29
唉，我都忘了我哪些地方用了这种写法，但是肯定不止目前发现的这个位置。真是郁闷了！
要是能通过什么编译 ...

我感觉这种有些“讨巧”的语法，还是少用心里更安一些呢

pasuka · 发表于 2017-8-6 17:22:43

ia_index的第1行数据为啥非要设置重复元素？

andy8496 · 发表于 2017-8-6 18:30:55

pasuka 发表于 2017-8-6 17:22
ia_index的第1行数据为啥非要设置重复元素？

这个是实际应用的需要。

[讨论] 一个数组运算的问题

美女勋章

热心勋章

星光勋章

新人勋章

贡献勋章

管理勋章

帅哥勋章

爱心勋章

规矩勋章

元老勋章

水王勋章