如何通過Lapack包爲Fortran乘以兩個矩陣?我在Ubuntu中使用gfortran編譯器。 我的代碼不工作是:使用lapack的dgemm(鏈接器錯誤)
program main
integer, parameter :: n = 10
double precision :: alpha = 1.0, beta = 0.0
real, dimension(10,10) :: a
do i1 = 1,n
do j1 = 1,n
a(i1,j1) = j1 + (i1-1)*n
end do
end do
call cpu_time(start)
call DGEMM('N', 'N', n, n, n, alpha, a, n, a, n, beta, a, n)
call cpu_time(end)
print *, end - start
end program main
我用:
gfortran 0.f90 -llapack
,它返回:
/tmp/ccPy78g5.o: In function `MAIN__':
0.f90:(.text+0x110): undefined reference to `dgemm_'
collect2: ld returned 1 exit status
你看過手冊嗎?什麼樣的矩陣? – Rook
我添加了我使用過的代碼並不起作用。 – MOON
確切的問題是什麼?一個簡單的「它不起作用」並不能幫助你。 – fuesika