fortran二维数组如何赋值给四维数组？？？

放耗子的绅

各位大佬们，请问二维数组如何赋值给四维数组呢？我是这样写的，但出错了：

[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

劳烦各位大佬们帮忙看看，具体该如何实现呢？谢谢啦！

vvt

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))

放耗子的绅

vvt 发表于 2019-11-14 19:32
Fortran是矢量化的语言，所以很多操作不需要循环
如果数组的大小是一样的，都是64的话。
real :: B(64,64,2 ...

非常感谢，问题已解决！！

GuihangShao

按照楼主代码，在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
    

   
结果截图：

D:\MarineCode\20181022Test\20191116forFcode

放耗子的绅

GuihangShao 发表于 2019-11-16 13:06
按照楼主代码，在vs2012+ivf2013环境下并没有报错。

代码：

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

necrohan

我感觉下面这样比较直观
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