/*
 * CINELERRA
 * Copyright (C) 2010 Adam Williams <broadcast at earthling dot net>
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 *
 */

#include "filexml.h"
#include "blur.h"
#include "blurwindow.h"
#include "bchash.h"
#include "keyframe.h"
#include "language.h"
#include "vframe.h"

#include <math.h>
#include <stdint.h>
#include <string.h>





#define MIN_RADIUS 2


BlurConfig::BlurConfig()
{
	vertical = 1;
	horizontal = 1;
	radius = 5;
	a_key = 0;
	a = r = g = b = 1;
}

int BlurConfig::equivalent(BlurConfig &that)
{
	return (vertical == that.vertical &&
		horizontal == that.horizontal &&
		radius == that.radius &&
		a_key == that.a_key &&
		a == that.a &&
		r == that.r &&
		g == that.g &&
		b == that.b);
}

void BlurConfig::copy_from(BlurConfig &that)
{
	vertical = that.vertical;
	horizontal = that.horizontal;
	radius = that.radius;
	a_key = that.a_key;
	a = that.a;
	r = that.r;
	g = that.g;
	b = that.b;
}

void BlurConfig::interpolate(BlurConfig &prev,
	BlurConfig &next,
	int64_t prev_frame,
	int64_t next_frame,
	int64_t current_frame)
{
	double next_scale = (double)(current_frame - prev_frame) / (next_frame - prev_frame);
	double prev_scale = (double)(next_frame - current_frame) / (next_frame - prev_frame);


//printf("BlurConfig::interpolate %d %d %d\n", prev_frame, next_frame, current_frame);
	this->vertical = prev.vertical;
	this->horizontal = prev.horizontal;
	this->radius = round(prev.radius * prev_scale + next.radius * next_scale);
	a_key = prev.a_key;
	a = prev.a;
	r = prev.r;
	g = prev.g;
	b = prev.b;
}






REGISTER_PLUGIN(BlurMain)








BlurMain::BlurMain(PluginServer *server)
 : PluginVClient(server)
{
	need_reconfigure = 1;
	engine = 0;
	overlayer = 0;
}

BlurMain::~BlurMain()
{
//printf("BlurMain::~BlurMain 1\n");

	if(engine)
	{
		for(int i = 0; i < (get_project_smp() + 1); i++)
			delete engine[i];
		delete [] engine;
	}

	if(overlayer) delete overlayer;
}

const char* BlurMain::plugin_title() { return N_("Blur"); }
int BlurMain::is_realtime() { return 1; }


NEW_WINDOW_MACRO(BlurMain, BlurWindow)

LOAD_CONFIGURATION_MACRO(BlurMain, BlurConfig)




int BlurMain::process_buffer(VFrame *frame,
	int64_t start_position,
	double frame_rate)
{
	int i;

	need_reconfigure |= load_configuration();

	read_frame(frame,
		0,
		start_position,
		frame_rate,
		0);

// Create temp based on alpha keying.
// Alpha keying needs 2x oversampling.

	if(config.a_key)
	{
		PluginVClient::new_temp(frame->get_w() * 2,
				frame->get_h() * 2,
				frame->get_color_model());
		if(!overlayer)
		{
			overlayer = new OverlayFrame(PluginClient::get_project_smp() + 1);
		}

		overlayer->overlay(PluginVClient::get_temp(),
			frame,
			0,
			0,
			frame->get_w(),
			frame->get_h(),
			0,
			0,
			PluginVClient::get_temp()->get_w(),
			PluginVClient::get_temp()->get_h(),
			1,        // 0 - 1
			TRANSFER_REPLACE,
			NEAREST_NEIGHBOR);
		input_frame = PluginVClient::get_temp();
	}
	else
	{
		PluginVClient::new_temp(frame->get_w(),
				frame->get_h(),
				frame->get_color_model());
		input_frame = frame;
	}


//printf("BlurMain::process_realtime 1 %d %d\n", need_reconfigure, config.radius);
	if(need_reconfigure)
	{
		if(!engine)
		{
			engine = new BlurEngine*[(get_project_smp() + 1)];
			for(i = 0; i < (get_project_smp() + 1); i++)
			{
				engine[i] = new BlurEngine(this);
				engine[i]->start();
			}
		}

		for(i = 0; i < (get_project_smp() + 1); i++)
		{
			engine[i]->reconfigure(&engine[i]->forward_constants, config.radius);
			engine[i]->reconfigure(&engine[i]->reverse_constants, config.radius);
		}
		need_reconfigure = 0;
	}


	if(config.radius < MIN_RADIUS ||
		(!config.vertical && !config.horizontal))
	{
// Data never processed
	}
	else
	{
// Process blur
// Need to blur vertically to a temp and
// horizontally to the output in 2 discrete passes.

		for(i = 0; i < get_project_smp() + 1; i++)
		{
			engine[i]->set_range(
				input_frame->get_h() * i / (get_project_smp() + 1),
				input_frame->get_h() * (i + 1) / (get_project_smp() + 1),
				input_frame->get_w() * i / (get_project_smp() + 1),
				input_frame->get_w() * (i + 1) / (get_project_smp() + 1));
		}

		for(i = 0; i < (get_project_smp() + 1); i++)
		{
			engine[i]->do_horizontal = 0;
			engine[i]->start_process_frame(input_frame);
		}

		for(i = 0; i < (get_project_smp() + 1); i++)
		{
			engine[i]->wait_process_frame();
		}

		for(i = 0; i < (get_project_smp() + 1); i++)
		{
			engine[i]->do_horizontal = 1;
			engine[i]->start_process_frame(input_frame);
		}

		for(i = 0; i < (get_project_smp() + 1); i++)
		{
			engine[i]->wait_process_frame();
		}
	}


// Downsample
	if(config.a_key)
	{
		overlayer->overlay(frame,
			PluginVClient::get_temp(),
			0,
			0,
			PluginVClient::get_temp()->get_w(),
			PluginVClient::get_temp()->get_h(),
			0,
			0,
			frame->get_w(),
			frame->get_h(),
			1,        // 0 - 1
			TRANSFER_REPLACE,
			NEAREST_NEIGHBOR);
	}

	return 0;
}


void BlurMain::update_gui()
{
	if(thread)
	{
		int reconfigure = load_configuration();
		if(reconfigure)
		{
			((BlurWindow*)thread->window)->lock_window("BlurMain::update_gui");
			((BlurWindow*)thread->window)->horizontal->update(config.horizontal);
			((BlurWindow*)thread->window)->vertical->update(config.vertical);
			((BlurWindow*)thread->window)->radius->update(config.radius);
			((BlurWindow*)thread->window)->radius_text->update((int64_t)config.radius);
			((BlurWindow*)thread->window)->a_key->update(config.a_key);
			((BlurWindow*)thread->window)->a->update(config.a);
			((BlurWindow*)thread->window)->r->update(config.r);
			((BlurWindow*)thread->window)->g->update(config.g);
			((BlurWindow*)thread->window)->b->update(config.b);
			((BlurWindow*)thread->window)->unlock_window();
		}
	}
}





void BlurMain::save_data(KeyFrame *keyframe)
{
	FileXML output;

// cause data to be stored directly in text
	output.set_shared_output(keyframe->get_data(), MESSAGESIZE);
	output.tag.set_title("BLUR");
	output.tag.set_property("VERTICAL", config.vertical);
	output.tag.set_property("HORIZONTAL", config.horizontal);
	output.tag.set_property("RADIUS", config.radius);
	output.tag.set_property("R", config.r);
	output.tag.set_property("G", config.g);
	output.tag.set_property("B", config.b);
	output.tag.set_property("A", config.a);
	output.tag.set_property("A_KEY", config.a_key);
	output.append_tag();
	output.tag.set_title("/BLUR");
	output.append_tag();
	output.append_newline();
	output.terminate_string();
}

void BlurMain::read_data(KeyFrame *keyframe)
{
	FileXML input;

	input.set_shared_input(keyframe->get_data(), strlen(keyframe->get_data()));

	int result = 0;

	while(!result)
	{
		result = input.read_tag();

		if(!result)
		{
			if(input.tag.title_is("BLUR"))
			{
				config.vertical = input.tag.get_property("VERTICAL", config.vertical);
				config.horizontal = input.tag.get_property("HORIZONTAL", config.horizontal);
				config.radius = input.tag.get_property("RADIUS", config.radius);
//printf("BlurMain::read_data 1 %d %d %s\n", get_source_position(), keyframe->position, keyframe->get_data());
				config.r = input.tag.get_property("R", config.r);
				config.g = input.tag.get_property("G", config.g);
				config.b = input.tag.get_property("B", config.b);
				config.a = input.tag.get_property("A", config.a);
				config.a_key = input.tag.get_property("A_KEY", config.a_key);
			}
		}
	}
}









BlurEngine::BlurEngine(BlurMain *plugin)
 : Thread(1, 0, 0)
{
	this->plugin = plugin;
	last_frame = 0;

// Strip size
	int size = plugin->get_input()->get_w() > plugin->get_input()->get_h() ?
		plugin->get_input()->get_w() : plugin->get_input()->get_h();
// Prepare for oversampling
	size *= 2;
	val_p = new pixel_f[size];
	val_m = new pixel_f[size];
	radius = new double[size];
	src = new pixel_f[size];
	dst = new pixel_f[size];

	set_synchronous(1);
	input_lock.lock();
	output_lock.lock();
}

BlurEngine::~BlurEngine()
{
	last_frame = 1;
	input_lock.unlock();
	join();
	delete [] val_p;
	delete [] val_m;
	delete [] src;
	delete [] dst;
	delete [] radius;
}

void BlurEngine::set_range(int start_y,
	int end_y,
	int start_x,
	int end_x)
{
	this->start_y = start_y;
	this->end_y = end_y;
	this->start_x = start_x;
	this->end_x = end_x;
}

int BlurEngine::start_process_frame(VFrame *frame)
{
	this->frame = frame;
	input_lock.unlock();
	return 0;
}

int BlurEngine::wait_process_frame()
{
	output_lock.lock();
	return 0;
}

void BlurEngine::run()
{
	int j, k;


	while(1)
	{
		input_lock.lock();
		if(last_frame)
		{
			output_lock.unlock();
			return;
		}





		color_model = frame->get_color_model();
		int w = frame->get_w();
		int h = frame->get_h();
// Force recalculation of filter
		prev_forward_radius = -65536;
		prev_reverse_radius = -65536;




#define BLUR(type, max, components) \
{ \
	type **input_rows = (type **)frame->get_rows(); \
	type **output_rows = (type **)frame->get_rows(); \
	type **current_input = input_rows; \
	type **current_output = output_rows; \
	vmax = max; \
 \
	if(!do_horizontal && plugin->config.vertical) \
	{ \
/* Vertical pass */ \
/* Render to temp if a horizontal pass comes next */ \
/*		if(plugin->config.horizontal) */ \
/*		{ */ \
/*			current_output = (type **)plugin->get_temp()->get_rows(); */ \
/*		} */ \
 \
		for(j = start_x; j < end_x; j++) \
		{ \
			bzero(val_p, sizeof(pixel_f) * h); \
			bzero(val_m, sizeof(pixel_f) * h); \
 \
			for(k = 0; k < h; k++) \
			{ \
				if(plugin->config.r) src[k].r = (double)current_input[k][j * components]; \
				if(plugin->config.g) src[k].g = (double)current_input[k][j * components + 1]; \
				if(plugin->config.b) src[k].b = (double)current_input[k][j * components + 2]; \
				if(components == 4) \
				{ \
					if(plugin->config.a || plugin->config.a_key) src[k].a = (double)current_input[k][j * components + 3]; \
				} \
			} \
 \
			if(components == 4) \
				blur_strip4(h); \
			else \
				blur_strip3(h); \
 \
			for(k = 0; k < h; k++) \
			{ \
				if(plugin->config.r) current_output[k][j * components] = (type)dst[k].r; \
				if(plugin->config.g) current_output[k][j * components + 1] = (type)dst[k].g; \
				if(plugin->config.b) current_output[k][j * components + 2] = (type)dst[k].b; \
				if(components == 4) \
					if(plugin->config.a || plugin->config.a_key) current_output[k][j * components + 3] = (type)dst[k].a; \
			} \
		} \
 \
/*		current_input = current_output; */ \
/*		current_output = output_rows; */ \
	} \
 \
 \
	if(do_horizontal && plugin->config.horizontal) \
	{ \
/* Vertical pass */ \
/*		if(plugin->config.vertical) */ \
/*		{ */ \
/*			current_input = (type **)plugin->get_temp()->get_rows(); */ \
/*		} */ \
 \
/* Horizontal pass */ \
		for(j = start_y; j < end_y; j++) \
		{ \
			bzero(val_p, sizeof(pixel_f) * w); \
			bzero(val_m, sizeof(pixel_f) * w); \
 \
			for(k = 0; k < w; k++) \
			{ \
				if(plugin->config.r) src[k].r = (double)current_input[j][k * components]; \
				if(plugin->config.g) src[k].g = (double)current_input[j][k * components + 1]; \
				if(plugin->config.b) src[k].b = (double)current_input[j][k * components + 2]; \
				if(components == 4) \
					if(plugin->config.a || plugin->config.a_key) src[k].a = (double)current_input[j][k * components + 3]; \
			} \
 \
 			if(components == 4) \
				blur_strip4(w); \
			else \
				blur_strip3(w); \
 \
			for(k = 0; k < w; k++) \
			{ \
				if(plugin->config.r) current_output[j][k * components] = (type)dst[k].r; \
				if(plugin->config.g) current_output[j][k * components + 1] = (type)dst[k].g; \
				if(plugin->config.b) current_output[j][k * components + 2] = (type)dst[k].b; \
				if(components == 4) \
				{ \
					if(plugin->config.a && !plugin->config.a_key) \
						current_output[j][k * components + 3] = (type)dst[k].a; \
					else if(plugin->config.a_key) \
						current_output[j][k * components + 3] = max; \
				} \
			} \
		} \
	} \
}



		switch(color_model)
		{
			case BC_RGB888:
			case BC_YUV888:
				BLUR(unsigned char, 0xff, 3);
				break;

			case BC_RGB_FLOAT:
				BLUR(float, 1.0, 3);
				break;

			case BC_RGBA8888:
			case BC_YUVA8888:
				BLUR(unsigned char, 0xff, 4);
				break;

			case BC_RGBA_FLOAT:
				BLUR(float, 1.0, 4);
				break;

			case BC_RGB161616:
			case BC_YUV161616:
				BLUR(uint16_t, 0xffff, 3);
				break;

			case BC_RGBA16161616:
			case BC_YUVA16161616:
				BLUR(uint16_t, 0xffff, 4);
				break;
		}

		output_lock.unlock();
	}
}

int BlurEngine::reconfigure(BlurConstants *constants, double radius)
{
// Blurring an oversampled temp
	if(plugin->config.a_key) radius *= 2;
	double std_dev = sqrt(-(double)(radius * radius) /
		(2 * log (1.0 / 255.0)));
	get_constants(constants, std_dev);
	return 0;
}

int BlurEngine::get_constants(BlurConstants *ptr, double std_dev)
{
	int i;
	double constants[8];
	double div;

	div = sqrt(2 * M_PI) * std_dev;
	constants[0] = -1.783 / std_dev;
	constants[1] = -1.723 / std_dev;
	constants[2] = 0.6318 / std_dev;
	constants[3] = 1.997  / std_dev;
	constants[4] = 1.6803 / div;
	constants[5] = 3.735 / div;
	constants[6] = -0.6803 / div;
	constants[7] = -0.2598 / div;

	ptr->n_p[0] = constants[4] + constants[6];
	ptr->n_p[1] = exp(constants[1]) *
				(constants[7] * sin(constants[3]) -
				(constants[6] + 2 * constants[4]) * cos(constants[3])) +
				exp(constants[0]) *
				(constants[5] * sin(constants[2]) -
				(2 * constants[6] + constants[4]) * cos(constants[2]));

	ptr->n_p[2] = 2 * exp(constants[0] + constants[1]) *
				((constants[4] + constants[6]) * cos(constants[3]) *
				cos(constants[2]) - constants[5] *
				cos(constants[3]) * sin(constants[2]) -
				constants[7] * cos(constants[2]) * sin(constants[3])) +
				constants[6] * exp(2 * constants[0]) +
				constants[4] * exp(2 * constants[1]);

	ptr->n_p[3] = exp(constants[1] + 2 * constants[0]) *
				(constants[7] * sin(constants[3]) -
				constants[6] * cos(constants[3])) +
				exp(constants[0] + 2 * constants[1]) *
				(constants[5] * sin(constants[2]) - constants[4] *
				cos(constants[2]));
	ptr->n_p[4] = 0.0;

	ptr->d_p[0] = 0.0;
	ptr->d_p[1] = -2 * exp(constants[1]) * cos(constants[3]) -
				2 * exp(constants[0]) * cos(constants[2]);

	ptr->d_p[2] = 4 * cos(constants[3]) * cos(constants[2]) *
				exp(constants[0] + constants[1]) +
				exp(2 * constants[1]) + exp (2 * constants[0]);

	ptr->d_p[3] = -2 * cos(constants[2]) * exp(constants[0] + 2 * constants[1]) -
				2 * cos(constants[3]) * exp(constants[1] + 2 * constants[0]);

	ptr->d_p[4] = exp(2 * constants[0] + 2 * constants[1]);

	for(i = 0; i < 5; i++) ptr->d_m[i] = ptr->d_p[i];

	ptr->n_m[0] = 0.0;
	for(i = 1; i <= 4; i++)
		ptr->n_m[i] = ptr->n_p[i] - ptr->d_p[i] * ptr->n_p[0];

	double sum_n_p, sum_n_m, sum_d;
	double a, b;

	sum_n_p = 0.0;
	sum_n_m = 0.0;
	sum_d = 0.0;
	for(i = 0; i < 5; i++)
	{
		sum_n_p += ptr->n_p[i];
		sum_n_m += ptr->n_m[i];
		sum_d += ptr->d_p[i];
	}

	a = sum_n_p / (1 + sum_d);
	b = sum_n_m / (1 + sum_d);

	for (i = 0; i < 5; i++)
	{
		ptr->bd_p[i] = ptr->d_p[i] * a;
		ptr->bd_m[i] = ptr->d_m[i] * b;
	}

	return 0;
}

#define BOUNDARY(x) if((x) > vmax) (x) = vmax; else if((x) < 0) (x) = 0;

int BlurEngine::transfer_pixels(pixel_f *src1,
	pixel_f *src2,
	pixel_f *src,
	double *radius,
	pixel_f *dest,
	int size)
{
	int i;
	double sum;

// printf("BlurEngine::transfer_pixels %d %d %d %d\n",
// plugin->config.r,
// plugin->config.g,
// plugin->config.b,
// plugin->config.a);

	for(i = 0; i < size; i++)
    {
		sum = src1[i].r + src2[i].r;
		BOUNDARY(sum);
		dest[i].r = sum;

		sum = src1[i].g + src2[i].g;
		BOUNDARY(sum);
		dest[i].g = sum;

		sum = src1[i].b + src2[i].b;
		BOUNDARY(sum);
		dest[i].b = sum;

		sum = src1[i].a + src2[i].a;
		BOUNDARY(sum);
		dest[i].a = sum;

// 		if(radius[i] < 2)
// 		{
// 			double scale = 2.0 - (radius[i] * radius[i] - 2.0);
// 			dest[i].r /= scale;
// 			dest[i].g /= scale;
// 			dest[i].b /= scale;
// 			dest[i].a /= scale;
// 		}
    }
	return 0;
}


int BlurEngine::multiply_alpha(pixel_f *row, int size)
{
//	int i;
//	double alpha;

// 	for(i = 0; i < size; i++)
// 	{
// 		alpha = (double)row[i].a / vmax;
// 		row[i].r *= alpha;
// 		row[i].g *= alpha;
// 		row[i].b *= alpha;
// 	}
	return 0;
}

int BlurEngine::separate_alpha(pixel_f *row, int size)
{
//	int i;
//	double alpha;
//	double result;

// 	for(i = 0; i < size; i++)
// 	{
// 		if(row[i].a > 0 && row[i].a < vmax)
// 		{
// 			alpha = (double)row[i].a / vmax;
// 			result = (double)row[i].r / alpha;
// 			row[i].r = (result > vmax ? vmax : result);
// 			result = (double)row[i].g / alpha;
// 			row[i].g = (result > vmax ? vmax : result);
// 			result = (double)row[i].b / alpha;
// 			row[i].b = (result > vmax ? vmax : result);
// 		}
// 	}
	return 0;
}

int BlurEngine::blur_strip3(int &size)
{
//	multiply_alpha(src, size);

	pixel_f *sp_p = src;
	pixel_f *sp_m = src + size - 1;
	pixel_f *vp = val_p;
	pixel_f *vm = val_m + size - 1;

	initial_p = sp_p[0];
	initial_m = sp_m[0];

	int l;
	for(int k = 0; k < size; k++)
	{
		terms = (k < 4) ? k : 4;

		radius[k] = plugin->config.radius;

		for(l = 0; l <= terms; l++)
		{
			if(plugin->config.r)
			{
				vp->r += forward_constants.n_p[l] * sp_p[-l].r - forward_constants.d_p[l] * vp[-l].r;
				vm->r += reverse_constants.n_m[l] * sp_m[l].r - reverse_constants.d_m[l] * vm[l].r;
			}
			if(plugin->config.g)
			{
				vp->g += forward_constants.n_p[l] * sp_p[-l].g - forward_constants.d_p[l] * vp[-l].g;
				vm->g += reverse_constants.n_m[l] * sp_m[l].g - reverse_constants.d_m[l] * vm[l].g;
			}
			if(plugin->config.b)
			{
				vp->b += forward_constants.n_p[l] * sp_p[-l].b - forward_constants.d_p[l] * vp[-l].b;
				vm->b += reverse_constants.n_m[l] * sp_m[l].b - reverse_constants.d_m[l] * vm[l].b;
			}
		}

		for( ; l <= 4; l++)
		{
			if(plugin->config.r)
			{
				vp->r += (forward_constants.n_p[l] - forward_constants.bd_p[l]) * initial_p.r;
				vm->r += (reverse_constants.n_m[l] - reverse_constants.bd_m[l]) * initial_m.r;
			}
			if(plugin->config.g)
			{
				vp->g += (forward_constants.n_p[l] - forward_constants.bd_p[l]) * initial_p.g;
				vm->g += (reverse_constants.n_m[l] - reverse_constants.bd_m[l]) * initial_m.g;
			}
			if(plugin->config.b)
			{
				vp->b += (forward_constants.n_p[l] - forward_constants.bd_p[l]) * initial_p.b;
				vm->b += (reverse_constants.n_m[l] - reverse_constants.bd_m[l]) * initial_m.b;
			}
		}
		sp_p++;
		sp_m--;
		vp++;
		vm--;
	}

	transfer_pixels(val_p, val_m, src, radius, dst, size);
//	separate_alpha(dst, size);
	return 0;
}


int BlurEngine::blur_strip4(int &size)
{

//	multiply_alpha(src, size);

	pixel_f *sp_p = src;
	pixel_f *sp_m = src + size - 1;
	pixel_f *vp = val_p;
	pixel_f *vm = val_m + size - 1;

	initial_p = sp_p[0];
	initial_m = sp_m[0];

	int l;
	for(int k = 0; k < size; k++)
	{
		if(plugin->config.a_key)
			radius[k] = (double)plugin->config.radius * src[k].a / vmax;
		else
			radius[k] = plugin->config.radius;
	}

	for(int k = 0; k < size; k++)
	{
		terms = (k < 4) ? k : 4;

		if(plugin->config.a_key)
		{
			if(radius[k] != prev_forward_radius)
			{
				prev_forward_radius = radius[k];
				if(radius[k] >= MIN_RADIUS / 2)
				{
					reconfigure(&forward_constants, radius[k]);
				}
				else
				{
					reconfigure(&forward_constants, (double)MIN_RADIUS / 2);
				}
			}

			if(radius[size - 1 - k] != prev_reverse_radius)
			{
//printf("BlurEngine::blur_strip4 %f\n", sp_m->a);
				prev_reverse_radius = radius[size - 1 - k];
				if(radius[size - 1 - k] >= MIN_RADIUS / 2)
				{
					reconfigure(&reverse_constants, radius[size - 1 - k]);
				}
				else
				{
					reconfigure(&reverse_constants, (double)MIN_RADIUS / 2);
				}
			}

// Force alpha to be copied regardless of alpha blur enabled
			vp->a = sp_p->a;
			vm->a = sp_m->a;
		}


		for(l = 0; l <= terms; l++)
		{
			if(plugin->config.r)
			{
				vp->r += forward_constants.n_p[l] * sp_p[-l].r - forward_constants.d_p[l] * vp[-l].r;
				vm->r += reverse_constants.n_m[l] * sp_m[l].r - reverse_constants.d_m[l] * vm[l].r;
			}
			if(plugin->config.g)
			{
				vp->g += forward_constants.n_p[l] * sp_p[-l].g - forward_constants.d_p[l] * vp[-l].g;
				vm->g += reverse_constants.n_m[l] * sp_m[l].g - reverse_constants.d_m[l] * vm[l].g;
			}
			if(plugin->config.b)
			{
				vp->b += forward_constants.n_p[l] * sp_p[-l].b - forward_constants.d_p[l] * vp[-l].b;
				vm->b += reverse_constants.n_m[l] * sp_m[l].b - reverse_constants.d_m[l] * vm[l].b;
			}
			if(plugin->config.a && !plugin->config.a_key)
			{
				vp->a += forward_constants.n_p[l] * sp_p[-l].a - forward_constants.d_p[l] * vp[-l].a;
				vm->a += reverse_constants.n_m[l] * sp_m[l].a - reverse_constants.d_m[l] * vm[l].a;
			}
		}

		for( ; l <= 4; l++)
		{
			if(plugin->config.r)
			{
				vp->r += (forward_constants.n_p[l] - forward_constants.bd_p[l]) * initial_p.r;
				vm->r += (reverse_constants.n_m[l] - reverse_constants.bd_m[l]) * initial_m.r;
			}
			if(plugin->config.g)
			{
				vp->g += (forward_constants.n_p[l] - forward_constants.bd_p[l]) * initial_p.g;
				vm->g += (reverse_constants.n_m[l] - reverse_constants.bd_m[l]) * initial_m.g;
			}
			if(plugin->config.b)
			{
				vp->b += (forward_constants.n_p[l] - forward_constants.bd_p[l]) * initial_p.b;
				vm->b += (reverse_constants.n_m[l] - reverse_constants.bd_m[l]) * initial_m.b;
			}
			if(plugin->config.a && !plugin->config.a_key)
			{
				vp->a += (forward_constants.n_p[l] - forward_constants.bd_p[l]) * initial_p.a;
				vm->a += (reverse_constants.n_m[l] - reverse_constants.bd_m[l]) * initial_m.a;
			}
		}

		sp_p++;
		sp_m--;
		vp++;
		vm--;
	}

	transfer_pixels(val_p, val_m, src, radius, dst, size);
//	separate_alpha(dst, size);
	return 0;
}