nonuniform_spectrum_sampling.h

Go to the documentation of this file.

#ifndef NONUNIFORM_SPECTRUM_SAMPLING_H
#define NONUNIFORM_SPECTRUM_SAMPLING_H
#include <boost/shared_ptr.hpp>
#include "spectrum_sampling.h"
#include "fp_exception.h"
#include <vector>

namespace FullPhysics {
/****************************************************************/
class NonuniformSpectrumSampling : public SpectrumSampling {
public:
  NonuniformSpectrumSampling(const SpectralDomain& Grid,
     const boost::shared_ptr<SpectrumSampling>& Interpolated_sampling);
  NonuniformSpectrumSampling(const SpectralDomain& Grid1,
     const SpectralDomain& Grid2,
     const SpectralDomain& Grid3,
     const boost::shared_ptr<SpectrumSampling>& Interpolated_sampling);

  virtual ~NonuniformSpectrumSampling() {}

  virtual SpectralDomain spectral_domain_interpolated(int Spec_index, 
                 const SpectralDomain& Lowres_grid, 
                 const DoubleWithUnit& Ils_half_width) const
  { return interpolated_sampling->
      spectral_domain(Spec_index, Lowres_grid, Ils_half_width); }
  virtual SpectralDomain spectral_domain(int spec_index,
                 const SpectralDomain& Lowres_grid, 
                 const DoubleWithUnit& Ils_half_width) const;

  virtual bool need_interpolation(int Spec_index) const
  { return (spec_domain[Spec_index].data().rows() > 0); }
  virtual void print(std::ostream& Os) const;
private:
  SpectralDomain sort_sd(const SpectralDomain& In) const;
  boost::shared_ptr<SpectrumSampling> interpolated_sampling;
  std::vector<SpectralDomain> spec_domain;
};
}
#endif