Dust.h

Go to the documentation of this file.

#ifndef DUST_H
#define DUST_H

#include <gsl/gsl_spline.h>

#include "mutils/CommandLine.h"
#include "mutils/util.h"

namespace manticore {

namespace dust {
  extern void printModels(std::vector<mutils::CommandLine::Option> &opts,
                          char shortName);
}


class Dust {
public:
  Dust(const std::string &model = "OH5", double rho = 100.0)
  { setModel(model, rho); }

  ~Dust() {
    if (spline_) { gsl_spline_free(spline_); }
  }

  Dust(double nu0, double kappa0, double beta)
  { setModel(nu0, kappa0, beta); }

  void setModel(const std::string &model, double rho = 100.0);

  void setModel(double nu0, double kappa0, double beta);

  double kappa(double nu, gsl_interp_accel *acc = nullptr) const;

  const std::string &name() const noexcept { return name_; }

  double rho() const noexcept { return rho_; }

protected:
  struct Opacity { double inu0, kappa0, beta; };

  std::string name_;

  double rho_;

  Opacity opacity_ = {0.0, 0.0, 0.0};

  gsl_spline *spline_ = nullptr;

  bool ln_ = false;
};

}  // namespace manticore

#endif  // DUST_H