Changeset 498

Timestamp:

02/22/05 13:31:54 (20 years ago)

Author:

kil064

Message:

return to convention I=(XX+YY)/2

File:

• trunk/src/SDPol.cc (modified) (5 diffs)

trunk/src/SDPol.cc

@@ -167 +167 @@
 // array of shape (nBeam,nIF,nPol,nChan)
 //
-// We use the scaling convention I=(XX+YY) 
-// 
+// We use the scaling convention I=(XX+YY)/2
+
 {
 
@@ -203 +203 @@
 //
    if (nPol==1) {
-      I = Float(2.0)*C1;
+      I = C1;
       return;
    }
@@ -211 +211 @@
    Array<Float> C2 = input2(start,end);          // Input : C1
 //
-   I = C1 + C2;
+   I = Float(0.5) * (C1 + C2);
    if (nPol <= 2) return;
 //
@@ -225 +225 @@
    end(polAxis) = 1;
    Array<Float> Q = output(start,end);           // Output : Q
-   Q = C1 - C2;
+   Q = Float(0.5) * (C1 - C2);
 //
    start(polAxis) = 2;
    end(polAxis) = 2;
    Array<Float> U = output(start,end);           // Output : U
-   U = Float(2.0)*C3;
+   U = C3;
 //
    start(polAxis) = 3;
    end(polAxis) = 3;
    Array<Float> V = output(start,end);           // Output : V
-   V = Float(2.0)*C4;
+   V = C4;
 }
 
@@ -315 +315 @@
                                                         Array<Float>& V, 
                                                         Bool doRR)
+//
+// We use the convention
+//  I = (RR+LL)/2
+//
 {
    if (doRR) {
-      return Float(0.5)*(I+V);
+      return I + V;
    } else {
-      return Float(0.5)*(I-V);
+      return I - V;
    }
 }