From f8d8ce4d35bfbba5391af2a3fbce2ad88f94dd69 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Juha=20Vehvil=C3=A4inen?= <jusu@users.sourceforge.net>
Date: Tue, 9 Jul 2002 12:36:46 +0000
Subject: PeRColate

svn path=/trunk/Framestein/; revision=31
---
 Plugins/rgbavg02.c | 799 +++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 799 insertions(+)
 create mode 100644 Plugins/rgbavg02.c

(limited to 'Plugins/rgbavg02.c')

diff --git a/Plugins/rgbavg02.c b/Plugins/rgbavg02.c
new file mode 100644
index 0000000..5fc773e
--- /dev/null
+++ b/Plugins/rgbavg02.c
@@ -0,0 +1,799 @@
+// 242.rgbavg02 -- does 2-source funky pixel averaging.
+//		by r. luke dubois (luke@music.columbia.edu),
+//			computer music center, columbia university, 2001.
+//
+//  Pd / Framestein port by Olaf Matthes <olaf.matthes@gmx.de>, June 2002
+//
+//	objects and source are provided without warranty of any kind, express or implied.
+//
+//  usage: rgbavg <red> <green> <blue> <mode>
+//
+
+#include <stdio.h>
+#include <string.h>
+#include <math.h>
+#include "plugin.h"
+
+#pragma warning( disable : 4761 )	// that's why it's funky !!
+
+void perform_effect(struct frame f1, struct frame f2, struct args a)
+{
+	printf("Using rgbavg as effect does nothing!\n");
+}
+
+void perform_copy(struct frame f1, struct frame f2, struct args a)
+{
+	short x,y, h, w;
+	long red, green, blue;
+	short alpha, mode;
+	long redpix, greenpix, bluepix, rp, gp, bp;
+	pixel16 *c1_16, *c2_16;
+	pixel24 *c1_24, *c2_24;
+	pixel32 *c1_32, *c2_32;
+	char *t;
+
+	// get params
+	if(!a.s) return;
+	red = atoi(a.s);
+	if(!(t = strstr(a.s, " "))) return;
+	green = atoi(t+1);
+	if(!(t = strstr(t+1, " "))) return;
+	blue = atoi(t+1);
+	if(!(t = strstr(t+1, " "))) return;
+	mode = atoi(t+1);	// mode=0 = bypass
+
+	printf("rgbavg: %d %d %d - %d\n", red, green, blue, mode);
+
+	w = f1.width<f2.width ? f1.width : f2.width;
+	h = f1.height<f2.height ? f1.height : f2.height;
+
+	// perform routines do different pixel averagings in the main loop
+	switch (mode)
+	{
+		case 1:	// output: color / (sum of other two colors)
+			switch (f1.pixelformat)
+			{
+				case 16:
+					for(y = 0; y < h; y++) 
+					{
+						c1_16 = scanline16(f1, y);
+						c2_16 = scanline16(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r16(c1_16[x]) + red;
+							gp = g16(c1_16[x]) + green;
+							bp = b16(c1_16[x]) + blue;
+							redpix = r16(c2_16[x]) + red;
+							greenpix = g16(c2_16[x]) + green;
+							bluepix = b16(c2_16[x]) + blue;
+							c2_16[x] = rgbtocolor16(klamp255(redpix/(gp+bp)), 
+								                    klamp255(greenpix/(rp+bp)), 
+											        klamp255(bluepix/(rp+gp)));
+						}
+					}
+					break;
+				case 24:
+					for(y = 0; y < h; y++) 
+					{
+						c1_24 = scanline24(f1, y);
+						c2_24 = scanline24(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r24(c1_24[x]) + red;
+							gp = g24(c1_24[x]) + green;
+							bp = b24(c1_24[x]) + blue;
+							redpix = r24(c2_24[x]) + red;
+							greenpix = g24(c2_24[x]) + green;
+							bluepix = b24(c2_24[x]) + blue;
+							c2_24[x] = rgbtocolor24(klamp255(redpix/(gp+bp)), 
+								                    klamp255(greenpix/(rp+bp)), 
+											        klamp255(bluepix/(rp+gp)));
+						}
+					}
+					break;
+				case 32:
+					for(y = 0; y < h; y++) 
+					{
+						c1_32 = scanline32(f1, y);
+						c2_32 = scanline32(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r32(c1_32[x]) + red;
+							gp = g32(c1_32[x]) + green;
+							bp = b32(c1_32[x]) + blue;
+							redpix = r32(c2_32[x]) + red;
+							greenpix = g32(c2_32[x]) + green;
+							bluepix = b32(c2_32[x]) + blue;
+							c2_32[x] = rgbtocolor32(klamp255(redpix/(gp+bp)), 
+								                    klamp255(greenpix/(rp+bp)), 
+											        klamp255(bluepix/(rp+gp)));
+						}
+					}
+					break;
+			}
+			break;
+		case 2:	// output: (sum of two colors) / third color
+			switch (f1.pixelformat)
+			{
+				case 16:
+					for(y = 0; y < h; y++) 
+					{
+						c1_16 = scanline16(f1, y);
+						c2_16 = scanline16(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r16(c1_16[x]) + red;
+							gp = g16(c1_16[x]) + green;
+							bp = b16(c1_16[x]) + blue;
+							redpix = r16(c2_16[x]) + red;
+							greenpix = g16(c2_16[x]) + green;
+							bluepix = b16(c2_16[x]) + blue;
+							c2_16[x] = rgbtocolor16(klamp255((gp+bp)/redpix), 
+								                    klamp255((rp+bp)/greenpix), 
+											        klamp255((rp+gp)/bluepix));
+						}
+					}
+					break;
+				case 24:
+					for(y = 0; y < h; y++) 
+					{
+						c1_24 = scanline24(f1, y);
+						c2_24 = scanline24(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r24(c1_24[x]) + red;
+							gp = g24(c1_24[x]) + green;
+							bp = b24(c1_24[x]) + blue;
+							redpix = r24(c2_24[x]) + red;
+							greenpix = g24(c2_24[x]) + green;
+							bluepix = b24(c2_24[x]) + blue;
+							c2_24[x] = rgbtocolor24(klamp255((gp+bp)/redpix), 
+								                    klamp255((rp+bp)/greenpix), 
+											        klamp255((rp+gp)/bluepix));
+						}
+					}
+					break;
+				case 32:
+					for(y = 0; y < h; y++) 
+					{
+						c1_32 = scanline32(f1, y);
+						c2_32 = scanline32(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r32(c1_32[x]) + red;
+							gp = g32(c1_32[x]) + green;
+							bp = b32(c1_32[x]) + blue;
+							redpix = r32(c2_32[x]) + red;
+							greenpix = g32(c2_32[x]) + green;
+							bluepix = b32(c2_32[x]) + blue;
+							c2_32[x] = rgbtocolor32(klamp255((gp+bp)/redpix), 
+								                    klamp255((rp+bp)/greenpix), 
+											        klamp255((rp+gp)/bluepix));
+						}
+					}
+					break;
+			}
+			break;
+		case 3:	// output: color / (sum of other two colors) + color(input)
+			switch (f1.pixelformat)
+			{
+				case 16:
+					for(y = 0; y < h; y++) 
+					{
+						c1_16 = scanline16(f1, y);
+						c2_16 = scanline16(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r16(c1_16[x]) + 1;
+							gp = g16(c1_16[x]) + 1;
+							bp = b16(c1_16[x]) + 1;
+							redpix = r16(c2_16[x]) + 1;
+							greenpix = g16(c2_16[x]) + 1;
+							bluepix = b16(c2_16[x]) + 1;
+							c2_16[x] = rgbtocolor16(klamp255(red + ((gp+bp)/redpix)), 
+								                    klamp255(green + ((rp+bp)/greenpix)),
+													klamp255(blue + ((rp+gp)/bluepix)));
+						}
+					}
+					break;
+				case 24:
+					for(y = 0; y < h; y++) 
+					{
+						c1_24 = scanline24(f1, y);
+						c2_24 = scanline24(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r24(c1_24[x]) + 1;
+							gp = g24(c1_24[x]) + 1;
+							bp = b24(c1_24[x]) + 1;
+							redpix = r24(c2_24[x]) + 1;
+							greenpix = g24(c2_24[x]) + 1;
+							bluepix = b24(c2_24[x]) + 1;
+							c2_24[x] = rgbtocolor24(klamp255(red + ((gp+bp)/redpix)), 
+								                    klamp255(green + ((rp+bp)/greenpix)),
+													klamp255(blue + ((rp+gp)/bluepix)));
+						}
+					}
+					break;
+				case 32:
+					for(y = 0; y < h; y++) 
+					{
+						c1_32 = scanline32(f1, y);
+						c2_32 = scanline32(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r32(c1_32[x]) + 1;
+							gp = g32(c1_32[x]) + 1;
+							bp = b32(c1_32[x]) + 1;
+							redpix = r32(c2_32[x]) + 1;
+							greenpix = g32(c2_32[x]) + 1;
+							bluepix = b32(c2_32[x]) + 1;
+							c2_32[x] = rgbtocolor32(klamp255(red + ((gp+bp)/redpix)), 
+								                    klamp255(green + ((rp+bp)/greenpix)),
+													klamp255(blue + ((rp+gp)/bluepix)));
+						}
+					}
+					break;
+			}
+			break;
+		case 4:
+			switch (f1.pixelformat)
+			{
+				case 16:
+					for(y = 0; y < h; y++) 
+					{
+						c1_16 = scanline16(f1, y);
+						c2_16 = scanline16(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r16(c1_16[x]) + 1;
+							gp = g16(c1_16[x]) + 1;
+							bp = b16(c1_16[x]) + 1;
+							redpix = r16(c2_16[x]) + 1;
+							greenpix = g16(c2_16[x]) + 1;
+							bluepix = b16(c2_16[x]) + 1;
+							c2_16[x] = rgbtocolor16(klamp255(((gp*green)+(bp*blue))/redpix), 
+								                    klamp255(((rp*red)+(bp*blue))/greenpix), 
+											        klamp255(((rp*red)+(gp*green))/bluepix));
+						}
+					}
+					break;
+				case 24:
+					for(y = 0; y < h; y++) 
+					{
+						c1_24 = scanline24(f1, y);
+						c2_24 = scanline24(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r24(c1_24[x]) + 1;
+							gp = g24(c1_24[x]) + 1;
+							bp = b24(c1_24[x]) + 1;
+							redpix = r24(c2_24[x]) + 1;
+							greenpix = g24(c2_24[x]) + 1;
+							bluepix = b24(c2_24[x]) + 1;
+							c2_24[x] = rgbtocolor24(klamp255(((gp*green)+(bp*blue))/redpix), 
+								                    klamp255(((rp*red)+(bp*blue))/greenpix), 
+											        klamp255(((rp*red)+(gp*green))/bluepix));
+						}
+					}
+					break;
+				case 32:
+					for(y = 0; y < h; y++) 
+					{
+						c1_32 = scanline32(f1, y);
+						c2_32 = scanline32(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r32(c1_32[x]) + 1;
+							gp = g32(c1_32[x]) + 1;
+							bp = b32(c1_32[x]) + 1;
+							redpix = r32(c2_32[x]) + 1;
+							greenpix = g32(c2_32[x]) + 1;
+							bluepix = b32(c2_32[x]) + 1;
+							c2_32[x] = rgbtocolor32(klamp255(((gp*green)+(bp*blue))/redpix), 
+								                    klamp255(((rp*red)+(bp*blue))/greenpix), 
+											        klamp255(((rp*red)+(gp*green))/bluepix));
+						}
+					}
+					break;
+			}
+			break;
+		case 5:
+			switch (f1.pixelformat)
+			{
+				case 16:
+					for(y = 0; y < h; y++) 
+					{
+						c1_16 = scanline16(f1, y);
+						c2_16 = scanline16(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r16(c1_16[x]) + 1;
+							gp = g16(c1_16[x]) + 1;
+							bp = b16(c1_16[x]) + 1;
+							redpix = r16(c2_16[x]) + 1;
+							greenpix = g16(c2_16[x]) + 1;
+							bluepix = b16(c2_16[x]) + 1;
+							c2_16[x] = rgbtocolor16(klamp255(((gp*green)+(bp*blue))/(redpix*red)), 
+								                    klamp255(((rp*red)+(bp*blue))/(greenpix*green)), 
+											        klamp255(((rp*red)+(gp*green))/(bluepix*blue)));
+						}
+					}
+					break;
+				case 24:
+					for(y = 0; y < h; y++) 
+					{
+						c1_24 = scanline24(f1, y);
+						c2_24 = scanline24(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r24(c1_24[x]) + 1;
+							gp = g24(c1_24[x]) + 1;
+							bp = b24(c1_24[x]) + 1;
+							redpix = r24(c2_24[x]) + 1;
+							greenpix = g24(c2_24[x]) + 1;
+							bluepix = b24(c2_24[x]) + 1;
+							c2_24[x] = rgbtocolor24(klamp255(((gp*green)+(bp*blue))/(redpix*red)), 
+								                    klamp255(((rp*red)+(bp*blue))/(greenpix*green)), 
+											        klamp255(((rp*red)+(gp*green))/(bluepix*blue)));
+						}
+					}
+					break;
+				case 32:
+					for(y = 0; y < h; y++) 
+					{
+						c1_32 = scanline32(f1, y);
+						c2_32 = scanline32(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r32(c1_32[x]) + 1;
+							gp = g32(c1_32[x]) + 1;
+							bp = b32(c1_32[x]) + 1;
+							redpix = r32(c2_32[x]) + 1;
+							greenpix = g32(c2_32[x]) + 1;
+							bluepix = b32(c2_32[x]) + 1;
+							c2_32[x] = rgbtocolor32(klamp255(((gp*green)+(bp*blue))/(redpix*red)), 
+								                    klamp255(((rp*red)+(bp*blue))/(greenpix*green)), 
+											        klamp255(((rp*red)+(gp*green))/(bluepix*blue)));
+						}
+					}
+					break;
+			}
+			break;
+		case 6:
+			switch (f1.pixelformat)
+			{
+				case 16:
+					for(y = 0; y < h; y++) 
+					{
+						c1_16 = scanline16(f1, y);
+						c2_16 = scanline16(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r16(c1_16[x]) + 1;
+							gp = g16(c1_16[x]) + 1;
+							bp = b16(c1_16[x]) + 1;
+							redpix = r16(c2_16[x]) + 1;
+							greenpix = g16(c2_16[x]) + 1;
+							bluepix = b16(c2_16[x]) + 1;
+							c2_16[x] = rgbtocolor16(klamp255((redpix*red)/(gp*green)+(bp*blue)), 
+								                    klamp255((greenpix*green)/(rp*red)+(bp*blue)), 
+											        klamp255((bluepix*blue)/(rp*red)+(gp*green)));
+						}
+					}
+					break;
+				case 24:
+					for(y = 0; y < h; y++) 
+					{
+						c1_24 = scanline24(f1, y);
+						c2_24 = scanline24(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r24(c1_24[x]) + 1;
+							gp = g24(c1_24[x]) + 1;
+							bp = b24(c1_24[x]) + 1;
+							redpix = r24(c2_24[x]) + 1;
+							greenpix = g24(c2_24[x]) + 1;
+							bluepix = b24(c2_24[x]) + 1;
+							c2_24[x] = rgbtocolor24(klamp255((redpix*red)/(gp*green)+(bp*blue)), 
+								                    klamp255((greenpix*green)/(rp*red)+(bp*blue)), 
+											        klamp255((bluepix*blue)/(rp*red)+(gp*green)));
+						}
+					}
+					break;
+				case 32:
+					for(y = 0; y < h; y++) 
+					{
+						c1_32 = scanline32(f1, y);
+						c2_32 = scanline32(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r32(c1_32[x]) + 1;
+							gp = g32(c1_32[x]) + 1;
+							bp = b32(c1_32[x]) + 1;
+							redpix = r32(c2_32[x]) + 1;
+							greenpix = g32(c2_32[x]) + 1;
+							bluepix = b32(c2_32[x]) + 1;
+							c2_32[x] = rgbtocolor32(klamp255((redpix*red)/(gp*green)+(bp*blue)), 
+								                    klamp255((greenpix*green)/(rp*red)+(bp*blue)), 
+											        klamp255((bluepix*blue)/(rp*red)+(gp*green)));
+						}
+					}
+					break;
+			}
+			break;
+		case 7:
+			switch (f1.pixelformat)
+			{
+				case 16:
+					for(y = 0; y < h; y++) 
+					{
+						c1_16 = scanline16(f1, y);
+						c2_16 = scanline16(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r16(c1_16[x]) + red;
+							gp = g16(c1_16[x]) + green;
+							bp = b16(c1_16[x]) + blue;
+							redpix = r16(c2_16[x]) + red;
+							greenpix = g16(c2_16[x]) + green;
+							bluepix = b16(c2_16[x]) + blue;
+							c2_16[x] = rgbtocolor16(klamp255((redpix/(gp+bp))*red), 
+								                    klamp255((greenpix/(rp+bp))*green), 
+											        klamp255((bluepix/(rp+gp))*blue));
+						}
+					}
+					break;
+				case 24:
+					for(y = 0; y < h; y++) 
+					{
+						c1_24 = scanline24(f1, y);
+						c2_24 = scanline24(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r24(c1_24[x]) + red;
+							gp = g24(c1_24[x]) + green;
+							bp = b24(c1_24[x]) + blue;
+							redpix = r24(c2_24[x]) + red;
+							greenpix = g24(c2_24[x]) + green;
+							bluepix = b24(c2_24[x]) + blue;
+							c2_24[x] = rgbtocolor24(klamp255((redpix/(gp+bp))*red), 
+								                    klamp255((greenpix/(rp+bp))*green), 
+											        klamp255((bluepix/(rp+gp))*blue));
+						}
+					}
+					break;
+				case 32:
+					for(y = 0; y < h; y++) 
+					{
+						c1_32 = scanline32(f1, y);
+						c2_32 = scanline32(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r32(c1_32[x]) + red;
+							gp = g32(c1_32[x]) + green;
+							bp = b32(c1_32[x]) + blue;
+							redpix = r32(c2_32[x]) + red;
+							greenpix = g32(c2_32[x]) + green;
+							bluepix = b32(c2_32[x]) + blue;
+							c2_32[x] = rgbtocolor32(klamp255((redpix/(gp+bp))*red), 
+								                    klamp255((greenpix/(rp+bp))*green), 
+											        klamp255((bluepix/(rp+gp))*blue));
+						}
+					}
+					break;
+			}
+			break;
+		case 8:
+			switch (f1.pixelformat)
+			{
+				case 16:
+					for(y = 0; y < h; y++) 
+					{
+						c1_16 = scanline16(f1, y);
+						c2_16 = scanline16(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r16(c1_16[x]) + red;
+							gp = g16(c1_16[x]) + green;
+							bp = b16(c1_16[x]) + blue;
+							redpix = r16(c2_16[x]) + red;
+							greenpix = g16(c2_16[x]) + green;
+							bluepix = b16(c2_16[x]) + blue;
+							c2_16[x] = rgbtocolor16(klamp255(((gp+bp)/redpix)*red), 
+								                    klamp255(((rp+bp)/greenpix)*green), 
+											        klamp255(((rp+gp)/bluepix)*blue));
+						}
+					}
+					break;
+				case 24:
+					for(y = 0; y < h; y++) 
+					{
+						c1_24 = scanline24(f1, y);
+						c2_24 = scanline24(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r24(c1_24[x]) + red;
+							gp = g24(c1_24[x]) + green;
+							bp = b24(c1_24[x]) + blue;
+							redpix = r24(c2_24[x]) + red;
+							greenpix = g24(c2_24[x]) + green;
+							bluepix = b24(c2_24[x]) + blue;
+							c2_24[x] = rgbtocolor24(klamp255(((gp+bp)/redpix)*red), 
+								                    klamp255(((rp+bp)/greenpix)*green), 
+											        klamp255(((rp+gp)/bluepix)*blue));
+						}
+					}
+					break;
+				case 32:
+					for(y = 0; y < h; y++) 
+					{
+						c1_32 = scanline32(f1, y);
+						c2_32 = scanline32(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r32(c1_32[x]) + red;
+							gp = g32(c1_32[x]) + green;
+							bp = b32(c1_32[x]) + blue;
+							redpix = r32(c2_32[x]) + red;
+							greenpix = g32(c2_32[x]) + green;
+							bluepix = b32(c2_32[x]) + blue;
+							c2_32[x] = rgbtocolor32(klamp255(((gp+bp)/redpix)*red), 
+								                    klamp255(((rp+bp)/greenpix)*green), 
+											        klamp255(((rp+gp)/bluepix)*blue));
+						}
+					}
+					break;
+			}
+			break;
+		case 9:
+			switch (f1.pixelformat)
+			{
+				case 16:
+					for(y = 0; y < h; y++) 
+					{
+						c1_16 = scanline16(f1, y);
+						c2_16 = scanline16(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r16(c1_16[x]) + 1;
+							gp = g16(c1_16[x]) + 1;
+							bp = b16(c1_16[x]) + 1;
+							redpix = r16(c2_16[x]) + 1;
+							greenpix = g16(c2_16[x]) + 1;
+							bluepix = b16(c2_16[x]) + 1;
+							c2_16[x] = rgbtocolor16(klamp255((red + ((gp+bp)/redpix))*red), 
+								                    klamp255((green + ((rp+bp)/greenpix))*green), 
+											        klamp255((blue + ((rp+gp)/bluepix))*blue));
+						}
+					}
+					break;
+				case 24:
+					for(y = 0; y < h; y++) 
+					{
+						c1_24 = scanline24(f1, y);
+						c2_24 = scanline24(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r24(c1_24[x]) + 1;
+							gp = g24(c1_24[x]) + 1;
+							bp = b24(c1_24[x]) + 1;
+							redpix = r24(c2_24[x]) + 1;
+							greenpix = g24(c2_24[x]) + 1;
+							bluepix = b24(c2_24[x]) + 1;
+							c2_24[x] = rgbtocolor24(klamp255((red + ((gp+bp)/redpix))*red), 
+								                    klamp255((green + ((rp+bp)/greenpix))*green), 
+											        klamp255((blue + ((rp+gp)/bluepix))*blue));
+						}
+					}
+					break;
+				case 32:
+					for(y = 0; y < h; y++) 
+					{
+						c1_32 = scanline32(f1, y);
+						c2_32 = scanline32(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r32(c1_32[x]) + 1;
+							gp = g32(c1_32[x]) + 1;
+							bp = b32(c1_32[x]) + 1;
+							redpix = r32(c2_32[x]) + 1;
+							greenpix = g32(c2_32[x]) + 1;
+							bluepix = b32(c2_32[x]) + 1;
+							c2_32[x] = rgbtocolor32(klamp255((red + ((gp+bp)/redpix))*red), 
+								                    klamp255((green + ((rp+bp)/greenpix))*green), 
+											        klamp255((blue + ((rp+gp)/bluepix))*blue));
+						}
+					}
+					break;
+			}
+			break;
+		case 10:
+			switch (f1.pixelformat)
+			{
+				case 16:
+					for(y = 0; y < h; y++) 
+					{
+						c1_16 = scanline16(f1, y);
+						c2_16 = scanline16(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r16(c1_16[x]) + 1;
+							gp = g16(c1_16[x]) + 1;
+							bp = b16(c1_16[x]) + 1;
+							redpix = r16(c2_16[x]) + 1;
+							greenpix = g16(c2_16[x]) + 1;
+							bluepix = b16(c2_16[x]) + 1;
+							c2_16[x] = rgbtocolor16(klamp255((((gp*green)+(bp*blue))/redpix)*red), 
+								                    klamp255((((rp*red)+(bp*blue))/greenpix)*green), 
+											        klamp255((((rp*red)+(gp*green))/bluepix)*blue));
+						}
+					}
+					break;
+				case 24:
+					for(y = 0; y < h; y++) 
+					{
+						c1_24 = scanline24(f1, y);
+						c2_24 = scanline24(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r24(c1_24[x]) + 1;
+							gp = g24(c1_24[x]) + 1;
+							bp = b24(c1_24[x]) + 1;
+							redpix = r24(c2_24[x]) + 1;
+							greenpix = g24(c2_24[x]) + 1;
+							bluepix = b24(c2_24[x]) + 1;
+							c2_24[x] = rgbtocolor24(klamp255((((gp*green)+(bp*blue))/redpix)*red), 
+								                    klamp255((((rp*red)+(bp*blue))/greenpix)*green), 
+											        klamp255((((rp*red)+(gp*green))/bluepix)*blue));
+						}
+					}
+					break;
+				case 32:
+					for(y = 0; y < h; y++) 
+					{
+						c1_32 = scanline32(f1, y);
+						c2_32 = scanline32(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r32(c1_32[x]) + 1;
+							gp = g32(c1_32[x]) + 1;
+							bp = b32(c1_32[x]) + 1;
+							redpix = r32(c2_32[x]) + 1;
+							greenpix = g32(c2_32[x]) + 1;
+							bluepix = b32(c2_32[x]) + 1;
+							c2_32[x] = rgbtocolor32(klamp255((((gp*green)+(bp*blue))/redpix)*red), 
+								                    klamp255((((rp*red)+(bp*blue))/greenpix)*green), 
+											        klamp255((((rp*red)+(gp*green))/bluepix)*blue));
+						}
+					}
+					break;
+			}
+			break;
+		case 11:
+			switch (f1.pixelformat)
+			{
+				case 16:
+					for(y = 0; y < h; y++) 
+					{
+						c1_16 = scanline16(f1, y);
+						c2_16 = scanline16(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r16(c1_16[x]) + 1;
+							gp = g16(c1_16[x]) + 1;
+							bp = b16(c1_16[x]) + 1;
+							redpix = r16(c2_16[x]) + 1;
+							greenpix = g16(c2_16[x]) + 1;
+							bluepix = b16(c2_16[x]) + 1;
+							c2_16[x] = rgbtocolor16(klamp255((((gp*green)+(bp*blue))/(redpix*red))*red), 
+								                    klamp255((((rp*red)+(bp*blue))/(greenpix*green))*green), 
+											        klamp255((((rp*red)+(gp*green))/(bluepix*blue))*blue));
+						}
+					}
+					break;
+				case 24:
+					for(y = 0; y < h; y++) 
+					{
+						c1_24 = scanline24(f1, y);
+						c2_24 = scanline24(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r24(c1_24[x]) + 1;
+							gp = g24(c1_24[x]) + 1;
+							bp = b24(c1_24[x]) + 1;
+							redpix = r24(c2_24[x]) + 1;
+							greenpix = g24(c2_24[x]) + 1;
+							bluepix = b24(c2_24[x]) + 1;
+							c2_24[x] = rgbtocolor24(klamp255((((gp*green)+(bp*blue))/(redpix*red))*red), 
+								                    klamp255((((rp*red)+(bp*blue))/(greenpix*green))*green), 
+											        klamp255((((rp*red)+(gp*green))/(bluepix*blue))*blue));
+						}
+					}
+					break;
+				case 32:
+					for(y = 0; y < h; y++) 
+					{
+						c1_32 = scanline32(f1, y);
+						c2_32 = scanline32(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r32(c1_32[x]) + 1;
+							gp = g32(c1_32[x]) + 1;
+							bp = b32(c1_32[x]) + 1;
+							redpix = r32(c2_32[x]) + 1;
+							greenpix = g32(c2_32[x]) + 1;
+							bluepix = b32(c2_32[x]) + 1;
+							c2_32[x] = rgbtocolor32(klamp255((((gp*green)+(bp*blue))/(redpix*red))*red), 
+								                    klamp255((((rp*red)+(bp*blue))/(greenpix*green))*green), 
+											        klamp255((((rp*red)+(gp*green))/(bluepix*blue))*blue));
+						}
+					}
+					break;
+			}
+			break;
+		case 12:
+			switch (f1.pixelformat)
+			{
+				case 16:
+					for(y = 0; y < h; y++) 
+					{
+						c1_16 = scanline16(f1, y);
+						c2_16 = scanline16(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r16(c1_16[x]) + 1;
+							gp = g16(c1_16[x]) + 1;
+							bp = b16(c1_16[x]) + 1;
+							redpix = r16(c2_16[x]) + 1;
+							greenpix = g16(c2_16[x]) + 1;
+							bluepix = b16(c2_16[x]) + 1;
+							c2_16[x] = rgbtocolor16(klamp255(((redpix*red)/(gp*green)+(bp*blue))*red), 
+								                    klamp255(((greenpix*green)/(rp*red)+(bp*blue))*green), 
+											        klamp255(((bluepix*blue)/(rp*red)+(gp*green))*blue));
+						}
+					}
+					break;
+				case 24:
+					for(y = 0; y < h; y++) 
+					{
+						c1_24 = scanline24(f1, y);
+						c2_24 = scanline24(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r24(c1_24[x]) + 1;
+							gp = g24(c1_24[x]) + 1;
+							bp = b24(c1_24[x]) + 1;
+							redpix = r24(c2_24[x]) + 1;
+							greenpix = g24(c2_24[x]) + 1;
+							bluepix = b24(c2_24[x]) + 1;
+							c2_24[x] = rgbtocolor24(klamp255(((redpix*red)/(gp*green)+(bp*blue))*red), 
+								                    klamp255(((greenpix*green)/(rp*red)+(bp*blue))*green), 
+											        klamp255(((bluepix*blue)/(rp*red)+(gp*green))*blue));
+						}
+					}
+					break;
+				case 32:
+					for(y = 0; y < h; y++) 
+					{
+						c1_32 = scanline32(f1, y);
+						c2_32 = scanline32(f2, y);
+						for(x = 0; x < w; x++)
+						{
+							rp = r32(c1_32[x]) + 1;
+							gp = g32(c1_32[x]) + 1;
+							bp = b32(c1_32[x]) + 1;
+							redpix = r32(c2_32[x]) + 1;
+							greenpix = g32(c2_32[x]) + 1;
+							bluepix = b32(c2_32[x]) + 1;
+							c2_32[x] = rgbtocolor32(klamp255(((redpix*red)/(gp*green)+(bp*blue))*red), 
+								                    klamp255(((greenpix*green)/(rp*red)+(bp*blue))*green), 
+											        klamp255(((bluepix*blue)/(rp*red)+(gp*green))*blue));
+						}
+					}
+					break;
+			}
+			break;
+	}
+
+}
\ No newline at end of file
-- 
cgit v1.2.1