如何在Matlab中分析MEX函数

时间:2012-06-27 05:59:01

标签: c++ matlab profiling mex

我有一个Mex函数(你可以从Matlab调用的c ++函数),我已经编写过了,我想用valgrind / kcachegrind对它进行分析。如果你直接运行一个c ++程序,我知道如何使用valgrind / kcachegrind,但如果我从Matlab调用c ++程序,有没有办法进行这种分析?

2 个答案:

答案 0 :(得分:10)

分析MEX文件很棘手,因为MEX文件是共享库。使用标准的'gprof'方法无法在Linux上完成 - gprof根本不这样做。我尝试使用sprof,但我得到“PLTREL not found error” - sprof也不能使用。有一篇帖子here,但没有人给出最终答案。

幸运的是,有一种方法可以在 Linux 上使用valgrind来实现。首先,我们需要编写“运行”代码来加载mex文件,提供mexFunction符号供我们调用,并设置MEX文件的参数。我选择使用推荐的方法使用MATLAB执行此操作 - 使用MATLAB engine。以下代码(另存为test.c)加载MEX文件并查找mexFunction符号,从先前保存为'input.mat'的文件中加载输入数据(可以使用save命令在MATLAB中完成),并调用mexFunction。

#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <errno.h>
#include <dlfcn.h>
#include "engine.h"

typedef void (*mexFunction_t)(int nargout, mxArray *pargout [ ], int nargin, const mxArray *pargin[]);

int main(int argc, const char *argv[])

{
  Engine *ep;
  char buff[1024];
  int i;

  /* matlab must be in the PATH! */
  if (!(ep = engOpen("matlab -nodisplay"))) {
    fprintf(stderr, "Can't start MATLAB engine\n");
    return -1;
  }
  engOutputBuffer(ep, buff, 1023);

  /* load the mex file */
  if(argc<2){
    fprintf(stderr, "Error. Give full path to the MEX file as input parameter.\n");
    return -1;
  }
  void *handle = dlopen(argv[1], RTLD_NOW);
  if(!handle){
    fprintf(stderr, "Error loading MEX file: %s\n", strerror(errno));
    return -1;
  }

  /* grab mexFunction handle */
  mexFunction_t mexfunction = (mexFunction_t)dlsym(handle, "mexFunction");
  if(!mexfunction){
    fprintf(stderr, "MEX file does not contain mexFunction\n");
    return -1;
  }

  /* load input data - for convenience do that using MATLAB engine */
  /* NOTE: parameters are MEX-file specific, so one has to modify this*/
  /* to fit particular needs */
  engEvalString(ep, "load input.mat");
  mxArray *arg1 = engGetVariable(ep, "Ain");
  mxArray *arg2 = engGetVariable(ep, "opts");
  mxArray *pargout[1] = {0};
  const mxArray *pargin[2] = {arg1, arg2};

  /* execute the mex function */
  mexfunction(1, pargout, 2, pargin);

  /* print the results using MATLAB engine */
  engPutVariable(ep, "result", pargout[0]);
  engEvalString(ep, "result");
  printf("%s\n", buff);

  /* cleanup */
  mxDestroyArray(pargout[0]);
  engEvalString(ep, "clear all;");
  dlclose(handle);
  engClose(ep);

  return 0;
}

MEX文件本身也应使用mex -g开关编译。上面的代码必须使用mex -g编译并使用engopts.sh作为编译参数。从MATLAB命令行类型

mex('-v', '-f', fullfile(matlabroot,...
    'bin','engopts.sh'),...
    'test.c');

或在标准的Linux终端运行

/path/to/matlab/bin/mex -g -f /path/to/matlab/bin/engopts.sh test.c

使用valgrind分析MEX文件需要从命令行运行'test'程序。在test和MEX文件都驻留的目录中键入命令:

PATH=$PATH:/path/to/matlab/bin/ LD_LIBRARY_PATH=/path/to/matlab/bin/glnxa64/:/path/to/matlab/sys/os/glnxa64/ valgrind --tool=callgrind ./test ./mex_file.mexa64

请注意,需要设置MATLAB的路径和正确的架构相关库路径! matlab可执行文件必须存在于PATH中,否则'test'将失败。

还有一个问题。 MATLAB引擎需要在系统上安装csh(你可以使用任何shell,csh只需要存在于/ bin中)。因此,如果您没有它,则必须安装它才能使其正常工作。

答案 1 :(得分:4)

您可以使用-D选项启动MATLAB,如MatlabCentral thread

所述

matlab -nojvm -nodesktop -nosplash -D"valgrind --error-limit=no --leak-check=yes --tool=memcheck -v --log-file=valgrind.log"

我会添加以确保您拥有最新版本的valgrind。当我尝试使用valgrind版本3.6调试我的MEX文件时,valgrind崩溃而不是报告内存错误。