--- /dev/null
+
+/*
+ * 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 _("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)
+{
+// register int i;
+// register 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)
+{
+// register int i;
+// register double alpha;
+// register 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;
+}
+
+
+
+
+
+