存儲多維陣列

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