c——eoverlap.docVIP

#include mex.h #include math.h // This mex function implements the equivalent of: // // function y = example_mat(x) // y = 2 * x + 1; // return // [y, z] = foo(a, b, c) // // nrhs = 3 // a : prhs[0] // b : prhs[1] // c : prhs[2] // // nlhs = 2 // y : plhs[0] // z : plhs[1] void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { // must have single double array input. if (nrhs != 3) { printf(in nrhs != 3\n); return; } if (mxGetClassID(prhs[0]) != mxDOUBLE_CLASS) { printf(input must be a double array\n); return; } if (mxGetClassID(prhs[1]) != mxDOUBLE_CLASS) { printf(input must be a double array\n); return; } // must have single output. if (nlhs != 4) { printf(out nlhs != 4\n); return; } // get size of x. int const num_dims1 = mxGetNumberOfDimensions(prhs[0]); int const *dims1 = mxGetDimensions(prhs[0]); int const num_dims2 = mxGetNumberOfDimensions(prhs[1]); int const *dims2 = mxGetDimensions(prhs[1]); int const num_dims3 = mxGetNumberOfDimensions(prhs[2]); int const *dims3 = mxGetDimensions(prhs[2]); if(dims1[0]9 || dims2[0]9 || dims1[0]!=dims2[0] || dims3[0]!=1){ printf(dims != 9\n); return; } // create new array of the same size as x. int *odim = new int[2]; odim[0]=dims2[1]; odim[1]=dims1[1]; plhs[0] = mxCreateNumericArray(2, odim, mxDOUBLE_CLASS, mxREAL); plhs[1] = mxCreateNumericArray(2, odim, mxDOUBLE_CLASS, mxREAL); plhs[2] = mxCreateNumericArray(2, odim, mxDOUBLE_CLASS, mxREAL); plhs[3] = mxCreateNumericArray(2, odim, mxDOUBLE_CLASS, mxREAL); // get pointers to beginning of x and y. double const *feat1 = (double *)mxGetData(prhs[0]); double const *feat2 = (double *)mxGetData(prhs[1]); double const *flag = (double *)mxGetData(prhs[2]); double *over_out = (double *)mxGetData(plhs[0]); double *mover_out = (double *)mxGetData(plhs[1]); double *desc_out = (double *)mxGetData(plhs[2]); double *mdesc_out = (doubl