/*
 * mexAdd.cpp - matlab külső futtatható kód példa
 *
 * Összead két bemeneti vektor (v1, v2)
 * a kimenetre adja az eredményt (result)
 *
 * The calling syntax is:
 *
 *		result = mexAdd(v1, v2)
*/

#include "mex.h"

void arrayAdd(double *x, double *y, double *z, mwSize n)
{
  mwSize i;
  
  for (i=0; i<n; i++) {
    z[i] = x[i] + y[i];
  }
}

/* The gateway function */
void mexFunction(int nlhs, mxArray *plhs[],
                 int nrhs, const mxArray *prhs[])
{
    if(nrhs != 2) {
        mexErrMsgIdAndTxt("MyToolbox:arrayProduct:nrhs",
                          "Two inputs required.");
    }

    if(nlhs != 1) {
        mexErrMsgIdAndTxt("MyToolbox:arrayProduct:nlhs",
                          "One output required.");
    }
    
    if( !mxIsDouble(prhs[0]) || 
         mxIsComplex(prhs[0])) {
        mexErrMsgIdAndTxt("MyToolbox:arrayProduct:notDouble",
                          "Input multiplier must be type double.");
    }

    if( !mxIsDouble(prhs[1]) || 
         mxIsComplex(prhs[1])) {
        mexErrMsgIdAndTxt("MyToolbox:arrayProduct:notDouble",
                          "Input matrix must be type double.");
    }
    
    if(mxGetM(prhs[0]) != 1) {
        mexErrMsgIdAndTxt("MyToolbox:arrayProduct:notRowVector",
                          "First input must be a row vector.");
    }
    
    if(mxGetM(prhs[1]) != 1) {
        mexErrMsgIdAndTxt("MyToolbox:arrayProduct:notRowVector",
                          "Second input must be a row vector.");
    }
    
    if(mxGetN(prhs[0]) != mxGetN(prhs[1])) {
        mexErrMsgIdAndTxt("MyToolbox:arrayProduct:notRowVector",
                          "Two input must have the same size.");
    }
    
    double *inVector1 = mxGetPr(prhs[0]);
    double *inVector2 = mxGetPr(prhs[1]);
    
    mwSize N;
    
    double *outVector;

    N = mxGetN(prhs[0]);
    
    plhs[0] = mxCreateDoubleMatrix(1,N,mxREAL);
    outVector = mxGetPr(plhs[0]);
    
    arrayAdd(inVector1,inVector2,outVector,N);

}