spectral_range.h

Go to the documentation of this file.

#ifndef SPECTRAL_RANGE_H
#define SPECTRAL_RANGE_H
#include "printable.h"
#include "unit.h"
#include "array_with_unit.h"
#include "array_ad.h"

namespace FullPhysics {
/****************************************************************/
class SpectralRange: public Printable<SpectralRange> {
public:
  SpectralRange(const ArrayWithUnit<double, 1>& Data)
    : data_(Data.value), units_(Data.units) {}
  SpectralRange(const ArrayAd<double, 1>& Data, 
                const Unit& U)
    : data_(Data), units_(U) {}
  SpectralRange(const blitz::Array<double, 1>& Data, 
                const Unit& U)
    : data_(Data), units_(U) {}
  SpectralRange(const blitz::Array<double, 1>& Data, 
                const Unit& U,
                const blitz::Array<double, 1>& Uncertainty)
    : data_(Data), units_(U), 
      uncertainty_(Uncertainty) {}
  SpectralRange(const ArrayAd<double, 1>& Data, 
                const Unit& U,
                const blitz::Array<double, 1>& Uncertainty)
    : data_(Data), units_(U), 
      uncertainty_(Uncertainty) {}

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  
  SpectralRange& operator=(const SpectralRange& V)
  { 
    data_.reference(V.data_ad().copy());
    units_ = V.units();
    uncertainty_.reference(V.uncertainty().copy());
    return *this;
  }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  const blitz::Array<double, 1>& data() const {return data_.value(); }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  blitz::Array<double, 1>& data() {return data_.value(); }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  const ArrayAd<double, 1>& data_ad() const {return data_; }
  

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  ArrayAd<double, 1>& data_ad() {return data_; }
  
//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  const blitz::Array<double, 1>& uncertainty() const {return uncertainty_;}

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  const SpectralRange clone() const { return SpectralRange(data_.copy(), units_, uncertainty_.copy()); }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  const Unit& units() const {return units_;}
  SpectralRange convert(const Unit& R) const;
  void print(std::ostream& Os) const { Os << "SpectralRange";}

  SpectralRange() {}

private:
  ArrayAd<double, 1> data_;
  Unit units_;
  blitz::Array<double, 1> uncertainty_;
};
}
#endif