Changeset 1189 for trunk/src/STPolLinear.cpp

Timestamp:

08/28/06 12:28:27 (18 years ago)

Author:

mar637

Message:

Fix for ticket #61 - changing of default polarisation. Done ticket #62 as far as possible. Tid data gets automatically set as 'circular'\nAlso add circular pol class which was missing so far.

File:

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

trunk/src/STPolLinear.cpp

@@ -49 +49 @@
         break;
     }
+  } else if (nspec() == 2) {
+    if ( index == 0 )
+      out = Vector<Float>(getSpectrum(0) + getSpectrum(1));
   }
   return out;
@@ -55 +58 @@
 Vector<Float> asap::STPolLinear::getLinPol( uInt index )
 {
+  if (nspec() != 4) {
+    throw(AipsError("You must have 4 linear polarizations to run this function"));
+  }
   if ( index < 0 || index >4 ) throw(AipsError("LinPol index out of range"));
   Vector<Float> out,q,u;
@@ -87 +93 @@
   //   We use the convention
   //    I = (RR+LL)  // definition changed
-
+  if (nspec() != 4) {
+    throw(AipsError("You must have 4 linear polarizations to run this function"));
+  }
   if ( index == 2 || index ==3  ) throw(AipsError("Re/Imag RL not implemented"));
   Vector<Float> I,V,out;
@@ -121 +129 @@
 {
   // phase isnt used, just ro keep interface the same for all pol operations
+  if (nspec() != 4) {
+    throw(AipsError("You must have 4 linear polarizations to run this function"));
+  }
   Matrix<Float>& specs = getSpectra();
   Vector<Float> I = specs.column(3);
@@ -135 +146 @@
 // C1 = XX; C2 = YY, C3 = Real(XY)
 //
+  if (nspec() != 4) {
+    throw(AipsError("You must have 4 linear polarizations to run this function"));
+  }
   Vector<Float> I,Q,U;
   I = getStokes(0);