※ 引述《kuang01 (我叫小劉)》之銘言:
: 想請問各位大大
: 在fortran 一維矩陣裡
: do i=1,5
: A(i)=2
: enddo
: 這樣會是 一colum vector 還是 row vector
: 因為我在matlab裡
: 如果把 do 換成for
: 這樣去寫 會出現 row vector
: 想請問 fortran 的一維跟matlab是否不太一樣
底下的人回答大致正確,因為對指標上來說,確實沒差別,
但是你有這種問題,甚至刻意的去區別到底一個一維陣
列是row還是column vector也不是壞事.
因為在記憶體的使用上,fortran確實是by column,而C
based的語言,例如C/C++,Matlab,Scilab,Python都是
by row.也就是說,當你往記憶體裡面填入一個二維陣
列時,fortran會讓同一個column的資料作連續的排列.
因此當你在記憶體裡面填值或讀值的時候,在fortran應
該要先作column vector,而Matlab,C/C++,Python等應該要
先做row vector,在記憶體的使用上才是連續的讀寫,不會
一直跳躍.
這個影響還不只是讀寫效率的問題,那還只是其次的,作
不好頂多效率比較低,但是當你要讓by column跟by row
的語言相互傳值時,問題才出現.
例如你想在C-based的語言裡面使用fortran寫成的library,
例如lapack.因為你會把C裡面的一塊記憶體傳到fortran中,
但因為C是將同一個row的資料在記憶體中作連續排列,跟
fortran對記憶體的使用方式剛好相反,所以當你把資料傳入
時,fortran在理解這個陣列時,行列都會顛倒了.所以當你要
把C裡面的值傳到fortran時,往往需要先轉秩才能傳入.反之
亦然.