Changeset 3894 for trunk/Ohana/src/opihi/lib.data/gaussian.c

Timestamp:

May 11, 2005, 9:35:40 PM (21 years ago)

Author:

eugene

Message:

cleaned gaussdeviate, added vsmooth

File:

• trunk/Ohana/src/opihi/lib.data/gaussian.c (modified) (4 diffs)

trunk/Ohana/src/opihi/lib.data/gaussian.c

@@ -1 +1 @@
 # include "data.h"
 
-static int Ngaussint;
+static int Ngaussint = 0;
 static double *gaussint;
 
@@ -24 +24 @@
   double mean, sigma;
  
+  /* no need to generate this if it already exists */
+  if (Ngaussint == Nbin) return;
+
   A = time(NULL);
   for (B = 0; A == time(NULL); B++);
@@ -45 +48 @@
           3.0*gaussian(x+dx3, mean, sigma)) * (dx1/8.0);
     val += df;
-    if (val > (i / (double) Ngaussint)) {
+    if (val > (i + 0.5) / (double) Ngaussint) {
       gaussint[i] = x + dx / 2.0;
       i++;
@@ -65 +68 @@
 }
  
+double int_gauss (int i) {
+  double y;
+  y = gaussint[i];
+  return (y);
+}
+