#include <synthesis/TransformMachines/VLACalcIlluminationConvFunc.h>
#include <coordinates/Coordinates/DirectionCoordinate.h>
#include <coordinates/Coordinates/SpectralCoordinate.h>
#include <coordinates/Coordinates/LinearCoordinate.h>
#include <coordinates/Coordinates/StokesCoordinate.h>
#include <synthesis/TransformMachines/Utils.h>
#include <synthesis/TransformMachines/SynthesisError.h>
#include <ms/MeasurementSets/MSColumns.h>
#include <lattices/LEL/LatticeExpr.h>
#include <lattices/LatticeMath/LatticeFFT.h>
#include <images/Images/ImageRegrid.h>
#include <images/Images/PagedImage.h>
#include <casa/Arrays/ArrayMath.h>
#include <casa/Utilities/Assert.h>
#include <casa/OS/File.h>
#include <casa/sstream.h>
#include <casa/OS/Timer.h>
using namespace casacore;
VLACalcIlluminationConvFunc::VLACalcIlluminationConvFunc():IlluminationConvFunc()
LogIO logIO(LogOrigin("VLACalcIlluminationConvFunc","ctor"));
ap.x0=-13.0; ap.y0=-13.0;
ap.band = BeamCalc_VLA_L;
IPosition shape(4,ap.nx,ap.ny,1,1);
ap.aperture = new TempImage<Complex>();
if (maximumCacheSize() > 0) ap.aperture->setMaximumCacheSize(maximumCacheSize());
ap.aperture->resize(shape);
CoordinateSystem VLACalcIlluminationConvFunc::makeUVCoords(CoordinateSystem& imageCoordSys,
CoordinateSystem FTCoords = imageCoordSys;
Int dirIndex=FTCoords.findCoordinate(Coordinate::LINEAR);
if (dirIndex >= 0) return FTCoords;
dirIndex=FTCoords.findCoordinate(Coordinate::DIRECTION);
DirectionCoordinate dc=imageCoordSys.directionCoordinate(dirIndex);
Vector<Bool> axes(2); axes=true;