Changeset 3928 in Sophya for trunk/AddOn/TAcq/svv2mtx2_1210.cc

Timestamp:

Dec 17, 2010, 4:24:06 PM (15 years ago)

Author:

cmv

Message:

gestion de la conjugaison complexe pour visi # de <E.conj(W)>, cmv 17/12/2010

File:

• trunk/AddOn/TAcq/svv2mtx2_1210.cc (modified) (6 diffs)

trunk/AddOn/TAcq/svv2mtx2_1210.cc

@@ -24 +24 @@
 cout<<"svv2mtx2_1012 [options] dir : lecture des fichiers acq de Pittsburgh Dec 2010"<<endl
     <<" dir : repertoire ou se trouvent les fichiers d'acq"<<endl
+    <<" -C : compute and store the complexe conjugated visi"<<endl
     <<" -D : visi is a duplicated one"<<endl
     <<" -o visi.ppf : nom du ficher ppf pour ecrire la visi temps-frequence"<<endl
@@ -40 +41 @@
   string outname = "";
   int numthread = -1, numrow = -1;
-  bool dupli = false;
+  bool dupli = false, doconj = false;
   double freq0 = 0.;
   int ifilmin=0, ifilmax=99999;
@@ -47 +48 @@
 
   char c;
-  while((c = getopt(narg,arg,"hDo:t:r:f:T:F:")) != -1) {
+  while((c = getopt(narg,arg,"hDCo:t:r:f:T:F:")) != -1) {
     switch (c) {
     case 'D' :
       dupli = true;
       break;
+    case 'C' :
+      doconj = true;
+      break;
     case 'o' :
       outname = optarg;
@@ -83 +87 @@
   cout<<"numrow="<<numrow<<endl;  if(numrow<0) return -2;
   cout<<"dupli="<<(int)dupli<<endl;
+  cout<<"doconj="<<(int)doconj<<endl;
   cout<<"outname="<<outname<<endl; if(outname.size()<=0) return -2;
   cout<<"indir="<<indir<<endl;
@@ -204 +209 @@
       Npaqsum(ntimefill) = nf * npaqsum;
       MVisi(ntimefill,i) /= double(Npaqsum(ntimefill));
+      if(doconj) MVisi(ntimefill,i) = conj(MVisi(ntimefill,i));
     }
     ntimefill++;
@@ -223 +229 @@
   MVisi.Info()["dir"] = indir;
   MVisi.Info()["dupli"] = (dupli) ? 1: 0;
+  MVisi.Info()["doconj"] = (doconj) ? 1: 0;
   MVisi.Info()["TUobs_0"] = tudeb;
   MVisi.Info()["TUobs_N"] = tufin;