Changeset 797 for branches/Release12

Timestamp:

12/19/05 15:32:36 (18 years ago)

Author:

phi196

Message:

Added swapPol & invertPhase

File:

• branches/Release12/src/SDMath.cc (modified) (1 diff)

branches/Release12/src/SDMath.cc

@@ -1384 +1384 @@
 }
 
+void SDMath::invertPhase(SDMemTable& in, Bool doAll)
+//
+// Takes complex conjugate of complex correlation
+// should work ok for linear and circular polarisation
+//
+{
+  if (in.nPol() != 4) {
+    throw(AipsError("You must have 4 polarizations to run this function"));
+  }
+
+  //SDHeader sh = in.getSDHeader();
+  //Instrument inst = SDAttr::convertInstrument(sh.antennaname, False);
+  //   SDAttr sdAtt;
+  // if (sdAtt.feedPolType(inst) != LINEAR) {
+  //   throw(AipsError("Only linear polarizations are supported"));
+  //}
+//
+   const Table& tabIn = in.table();
+   ArrayColumn<Float> specCol(tabIn,"SPECTRA");
+   IPosition start(asap::nAxes,0);
+   IPosition end(asap::nAxes);
+
+// Set cursor slice. Assumes shape the same for all rows
+
+   setCursorSlice (start, end, doAll, in);
+
+   IPosition start4(start);
+   start4(asap::PolAxis) = 3;                 // Imag (XY)
+   IPosition end4(end);
+   end4(asap::PolAxis) = 3;
+//
+   uInt nRow = in.nRow();
+   Array<Float> data;
+   for (uInt i=0; i<nRow;++i) {
+      specCol.get(i,data);
+      IPosition shape = data.shape();
+
+      // Get polarization slice references
+      Array<Float> C4 = data(start4,end4);
+
+      // Invert
+      C4 = -C4
+
+      // Put
+      specCol.put(i,data);
+   }
+}
+
+void SDMath::swapPol(SDMemTable& in, Bool doAll)
+//
+// phase in degrees
+// assumes linear correlations
+//
+{
+   if (in.nPol() != 4 && in.nPol() != 2) {
+      throw(AipsError("You must have 2 or 4 polarizations to run this function"));
+   }
+//
+//   SDHeader sh = in.getSDHeader();
+//   Instrument inst = SDAttr::convertInstrument(sh.antennaname, False);
+//   SDAttr sdAtt;
+   //  if (sdAtt.feedPolType(inst) != LINEAR) {
+   //   throw(AipsError("Only linear polarizations are supported"));
+   //  }
+//
+   const Table& tabIn = in.table();
+   ArrayColumn<Float> specCol(tabIn,"SPECTRA");
+   ROArrayColumn<Float> stokesCol(tabIn,"STOKES");
+   IPosition start(asap::nAxes,0);
+   IPosition end(asap::nAxes);
+
+// Set cursor slice. Assumes shape the same for all rows
+
+   setCursorSlice (start, end, doAll, in);
+//
+   IPosition start1(start);
+   start1(asap::PolAxis) = 0;                // C1 (XX)
+   IPosition end1(end);
+   end1(asap::PolAxis) = 0;
+//
+   IPosition start2(start);
+   start2(asap::PolAxis) = 1;                 // C2 (YY)
+   IPosition end2(end);
+   end2(asap::PolAxis) = 1;
+//
+   uInt nRow = in.nRow();
+   Array<Float> data, stokes;
+   for (uInt i=0; i<nRow;++i) {
+      specCol.get(i,data);
+      stokesCol.get(i,stokes);
+      IPosition shape = data.shape();
+
+// Get linear polarization slice references
+
+      Array<Float> C1 = data(start1,end1);
+      Array<Float> C2 = data(start2,end2);
+
+// Swap
+
+      Array<Float> tmp = C1;
+      C1 = C2;
+      C2 = C1;
+
+// Put
+
+      specCol.put(i,data);
+   }
+}
+
 // 'private' functions