Changeset 40313 for trunk/Ohana/src/opihi/cmd.data/periodogram-fm.c

Timestamp:

Jan 18, 2018, 11:28:51 AM (8 years ago)

Author:

eugene

Message:

adding option for optimal sampling

File:

• trunk/Ohana/src/opihi/cmd.data/periodogram-fm.c (modified) (3 diffs)

trunk/Ohana/src/opihi/cmd.data/periodogram-fm.c

@@ -8 +8 @@
   if ((N = get_argument (argc, argv, "-v"))) {
     VERBOSE = TRUE;
+    remove_argument (N, &argc, argv);
+  }
+
+  int OPTIMAL = FALSE;
+  if ((N = get_argument (argc, argv, "-optimal"))) {
+    OPTIMAL = TRUE;
     remove_argument (N, &argc, argv);
   }
@@ -82 +88 @@
   // XXX consider choice of period or frequency step
   double P = minP;
-  for (Np = 0; Np < Nperiods; Np++) {
+  for (Np = 0; (Np < Nperiods) && (P <= maxP); Np++) {
     double w = 2*M_PI/P;
     
@@ -154 +160 @@
     if (VERBOSE) gprint (GP_ERR, "tau: %f, P: %f, ratio: %f, dTime: %f, nextP: %f\n", tau, P, dP, dTime, P*dP);
 
+    if (OPTIMAL) {
+      dP = 0.4 * SQ(P) / dTime;
+      P += dP;
+      continue;
+    }
+
     if (LINEAR) {
       P += dP;