Changeset 3934 in Sophya for trunk/Cosmo/RadioBeam/mdish.cc

Timestamp:

Dec 24, 2010, 11:03:28 AM (15 years ago)

Author:

ansari

Message:

optimisation calcul reponse multi-dish ds plan (u,v), Reza 24/12/2010

File:

• trunk/Cosmo/RadioBeam/mdish.cc (modified) (3 diffs)

trunk/Cosmo/RadioBeam/mdish.cc

@@ -309 +309 @@
 
   for(size_t i=0; i<dishes_.size(); i++) {
-    for(size_t j=0; j<dishes_.size(); j++) {
+    for(size_t j=i; j<dishes_.size(); j++) {
       double kx0 = DeuxPI*(dishes_[i].X-dishes_[j].X)/lambda_;
       double ky0 = DeuxPI*(dishes_[i].Y-dishes_[j].Y)/lambda_;
@@ -322 +322 @@
           if ((ix>0)&&(jy>0)) {
             sumw += AddToHisto(kx0, ky0, x, y, rd(x,y), fgfh);
+            if (j!=i) sumw += AddToHisto(-kx0, -ky0, x, y, rd(x,y), fgfh);
           }
           else {
-            if ((ix==0)&&(jy==0)) 
+            if ((ix==0)&&(jy==0)) {
               sumw += h2w_.Add(kx0, ky0, rd(0.,0.), fgfh);
+              if (j!=i) sumw += h2w_.Add(-kx0, -ky0, rd(0.,0.), fgfh);
+            }
             else {
               double w = rd(x,y);
@@ -331 +334 @@
                 sumw += h2w_.Add(kx0, ky0+y, w, fgfh);
                 sumw += h2w_.Add(kx0, ky0-y, w, fgfh);
+                if (j!=i) {
+                  sumw += h2w_.Add(-kx0, -ky0+y, w, fgfh);
+                  sumw += h2w_.Add(-kx0, -ky0-y, w, fgfh);
+                }
               }
               else {
                 sumw += h2w_.Add(kx0+x, ky0, w, fgfh);
                 sumw += h2w_.Add(kx0-x, ky0, w, fgfh);
+                if (j!=i) {
+                  sumw += h2w_.Add(-kx0+x, -ky0, w, fgfh);
+                  sumw += h2w_.Add(-kx0-x, -ky0, w, fgfh);
+                }               
               }
             }