存储多维阵列

Question

前言

的Fortran程序，我正在写应处理取决于ndims 1D，2D和3D的问题，它可以是1，2或3并从输入文件中读取。

在这些情况下的量/所关注IES可以（一个可以被命名为phi）存储在数组

秩 dims（ ALLOCATABLE(:)或 ALLOCATABLE(:,:)或 ALLOCATABLE(:,:,:)）的
1. ，
2. 或在等级阵列3（ALLOCATABLE(:,:,:)，第三维在2D中等于1或者第二和第三维在1D中等于1）;

这两种情况在this answer都有很好的解释。第一种方法对我来说似乎更加优雅，但在下面我假设第二种方法，这绝对简单。

这些量必须由几个子程序上进行操作（如mysub）沿ndims尺寸（沿“铅笔”应该给图形的想法），所以我应该叫像

SELECT CASE (ndims) 

! 3D case 
CASE (3) 
    DO j = ... 
    DO k = ... 
     CALL mysub(phi(:,j,k)) 
    END DO 
    END DO 
    DO i = ... 
    DO k = ... 
     CALL mysub(phi(i,:,k)) 
    END DO 
    END DO 
    DO i = ... 
    DO j = ... 
     CALL mysub(phi(i,j,:)) 
    END DO 
    END DO 

! 2D case 
CASE (2) 
    DO j = ... 
    DO k = ... 
     CALL mysub(phi(:,j,1)) 
    END DO 
    END DO 
    DO i = ... 
    DO k = ... 
     CALL mysub(phi(i,:,1)) 
    END DO 
    END DO 

! 1D case 
CASE (1) 
    DO j = ... 
    DO k = ... 
     CALL mysub(phi(:,1,1)) 
    END DO 
    END DO 
END SELECT 

实际问题

任何人都可以建议我（或帮助我设计！）一种不同的存储方式phi（可能涉及派生数据类型？），以便我可以折叠前面的代码如下？

DO id = 1, ndims 
    CALL mysub2(phi,id) 
END DO 

（这里​​作用mysub的地方。）

所以，问题是我应该如何存储披，这样我就可以用第二个替代第一码？

也许我可以回到前言，并决定按照1点，在这种情况下会更容易编写一个通用的接口。我认为，这只是一种“隐藏”SELECT CASE会做什么的方法。哪两个（SELECT CASE /通用INTERFACE）会更有效率？

这些是面对这个问题只有两个办法？

Answer 1

也许是我误会了，不过我想答案的具体问题是不能对存储或披声明所有的任何变化。

在原始代码中，三维数据切片中沿着第一维度被处理，则第二，然后第三个（从用于存储数据的阵列的等级区分的数据的等级）。二维数据沿着第一个处理，然后第二个和一维数据沿着第一个处理。

ID为1将在数据的维数

所以，请考虑以下实施​​：

SUBROUTINE mysub2(phi, id) 
    TYPE(pink_elephant), INTENT(IN) :: phi(:,:,:) 
    INTEGER, INTENT(IN) :: id 

    INTEGER :: i, j, k 

    SELECT CASE (id) 
    CASE (1) 
    DO j = ... 
     DO k = ... 
     CALL mysub(phi(:,j,k)) 
     END DO 
    END DO 

    CASE (2) 
    DO i = ... 
     DO k = ... 
     CALL mysub(phi(i,:,k)) 
     END DO 
    END DO 

    CASE (3) 
    DO i = ... 
     DO j = ... 
     CALL mysub(phi(i,j,:)) 
     END DO 
    END DO 

    END SELECT 
END SUBROUTINE mysub2 

~~

通用接口，可以随时在“编译时”来解决 - 可以通过查看代码中的声明来确定将由特定CALL语句或函数引用调用的特定过程（不是类型绑定）或绑定（类型绑定）。

如果您遇到了“运行时间”信息将影响过程选择的情况，那么必须有一些其他可执行机制，除了通用解析之外或附加解决方案 - 即使因此声明后，选择的情况下，动态调度，等，等，等

询问通用的分辨率是否为可执行的决定更有效的是没有什么特殊意义 - 它们是不同的东西。

Answer 2

你可能想是这样的：

program test 

    integer :: j,ndims 
    integer :: n ! rank of each dimension, could also be read from input an allocated separately 

    type arr 
     real(8) :: x(n) ! one array for each dimension 
    end type 

    type(arr),allocatable :: phi 

    read(*,*) ndims 
    allocate(phi(ndims)) 

    do j=1,ndims 
     call mysub(phi(j)%x) ! acts on the array in dimension j 
    end do 

contains 

    subroutine mysub(x) 
    ... 
    end subroutine 

end program