solve.cc File Reference

SED fitting routines. More...

#include <algorithm>
#include <ctime>
#include <thread>
#include <gsl/gsl_multifit_nlin.h>
#include "Graybody.h"
#include "PACS.h"
#include "SPIRE.h"

Include dependency graph for solve.cc:

Go to the source code of this file.

Classes
struct	manticore::GrayData
	Observation fitting data structure. More...

Namespaces
	manticore
	Package namespace.

Functions
std::string	manticore::asctime ()
	Current time a la std::asctime(). More...
int	manticore::df_gray (const gsl_vector *x, void *p, gsl_matrix *J)
	Single-temperature observation fitting function Jacobian. More...
int	manticore::df_gray2 (const gsl_vector *x, void *p, gsl_matrix *J)
	Two-temperature observation fitting function Jacobian. More...
int	manticore::f_gray (const gsl_vector *x, void *p, gsl_vector *f)
	Single-temperature observation fitting function. More...
int	manticore::f_gray2 (const gsl_vector *x, void *p, gsl_vector *f)
	Two-temperature observation fitting function. More...
int	manticore::fdf_gray (const gsl_vector *x, void *p, gsl_vector *f, gsl_matrix *J)
	Single-temperature observation fitting function plus Jacobian. More...
int	manticore::fdf_gray2 (const gsl_vector *x, void *p, gsl_vector *f, gsl_matrix *J)
	Two-temperature observation fitting function plus Jacobian. More...
mapDataType	manticore::getChi2 (const mu::CommandLine &cli)
	Read 2-param reduced \(\chi^2\) map. More...
void	manticore::makeModelMap (std::vector< mapDataType > &bandData, std::vector< mapDataType > &bandError, const GrayData &gdata)
	Replace input observations with theoretical model. More...
void	manticore::solve (std::vector< mapDataType > &outMap, std::vector< CCfits::ExtHDU * > &outHDU, CCfits::FITS *outFITS, const std::vector< mapDataType > &bandData, const std::vector< mapDataType > &bandError, const std::vector< CCfits::PHDU * > &bandHDU, const mu::CommandLine &cli)
	Performs least-squares fitting of band SEDs. More...
void	manticore::solveStage1 (const std::vector< mapDataType > &bandData, const std::vector< mapDataType > &bandError, const mapDataType &chi2Map, const mu::CommandLine &cli, const GrayData &gdata0, unsigned id, unsigned nThreads)
	Stage 1 solver. More...
void	manticore::solveStage2 (const std::vector< mapDataType > &bandData, const std::vector< mapDataType > &bandError, const mapDataType &chi2Map, const mu::CommandLine &cli, const GrayData &gdata0, unsigned id, unsigned nThreads)
	Stage 2 solver. More...

Variables
const double	manticore::imu_g = 1.0/mu_g
	Inverse mean molecular weight. More...
const double	manticore::mu_g = (2.75*m_H)
	Mean molecular weight. More...

Detailed Description

SED fitting routines.

Warning

This file assumes GSL version 2.1 or earlier. With v2.2, the non-linear least-squares API was completely rewritten!!

Author:  Kevin P. Rauch  <rauch.nosp@m.@ast.nosp@m.ro.um.nosp@m.d.ed.nosp@m.u>

Definition in file solve.cc.