casa6

#ifndef SYNTHESIS_CPPSOLVERS_H

#define SYNTHESIS_CPPSOLVERS_H

​

#include <ms/MeasurementSets/MeasurementSet.h>

#include <casa/Arrays/Matrix.h>

#include <casa/Arrays/IPosition.h>

#include <images/Images/ImageInterface.h>

#include <images/Images/PagedImage.h>

#include <images/Images/TempImage.h>

​

#include <scimath/Mathematics/FFTServer.h>

#include <scimath/Functionals/Gaussian2D.h>

​

#include <gsl/gsl_sf_bessel.h>

#include <gsl/gsl_matrix.h>

#include <gsl/gsl_multimin.h>

#include <gsl/gsl_min.h>

#include <gsl/gsl_cblas.h>

#include <gsl/gsl_blas.h>

​

#include <Eigen/Core>

#include <chrono>

​

namespace casa { //# NAMESPACE CASA - BEGIN

​

class ParamObj

{

private:

  int nX;

  int nY;

  unsigned int AspLen;

  casacore::Matrix<casacore::Float> itsMatDirty;

  casacore::Matrix<casacore::Complex> itsPsfFT;

  std::vector<casacore::IPosition> center;

  //genie

  //Eigen::Matrix<float, Eigen::Dynamic, Eigen::Dynamic, Eigen::ColMajor> newResidual;

  casacore::Matrix<casacore::Float> newResidual;

  casacore::Matrix<casacore::Float> AspConvPsf;

  casacore::Matrix<casacore::Float> dAspConvPsf;

  casacore::FFTServer<casacore::Float,casacore::Complex> fft;

  casacore::Matrix<casacore::Float> Asp;

  casacore::Matrix<casacore::Float> dAsp;

​

public:

  ParamObj(const casacore::Matrix<casacore::Float>& dirty,

    const casacore::Matrix<casacore::Complex>& psf,

    const std::vector<casacore::IPosition>& positionOptimum) :

    itsMatDirty(dirty),

    itsPsfFT(psf),

    center(positionOptimum)/*,

    newResidual(dirty)*/

  {

    nX = itsMatDirty.shape()(0);

    nY = itsMatDirty.shape()(1);

    AspLen = center.size();

    //genie

    //newResidual = Eigen::MatrixXf::Zero(nX, nY);

    newResidual.resize(nX, nY);

    AspConvPsf.resize(nX, nY);

    dAspConvPsf.resize(nX, nY);

    fft = casacore::FFTServer<casacore::Float,casacore::Complex>(itsMatDirty.shape());

    Asp.resize(nX, nY);

    dAsp.resize(nX, nY);