Changeset 79 for trunk/src/SDContainer.cc

Timestamp:

08/10/04 16:26:33 (20 years ago)

Author:

mar637

Message:

Added support of source direction.

File:

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

trunk/src/SDContainer.cc

@@ -66 +66 @@
     flags_(IPosition(4,nBeam,nIF,nPol,nChan)),
     tsys_(IPosition(4,nBeam,nIF,nPol,nChan)),
-    freqidx_(nIF) {
+    freqidx_(nIF),
+    direction_(IPosition(2,nBeam,2)) {
   uChar x = 0;
   flags_ = ~x;
@@ -79 +80 @@
     flags_(shp),
     tsys_(shp),
-    freqidx_(shp(1)){
+    freqidx_(shp(1)) {
+  IPosition ip(2,shp(0),2);
+  direction_.resize(ip);
   uChar x = 0;
   flags_ = ~x;
@@ -96 +99 @@
   tsys_.resize(shp);
   freqidx_.resize(shp(1));
+  IPosition ip(2,shp(0),2);
+  direction_.resize(ip);
 }
 
@@ -284 +289 @@
 }
 
+Array<Double> SDContainer::getDirection(uInt whichBeam) const {
+  Vector<Double> direct(2);
+  ArrayAccessor<Double, Axis<0> > i0(direction_);
+  i0.reset(i0.begin(whichBeam));
+  ArrayAccessor<Double, Axis<0> > o0(direct);
+  ArrayAccessor<Double, Axis<1> > i1(i0);
+  while (i1 != i1.end()) {
+    *o0 = *i1;
+    i1++;
+    o0++;
+  }  
+  return direct.copy();
+}
+
+
 Bool SDContainer::setFrequencyMap(uInt freqslot, uInt whichIF) {
   freqidx_[whichIF] = freqslot;
@@ -292 +312 @@
   freqidx_.resize();
   freqidx_ = freqs;
+  return True;
+}
+
+Bool SDContainer::setDirection(const Vector<Double>& point, uInt whichBeam) {
+  if (point.nelements() != 2) return False;
+  ArrayAccessor<Double, Axis<0> > aa0(direction_);
+  aa0.reset(aa0.begin(whichBeam));
+  ArrayAccessor<Double, Axis<0> > jj(point);
+  for (ArrayAccessor<Double, Axis<1> > i(aa0);i != i.end(); ++i) {
+    
+    (*i) = (*jj);
+    jj++;
+  }
+  return True;
+}
+
+Bool SDContainer::putDirection(const Array<Double>& dir) {
+  direction_.resize();
+  direction_ = dir;
   return True;
 }