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