Changeset 1963 in Sophya for trunk/ArchTOIPipe/Processors

Timestamp:

Apr 11, 2002, 6:12:17 PM (23 years ago)

Author:

aubourg

Message:

new output for wiener, unconnected outs ok in processor

File:

• trunk/ArchTOIPipe/Processors/wienerdecor.cc (modified) (5 diffs)

trunk/ArchTOIPipe/Processors/wienerdecor.cc

@@ -16 +16 @@
   declareInput("probe");
   declareOutput("signal");
+  declareOutput("noiseestim");
   name="WienerDecorrelator";
   setNeededHistory(nsamples+lcorr+1);
@@ -33 +34 @@
   double* y = new double[lcorr+1];  // corr vector
   double* window = new double[lcorr];
-  double* filter = w+1;
+  double* filter = new double[lcorr];
+  for (int i=0; i<lcorr; i++) filter[i]=0;
   
   int sn = snb;
@@ -57 +59 @@
       }
       dtoeplz(r,w,y,lcorr);
-      cout << "Wiener filter : \n ";
+      if (!isnan(w[1])) {
+        for (int i=0; i<lcorr; i++) {
+          filter[i] = w[i+1];
+        }
+      } else {
+        cout << "Bad inversion, keeping previous filter\n";
+      }
+      cout << "Wiener filter : " << sn << "\n ";
       for (int i=0; i<lcorr; i++) {
         cout << filter[i] << " ";
@@ -71 +80 @@
       }
       putData(0, sn, getData(0, sn) - outSig);
+      putData(1, sn, outSig);
     }
     sn++;
@@ -78 +88 @@
   delete[] w;
   delete[] r;
+  delete[] filter;
 }