回 帖 发 新 帖 刷新版面

主题:大牛们,帮帮忙!谢谢,  调用mkl blas出错

如下代码,调用DGEMM时,提示Program Exception - access violation错误,请问,怎么回事?
(IVF 11.1.038自带的mkl)
谢谢!!!
program test_blas
implicit none
  real(kind=8):: A(200, 200), B(200, 200), C(200, 200)
  real(kind = 8) :: time_begin, time_end
  
  CALL RANDOM(A)
  CALL RANDOM(B)

  CALL CPU_TIME(time_begin)
    C=matmul(A, B)
  CALL CPU_TIME(time_end)
  WRITE(*,*)"consumed CPU_time(s):", time_end - time_begin

  CALL CPU_TIME(time_begin)
  CALL dgemm('N', 'N', 2000, 2000, 2000, 1.0, A, 2000, B, 2000, 0.0, C, 2000)
  CALL CPU_TIME(time_end)
  WRITE(*,*)"consumed CPU_time(s):", time_end - time_begin
end program


错误提示:

consumed CPU_time(s):  0.109375000000000

forrtl: severe (157): Program Exception - access violation
Image              PC        Routine            Line        Source
mkl_core.dll       005E4D61  Unknown               Unknown  Unknown

回复列表 (共2个回复)

沙发

大哥们,帮帮忙

板凳

你是用这个帖子的程序算的吧?
http://bbs.pfan.cn/post-320160.html

既然你把矩阵的大小A, B改成200, 那么下面调用的函数2000也要相应作出点修改啊.

我来回复

您尚未登录,请登录后再回复。点此登录或注册