[請益] FORTRAN do loop 在某問題上的寫法

作者: py94g (...)   2012-05-07 02:30:28
大家好,我想請教一個對我來說蠻困難的問題
我要讀一個檔案,檔名testyyy.dat,檔案內容是3行10列的數值。
檔案內容如下:
1 11 -0.1
2 12 -0.2
3 13 -0.3
4 14 -0.4
5 15 -0.5
6 16 -0.4
7 17 -0.3
8 18 -0.2
9 19 -0.1
10 20 -0.1
以下程式可以把檔案讀出並列在螢幕上:
Dimension POINT(3, N) ! 設定矩陣維度,3行N列,N=10。
10 FORMAT(1X, 3(3F10.5)) ! 輸出的檔案格式 Raw data。
C 讀檔案,把檔案中的數值填入建立POINT(3,N)的矩陣中。
OPEN(323, FILE='testyyy.dat') ! 把要讀入的檔案檔名+路徑設定好。
READ(323,*) POINT ! 讀入。
Do i = 1,N ! 從第一個數據開始讀到最後一個。
WRITE(*,10) POINT(1,i), POINT(2,i), POINT(3,i)
enddo ! 讀檔案結束。
接著,我要把內容分成兩個部分並輸出。
1 11 -0.1
2 12 -0.2
3 13 -0.3
4 14 -0.4
5 15 -0.5
輸出檔案a.dat
6 16 -0.4
7 17 -0.3
8 18 -0.2
9 19 -0.1
10 20 -0.1
輸出檔案b.dat
Do k = 1,5 ! 第一組1~5列
WRITE(*,10) POINT(1,k), POINT(2,k), POINT(3,k)
enddo
WRITE(*,*)
Do j = 6,10 ! 第二組6~10列
WRITE(*,10) POINT(1,j), POINT(2,j), POINT(3,j)
enddo
真正的問題是我的數據是3行1000000列,
依序由上到下分成1000組,每一組挑出一個極小值(這簡單可以自己做)。
如果依照上面的寫法是不行的,想了好幾天只好來這裡求助高手給我點提示。
感恩
作者: janhon (海牛)   2012-05-07 11:02:00
大概矩陣過大 記憶體或者編譯器不支援
作者: terryys (terry)   2012-05-07 11:40:00
你的問題是不知道怎麼樣把1M個資料分成1000個loop寫出來嗎另外這個矩陣其實不大,就算是real(8)也才24MB而已
作者: py94g (...)   2012-05-07 13:38:00
我用10列3行來練習,寫一組回圈分出兩組,結果怎麼寫都不對分別寫兩個loop j和k就OK,那分1000組豈不是寫1000次loop呢所以想請問大家是否有更好的寫法不必寫1000次loop分出1000組
作者: terryys (terry)   2012-05-07 17:24:00
其實雙層loop就可以解決你的問題了,不用寫1k個loop出來

Links booklink

Contact Us: admin [ a t ] ucptt.com