Fortran Coder

查看: 40|回复: 6

[数学库] fortran代码中dgemm函数调用的疑问

[复制链接]

4

帖子

1

主题

0

精华

新人

F 币
23 元
贡献
10 点
发表于 6 天前 | 显示全部楼层 |阅读模式
代码如下:
[Fortran] 纯文本查看 复制代码
        program TestMKLDgemm
        implicit none
        real*8  :: alpha, beta
        integer :: i, j
        integer, parameter :: m = 5, k = 5, n = 5
        real*8 :: A(m,k), B(k,n), C(m,n)
        !// C = alpha*A*B + beta*C
        alpha = 1.d0
        beta  = 0.d0

        A=5.d0

        B=2.d0

        C = 1.d0

        call DGEMM('N', 'N', M, N, K, alpha, A, M, B, K, beta, C, M)
        print*, B
        PRINT*
        print*, A
        PRINT*
         print*, C
         READ*
        end program TestMKLDgemm

生成解决方案没有问题,但是得出的结果如图,根本不对啊,请教是什么问题呀
aaaa.png
回复

使用道具 举报

4

帖子

1

主题

0

精华

新人

F 币
23 元
贡献
10 点
 楼主| 发表于 6 天前 | 显示全部楼层
添加的库是mkl_blas95_ilp64.lib, mkl_intel_lp64_dll.lib, mkl_intel_thread_dll.lib, mkl_core_dll.lib, libiomp5md.lib

1408

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
646 元
贡献
500 点

新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 6 天前 | 显示全部楼层
哪儿不对?
【2】【5】不就是 【50】 吗?2*5+2*5+2*5+2*5+2*5 = 50

4

帖子

1

主题

0

精华

新人

F 币
23 元
贡献
10 点
 楼主| 发表于 6 天前 | 显示全部楼层
fcode 发表于 2019-7-11 17:46
哪儿不对?
【2】【5】不就是 【50】 吗?2*5+2*5+2*5+2*5+2*5 = 50

对不起对不起,脑子懵了,还以为是10。。。

1408

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
646 元
贡献
500 点

新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 6 天前 | 显示全部楼层
建议你用95接口,调用更简单

use blas95
call GEMM(a, b, c, 'N','N', alpha, beta)

4

帖子

1

主题

0

精华

新人

F 币
23 元
贡献
10 点
 楼主| 发表于 6 天前 | 显示全部楼层
fcode 发表于 2019-7-11 18:23
建议你用95接口,调用更简单

use blas95

老师我想请教下,95接口的使用条件是什么呢?是跟代码的语法是否是for95有关系吗?

1408

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
646 元
贡献
500 点

新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 5 天前 | 显示全部楼层
条件就是你的编译器提供95接口,很显然ivf是满足这个条件的。

代码的语法是不是95没关系。而且,只要编译器支持,一段代码可以77,90,95,2003,2008,2015混合着用(这些规范之间本来也有很大的交集)。没有人能明确的指出一段代码到底是什么标准。
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|群聊|Fcode

GMT+8, 2019-7-17 06:54

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表