Блог

У меня есть программа на C, которая умножает скалярное число на массив. Я хотел бы создать mex-файл для вызова этой функции в MATLAB. Я перешел по этой ссылке, но у меня возникли некоторые проблемы с компиляцией, и мне пришлось изменить часть кода.

Функция C, которую я реализовал, является:

 #include "mex.h"
#include <stdio.h>
#include <stdlib.h>


void arrayProduct(double x, double *y, double *z, int n){
    int i;

    for (i=0;i<n;i  ){
        z[i]=x * y[i];
    }
}

void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]){

double multiplier;
double *inMatrix;
mwSize ncols;
double *outMatrix;

multiplier = mxGetScalar(prhs[0]);
#if MX_HAS_INTERLEAVED_COMPLEX
    inMatrix=mxGetDoubles(prhs[1]);
#else
    inMatrix = mxGetPr(prhs[1]);
#endif
ncols=mxGetN(prhs[1]);

plhs[0]=mxCreateDoubleMatrix(1,ncols,mxREAL);

#if MX_HAS_INTERLEAVED_COMPLEX
    outMatrix=mxGetDoubles(plhs[0]);
#else
    outMatrix=mxGetPr(prhs[0]);
#endif
arrayProduct(multiplier,inMatrix,outMatrix,ncols);

if(nrhs != 2){
    mexErrMsgIdAndTxt("MyToolbox:arrayProduct:nrhs","Two inputs required.");
}
if(nlhs != 1){
    mexErrMsgIdAndTxt("MyToolbox:arrayProduct:nlhs","One input required.");
}
if( !mxIsDouble(prhs[0]) || mxIsComplex(prhs[0]) || mxGetNumberOfElements(prhs[0]) != 1){
    mexErrMsgIdAndTxt("MyToolbox:arrayProduct:notScalar","Input multiplier must be scalar.");
}


if(!mxIsDouble(prhs[1]) || mxIsComplex(prhs[1])){
    mexErrMsgIdAndTxt("MyToolbox:arrayProduct:notDouble","Input matrix must be type double.");
}
 
if(mxGetM(prhs[1]) != 1){
    mexErrMsgIdAndTxt("MyToolbox:arrayProduct:notRowVector","Input must be a row vector.");
}
}
 

Прямо сейчас я могу создать mex-файл (.mexmaci64, я использую macOS), но когда я вызываю функцию arrayProduct(s,A), где s = 5 и A = [1.5,2,9], я получаю это сообщение:

Ошибка при использовании arrayProduct Требуется один ввод.

Как я могу исправить эту ошибку, чтобы сгенерировать правильный mex-файл или правильно запустить функцию?

Любая помощь приветствуется.