comparison src/delta.c @ 239:b015187746fa

Add normalisation capability to xtract_lnorm()

  */
 
 /* xtract_delta.c: defines functions that extract a feature as a single value from more than one input vector */
 
 #include <math.h>
-
 #include "../xtract/libxtract.h"
 
 int xtract_flux(const double *data, const int N, const void *argv , double *result)
 {
 

 
 int xtract_lnorm(const double *data, const int N, const void *argv , double *result)
 {
 
     int n,
-        type;
+        type,
+        normalise,
+        k = 0;
 
     double order;
 
     order = *(double *)argv;
     type = *((double *)argv+1);
+    normalise = (int)*((double *)argv+2);
 
     order = order > 0 ? order : 2.0;
 
     *result = 0.0;
-
+    
     switch(type)
     {
 
     case XTRACT_POSITIVE_SLOPE:
         for(n = 0; n < N; n++)
         {
             if(data[n] > 0)
+            {
                 *result += pow(data[n], order);
+                ++k;
+            }
         }
         break;
     default:
         for(n = 0; n < N; n++)
-            *result += pow(data[n], order);
+        {
+            *result += pow(fabs(data[n]), order);
+            ++k;
+        }
         break;
 
     }
 
     *result = pow(*result, 1.0 / order);
+    
+    if (k == 0)
+    {
+        return XTRACT_NO_RESULT;
+    }
+    
+    if (normalise == 1)
+    {
+        *result = log(1 + *result);
+    }
 
     return XTRACT_SUCCESS;
 
 }