Matlab mex-файл с mexCallMATLAB почти в 300 раз медленнее соответствующего m-файла

Я начал реализовывать несколько M-файлов на C++, чтобы сократить время выполнения. M-файлы создают N-мерные точки и оценивают значения функций в этих точках. Функции определяются пользователем и передаются в m-файлы и mex-файлы в качестве дескрипторов функций. Mex-файлы используют mexCallMATLAB с feval для поиска значений функций.

Я построил приведенный ниже пример, где дескриптор функции fn, построенный в командной строке Matlab, передается matlabcallingmatlab.м и mexcallingmatlab.процедуры ЧГК. С недавно открытым Matlab, mexcallingmatlab оценивает эту функцию 200000 в 241.5 секунд, в то время как matlabcallingmatlab оценивает его в 0.81522 секунд, поэтому 296 раз замедляется с реализацией mex. Эти времена являются результатами вторых запусков, поскольку первые запуски кажутся больше, вероятно, из-за некоторых накладных расходов, связанных с первой загрузкой программы и т. д.

Я потратил много дней на поиск в интернете по этой проблеме и попробовал предложения по этому поводу. Я пробовал разные флаги компиляции mex для оптимизации mex, но почти не было разницы в производительности. Предыдущий пост в Stackoverflow заявил, что обновление Matlab было решением, но я использую, вероятно, последнюю версию MATLAB: 8.1.0.604 (R2013a) на Mac OS X версия: 10.8.4. Я скомпилировал файл mex с флагом largeArrayDims и без него, но это тоже не имело никакого значения. Некоторые предположили, что содержимое дескриптора функции может быть непосредственно кодируется в файле cpp, но это невозможно, так как я хотел бы предоставить этот код любому пользователю с любым типом функции с векторным входом и выходом реального числа.

насколько я узнал, файлы mex должны проходить через функцию feval для использования дескриптора функции, тогда как m-файлы могут напрямую вызывать дескрипторы функций при условии, что версия Matlab новее, чем некоторая версия.

любая помощь была бы весьма признательна.

простая созданная ручка функции в командной строке Matlab:

fn = @(x) x'*x 

matlabcallingmatlab.м :

function matlabcallingmatlab( fn )
x = zeros(2,1); 
for i = 0 : 199999
    x(2) = i; 
    f = fn( x ); 
end

mexcallingmatlab.cpp:

#include "mex.h"
#include <cstring>

void mexFunction( int nlhs, mxArray *plhs[],
                  int nrhs, const mxArray *prhs[] )
{
    mxArray *lhs[1], *rhs[2]; //parameters to be passed to feval
    double f, *xptr, x[] = {0.0, 0.0}; // x: input to f and f=f(x)
    int n = 2, nbytes = n * sizeof(double);  // n: dimension of input x to f

    // prhs[0] is the function handle as first argument to feval
    rhs[0] = const_cast<mxArray *>( prhs[0] );

    // rhs[1] contains input x to the function
    rhs[1] = mxCreateDoubleMatrix( n, 1, mxREAL);
    xptr = mxGetPr( rhs[1] );

    for (int i = 0; i < 200000; ++i)
    {
        x[1] = double(i);   // change input 
        memcpy( xptr, x, nbytes );  // now rhs[1] has new x
        mexCallMATLAB(1, lhs, 2, rhs, "feval");
        f = *mxGetPr( lhs[0] );
    }
}

компиляция файла mex:

>> mex -v -largeArrayDims mexcallingmatlab.cpp