[問題] array 當argument 的問題

[yhd (行的正坐的直，不換ID)]

舉個簡單的例子，二個向量交換

program main
    implicit none

    real :: a(2), b(2)

    a(1) = 1.0E0
    a(2) = 2.0E0
    b(1) = 3.0E0
    b(2) = 4.0E0
    call swap(a, b)

    write(*, *) "a(1)=", a(1), "a(2)=", a(2)
    write(*, *) "b(1)=", b(1), "b(2)=", b(2)
end program main

SUBROUTINE swap(a, b)
    REAL, DIMENSION(:)       :: a, b
    REAL, DIMENSION(SIZE(a)) :: work
    work = a
    a = b
    b = work
END SUBROUTINE swap

可以編譯成功，但執行時卻會發生segmentation fault，
把swap 裡的
work = a
a = b
b = work
註解掉，然後加了
write(*,*) size(a)

看size(a) 的值，
可以正常執行，但size(a) 的結果卻是0

我用的compiler 是ifort 11.1

請問有誰知道這是怎麼回事嗎？謝謝 :-)

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 140.109.176.53
→ yhd:我發現問題似乎是a, b 這二個array 的傳遞失敗，如果在swap 裡
→ yhd:明確指定大小，i.e., REAL, DIMENSION(2) :: a, b
→ yhd:但這樣就失去希望work 為一個automatic array 的用意了
→ awer89:不過主程式宣告a b 2陣列的大小  不也限制了work的大小了?
→ yhd:沒有吧，work 會隨著輸入的a 的大小更改
→ awer89:不過這樣要再改一次主程式不是嗎?a b大小都宣告了
→ yhd:是呀，重點在這樣寫subroutine 就不用再傳二個向量的長度這個
→ yhd:引數了。我是看mkl 的blas95 interfaces 這麼寫，想如法炮製
推 charlesdc:你把副程式包成模組的方式使用就沒問題了:)
推 motoman:應該是沒寫interface吧...不過包成module就不用寫了
→ yhd:可以了，謝謝大家，看來我對interface 不熟，來好好研究一下XD