上一节文章提及了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分别为矩阵第一维度和第二维度大小
MATLAB+C+Openmp混合编程(2)-MATLAB矩阵索引规则
于 2022-02-01 22:04:21 首次发布