casa6
Public
Source
casacore::Bool convertGaussWidth = k==PCFSpectralElement::WIDTH && type == SpectralElement::GAUSSIAN;//# SpectralFit2.cc: Least Squares fitting of spectral elements: templated part//# Copyright (C) 2001,2002,2004//# 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: SpectralFit2.tcc 19935 2007-02-27 05:07:40Z Malte.Marquarding $//# Includes#include <components/SpectralComponents/SpectralFit.h>#include <components/SpectralComponents/CompiledSpectralElement.h>#include <components/SpectralComponents/GaussianSpectralElement.h>#include <components/SpectralComponents/LogTransformedPolynomialSpectralElement.h>#include <components/SpectralComponents/LorentzianSpectralElement.h>#include <components/SpectralComponents/PolynomialSpectralElement.h>#include <components/SpectralComponents/PowerLogPolynomialSpectralElement.h>#include <casacore/scimath/Fitting/NonLinearFitLM.h>#include <casacore/scimath/Functionals/CompiledFunction.h>#include <casacore/scimath/Functionals/CompoundFunction.h>#include <casacore/scimath/Functionals/CompoundParam.h>#include <casacore/scimath/Functionals/Gaussian1D.h>#include <casacore/scimath/Functionals/Lorentzian1D.h>#include <casacore/scimath/Functionals/Polynomial.h>#include <casacore/scimath/Functionals/PowerLogarithmicPolynomial.h>namespace casa { //# NAMESPACE CASA - BEGIN//# Templated member functionstemplate <class MT>bool SpectralFit::fit(const casacore::Vector<MT> &y, const casacore::Vector<MT> &x, const casacore::Vector<casacore::Bool> *mask) { casacore::Vector<MT> sigma(x.nelements()); sigma = 1.0; return fit(sigma, y, x, mask);}template <class MT>bool SpectralFit::fit( const casacore::Vector<MT> &sigma, const casacore::Vector<MT> &y, const casacore::Vector<MT> &x, const casacore::Vector<casacore::Bool> *mask) {