[問題] 傳入副程式陣列維度不同? 懸賞1000P...

[silverair (那真是太諷刺了紹安)]

最近在把fortran90改成C
但本身不是學f90，網路上資料能找到的也有限...

最近碰到一個奇怪的問題，
CALL GATEOPEN(QU, HO, INGT(IST2, 1), GT(1, 1, IST2), MDC_GT);

這邊的GT應該是一個三維陣列中取出來的REAL值



但是在GATEOPEN這個subroutine中

SUBROUTINE GATEOPEN( T,F,N,FT,MAXD )
        .....
        INTEGER, INTENT(IN)   :: N,MAXD
        REAL(SP), INTENT(IN)  :: T,FT(MAXD,2)
        ...

我實在是看不懂為何這邊會變成FT(MAXD,2)
後面使用也是變成二維陣列的形式

        DO I = 1, N
                IF( T <= FT(I,1) ) THEN
                        F = FT(I,2)
                        RETURN
                END IF
        END DO


求解...懸賞1000P
感謝各位大神。

--
我的意中人是個蓋世英雄，有一天他會踩著七色的雲彩來娶我

                                             我猜中了前頭，可是我猜不著這結局。

--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.124.44.60 (臺灣)
※ 文章網址: https://www.ptt.cc/bbs/Fortran/M.1574677143.A.D89.html
推 fragmentwing: 老實說我看不太懂你的程式碼
這是project在用的軟體的其中一部分，我的任務是轉成C
所以基本上這個Code應該是沒問題
推 fragmentwing: 1不過我記得是有很奇耙的技巧
→ fragmentwing: 會故意少輸入東西來跑
推 fragmentwing: 但問題是 我記得那個案例是用項在順序上跑的差異
→ fragmentwing: 來進行操作 而不是針對同一項內的陣列進行操作
→ fragmentwing: 還有 你的N怎麼從陣列變成整數了
→ fragmentwing: 這不是不會給過嗎?
→ fragmentwing: 它不是對到INGT(IST2, 1)這一項嗎?
我的理解是INGT(IST2, 1)
應該等於C語言的INGT[IST2][0]吧

INGT是一個儲存int的二維陣列
已經在其他地方宣告過了@@
所以進到副程式應該是整數沒錯啊

GT則是一個儲存real的三維陣列，
GT(1, 1, IST2)取出來應該是一個real數，可是進到副程式他怎麼宣告成陣列
推 fragmentwing: 阿哩 好像找到了 我拍一下
推 fragmentwing: http://i.imgur.com/8oUs7kH.jpg
→ fragmentwing: http://i.imgur.com/3lTyH84.jpg
→ fragmentwing: http://i.imgur.com/6bhNAPF.jpg
→ fragmentwing: 好像是這個?雖然書上說很常用，不過我沒用過，這樣
→ fragmentwing: 交接起來會有麻煩吧w
→ fragmentwing: 我都手動寫程式碼來降維的
感謝大大
推 youlieme: fortran傳的只是reference 裡面的資料排列如果清楚
→ youlieme: 陣列不需要一樣大 他會從GT(1,1,ist2)這個位址開始
→ youlieme: 往後取MAXD*2的記憶體當作FT
→ youlieme: 我不知道你的GT是什麼資料 但看起來會是成對或可能
→ youlieme: 是兩兩一組 subroutine只判斷奇數值是否=>T而已
好的我研究一下，看起來可能是這個問題
→ blc: fortran的陣列在記憶體中是連續的，所以可以從陣列第一個
→ blc: 元素跟陣列大小得到整個陣列。就像C的指標那樣。
→ blc: 優點是會快很多(以前)，現在編譯器最佳化不知道
原來如此，感謝大大
※ 編輯: silverair (220.141.70.95 臺灣), 11/26/2019 22:49:06
推 fragmentwing: 收到推
→ fragmentwing: 另外我果然還需要精進阿......
→ blc: 晚了，我也有收到，感謝。