mxsmdvquad.c

#include "smdvquad.h"

void mexFunction(int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]) {
    // y = smdv_quad(Q,x)
    //check_smdv_quad_input(nrhs, prhs);

    const mxArray *Q, *x;
    double L;
    double *output;

    x = prhs[1];
    Q = prhs[0];

    const double *xVec = mxGetPr(x);
    smdv_quad(Q, &xVec, &L);

    plhs[0] = mxCreateDoubleMatrix(1, 1, mxREAL);
    output = mxGetPr(plhs[0]);
    output[0] = L;
}