comparison toolboxes/FullBNT-1.0.7/KPMtools/normaliseC.c @ 0:e9a9cd732c1e tip

first hg version after svn

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