#include <limits>
#include <stdexcept> /* std::domain_error(-) */
#include "igs_motion_wind_pixel.h"

namespace
{
template <class T>
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>
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;
	}
}
template <class T>
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);
	}
}
//-----------------------------------------------------------
template <class T, class R>
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();
	const double div_val = static_cast<double>(max_val);
	const double mul_val = div_val + 0.999999;

	const unsigned int ref_max_val = std::numeric_limits<R>::max();
	const double ref_div_val = static_cast<double>(ref_max_val);

	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 T, class R>
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();
	const double div_val = static_cast<double>(max_val);
	const double mul_val = div_val + 0.999999;

	const unsigned int ref_max_val = std::numeric_limits<R>::max();
	const double ref_div_val = static_cast<double>(ref_max_val);

	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);
	const int refbi = static_cast<int>(rb);
	if (std::numeric_limits<unsigned char>::digits == bi) {
		if (std::numeric_limits<unsigned short>::digits == refbi) {
			if ((0 == direction) || (2 == direction)) {
				change_template_hori_(
					in, ww, hh, cc, reinterpret_cast<const unsigned short *>(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, reinterpret_cast<const unsigned short *>(ref), rw, 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) {
		if (std::numeric_limits<unsigned short>::digits == refbi) {
			if ((0 == direction) || (2 == direction)) {
				change_template_hori_(
					reinterpret_cast<unsigned short *>(in), ww, hh, cc, reinterpret_cast<const unsigned short *>(ref), rw, rh, rc
					//,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 unsigned short *>(ref), rw, 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_(
					reinterpret_cast<unsigned short *>(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_(
					reinterpret_cast<unsigned short *>(in), ww, hh, cc, ref, rw, rh, rc
					//,rz
					,
					ref_mode, pix_rend, (1 == direction) ? false : true);
			}
		}
	} else {
		throw std::domain_error("Bad bits,Not uchar/ushort");
	}

	pix_rend.clear();
}