将参数从 Fortran 程序传递给 Fortran 函数

您可以从 Fortran MEX 文件中调用 LAPACK 和 BLAS 函数。以下示例使用两个矩阵，并通过调用 BLAS 例程 dgemm 将这两个矩阵相乘。要运行该示例，请将代码复制到编辑器中并将文件命名为 calldgemm.F。

#include "fintrf.h"

subroutine mexFunction(nlhs, plhs, nrhs, prhs)

mwPointer plhs(*), prhs(*)

integer nlhs, nrhs

mwPointer mxcreatedoublematrix

mwPointer mxgetpr

mwPointer A, B, C

mwSize mxgetm, mxgetn

mwSignedIndex m, n, p

mwSize numel

double precision one, zero, ar, br

character ch1, ch2

ch1 = 'N'

ch2 = 'N'

one = 1.0

zero = 0.0

A = mxgetpr(prhs(1))

B = mxgetpr(prhs(2))

m = mxgetm(prhs(1))

p = mxgetn(prhs(1))

n = mxgetn(prhs(2))

plhs(1) = mxcreatedoublematrix(m, n, 0.0)

C = mxgetpr(plhs(1))

numel = 1

call mxcopyptrtoreal8(A, ar, numel)

call mxcopyptrtoreal8(B, br, numel)

call dgemm(ch1, ch2, m, n, p, one, %val(A), m,

+ %val(B), p, zero, %val(C), m)

return

end

链接到 BLAS 库，其中包含 dgemm 函数。

mex -v -R2017b calldgemm.F -lmwblas