Daniel@0: /* C mex version of normalise.m in misc directory  */
Daniel@0: 
Daniel@0: #include "mex.h"
Daniel@0: 
Daniel@0: void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
Daniel@0: {
Daniel@0:   double *T, *sum_ptr, sum;
Daniel@0:   int i, N;  
Daniel@0: 
Daniel@0:   plhs[0] = mxDuplicateArray(prhs[0]);
Daniel@0:   T = mxGetPr(plhs[0]);
Daniel@0:   if(mxIsSparse(plhs[0])) N = mxGetJc(plhs[0])[mxGetN(plhs[0])];
Daniel@0:   else N = mxGetNumberOfElements(plhs[0]);
Daniel@0: 
Daniel@0:   plhs[1] = mxCreateDoubleMatrix(1, 1, mxREAL);
Daniel@0:   sum_ptr = mxGetPr(plhs[1]);
Daniel@0: 
Daniel@0:   sum = 0;
Daniel@0:   for (i = 0; i < N; i++) {
Daniel@0:     sum += *T++;
Daniel@0:   }
Daniel@0:   T = mxGetPr(plhs[0]);
Daniel@0:   if (sum > 0) {
Daniel@0:     for (i = 0; i < N; i++) {
Daniel@0:       *T++ /= sum;
Daniel@0:     }
Daniel@0:   } 
Daniel@0:   *sum_ptr = sum;
Daniel@0: }