放耗子的绅 发表于 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

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

vvt 发表于 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(,shape(b))



放耗子的绅 发表于 2019-11-15 10:47:32

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

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

GuihangShao 发表于 2019-11-16 13:06:04

本帖最后由 GuihangShao 于 2019-11-16 13:11 编辑

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

代码：


    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

放耗子的绅 发表于 2019-11-20 11:24:27

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

代码：


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

necrohan 发表于 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
页: [1]
查看完整版本: fortran二维数组如何赋值给四维数组？？？