Fortran Coder

标题: fortran二维数组如何赋值给四维数组？？？ [打印本页]

作者: 放耗子的绅 时间: 2019-11-14 17:40
标题: fortran二维数组如何赋值给四维数组？？？
各位大佬们，请问二维数组如何赋值给四维数组呢？我是这样写的，但出错了：

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

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
本帖最后由 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))

作者: 放耗子的绅 时间: 2019-11-15 10:47

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

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

作者: GuihangShao 时间: 2019-11-16 13:06
本帖最后由 GuihangShao 于 2019-11-16 13:11 编辑

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

代码：

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode



    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

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

代码：

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

作者: necrohan 时间: 2019-11-26 11:03
我感觉下面这样比较直观
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

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