#include <limits></limits>

#include <stdexcept> /* std::domain_error(-) */</stdexcept>

#include "igs_motion_wind_pixel.h"

namespace {

template <class t=""></class>

const T *ref_pixel_(const T *ref, const int hh, const int ww, const int cc,

                    int xx, int yy) {

  while (xx < 0) {

    xx += ww;

  }

  while (ww <= xx) {

    xx -= ww;

  }

  while (yy < 0) {

    yy += hh;

  }

  while (hh <= yy) {

    yy -= hh;

  }

  return ref + cc * ww * yy + cc * xx;

}

template <class t=""></class>

void pixel_pop_(const T *in, const double div_val, const int cc, double *out) {

  /* 正規化し取り出し */

  for (int zz = 0; zz < cc; ++zz) {

    out[zz] = static_cast<double>(in[zz]) / div_val;</double>

  }

}

template <class t=""></class>

void pixel_push_(const double *in, const double mul_val, const int cc, T *out) {

  /* 値が変化したら格納する */

  for (int zz = 0; zz < cc; ++zz) {

    out[zz] = static_cast<t>(in[zz] * mul_val);</t>

  }

}

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

template <class class="" r="" t,=""></class>

void change_template_vert_(

    T *in, const int ww, const int hh, const int cc, const R *ref, const int rw,

    const int rh, const int rc

    //, const int rz

    ,

    const int ref_mode /* 0=R,1=G,2=B,3=A,4=Luminance,5=Nothing */

    ,

    igs::motion_wind::pixel &pix_rend, const bool reverse_sw) {

  const unsigned int max_val = std::numeric_limits<t>::max();</t>

  const double div_val       = static_cast<double>(max_val);</double>

  const double mul_val       = div_val + 0.999999;

  const unsigned int ref_max_val = std::numeric_limits<r>::max();</r>

  const double ref_div_val       = static_cast<double>(ref_max_val);</double>

  const int refprev = (reverse_sw) ? 1 : -1;

  for (int xx = 0; xx < ww; ++xx) {

    for (int yy = 0; yy < hh; ++yy) {

      T *next = in;

      if (reverse_sw) { /* 上から */

        next += (hh - 1 - yy) * ww * cc + xx * cc;

      } else { /* 下から */

        next += yy * ww * cc + xx * cc;

      }

      double pixel_in[igs::image::rgba::siz];

      pixel_pop_(next, div_val, cc, pixel_in);

      /****double pixel_ref[igs::image::rgba::siz];

double *rp_p = 0;

if (0 != ref) {

      pixel_pop_(

              ref_pixel_(ref,rh,rw,rc,xx,yy+refprev)

              , ref_div_val, rc, pixel_ref);

      rp_p = pixel_ref;

}

if (pix_rend.change(xx <= 0, rz, rp_p, cc, pixel_in)) {

      pixel_push_(pixel_in,mul_val,cc,next);

} ***/

      double ref_val = -1.0;

      if (0 != ref && 0 <= ref_mode && ref_mode <= 4) {

        ref_val =

            igs::color::ref_value(ref_pixel_(ref, rh, rw, rc, xx, yy + refprev),

                                  rc, ref_max_val, ref_mode);

      }

      if (pix_rend.change(xx <= 0, ref_val, cc, pixel_in)) {

        pixel_push_(pixel_in, mul_val, cc, next);

      }

    }

  }

}

template <class class="" r="" t,=""></class>

void change_template_hori_(

    T *in, const int ww, const int hh, const int cc, const R *ref, const int rw,

    const int rh, const int rc

    //, const int rz

    ,

    const int ref_mode /* 0=R,1=G,2=B,3=A,4=Luminance,5=Nothing */

    ,

    igs::motion_wind::pixel &pix_rend, const bool reverse_sw) {

  const unsigned int max_val = std::numeric_limits<t>::max();</t>

  const double div_val       = static_cast<double>(max_val);</double>

  const double mul_val       = div_val + 0.999999;

  const unsigned int ref_max_val = std::numeric_limits<r>::max();</r>

  const double ref_div_val       = static_cast<double>(ref_max_val);</double>

  const int refprev = (reverse_sw) ? 1 : -1;

  for (int yy = 0; yy < hh; ++yy) {

    T *next = in; /* 左から */

    if (reverse_sw) {

      next += (ww - 1) * cc;

    } /* 右から */

    for (int xx = 0; xx < ww; ++xx) {

      double pixel_in[igs::image::rgba::siz];

      pixel_pop_(next, div_val, cc, pixel_in);

      /***double pixel_ref[igs::image::rgba::siz];

double *rp_p = 0;

if (0 != ref) {

      pixel_pop_(

              ref_pixel_(ref,rh,rw,rc,xx+refprev,yy)

              , ref_div_val, rc, pixel_ref);

      rp_p = pixel_ref;

}

if (pix_rend.change(xx <= 0, rz, rp_p, cc, pixel_in)) {

      pixel_push_(pixel_in,mul_val,cc,next);

}***/

      double ref_val = -1.0;

      if (0 != ref && 0 <= ref_mode && ref_mode <= 4) {

        ref_val =

            igs::color::ref_value(ref_pixel_(ref, rh, rw, rc, xx, yy + refprev),

                                  rc, ref_max_val, ref_mode);

      }

      if (pix_rend.change(xx <= 0, ref_val, cc, pixel_in)) {

        pixel_push_(pixel_in, mul_val, cc, next);

      }

      if (reverse_sw) {

        next -= cc;

      } /* 右から左 */

      else {

        next += cc;

      } /* 左から右 */

    }

    in += ww * cc; /* 下から上 */

  }

}

}

/*-------------------------------------------------------

        Parameter

        方向			left-bottom is origin.

                0	:   0 : left to right

                1	:  90 : bottom to top

                2	: 180 : right to left

                3	: 270 : top to bottom

        暗い風へ変更するスイッチ

                true	: blow dark

                false	: blow light

        アルファにも風吹くスイッチ

                true	: blow alpha

                false	: not blow alpha

        風長さ(ランダム初期種)		1(0...)

        風長さ(最小値)			0(0...)

        風長さ(最大値)			18(0...)

        風長さ(最大最小のかたより値)	0.0(-1.0...0.0...1.0)

        風長さ(参照スイッチ)		false(true/false)

        風曲線(ランダム初期種)		1(0...)

        風曲線(最小値)			0(0...)

        風曲線(最大値)			1(0...)

        風曲線(最大最小のかたより値)	0.0(-1.0...0.0...1.0)

        風曲線(参照スイッチ)		false(true/false)

        風高さ(ランダム初期種)		1(0...)

        風高さ(最小値)			0(0...)

        風高さ(最大値)			1(0...)

        風高さ(最大最小のかたより値)	0.0(-1.0...0.0...1.0)

        風高さ(参照スイッチ)		false(true/false)

-------------------------------------------------------*/

#include "igs_motion_wind.h"

void igs::motion_wind::change(

    unsigned char *in

    ,

    const int hh, const int ww, const int cc, const int bb

    ,

    const unsigned char *ref

    ,

    const int rh, const int rw, const int rc, const int rb

    //, const int rz

    ,

    const int ref_mode /* 0=R,1=G,2=B,3=A,4=Luminance,5=Nothing */

    ,

    const int direction  // 0(0(LtoR),1(BtoT),2(RtoL),3(TtoB))

    ,

    const bool blow_dark_sw  // false(false,true)

    ,

    const bool blow_alpha_sw  // true(false,true)

    ,

    const unsigned long length_random_seed  // 0(0...)

    ,

    const double length_min  // 0.0...

    ,

    const double length_max  // 18.0...

    ,

    const double length_bias  // 1.0(0.0<...1.0...)

    ,

    const bool length_ref_sw  // false(false,true)

    ,

    const unsigned long force_random_seed  // 0(0...)

    ,

    const double force_min  // 0.0...

    ,

    const double force_max  // 1.0...

    ,

    const double force_bias  // 1.0(0.0<...1.0...)

    ,

    const bool force_ref_sw  // false(false,true)

    ,

    const unsigned long density_random_seed  // 0(0...)

    ,

    const double density_min  // 0.0...

    ,

    const double density_max  // 1.0...

    ,

    const double density_bias  // 1.0(0.0<...1.0...)

    ,

    const bool density_ref_sw  // false(false,true)

    ) {

  /* ユーザーによるノーアクションの指定、つまりゼロ長

  の時は処理せず終る */

  if ((length_max <= 0.0) && (length_min <= 0.0)) {

    return;

  }

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

  igs::motion_wind::pixel pix_rend(

      blow_dark_sw, blow_alpha_sw, length_random_seed, length_min, length_max,

      length_bias, length_ref_sw, force_random_seed, force_min, force_max,

      force_bias, force_ref_sw, density_random_seed, density_min, density_max,

      density_bias, density_ref_sw);

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

  const int bi    = static_cast<int>(bb);</int>

  const int refbi = static_cast<int>(rb);</int>

  if (std::numeric_limits<unsigned char="">::digits == bi) {</unsigned>

    if (std::numeric_limits<unsigned short="">::digits == refbi) {</unsigned>

      if ((0 == direction) || (2 == direction)) {

        change_template_hori_(

            in, ww, hh, cc, reinterpret_cast<const *="" short="" unsigned="">(ref), rw,</const>

            rh, rc

            //,rz

            ,

            ref_mode, pix_rend, (0 == direction) ? false : true);

      } else if ((1 == direction) || (3 == direction)) {

        change_template_vert_(

            in, ww, hh, cc, reinterpret_cast<const *="" short="" unsigned="">(ref), rw,</const>

            rh, rc

            //,rz

            ,

            ref_mode, pix_rend, (1 == direction) ? false : true);

      }

    } else { /* refbi is 8(uchar) or other(0...) */

      if ((0 == direction) || (2 == direction)) {

        change_template_hori_(in, ww, hh, cc, ref, rw, rh, rc

                              //,rz

                              ,

                              ref_mode, pix_rend,

                              (0 == direction) ? false : true);

      } else if ((1 == direction) || (3 == direction)) {

        change_template_vert_(in, ww, hh, cc, ref, rw, rh, rc

                              //,rz

                              ,

                              ref_mode, pix_rend,

                              (1 == direction) ? false : true);

      }

    }

  } else if (std::numeric_limits<unsigned short="">::digits == bi) {</unsigned>

    if (std::numeric_limits<unsigned short="">::digits == refbi) {</unsigned>

      if ((0 == direction) || (2 == direction)) {

        change_template_hori_(

            reinterpret_cast<unsigned *="" short="">(in), ww, hh, cc,</unsigned>

            reinterpret_cast<const *="" short="" unsigned="">(ref), rw, rh, rc</const>

            //,rz

            ,

            ref_mode, pix_rend, (0 == direction) ? false : true);

      } else if ((1 == direction) || (3 == direction)) {

        change_template_vert_(

            reinterpret_cast<unsigned *="" short="">(in), ww, hh, cc,</unsigned>

            reinterpret_cast<const *="" short="" unsigned="">(ref), rw, rh, rc</const>

            //,rz

            ,

            ref_mode, pix_rend, (1 == direction) ? false : true);

      }

    } else { /* refbi is 8(uchar) or other(0...) */

      if ((0 == direction) || (2 == direction)) {

        change_template_hori_(

            reinterpret_cast<unsigned *="" short="">(in), ww, hh, cc, ref, rw, rh, rc</unsigned>

            //,rz

            ,

            ref_mode, pix_rend, (0 == direction) ? false : true);

      } else if ((1 == direction) || (3 == direction)) {

        change_template_vert_(

            reinterpret_cast<unsigned *="" short="">(in), ww, hh, cc, ref, rw, rh, rc</unsigned>

            //,rz

            ,

            ref_mode, pix_rend, (1 == direction) ? false : true);

      }

    }

  } else {

    throw std::domain_error("Bad bits,Not uchar/ushort");

  }

  pix_rend.clear();

}