#pragma once

#ifndef TCG_PIXEL_OPS_H
#define TCG_PIXEL_OPS_H

// tcg includes
#include "base.h"

namespace tcg {

//**************************************************************************
//    Pixel Categories
//**************************************************************************

struct grayscale_pixel_tag {};
struct rgb_pixel_tag {};
struct rgbm_pixel_tag : public rgb_pixel_tag {};
struct indexed_pixel_tag {};

//**************************************************************************
//    Pixel Traits Types
//**************************************************************************

template <typename Pix, typename pixel_category = typename Pix::pixel_category>
struct pixel_traits_types;

template <typename Pix>
struct pixel_traits_types<Pix, grayscale_pixel_tag> {
  typedef Pix pixel_type;
  typedef grayscale_pixel_tag pixel_category;
  typedef typename Pix::channel_type channel_type;

  enum { channels_count = 1 };
};

template <typename Pix>
struct pixel_traits_types<Pix, rgb_pixel_tag> {
  typedef Pix pixel_type;
  typedef rgb_pixel_tag pixel_category;
  typedef typename Pix::channel_type channel_type;

  enum { channels_count = 3 };
};

template <typename Pix>
struct pixel_traits_types<Pix, rgbm_pixel_tag>
    : public pixel_traits_types<Pix, rgb_pixel_tag> {
  typedef Pix pixel_type;
  typedef rgbm_pixel_tag pixel_category;
  typedef typename Pix::channel_type channel_type;

  enum { channels_count = 4 };
};

template <typename Pix>
struct pixel_traits_types<Pix, indexed_pixel_tag> {
  typedef Pix pixel_type;
  typedef indexed_pixel_tag pixel_category;
  typedef typename Pix::index_type index_type;
};

//**************************************************************************
//    Pixel Traits
//**************************************************************************

template <typename Pix, typename pixel_category = typename Pix::pixel_category>
struct pixel_traits;

template <typename Pix>
struct pixel_traits<Pix, grayscale_pixel_tag>
    : public pixel_traits_types<Pix, grayscale_pixel_tag> {
  typedef pixel_traits_types<Pix, grayscale_pixel_tag> tr;

public:
  static typename tr::channel_type max_channel_value();

  static typename tr::channel_type l(const typename tr::pixel_type &pix);
  static typename tr::channel_type &l(typename tr::pixel_type &pix);
};

template <typename Pix>
struct pixel_traits<Pix, rgb_pixel_tag> {
  typedef pixel_traits_types<Pix, rgb_pixel_tag> tr;

public:
  static typename tr::channel_type max_channel_value();

  static typename tr::channel_type r(const typename tr::pixel_type &pix);
  static typename tr::channel_type &r(typename tr::pixel_type &pix);

  static typename tr::channel_type g(const typename tr::pixel_type &pix);
  static typename tr::channel_type &g(typename tr::pixel_type &pix);

  static typename tr::channel_type b(const typename tr::pixel_type &pix);
  static typename tr::channel_type &b(typename tr::pixel_type &pix);
};

template <typename Pix>
struct pixel_traits<Pix, rgbm_pixel_tag>
    : public pixel_traits<Pix, rgb_pixel_tag> {
  typedef pixel_traits_types<Pix, rgbm_pixel_tag> tr;

public:
  static typename tr::channel_type m(const typename tr::pixel_type &pix);
  static typename tr::channel_type &m(typename tr::pixel_type &pix);
};

template <typename Pix>
struct pixel_traits<Pix, indexed_pixel_tag> {
  typedef pixel_traits_types<Pix, indexed_pixel_tag> tr;
  typedef typename tr::pixel_type pixel_type;
  typedef typename tr::index_type index_type;

public:
  static typename tr::index_type index(const typename tr::pixel_type &pix);
  static typename tr::index_type &index(typename tr::pixel_type &pix);
};

namespace pixel_ops {

//**************************************************************************
//    Pixel Functions
//**************************************************************************

template <typename PixIn, typename PixOut>
inline PixOut _cast(const PixIn &pix, grayscale_pixel_tag) {
  return PixOut(pixel_traits<PixIn>::l(pix));
}

template <typename PixIn, typename PixOut>
inline PixOut _cast(const PixIn &pix, rgb_pixel_tag) {
  return PixOut(pixel_traits<PixIn>::r(pix), pixel_traits<PixIn>::g(pix),
                pixel_traits<PixIn>::b(pix));
}

template <typename PixIn, typename PixOut>
inline PixOut _cast(const PixIn &pix, rgbm_pixel_tag) {
  return PixOut(pixel_traits<PixIn>::r(pix), pixel_traits<PixIn>::g(pix),
                pixel_traits<PixIn>::b(pix), pixel_traits<PixIn>::m(pix));
}

template <typename PixIn, typename PixOut>
inline PixOut cast(const PixIn &pix) {
  return _cast<PixIn, PixOut>(pix,
                              typename pixel_traits<PixIn>::pixel_category());
}

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

template <typename PixIn, typename PixOut>
inline void _assign(PixOut &pixout, const PixIn &pixin, grayscale_pixel_tag) {
  pixel_traits<PixOut>::l(pixout) = pixel_traits<PixIn>::l(pixin);
}

template <typename PixIn, typename PixOut>
inline void _assign(PixOut &pixout, const PixIn &pixin, rgb_pixel_tag) {
  pixel_traits<PixOut>::r(pixout) = pixel_traits<PixIn>::r(pixin);
  pixel_traits<PixOut>::g(pixout) = pixel_traits<PixIn>::g(pixin);
  pixel_traits<PixOut>::b(pixout) = pixel_traits<PixIn>::b(pixin);
}

template <typename PixIn, typename PixOut>
inline void _assign(PixOut &pixout, const PixIn &pixin, rgbm_pixel_tag) {
  pixel_traits<PixOut>::r(pixout) = pixel_traits<PixIn>::r(pixin);
  pixel_traits<PixOut>::g(pixout) = pixel_traits<PixIn>::g(pixin);
  pixel_traits<PixOut>::b(pixout) = pixel_traits<PixIn>::b(pixin);
  pixel_traits<PixOut>::m(pixout) = pixel_traits<PixIn>::m(pixin);
}

template <typename PixIn, typename PixOut>
inline void assign(PixOut &pixout, const PixIn &pixin) {
  _assign<PixIn, PixOut>(pixout, pixin,
                         typename pixel_traits<PixIn>::pixel_category());
}

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

template <typename Pix1, typename Pix2>
inline Pix1 _sum(const Pix1 &a, const Pix2 &b, grayscale_pixel_tag) {
  return Pix1(pixel_traits<Pix1>::l(a) + pixel_traits<Pix2>::l(b));
}

template <typename Pix1, typename Pix2>
inline Pix1 _sum(const Pix1 &a, const Pix2 &b, rgb_pixel_tag) {
  return Pix1(pixel_traits<Pix1>::r(a) + pixel_traits<Pix2>::r(b),
              pixel_traits<Pix1>::g(a) + pixel_traits<Pix2>::g(b),
              pixel_traits<Pix1>::b(a) + pixel_traits<Pix2>::b(b));
}

template <typename Pix1, typename Pix2>
inline Pix1 _sum(const Pix1 &a, const Pix2 &b, rgbm_pixel_tag) {
  return Pix1(pixel_traits<Pix1>::r(a) + pixel_traits<Pix2>::r(b),
              pixel_traits<Pix1>::g(a) + pixel_traits<Pix2>::g(b),
              pixel_traits<Pix1>::b(a) + pixel_traits<Pix2>::b(b),
              pixel_traits<Pix1>::m(a) + pixel_traits<Pix2>::m(b));
}

template <typename Pix1, typename Pix2>
inline Pix1 operator+(const Pix1 &a, const Pix2 &b) {
  return _sum(a, b, typename pixel_traits<Pix2>::pixel_category());
}

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

template <typename Pix1, typename Pix2>
inline Pix1 _sub(const Pix1 &a, const Pix2 &b, grayscale_pixel_tag) {
  return Pix1(pixel_traits<Pix1>::l(a) - pixel_traits<Pix2>::l(b));
}

template <typename Pix1, typename Pix2>
inline Pix1 _sub(const Pix1 &a, const Pix2 &b, rgb_pixel_tag) {
  return Pix1(pixel_traits<Pix1>::r(a) - pixel_traits<Pix2>::r(b),
              pixel_traits<Pix1>::g(a) - pixel_traits<Pix2>::g(b),
              pixel_traits<Pix1>::b(a) - pixel_traits<Pix2>::b(b));
}

template <typename Pix1, typename Pix2>
inline Pix1 _sub(const Pix1 &a, const Pix2 &b, rgbm_pixel_tag) {
  return Pix1(pixel_traits<Pix1>::r(a) - pixel_traits<Pix2>::r(b),
              pixel_traits<Pix1>::g(a) - pixel_traits<Pix2>::g(b),
              pixel_traits<Pix1>::b(a) - pixel_traits<Pix2>::b(b),
              pixel_traits<Pix1>::m(a) - pixel_traits<Pix2>::m(b));
}

template <typename Pix1, typename Pix2>
inline Pix1 operator-(const Pix1 &a, const Pix2 &b) {
  return _sub(a, b, typename pixel_traits<Pix2>::pixel_category());
}

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

template <typename Pix, typename Scalar>
inline Pix _mult(const Pix &a, Scalar k, grayscale_pixel_tag) {
  return Pix(pixel_traits<Pix>::l(a) * k);
}

template <typename Pix, typename Scalar>
inline Pix _mult(const Pix &a, Scalar k, rgb_pixel_tag) {
  return Pix(pixel_traits<Pix>::r(a) * k, pixel_traits<Pix>::g(a) * k,
             pixel_traits<Pix>::b(a) * k);
}

template <typename Pix, typename Scalar>
inline Pix _mult(const Pix &a, Scalar k, rgbm_pixel_tag) {
  return Pix(pixel_traits<Pix>::r(a) * k, pixel_traits<Pix>::g(a) * k,
             pixel_traits<Pix>::b(a) * k, pixel_traits<Pix>::m(a) * k);
}

template <typename Pix, typename Scalar>
inline Pix operator*(Scalar k, const Pix &pix) {
  return _mult(pix, k, typename pixel_traits<Pix>::pixel_category());
}

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

template <typename Pix, typename Scalar>
inline Pix _div(const Pix &a, Scalar k, grayscale_pixel_tag) {
  return Pix(pixel_traits<Pix>::l(a) / k);
}

template <typename Pix, typename Scalar>
inline Pix _div(const Pix &a, Scalar k, rgb_pixel_tag) {
  return Pix(pixel_traits<Pix>::r(a) / k, pixel_traits<Pix>::g(a) / k,
             pixel_traits<Pix>::b(a) / k);
}

template <typename Pix, typename Scalar>
inline Pix _div(const Pix &a, Scalar k, rgbm_pixel_tag) {
  return Pix(pixel_traits<Pix>::r(a) / k, pixel_traits<Pix>::g(a) / k,
             pixel_traits<Pix>::b(a) / k, pixel_traits<Pix>::m(a) / k);
}

template <typename Pix, typename Scalar>
inline Pix operator/(const Pix &pix, Scalar k) {
  return _div(pix, k, typename pixel_traits<Pix>::pixel_category());
}

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

template <typename Pix, typename Scalar>
inline void premultiply(Pix &pix, Scalar = 0) {
  Scalar factor =
      pixel_traits<Pix>::m(pix) / Scalar(pixel_traits<Pix>::max_channel_value);

  pixel_traits<Pix>::r(pix) = pixel_traits<Pix>::r(pix) * factor;
  pixel_traits<Pix>::g(pix) = pixel_traits<Pix>::g(pix) * factor;
  pixel_traits<Pix>::b(pix) = pixel_traits<Pix>::b(pix) * factor;
}

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

template <typename Pix, typename Scalar>
inline void depremultiply(Pix &pix, Scalar = 0) {
  if (pixel_traits<Pix>::channel_type m = pixel_traits<Pix>::m(pix)) {
    Scalar factor = pixel_traits<Pix>::max_channel_value / Scalar(m);

    pixel_traits<Pix>::r(pix) = pixel_traits<Pix>::r(pix) * factor;
    pixel_traits<Pix>::g(pix) = pixel_traits<Pix>::g(pix) * factor;
    pixel_traits<Pix>::b(pix) = pixel_traits<Pix>::b(pix) * factor;
  }
}

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

template <typename Pix, typename Scalar>
inline Pix _blend(const Pix &p0, const Pix &p1, Scalar t, grayscale_pixel_tag) {
  return Pix((1 - t) * pixel_traits<Pix>::l(p0) + t * pixel_traits<Pix>::l(p1));
}

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

template <typename Pix, typename Scalar>
inline Pix _blend(const Pix &p0, const Pix &p1, Scalar t, rgb_pixel_tag) {
  Scalar one_t = 1 - t;
  return Pix(one_t * pixel_traits<Pix>::r(p0) + t * pixel_traits<Pix>::r(p1),
             one_t * pixel_traits<Pix>::g(p0) + t * pixel_traits<Pix>::g(p1),
             one_t * pixel_traits<Pix>::m(p0) + t * pixel_traits<Pix>::b(p1));
}

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

template <typename Pix, typename Scalar>
inline Pix _blend(const Pix &p0, const Pix &p1, Scalar t, rgbm_pixel_tag) {
  Scalar one_t = 1 - t;
  return Pix(one_t * pixel_traits<Pix>::r(p0) + t * pixel_traits<Pix>::r(p1),
             one_t * pixel_traits<Pix>::g(p0) + t * pixel_traits<Pix>::g(p1),
             one_t * pixel_traits<Pix>::b(p0) + t * pixel_traits<Pix>::b(p1),
             one_t * pixel_traits<Pix>::m(p0) + t * pixel_traits<Pix>::m(p1));
}

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

template <typename Pix, typename Scalar>
inline Pix blend(const Pix &p0, const Pix &p1, Scalar t) {
  return _blend(p0, p1, t, pixel_traits<Pix>::pixel_category());
}

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

template <typename Pix, typename Scalar>
inline void over_premult(Pix &down, const Pix &up, Scalar = 0) {
  Scalar t =
      (1 -
       pixel_traits<Pix>::m(up) / Scalar(pixel_traits<Pix>::max_channel_value));

  pixel_traits<Pix>::r(down) =
      pixel_traits<Pix>::r(up) + t * pixel_traits<Pix>::r(down);
  pixel_traits<Pix>::g(down) =
      pixel_traits<Pix>::g(up) + t * pixel_traits<Pix>::g(down);
  pixel_traits<Pix>::b(down) =
      pixel_traits<Pix>::b(up) + t * pixel_traits<Pix>::b(down);
  pixel_traits<Pix>::m(down) =
      pixel_traits<Pix>::m(up) + t * pixel_traits<Pix>::m(down);
}

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

template <typename Pix, typename Scalar>
inline void over(Pix &down, const Pix &up, Scalar = 0) {
  Scalar t = (1 -
              pixel_traits<Pix>::m(up) /
                  Scalar(pixel_traits<Pix>::max_channel_value)) *
             pixel_traits<Pix>::m(down);
  Scalar m = pixel_traits<Pix>::m(up) + t;

  Scalar up_fac = pixel_traits<Pix>::m(up) / m;
  Scalar dn_fac = t / m;

  pixel_traits<Pix>::r(down) =
      up_fac * pixel_traits<Pix>::r(up) + dn_fac * pixel_traits<Pix>::r(down);
  pixel_traits<Pix>::g(down) =
      up_fac * pixel_traits<Pix>::g(up) + dn_fac * pixel_traits<Pix>::g(down);
  pixel_traits<Pix>::b(down) =
      up_fac * pixel_traits<Pix>::b(up) + dn_fac * pixel_traits<Pix>::b(down);
  pixel_traits<Pix>::m(down) = m;
}
}
}  // namespace tcg::pixel_ops

#endif  // TCG_PIXEL_OPS_H