- Timestamp:
- 02/22/05 13:31:54 (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/SDPol.cc
r492 r498 167 167 // array of shape (nBeam,nIF,nPol,nChan) 168 168 // 169 // We use the scaling convention I=(XX+YY) 170 // 169 // We use the scaling convention I=(XX+YY)/2 170 171 171 { 172 172 … … 203 203 // 204 204 if (nPol==1) { 205 I = Float(2.0)*C1;205 I = C1; 206 206 return; 207 207 } … … 211 211 Array<Float> C2 = input2(start,end); // Input : C1 212 212 // 213 I = C1 + C2;213 I = Float(0.5) * (C1 + C2); 214 214 if (nPol <= 2) return; 215 215 // … … 225 225 end(polAxis) = 1; 226 226 Array<Float> Q = output(start,end); // Output : Q 227 Q = C1 - C2;227 Q = Float(0.5) * (C1 - C2); 228 228 // 229 229 start(polAxis) = 2; 230 230 end(polAxis) = 2; 231 231 Array<Float> U = output(start,end); // Output : U 232 U = Float(2.0)*C3;232 U = C3; 233 233 // 234 234 start(polAxis) = 3; 235 235 end(polAxis) = 3; 236 236 Array<Float> V = output(start,end); // Output : V 237 V = Float(2.0)*C4;237 V = C4; 238 238 } 239 239 … … 315 315 Array<Float>& V, 316 316 Bool doRR) 317 // 318 // We use the convention 319 // I = (RR+LL)/2 320 // 317 321 { 318 322 if (doRR) { 319 return Float(0.5)*(I+V);323 return I + V; 320 324 } else { 321 return Float(0.5)*(I-V);325 return I - V; 322 326 } 323 327 }
Note:
See TracChangeset
for help on using the changeset viewer.