casa6

​

#include <msvis/MSVis/VisibilityIterator2.h>

#include <msvis/MSVis/ViImplementation2.h>

#include <msvis/MSVis/VisBufferImpl2.h>

#include <casacore/measures/Measures/MDirection.h>

#include <casacore/measures/Measures/MEpoch.h>

#include <casacore/casa/Quanta/MVTime.h>

#include <casacore/ms/MSOper/MSDerivedValues.h>

#include <casacore/casa/Arrays.h>

​

#include <algorithm>

using namespace std;

​

using namespace casacore;

namespace casa {

​

namespace vi {

​

// #pragma message "Change to pure abstract method before checkin."

// Vector<MPosition>

// ViImplementation2::antennaPositions () const

// {

//     return Vector<MPosition> ();

// }

​

​

VisBuffer2 *

ViImplementation2::createAttachedVisBuffer(VisBufferOptions options)

{

    return new VisBufferImpl2(this, options);

}

​

VisBuffer2 *

ViImplementation2::getVisBuffer(const VisibilityIterator2 * vi) const

{

    VisBuffer2 *result = getVisBuffer();

    ThrowIf(result == nullptr, "VI Implementation has no VisBuffer.");

    result->associateWithVi2(vi);

    return result;

}

​

void

ViImplementation2::azel0Calculate (Double time, MSDerivedValues & msd,

        MDirection & azel0, const MEpoch & mEpoch0)

{

    // Refactored into a static method to allow VisBufferAsync to use

​

    MEpoch mEpoch = mEpoch0;

​

    mEpoch.set (MVEpoch (Quantity (time, "s")));

​

    msd.setEpoch (mEpoch);

​

    msd.setAntenna (-1);

​

    azel0 = msd.azel ();

​

    if (aips_debug) {

        cout << "At time: " << MVTime (mEpoch.getValue ()) <<

             " AzEl = " << azel0.getAngle ("deg") << endl;

    }

}

​