casa5

//# PolAverageTVI_GTest:   test of polarization averaging TVIs

//# Copyright (C) 1995,1999,2000,2001,2016

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

//#

//# This program 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 2 of the License, or (at your option)

//# any later version.

//#

//# This program 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.

//#

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

//# with this program; if not, write to the Free Software Foundation, Inc.,

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

//#

//# Correspondence concerning AIPS++ should be addressed 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$

​

#include <mstransform/TVI/PolAverageTVI.h>

​

#include <casacore/casa/aips.h>

#include <casacore/casa/OS/EnvVar.h>

#include <casacore/casa/OS/File.h>

#include <casacore/casa/OS/RegularFile.h>

#include <casacore/casa/OS/SymLink.h>

#include <casacore/casa/OS/Directory.h>

#include <casacore/casa/OS/DirectoryIterator.h>

#include <casacore/casa/Exceptions/Error.h>

#include <casacore/casa/iostream.h>

#include <casacore/casa/Arrays/ArrayMath.h>

#include <casacore/casa/iomanip.h>

#include <casacore/casa/Containers/Record.h>

​

#include <msvis/MSVis/VisibilityIteratorImpl2.h>

#include <msvis/MSVis/LayeredVi2Factory.h>

​

#include <gtest/gtest.h>

#include <memory>

#include <typeinfo>

​

using namespace std;

using namespace casa;

using namespace casacore;

using namespace casa::vi;

​

namespace {

string GetCasaDataPath() {

  if (casacore::EnvironmentVariable::isDefined("CASAPATH")) {

    string casapath = casacore::EnvironmentVariable::get("CASAPATH");

    size_t endindex = casapath.find(" ");

    if (endindex != string::npos) {

      string casaroot = casapath.substr(0, endindex);

      cout << "casaroot = " << casaroot << endl;

      return (casaroot + "/data/");

    } else {

      cout << "hit npos" << endl;

      return "/data/";

    }

  } else {

    cout << "CASAPATH is not defined" << endl;

    return "";

  }

}

​

template<class T>

struct VerboseDeleterForNew {

  void operator()(T *p) {

    cout << "Destructing " << typeid(p).name() << endl;

    delete p;

  }

};

​

template<class T>

struct Filler {

  static void FillArrayReference(MeasurementSet const &ms,

      String const &columnName, Vector<uInt> const &rowIds, Array<T> &data) {

    //cout << "Start " << __func__ << endl;

    ArrayColumn<T> col(ms, columnName);

    RefRows rows(rowIds);

    col.getColumnCells(rows, data);

  }

  static void FillScalarReference(MeasurementSet const &ms,