casa5

//# BriggsCubeWeight.cc: Implementation for Briggs weighting for cubes

//# Copyright (C) 2018-2019

//# Associated Universities, Inc. Washington DC, USA.

//#

//# This library is free software; you can redistribute it and/or modify it

//# under the terms of the GNU General Public License as published by

//# the Free Software Foundation; either version 3 of the License, or (at your

//# option) any later version.

//#

//# This library is distributed in the hope that it will be useful, but WITHOUT

//# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or

//# FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public

//# License for more details.

//#

//# https://www.gnu.org/licenses/

//#

//# You should have received a copy of the GNU  General Public License

//# along with this library; if not, write to the Free Software Foundation,

//# Inc., 675 Massachusetts Ave, Cambridge, MA 02139, USA.

//#

//# Queries concerning CASA should be submitted at

//#        https://help.nrao.edu

//#

//#        Postal address: CASA Project Manager 

//#                        National Radio Astronomy Observatory

//#                        520 Edgemont Road

//#                        Charlottesville, VA 22903-2475 USA

//#

//#

#include <casa/Arrays/ArrayMath.h>

#include <casa/Arrays/Matrix.h>

#include <casa/Arrays/Vector.h>

#include <casa/Containers/Record.h>

#include<msvis/MSVis/VisImagingWeight.h>

#include <msvis/MSVis/VisBuffer2.h>

#include <msvis/MSVis/VisibilityIterator2.h>

#include<synthesis/TransformMachines2/FTMachine.h>

#include<synthesis/TransformMachines2/BriggsCubeWeightor.h>

​

​

namespace casa{//# CASA namespace

namespace refim {//# namespace refactor imaging

using namespace casacore;

using namespace casa;

using namespace casacore;

using namespace casa::refim;

using namespace casacore;

using namespace casa::vi;

​

  BriggsCubeWeightor::BriggsCubeWeightor(): grids_p(0), ft_p(0), f2_p(0), d2_p(0), uscale_p(0), vscale_p(0), uorigin_p(0),vorigin_p(0), nx_p(0), ny_p(0), rmode_p(""), noise_p(0.0), robust_p(2), superUniformBox_p(0), multiField_p(False),initialized_p(False), refFreq_p(-1.0), freqInterpMethod_p(InterpolateArray1D<Double, Complex>::nearestNeighbour) {

    multiFieldMap_p.clear();

  }

  BriggsCubeWeightor::BriggsCubeWeightor( const String& rmode, const Quantity& noise, const Double robust, const Int superUniformBox, const Bool multiField)  : grids_p(0), ft_p(0), f2_p(0), d2_p(0), uscale_p(0), vscale_p(0), uorigin_p(0),vorigin_p(0), nx_p(0), ny_p(0), initialized_p(False), refFreq_p(-1.0),freqInterpMethod_p(InterpolateArray1D<Double, Complex>::nearestNeighbour) {

​

    rmode_p=rmode;

    noise_p=noise;

    robust_p=robust;

    superUniformBox_p=superUniformBox;

    multiField_p=multiField;

    multiFieldMap_p.clear();

  }

​

​

  BriggsCubeWeightor::BriggsCubeWeightor(vi::VisibilityIterator2& vi,

                       const String& rmode, const Quantity& noise,

                       const Double robust,

                                         const ImageInterface<Complex>& templateimage, const RecordInterface& inrec,

                     const Int superUniformBox, const Bool multiField){

    rmode_p=rmode;

    noise_p=noise;

    robust_p=robust;

    superUniformBox_p=superUniformBox;

    multiField_p=multiField;

    initialized_p=False;

    refFreq_p=-1.0;

    init(vi, templateimage,inrec);

​

​

​

​

  }

 void BriggsCubeWeightor::init(vi::VisibilityIterator2& vi,

                   const ImageInterface<Complex>& templateimage, const RecordInterface& inRec)

 {

  LogIO os(LogOrigin("BriggsCubeWeightor", "constructor", WHERE));