Changeset 498 for trunk/src


Ignore:
Timestamp:
02/22/05 13:31:54 (20 years ago)
Author:
kil064
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/SDPol.cc

    r492 r498  
    167167// array of shape (nBeam,nIF,nPol,nChan)
    168168//
    169 // We use the scaling convention I=(XX+YY) 
    170 //
     169// We use the scaling convention I=(XX+YY)/2
     170
    171171{
    172172
     
    203203//
    204204   if (nPol==1) {
    205       I = Float(2.0)*C1;
     205      I = C1;
    206206      return;
    207207   }
     
    211211   Array<Float> C2 = input2(start,end);          // Input : C1
    212212//
    213    I = C1 + C2;
     213   I = Float(0.5) * (C1 + C2);
    214214   if (nPol <= 2) return;
    215215//
     
    225225   end(polAxis) = 1;
    226226   Array<Float> Q = output(start,end);           // Output : Q
    227    Q = C1 - C2;
     227   Q = Float(0.5) * (C1 - C2);
    228228//
    229229   start(polAxis) = 2;
    230230   end(polAxis) = 2;
    231231   Array<Float> U = output(start,end);           // Output : U
    232    U = Float(2.0)*C3;
     232   U = C3;
    233233//
    234234   start(polAxis) = 3;
    235235   end(polAxis) = 3;
    236236   Array<Float> V = output(start,end);           // Output : V
    237    V = Float(2.0)*C4;
     237   V = C4;
    238238}
    239239
     
    315315                                                        Array<Float>& V,
    316316                                                        Bool doRR)
     317//
     318// We use the convention
     319//  I = (RR+LL)/2
     320//
    317321{
    318322   if (doRR) {
    319       return Float(0.5)*(I+V);
     323      return I + V;
    320324   } else {
    321       return Float(0.5)*(I-V);
     325      return I - V;
    322326   }
    323327}
Note: See TracChangeset for help on using the changeset viewer.