Fortran Coder

标题: Fortran中传递数组，将三维数组传递到二维数组 [打印本页]

---

作者: 伟大无敌sl    时间: 2022-8-3 23:54
标题: Fortran中传递数组，将三维数组传递到二维数组
Fortran中传递数组，将三维数组传递到二维数组

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

program ex0819

implicit none

integer,parameter::dim1=2

integer,parameter::dim2=2

integer,parameter::dim3=2

integer::a(dim1,dim2,dim3)

a(:,:,1)=1

a(:,:,2)=2

call getarray1(a(:,:,1),dim1,dim2)

call getarray2(a(1,1,2),dim1)

stop

end program ex0819



subroutine getarray1(a,dim1,dim2)

implicit none

integer::dim1,dim2

integer::a(dim1,dim2)

write(*,*)a

return

end subroutine getarray1



subroutine getarray2(a,dim1)

implicit none

integer::dim1,dim2

integer::a(dim1,*)！最后一维可以不赋值数组大小

integer::i

write(*,*)(a(:,i),i=1,2)

return

end subroutine getarray2

第一个子程序
为什么出来的结果会是1，1，1，1，第二个子程序为什么出来的结果会是2，2，2，2？？？

（彭国伦Fortran95程序设计，p182页例子)

---

作者: 青衣巷    时间: 2022-8-4 09:01


a  :  2x2x2 三维数组。
a(:,:,1)   :   是 2x2 的二维数组，是 a 的前半部分。内容全部等于1
a(:,:,2)   :   也是 2x2 的二维数组，是 a 的后半部分。内容全部等于2

call getarray1(a(:,:,1),dim1,dim2)  这种是推荐使用的方式。没什么好说的。
call getarray2(a(1,1,2),dim1)  这种是不建议的使用方式，容易出错。不易于理解。

a 在内存里是这样存储的
a(1,1,1),a(2,1,1),a(1,2,1),a(2,2,1),a(1,1,2),a(2,1,2),a(1,2,2),a(2,2,2)
1       ,1       ,1       ,1       ,2       ,2       ,2       ,2

call getarray2(a(1,1,2),dim1) 把a(1,1,2) 传递给getarray2函数。
getarray2看到的部分是蓝色的部分：
a(1,1,1),a(2,1,1),a(1,2,1),a(2,2,1),a(1,1,2),a(2,1,2),a(1,2,2),a(2,2,2)
1       ,1       ,1       ,1       ,2       ,2       ,2       ,2      
因此对于getarray2来说，其值是4个2。

这种方式非常不好，不推荐。

参考：http://passarray.w.fcode.cn/

现代Fortran代码，传递数组全部用假定形状+module就好了。

---

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)

Powered by Discuz! X3.2