在matlab里面将计算复杂度高的函数写成c++是非常实用的,尤其是对于我这种对c++编译不是很理解的小白。以前每次百度,出来的都是一些非常简单的例子,根本不知道如何入手,发现最直接的方法就是看别人在github上分享的代码和matlab 的document. 所以我以一个不是那么简单的例子入手,看一些整个流程。
主要分为两部分(在matlab命令窗口(command window)下输入 doc mex 可以找到你需要的所有说明文件):
1.写(输入传入和传出,用matlab给出的接口函数,将数据传进cpp file并返回给matlab,数据一旦传入c file, 就得严格按照c++的规则进行编程);2.编译成matlab可以调用的函数.
下面是我从程序成截出来的主要部分:
在matlab里面调用
注意到,FV.faces 是一个N*3的矩阵,但是c++不支持,所以拆成了数列3个N*1的数列,同理FV.vertices. 我们这里有八个参数传入了polygan2voxel_double函数,接下来我们看一下c file是如何用mexFunction接入的。下图135行 “mex.h”包含我们需要借口函数,可以看到,mexFunction只有四个参数:
nlhs: 输出参数的个数(number),我们这里是1, 不需要我们在matlab里面给,它会自己计算,只是用来做check的函数(如下图的152-157行)
plhs:输出参数的指针(point),这里就是在c file中指向Volume的指针(不要害怕指针,我也挺不喜欢这个,就当c++的数组用),下图177-178告诉我们是如何创建的输出参数的。 “mexCreate**” 是创建函数,“mexGet**”是取值函数。去MATLAB的文档里面看细节,然后找你需要的函数。
nrhs:输入参数的个数, 我们这里是8, 不需要我们在matlab里面给,它会自己计算,只是用来做check的函数(如下图的152-157行)
prhs:输入参数的指针,在我们这里就是指向输入的八个参数,下图159-167告诉我们如何取出来
(notes: 这些参数不是缩写,所以不好记,但是还是要记住,mexFunction参数依次是输出(个数,指针)-->输入(个数,指针))
完成这些以后,就像写c++一样正常写就好了,然后就是编译 “mex polygan2voxel_double.c” (下面这些在一个名叫polygan2voxel_double.c文件里) 编译的时候可能会有很多问题,“help mex” 看一下可选参数。 我之前就碰到过gcc 版本问题,就是matlab mex默认版本和系统的gcc版本不匹配,后来是通过可选参数指定gcc 版本解决的。