#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, rh, rc</const>

					//,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, 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_(

					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, reinterpret_cast<const *="" short="" unsigned="">(ref), rw, rh, rc</const></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, reinterpret_cast<const *="" short="" unsigned="">(ref), rw, rh, rc</const></unsigned>

					//,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();

}