看不懂别人用blas和lapack的某种写法……

liudy02 · 发表于 2020-3-13 09:02:05

看到人家写的程序中有如下语句：
call dgemm('T', 'N', nfeat2(itype), num(itype), nfeat0(itype), 1.d0, pv(1,1,itype), nfeat0m, feat_type(1,1,itype), nfeat0m, 0.d0, feat2_type(1,1,itype), nfeat2m)

完全不理解这里边像pv(1,1,itype)这种东西是啥子用法……
是和c里边一样只是作为一个地址指针用的么？

liudy02 · 发表于 2020-3-13 10:53:59

kyra 发表于 2020-3-13 09:10
如果pv是一个函数，则表示调用pv函数，其中 1 , 1 , itype 是函数的参数。并把 pv 函数的返回值传递给dgemm ...

pv是一个矩阵无疑
但是不可能是通常的指其某个元素的意义
因为从dgemm这个函数的参数列表来说，这个地方应该是一个矩阵

liudy02 · 发表于 2020-3-13 22:29:01

li913 发表于 2020-3-13 15:59
查mkl帮助可知，此参数应为二维数组。这里给出的用法比较老旧，不符合现在的标准，但由于fortran默认地址传 ...

嗯，后来想明白了，这用法和C对数组的处理类似
觉得这个方法很不fortran，不是个好办法……

[子程序] 看不懂别人用blas和lapack的某种写法……