Source
636
636
Complex cOne(1.0);
637
637
638
638
for (Int irow=0;irow<vb.nRows();++irow) {
639
639
if (vb.flagRow()(irow)) {
640
640
// Row flagged, make sure cube also flagged, weight/data zeroed
641
641
c(Slice(),Slice(),Slice(irow,1,1))=0.0f;
642
642
w(Slice(),Slice(),Slice(irow,1,1))=0.0f;
643
643
fl(Slice(),Slice(),Slice(irow,1,1))=True;
644
644
}
645
645
else {
646
-
Bool *flp=&fl(0,0,irow);
646
+
// Bool *flp=&fl(0,0,irow);
647
647
Float *wtp=&w(0,0,irow);
648
648
Complex *cvp=&c(0,0,irow);
649
649
Complex *mvp=&m(0,0,irow);
650
650
651
651
for (Int ichan=0;ichan<vb.nChannels();++ichan) {
652
652
Complex Imod(0.0);
653
653
Float Iamp2(1.0);
654
654
if (!fl(0,ichan,irow) && !fl(nCorr-1,ichan,irow)) {
655
655
Imod=(m(0,ichan,irow) + m(nCorr-1,ichan,irow))/2.0f;
656
656
Iamp2=real(Imod*conj(Imod)); // squared model amp (for weight adjust)