MultiTermFT.h

Go to the documentation of this file.

//# MultiTermFT.h: Definition for MultiTermFT
//# Copyright (C) 1996,1997,1998,1999,2000,2002
//# 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 Library General Public License as published by
//# the Free Software Foundation; either version 2 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 Library General Public
//# License for more details.
//#
//# You should have received a copy of the GNU Library General Public License
//# along with this library; if not, write to the Free Software Foundation,
//# Inc., 675 Massachusetts Ave, Cambridge, MA 02139, USA.
//#
//# Correspondence concerning AIPS++ should be adressed as follows:
//#        Internet email: aips2-request@nrao.edu.
//#        Postal address: AIPS++ Project Office
//#                        National Radio Astronomy Observatory
//#                        520 Edgemont Road
//#                        Charlottesville, VA 22903-2475 USA
//#
//#
//# $Id$

#ifndef SYNTHESIS_MULTITERMFT_H
#define SYNTHESIS_MULTITERMFT_H

#include <synthesis/TransformMachines/FTMachine.h>
#include <casa/Arrays/Matrix.h>
#include <scimath/Mathematics/FFTServer.h>
#include <msvis/MSVis/VisBuffer.h>
#include <images/Images/ImageInterface.h>
#include <images/Images/ImageInterface.h>
#include <casa/Containers/Block.h>
#include <casa/Arrays/Array.h>
#include <casa/Arrays/Vector.h>
#include <casa/Arrays/Matrix.h>
#include <scimath/Mathematics/ConvolveGridder.h>
#include <lattices/Lattices/LatticeCache.h>
#include <lattices/Lattices/ArrayLattice.h>
//#include <synthesis/MeasurementComponents/SynthesisPeek.h>
#include <casa/OS/Timer.h>

namespace casa { //# NAMESPACE CASA - BEGIN

class UVWMachine;

class MultiTermFT : public FTMachine {
public:

  // Construct using an existing FT-Machine 
  MultiTermFT(FTMachine *subftm, String subFTMname, Int nterms=1, Double reffreq=0.0);

  // Construct from a Record containing the MultiTermFT state
  MultiTermFT(const RecordInterface& stateRec);

  // Copy constructor. 
  // This first calls the default "=" operator, and then instantiates objects for member pointers.
  MultiTermFT(const MultiTermFT &other);

  // Assignment operator --- leave it as the default
  MultiTermFT &operator=(const MultiTermFT &other);

  // Destructor
  ~MultiTermFT();

  // Called at the start of de-gridding : subftm->initializeToVis()
  // Note : Pre-de-gridding model-image divisions by PBs will go here.
  void initializeToVis(ImageInterface<Complex>& image, const VisBuffer& vb);

  // Called at the end of de-gridding : subftm->finalizeToVis()
  void finalizeToVis();

  // Called at the start of gridding : subftm->initializeToSky()
  void initializeToSky(ImageInterface<Complex>& image,  Matrix<Float>& weight,  const VisBuffer& vb);

  // Called at the end of gridding : subftm->finalizeToSky()
  void finalizeToSky();

  // Do the degridding via subftm->get() and modify model-visibilities by Taylor-weights
  void get(VisBuffer& vb, Int row=-1);

  // Modify imaging weights with Taylor-weights and do gridding via subftm->put()
  void put(VisBuffer& vb, Int row=-1, Bool dopsf=False,
           FTMachine::Type type=FTMachine::OBSERVED);

  // Have a const version for compatibility with other FTMs.. Throw an exception if called.
  void put(const VisBuffer& /*vb*/, Int /*row=-1*/, Bool /*dopsf=False*/,
           FTMachine::Type /*type=FTMachine::OBSERVED*/)
  {throw(AipsError("MultiTermFT::put called with a const vb. This FTM needs to modify the vb."));};

  // Calculate residual visibilities if possible.
  // The purpose is to allow rGridFT to make this multi-threaded
  virtual void ComputeResiduals(VisBuffer&vb, Bool useCorrected); 

  // Make an image : subftm->makeImage()
  void makeImage(FTMachine::Type type,
                 VisSet& vs,
                 ImageInterface<Complex>& image,
                 Matrix<Float>& weight);
  
  // Get the final image: do the Fourier transform grid-correct, then 
  // optionally normalize by the summed weights
  // Note : Post-gridding residual-image divisions by PBs will go here.
  //           For now, it just calls subftm->getImage()
  ImageInterface<Complex>& getImage(Matrix<Float>&, Bool normalize=True);

  // Place-holder for possible use with AWProject and AWProjectWB FTMs
  virtual void normalizeImage(Lattice<Complex>& /*skyImage*/,
                              const Matrix<Double>& /*sumOfWts*/,
                              Lattice<Float>& /*sensitivityImage*/,
                              Bool /*fftNorm*/)
    {throw(AipsError("MultiTermFT::normalizeImage() is not implemented"));}

  // Get the final weights image - this will hold PB2
  void getWeightImage(ImageInterface<Float>&, Matrix<Float>&);

  // Save and restore the MultiTermFT to and from a record
  virtual Bool toRecord(String& error, RecordInterface& outRec, Bool withImage=False, const String diskimage="");
  virtual Bool fromRecord(String& error, const RecordInterface& inRec);

  // Various small inline functions
  virtual Bool isFourier() {return True;}
  virtual void setNoPadding(Bool nopad){subftm_p->setNoPadding(nopad);};
  virtual String name() const {return machineName_p;};
  virtual void setMiscInfo(const Int qualifier){thisterm_p=qualifier;};

protected:

  // Multiply Imaging weights by Taylor-function weights - during "put"
  Bool modifyVisWeights(VisBuffer& vb);
  // Multiply model visibilities by Taylor-function weights - during "get"
  Bool modifyModelVis(VisBuffer &vb);
  // Restore vb.imagingweights to the original
  Bool restoreImagingWeights(VisBuffer &vb);
  // have to call the initmaps of subftm
  virtual void initMaps(const VisBuffer& vb);
  CountedPtr<FTMachine> subftm_p; 
  String  subFTMname_p;
  Int nterms_p;
  Int thisterm_p;
  Double reffreq_p;

  Matrix<Float> imweights_p;
  Double sumwt_p;
  String machineName_p;

  Bool dbg_p,dotime_p;
  Timer tmr_p;
  Double time_get, time_put, time_res;
};

} //# NAMESPACE CASA - END

#endif