Fortran Coder

查看: 87|回复: 5

[数值问题] fortran二维数组如何赋值给四维数组???

[复制链接]

6

帖子

2

主题

0

精华

入门

F 币
50 元
贡献
22 点
发表于 2019-11-14 17:40:21 | 显示全部楼层 |阅读模式
各位大佬们,请问二维数组如何赋值给四维数组呢?我是这样写的,但出错了:
[Fortran] 纯文本查看 复制代码
do i=1,64
	do j=1,64
					
		B(i,j,1,1)=A11(i,j)
		B(i,j,1,2)=A12(i,j)
		B(i,j,2,1)=A21(i,j)
		B(i,j,2,2)=A22(i,j)
					
	enddo
enddo


劳烦各位大佬们帮忙看看,具体该如何实现呢?谢谢啦!
回复

使用道具 举报

745

帖子

0

主题

0

精华

宗师

F 币
913 元
贡献
506 点

规矩勋章元老勋章新人勋章水王勋章

QQ
发表于 2019-11-14 19:32:47 | 显示全部楼层
本帖最后由 vvt 于 2019-11-14 19:33 编辑

Fortran是矢量化的语言,所以很多操作不需要循环
如果数组的大小是一样的,都是64的话。
real :: B(64,64,2,2) , A11(64,64) , A12(64,64) , A21(64,64) , A22(64,64)

这么一句话就可以了。

B = reshape([A11,A12,A21,A22],shape(b))



6

帖子

2

主题

0

精华

入门

F 币
50 元
贡献
22 点
 楼主| 发表于 2019-11-15 10:47:32 | 显示全部楼层
vvt 发表于 2019-11-14 19:32
Fortran是矢量化的语言,所以很多操作不需要循环
如果数组的大小是一样的,都是64的话。
real :: B(64,64,2 ...

非常感谢,问题已解决!!

1

帖子

0

主题

0

精华

新人

F 币
10 元
贡献
3 点
发表于 2019-11-16 13:06:04 | 显示全部楼层
本帖最后由 GuihangShao 于 2019-11-16 13:11 编辑

按照楼主代码,在vs2012+ivf2013环境下并没有报错。

代码:

[Fortran] 纯文本查看 复制代码
    program ex_fcode
    implicit none
    integer :: i,j
    real*8 A11(64,64),A12(64,64),A21(64,64),A22(64,64),B(64,64,2,2)
    
    A11=0.0d0
    A12=0.0d0
    A21=0.0d0
    A22=0.0d0
    
    do i=1,64
        do j=1,64
            B(i,j,1,1)=A11(i,j)
            B(i,j,1,2)=A12(i,j)
            B(i,j,2,1)=A21(i,j)
            B(i,j,2,2)=A22(i,j)
        enddo
    enddo
    
    pause
    endprogram ex_fcode
    

   
结果截图:
结果截图.jpg
D:\MarineCode\20181022Test\20191116forFcode

6

帖子

2

主题

0

精华

入门

F 币
50 元
贡献
22 点
 楼主| 发表于 2019-11-20 11:24:27 | 显示全部楼层
GuihangShao 发表于 2019-11-16 13:06
按照楼主代码,在vs2012+ivf2013环境下并没有报错。

代码:

恩恩,非常感谢您,我后面也发现这样是可以的!!

43

帖子

1

主题

0

精华

熟手

F 币
277 元
贡献
129 点

规矩勋章

发表于 2019-11-26 11:03:52 | 显示全部楼层
我感觉下面这样比较直观
real*8 A11(64,64),A12(64,64),A21(64,64),A22(64,64),B(64,64,2,2)
B(:,:,1,1)=A11
B(:,:,1,2)=A12
B(:,:,2,1)=A21
B(:,:,2,2)=A22
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|群聊|Fcode

GMT+8, 2019-12-14 05:38

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表