manticore.h

Go to the documentation of this file.

#ifndef MANTICORE_H
#define MANTICORE_H

#include <CCfits/CCfits>
#include <gsl/gsl_spline.h>

#include "mutils/CommandLine.h"

namespace manticore {

using mapDataType = std::valarray<float>;

constexpr int mapDataCode = FLOAT_IMG;

constexpr long modelMapLen = 961L;

using tableType = std::vector<std::pair<double,double>>;

constexpr double c_light = 2.99792458e10;

constexpr double h_Planck = 6.626070040e-27;

constexpr double k_Boltzman = 1.38064852e-16;

constexpr double m_H = 1.6737237e-24;

extern mapDataType
  getError(const mapDataType &data, CCfits::PHDU &hdu,
           const mutils::CommandLine &cli);

extern gsl_spline
 *initSpline(const tableType &table, bool ln = false, unsigned stride = 1);

extern void
  process(const mutils::CommandLine &cli);

extern void
  setKeys(CCfits::ExtHDU &hdu, const CCfits::PHDU &hdu0,
          const mutils::CommandLine &cli, const std::vector<long> &ll);

extern void
  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 std::vector<long> &axis,
        const mutils::CommandLine &cli,
        const std::vector<long> &ll,
        const std::vector<long> &ur,
        const std::vector<long> &st);

extern int
  getBand(const CCfits::HDU &hdu);

}  // namespace manticore

#endif  // MANTICORE_H