casa6

#include <casacore/casa/Arrays/ArrayMath.h>

#include <casacore/casa/Arrays/ArrayPartMath.h>

#include <casacore/casa/BasicMath/Functors.h>

#include <casacore/casa/Utilities/CountedPtr.h>

#include <stdcasa/UtilJ.h>

#include <msvis/MSVis/AveragingTvi2.h>

#include <msvis/MSVis/AveragingVi2Factory.h>

#include <msvis/MSVis/MsRows.h>

#include <msvis/MSVis/VisBufferComponents2.h>

#include <msvis/MSVis/VisBuffer2.h>

#include <msvis/MSVis/VisBufferImpl2.h>

#include <msvis/MSVis/Vbi2MsRow.h>

#include <msvis/MSVis/VisibilityIterator2.h>

#include <tuple>

#include <set>

​

using std::set;

​

using namespace casacore;

namespace casa {

​

namespace vi {

​

namespace avg {

​

using casa::ms::MsRow;

​

​

///////////////////////////////////////////////////////////

//

// VbAvg: Averaging VisBuffer

//

/*

​

Definition: A baseline sample (baseline for short) is a table row

with a particular (antenna1, antenna2) pair at a particular time.

​

​

Averaging does not cross chunk boundaries of the input VI so the

input VI determines what values are averaged together.  For example,

if the input VI is allows data from multiple scans into the same chunk

then some averaging across scans can occur; in this case the scan number

of the result will be the scan number of the first baseline instance

seen.

​

Row-level value treatment

=========================

​

The average is approximately computed on a per baseline basis:

​

averaged_baseline (antenna1, antenna2) =

    sumOverAveragingInterval (baseline (antenna1, antenna2)) /

    nBaselinesFoundInInterval

​

How row-level values are computed

---------------------------------

​

Time - Set to time of first baseline making up the average plus

       half of the averaging interval.

Antenna1 - Copied from instance of baseline

Antenna2 - Copied from instance of baseline

Feed1 - Copied from instance of baseline

Feed2 - Copied from instance of baseline