Graybody.cc

Go to the documentation of this file.

#include "Detector.h"
#include "Graybody.h"

namespace mu = mutils;

namespace manticore {

namespace graybody {

extern "C" {

struct FluxData {
  double T,      
         Sigma;  

  const Detector *d;  

  const Graybody *g;  

  double (Graybody::*I)(double, double, double) const;
};

double flux(double nu, void *fd) {
  auto &p = *static_cast<const FluxData*>(fd);
  return (p.g->*(p.I))(nu, p.T, p.Sigma)*(p.d ? p.d->response(nu) : 1.0);
}

}  // extern "C"

}  // namescape graybody


double Graybody::Itot(double (Graybody::*I)(double, double, double) const,
                      double T, double Sigma, const Detector *detect,
                      double sr) const
{
  std::pair<double, double> freqRange =
    (detect ? detect->freqRange() :
              std::pair<double,double>{0.0,12.0*k_Boltzman*T/h_Planck});

  double result, abserr;
  graybody::FluxData fd = {T, Sigma, detect, this, I};
  gsl_function f = {graybody::flux, &fd};
  gsl_integration_qag(&f, freqRange.first, freqRange.second,
                      0.0, err_, subLimit, intKey,
                      integ_, &result, &abserr);
  return result*sr;
}

}  // namespace manticore