matlab和fortran混合编程

  1 #include "fintrf.h" !必须有的头文件,里面有mxGetPr, mxGetM, mxGetN,mxCreateDoubleMatrix等函数的申明 http://yuehongliang2000.lofter.com/post/1d001431_ab72071
  2 
  3 
  4       subroutine mexFunction(OutSum,OutVar,InSum,InVar)!函数接口名称必须为mexFunction,
  5 
  6 
  7     !OutSum:输出参数个数
  8 
  9 
 10     !OutVar:输出参数数组指针
 11 
 12 
 13     !InSum:输入参数个数
 14 
 15 
 16     !InVar:输入参数数组指针
 17 
 18 
 19     !参数顺序不能随意更改
 20       Integer InSum,OutSum
 21 
 22 
 23       mwPointer InVar(*),OutVar(*)                           !mwPointer专门用于表示指针变量,这个不能随意用Integer代替
 24 
 25 
 26       mwPointer mxGetPr, mxGetM, mxGetN,mxCreateDoubleMatrix !这个对返回指针函数的再次申明,
 27       integer, parameter :: fp = selected_real_kind(15,300)
 28       Integer , parameter :: myINT  = SELECTED_INT_KIND( 8 )
 29 
 30 
 31       Real(fp),Allocatable::x(:,:),y(:,:)
 32       real(fp)::z
 33 
 34 
 35       Integer m,n
 36 
 37 
 38     
 39 
 40 
 41       If(InSum/=2)Then
 42 
 43 
 44 call mexErrMsgIdAndTxt('MATLAB:InputTooBig','输入参数个数必须为2个')
 45 
 46 
 47       Return
 48 
 49 
 50       EndIf
 51 
 52 
 53     
 54 
 55 
 56       m=mxGetM(InVar(1))!获取第1个输入参数的行数
 57 
 58 
 59       n=mxGetN(InVar(1))!获取第1个输入参数的列数
 60 
 61 
 62       Allocate(x(m,n),y(m,n))
 63 
 64 
 65       Call mxCopyPtrToReal8(mxGetPr(InVar(1)),x,m*n)!将第1个参数数组赋值给x变量
 66 
 67 
 68       Call mxCopyPtrToReal8(mxGetPr(InVar(2)),z,1)!将第2个整数变量赋值给z
 69 
 70 
 71 
 72 
 73       Call Add(x,y,z,m,n)!调用内部函数
 74 
 75 
 76     
 77 
 78 
 79       OutVar(1)=mxCreateDoubleMatrix(m,n,0)!给返回参数分配内存
 80 
 81 
 82       Call mxCopyReal8ToPtr(y,mxGetPr(OutVar(1)),m*n)!将返回参数赋值给分配的内存
 83 
 84 
 85     
 86 
 87 
 88       DeAllocate(x,y)!释放临时分配的内存
 89 
 90 
 91       Return
 92 
 93 
 94       End SubRoutine
 95 
 96 
 97 
 98 
 99     !  SubRoutine Add(x,y,z,m,n)
100     !
101     !  Integer,Intent(In)::m,n
102     !  integer, parameter :: fp = selected_real_kind(15,300)
103     !
104     !  Real(fp),Intent(In)::x(m,n),z
105     !
106     !  Real(fp),Intent(Out)::y(m,n)
107     !
108     !  Integer i,j
109     !
110     !  Do i=1,m
111     !
112     !      Do j=1,n
113     !
114     !            y(i,j)=x(i,j)+z
115     !
116     !      EndDo
117     !
118     !  EndDo
119     !
120     !  Return
121     !
122     !End SubRoutine
123     
124         SubRoutine Add(x,y,z,m,n)
125 
126 
127         Integer,Intent(In)::m,n
128         integer, parameter :: fp = selected_real_kind(15,300)
129 
130 
131         Real(fp),Intent(In)::x(m,n),z
132 
133 
134         Real(fp),Intent(Out)::y(m,n)
135 
136 
137         Integer i,j
138 
139 
140         Do i=1,m
141 
142 
143         Do j=1,n
144 
145 
146         y(i,j)=x(i,j)+z
147 
148 
149         EndDo
150 
151 
152         EndDo
153 
154 
155         Return
156 
157 
158         End SubRoutine