|
上一节文章提及了mex文件的编译指令，主要是要确定openmp文件的路径，难度不大，其实高阶的混编是MATLAB+CUDA+C/C++,但是这个部分设计的问题过于复杂且多数情况涉及异构编程，题主也不是很清晰，故暂时不涉及CUDA相关的内容
MATLAB矩阵在C中的索引方式（基础）
想要搞定混编，必须明确多个编程语言的数组索引规律，这里直接用MATLAB和C的索引方式进行分析，Python和C混编原理完全相同
以MATLAB二维矩阵为例：
A=[1,2;3,4]
A(2)=2，A(3)=3
可以看出索引是逐列索引的
以三维矩阵为例
A=[1,2;3,4]
B=[5,6;7,8]
C=zeros(:,:,2)
C(:,:,1)=A
C(:,:,2)=B
C(1)=1;C(2)=3;C(3)=2;C(4)=4;
C(5)=5;C(6)=7;C(7)=6;C(8)=8;
对于上面2x2x2的三维矩阵C有如下索引公式
C(i,j,k)=C((k-1)MN+(j-1)M+i)
M为矩阵行数
N为矩阵列数
Q为矩阵第三维度大小
例如
C(2,2,2)=C(4+2+2)=C(8)
C(1,1,2)=C(4+0+1)=C(5)
C(2,1,2)=C(4+0+2)=C(6)
本节非常基础但是当MATLAB矩阵传递给C后，不论矩阵多少维度，多复杂都遵循如下公式
C_matlab(i,j,k)=C_c((k-1)MN+(j-1)M+i)
M和N分别为矩阵第一维度和第二维度大小
下一节将给出基础的案例代码
|