spectrum.h

#ifdef _MSC_VER
#pragma once
#endif

#ifndef _SPICA_SPECTRUM_H_
#define _SPICA_SPECTRUM_H_

#include <iostream>
#include <vector>

#include "common.h"

namespace spica {

    class SPICA_EXPORTS RGBSpectrum {
    public:
        RGBSpectrum();
        explicit RGBSpectrum(double l);
        
        RGBSpectrum(double red, double green, double blue);

        RGBSpectrum(const std::vector<double>& nm,
                    const std::vector<double>& values);

        RGBSpectrum(const RGBSpectrum& RGBSpectrum);
        ~RGBSpectrum();

        RGBSpectrum& operator=(const RGBSpectrum& c);
        bool operator==(const RGBSpectrum& c) const;
        bool operator!=(const RGBSpectrum& c) const;
        RGBSpectrum& operator+=(const RGBSpectrum& c);
        RGBSpectrum& operator+=(double x);
        RGBSpectrum& operator-=(const RGBSpectrum& c);
        RGBSpectrum& operator-=(double x);
        RGBSpectrum& operator*=(const RGBSpectrum& c);
        RGBSpectrum& operator*=(double s);
        RGBSpectrum& operator/=(const RGBSpectrum& c);
        RGBSpectrum& operator/=(double s);
        RGBSpectrum operator-() const;

        static RGBSpectrum fromXYZ(double x, double y, double z);

        static RGBSpectrum minimum(const RGBSpectrum& c1, const RGBSpectrum& c2);
        static RGBSpectrum maximum(const RGBSpectrum& c1, const RGBSpectrum& c2);

        static RGBSpectrum sqrt(const RGBSpectrum& c);
        static RGBSpectrum exp(const RGBSpectrum& c);
        static RGBSpectrum log(const RGBSpectrum& c);
        static RGBSpectrum clamp(const RGBSpectrum& c,
                                 const RGBSpectrum& lo = RGBSpectrum(0.0, 0.0, 0.0),
                                 const RGBSpectrum& hi = RGBSpectrum(INFTY, INFTY, INFTY));

        bool isBlack() const;
        bool isInf() const;
        bool isNaN() const;
        bool isValid() const;
        double dot(const RGBSpectrum& RGBSpectrum) const;
        double norm() const;
        double squaredNorm() const;
        double gray() const;

        inline double red()   const { return r_; }
        inline double green() const { return g_; }
        inline double blue()  const { return b_; }

        inline double operator[](int i) const {
            Assertion(i >= 0 && i <= 2, "Index out of range!!");
            if (i == 0) return r_;
            if (i == 1) return g_;
            if (i == 2) return b_;
            return 0.0;
        }

        inline double& ref(int i) {
            Assertion(i >= 0 && i <= 2, "Index out of range!!");
            if (i == 0) return r_;
            if (i == 1) return g_;
            return b_;            
        }

        inline RGBSpectrum toRGB() const {
            return { r_, g_, b_ };
        }

        std::string toString() const;

        // Public fields
        static const int channels = 3;

    private:
        double r_, g_, b_;
    };

    using Spectrum = RGBSpectrum;

}  // namespace spica

SPICA_EXPORTS spica::RGBSpectrum operator+(const spica::RGBSpectrum& c1, const spica::RGBSpectrum& c2);
SPICA_EXPORTS spica::RGBSpectrum operator+(double x, const spica::RGBSpectrum& c);
SPICA_EXPORTS spica::RGBSpectrum operator+(const spica::RGBSpectrum& c, double x);
SPICA_EXPORTS spica::RGBSpectrum operator-(const spica::RGBSpectrum& c1, const spica::RGBSpectrum& c2);
SPICA_EXPORTS spica::RGBSpectrum operator-(double x, const spica::RGBSpectrum& c);
SPICA_EXPORTS spica::RGBSpectrum operator-(const spica::RGBSpectrum& c, double x);
SPICA_EXPORTS spica::RGBSpectrum operator*(const spica::RGBSpectrum& c1, const spica::RGBSpectrum& c2);
SPICA_EXPORTS spica::RGBSpectrum operator*(const spica::RGBSpectrum& c1, double s);
SPICA_EXPORTS spica::RGBSpectrum operator*(double s, const spica::RGBSpectrum& c2);
SPICA_EXPORTS spica::RGBSpectrum operator/(const spica::RGBSpectrum& c1, const spica::RGBSpectrum& c2);
SPICA_EXPORTS spica::RGBSpectrum operator/(const spica::RGBSpectrum& c1, double s);
SPICA_EXPORTS std::ostream& operator<<(std::ostream& os, const spica::RGBSpectrum& c);

#endif  // SPICA_RGBSpectrum_H_