Changeset 1007 for trunk/src/STPolLinear.cpp

Timestamp:

04/07/06 13:29:08 (18 years ago)

Author:

mar637

Message:

Fixed computiaons of stokes and back after discussions with Jim Caswell.

File:

• trunk/src/STPolLinear.cpp (modified) (3 diffs)

trunk/src/STPolLinear.cpp

@@ -40 +40 @@
         break;
       case 1:
-        out = Vector<Float>(q * cos(phase) - getSpectrum(2) * sin(phase));
+        out = Vector<Float>(q * cos(phase) - Float(2.0)*getSpectrum(2) * sin(phase));
         break;
       case 2:
-        out = Vector<Float>(q * sin(phase) + getSpectrum(2) * cos(phase));
+        out = Vector<Float>(q * sin(phase) + Float(2.0)*getSpectrum(2) * cos(phase));
         break;
       case 3:
-        out = getFeedHand() * Vector<Float>(getSpectrum(3));
+        out = getFeedHand() * Float(2.0) * Vector<Float>(getSpectrum(3));
         break;
     }
@@ -94 +94 @@
   switch(index) {
   case 0:
-    out = (I + V);
+    out = (I + V)/Float(2.0);
     break;
   case 1:
-    out = (I - V);
+    out = (I - V)/Float(2.0);
     break;
   default:
@@ -146 +146 @@
   Q = Q2;
   Matrix<Float>& specs = getSpectra();
-  specs.column(0) = I+Q;
-  specs.column(1) = I-Q;
-  specs.column(2) = U;
+  specs.column(0) = (I+Q)/Float(2.0);
+  specs.column(1) = (I-Q)/Float(2.0);
+  specs.column(2) = U/Float(2.0);
 
 }