casa6

//# Flux.cc:

//# Copyright (C) 1998,1999,2000,2001,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 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: Flux.tcc 21292 2012-11-28 14:58:19Z gervandiepen $

​

#include <components/ComponentModels/Flux.h>

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

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

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

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

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

#include <casacore/casa/Logging/LogIO.h>

#include <casacore/casa/Logging/LogOrigin.h>

#include <casacore/casa/BasicSL/Complex.h>

#include <casacore/casa/BasicSL/Constants.h>

#include <casacore/measures/Measures/Stokes.h>

#include <casacore/casa/Quanta/Quantum.h>

#include <casacore/casa/Quanta/QuantumHolder.h>

#include <casacore/casa/Utilities/Assert.h>

#include <casacore/casa/Utilities/DataType.h>

#include <casacore/casa/BasicSL/String.h>

#include <iostream>

​

namespace casa { //# NAMESPACE CASA - BEGIN

​

template<class T> const std::vector<casacore::Unit> FluxRep<T>::_allowedUnits {"Jy", "Jy.km/s", "K.rad2"};

​

template<class T> FluxRep<T>::

FluxRep()

  :itsVal(4, typename casacore::NumericTraits<T>::ConjugateType(0,0)),

   itsPol(ComponentType::STOKES),

   itsUnit("Jy"),

   itsErr(4, typename casacore::NumericTraits<T>::ConjugateType(0,0))

{

  itsVal(0) = 1;

}

​

template<class T> FluxRep<T>::

FluxRep(T i) 

  :itsVal(4, typename casacore::NumericTraits<T>::ConjugateType(0,0)),

   itsPol(ComponentType::STOKES),

   itsUnit("Jy"),

   itsErr(4, typename casacore::NumericTraits<T>::ConjugateType(0,0))

{

  itsVal(0) = i;

}

​

template<class T> FluxRep<T>::

FluxRep(T i, T q, T u, T v)

  :itsVal(4, typename casacore::NumericTraits<T>::ConjugateType(0,0)),

   itsPol(ComponentType::STOKES),

   itsUnit("Jy"),

   itsErr(4, typename casacore::NumericTraits<T>::ConjugateType(0,0))

{

  itsVal(0) = i;

  itsVal(1) = q;

  itsVal(2) = u;

  itsVal(3) = v;

}

​

template<class T> FluxRep<T>::

FluxRep(typename casacore::NumericTraits<T>::ConjugateType xx,

    typename casacore::NumericTraits<T>::ConjugateType xy,

    typename casacore::NumericTraits<T>::ConjugateType yx,

    typename casacore::NumericTraits<T>::ConjugateType yy,

    ComponentType::Polarisation pol)

  :itsVal(4),

   itsPol(pol),

   itsUnit("Jy"),