+/*
+ * Copyright (C) 2000 Red Hat, Inc
+ * mediaLib integration Copyright (c) 2001-2007 Sun Microsystems, Inc.
+ * All rights reserved. (Brian Cameron, Dmitriy Demin, James Cheng,
+ * Padraig O'Briain)
+ *
+ * This library is free software; you can redistribute it and/or
+ * modify it under the terms of the GNU Lesser General Public
+ * License as published by the Free Software Foundation; either
+ * version 2 of the License, or (at your option) any later version.
+ *
+ * This library 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
+ * Lesser General Public License for more details.
+ *
+ * You should have received a copy of the GNU Lesser General Public
+ * License along with this library; if not, write to the
+ * Free Software Foundation, Inc., 59 Temple Place - Suite 330,
+ * Boston, MA 02111-1307, USA.
+ */
+#include "config.h"
#include <math.h>
#include <glib.h>
-#include "config.h"
#include "pixops.h"
#include "pixops-internal.h"
#define SUBSAMPLE_MASK ((1 << SUBSAMPLE_BITS)-1)
#define SCALE_SHIFT 16
+static void
+_pixops_scale_real (guchar *dest_buf,
+ int render_x0,
+ int render_y0,
+ int render_x1,
+ int render_y1,
+ int dest_rowstride,
+ int dest_channels,
+ gboolean dest_has_alpha,
+ const guchar *src_buf,
+ int src_width,
+ int src_height,
+ int src_rowstride,
+ int src_channels,
+ gboolean src_has_alpha,
+ double scale_x,
+ double scale_y,
+ PixopsInterpType interp_type);
+
typedef struct _PixopsFilter PixopsFilter;
+typedef struct _PixopsFilterDimension PixopsFilterDimension;
+
+struct _PixopsFilterDimension
+{
+ int n;
+ double offset;
+ double *weights;
+};
struct _PixopsFilter
{
- int *weights;
- int n_x;
- int n_y;
- double x_offset;
- double y_offset;
+ PixopsFilterDimension x;
+ PixopsFilterDimension y;
+ double overall_alpha;
};
typedef guchar *(*PixopsLineFunc) (int *weights, int n_x, int n_y,
- guchar *dest, int dest_x, guchar *dest_end, int dest_channels, int dest_has_alpha,
- guchar **src, int src_channels, gboolean src_has_alpha,
- int x_init, int x_step, int src_width,
- int check_size, guint32 color1, guint32 color2);
+ guchar *dest, int dest_x, guchar *dest_end,
+ int dest_channels, int dest_has_alpha,
+ guchar **src, int src_channels,
+ gboolean src_has_alpha, int x_init,
+ int x_step, int src_width, int check_size,
+ guint32 color1, guint32 color2);
+typedef void (*PixopsPixelFunc) (guchar *dest, int dest_x, int dest_channels,
+ int dest_has_alpha, int src_has_alpha,
+ int check_size, guint32 color1,
+ guint32 color2,
+ guint r, guint g, guint b, guint a);
+
+#ifdef USE_MEDIALIB
+#include <stdlib.h>
+#include <dlfcn.h>
+#include <mlib_image.h>
+
+#ifdef HAVE_STRINGS_H
+#include <strings.h>
+#endif
+
+#ifdef HAVE_STRING_H
+#include <string.h>
+#endif
+
+#if defined(HAVE_SYS_SYSTEMINFO_H)
+#include <sys/systeminfo.h>
+#elif defined(HAVE_SYS_SYSINFO_H)
+#include <sys/sysinfo.h>
+#endif
+
+static void pixops_medialib_composite (guchar *dest_buf,
+ int dest_width,
+ int dest_height,
+ int dest_rowstride,
+ int dest_channels,
+ int dest_has_alpha,
+ const guchar *src_buf,
+ int src_width,
+ int src_height,
+ int src_rowstride,
+ int src_channels,
+ int src_has_alpha,
+ int dest_x,
+ int dest_y,
+ int dest_region_width,
+ int dest_region_height,
+ double offset_x,
+ double offset_y,
+ double scale_x,
+ double scale_y,
+ PixopsInterpType interp_type,
+ int overall_alpha);
+
+static void pixops_medialib_scale (guchar *dest_buf,
+ int dest_width,
+ int dest_height,
+ int dest_rowstride,
+ int dest_channels,
+ int dest_has_alpha,
+ const guchar *src_buf,
+ int src_width,
+ int src_height,
+ int src_rowstride,
+ int src_channels,
+ int src_has_alpha,
+ int dest_x,
+ int dest_y,
+ int dest_region_width,
+ int dest_region_height,
+ double offset_x,
+ double offset_y,
+ double scale_x,
+ double scale_y,
+ PixopsInterpType interp_type);
+
+typedef struct _mlInterp mlInterp;
+
+struct _mlInterp
+{
+ double tx;
+ double ty;
+ PixopsFilter po_filter;
+ void *interp_table;
+};
+
+static gboolean medialib_initialized = FALSE;
+static gboolean use_medialib = TRUE;
+
+/*
+ * Sun mediaLib(tm) support.
+ *
+ * http://www.sun.com/processors/vis/mlib.html
+ *
+ */
+static void
+_pixops_use_medialib ()
+{
+ char *mlib_version_string;
+ char sys_info[257];
+ long count;
+
+ medialib_initialized = TRUE;
+
+ if (getenv ("GDK_DISABLE_MEDIALIB"))
+ {
+ use_medialib = FALSE;
+ return;
+ }
-typedef void (*PixopsPixelFunc) (guchar *dest, int dest_x, int dest_channels, int dest_has_alpha,
- int src_has_alpha, int check_size, guint32 color1,
- guint32 color2,
- int r, int g, int b, int a);
+ /*
+ * The imaging functions we want to use were added in mediaLib version 2.
+ * So turn off mediaLib support if the user has an older version.
+ * mlib_version returns a string in this format:
+ *
+ * mediaLib:0210:20011101:v8plusa
+ * ^^^^^^^^ ^^^^ ^^^^^^^^ ^^^^^^^
+ * libname vers build ISALIST identifier
+ * date (in this case sparcv8plus+vis)
+ *
+ * The first 2 digits of the version are the major version. The 3rd digit
+ * is the minor version, and the 4th digit is the micro version. So the
+ * above string corresponds to version 2.1.0. In the following test we only
+ * care about the major version.
+ */
+ mlib_version_string = mlib_version ();
+
+ count = sysinfo (SI_ARCHITECTURE, &sys_info[0], 257);
+
+ if (count != -1)
+ {
+ if (strcmp (sys_info, "i386") == 0)
+ {
+ char *mlib_target_isa = &mlib_version_string[23];
+
+ /*
+ * For x86 processors mediaLib generic C implementation
+ * does not give any performance advantage so disable it
+ */
+ if (strncmp (mlib_target_isa, "sse", 3) != 0)
+ {
+ use_medialib = FALSE;
+ return;
+ }
+
+ /*
+ * For x86 processors use of libumem conflicts with
+ * mediaLib, so avoid using it.
+ */
+ if ((dlsym (RTLD_DEFAULT, "umem_alloc") != NULL) ||
+ (dlsym (RTLD_PROBE, "umem_alloc") != NULL) ||
+ (dlsym (RTLD_NEXT, "umem_alloc") != NULL) ||
+ (dlsym (RTLD_SELF, "umem_alloc") != NULL))
+ {
+ use_medialib = FALSE;
+ return;
+ }
+ }
+ }
+ else
+ {
+ /* Failed to get system architecture, disable mediaLib anyway */
+ use_medialib = FALSE;
+ return;
+ }
+}
+#endif
static int
get_check_shift (int check_size)
double scale_x,
double scale_y)
{
- int i, j;
+ int i;
int x;
int x_step = (1 << SCALE_SHIFT) / scale_x;
int y_step = (1 << SCALE_SHIFT) / scale_y;
-
-#define INNER_LOOP(SRC_CHANNELS,DEST_CHANNELS) \
- for (j=0; j < (render_x1 - render_x0); j++) \
- { \
- const guchar *p = src + (x >> SCALE_SHIFT) * SRC_CHANNELS; \
- \
- dest[0] = p[0]; \
- dest[1] = p[1]; \
- dest[2] = p[2]; \
- \
- if (DEST_CHANNELS == 4) \
- { \
- if (SRC_CHANNELS == 4) \
- dest[3] = p[3]; \
- else \
- dest[3] = 0xff; \
- } \
- \
- dest += DEST_CHANNELS; \
- x += x_step; \
- }
+ int xmax, xstart, xstop, x_pos, y_pos;
+ const guchar *p;
+
+#define INNER_LOOP(SRC_CHANNELS,DEST_CHANNELS,ASSIGN_PIXEL) \
+ xmax = x + (render_x1 - render_x0) * x_step; \
+ xstart = MIN (0, xmax); \
+ xstop = MIN (src_width << SCALE_SHIFT, xmax); \
+ p = src + (CLAMP (x, xstart, xstop) >> SCALE_SHIFT) * SRC_CHANNELS; \
+ while (x < xstart) \
+ { \
+ ASSIGN_PIXEL; \
+ dest += DEST_CHANNELS; \
+ x += x_step; \
+ } \
+ while (x < xstop) \
+ { \
+ p = src + (x >> SCALE_SHIFT) * SRC_CHANNELS; \
+ ASSIGN_PIXEL; \
+ dest += DEST_CHANNELS; \
+ x += x_step; \
+ } \
+ x_pos = x >> SCALE_SHIFT; \
+ p = src + CLAMP (x_pos, 0, src_width - 1) * SRC_CHANNELS; \
+ while (x < xmax) \
+ { \
+ ASSIGN_PIXEL; \
+ dest += DEST_CHANNELS; \
+ x += x_step; \
+ }
for (i = 0; i < (render_y1 - render_y0); i++)
{
- const guchar *src = src_buf + (((i + render_y0) * y_step + y_step / 2) >> SCALE_SHIFT) * src_rowstride;
- guchar *dest = dest_buf + i * dest_rowstride;
+ const guchar *src;
+ guchar *dest;
+ y_pos = ((i + render_y0) * y_step + y_step / 2) >> SCALE_SHIFT;
+ y_pos = CLAMP (y_pos, 0, src_height - 1);
+ src = src_buf + y_pos * src_rowstride;
+ dest = dest_buf + i * dest_rowstride;
x = render_x0 * x_step + x_step / 2;
{
if (dest_channels == 3)
{
- INNER_LOOP (3, 3);
+ INNER_LOOP (3, 3, dest[0]=p[0];dest[1]=p[1];dest[2]=p[2]);
}
else
{
- INNER_LOOP (3, 4);
+ INNER_LOOP (3, 4, dest[0]=p[0];dest[1]=p[1];dest[2]=p[2];dest[3]=0xff);
}
}
else if (src_channels == 4)
{
if (dest_channels == 3)
{
- INNER_LOOP (4, 3);
+ INNER_LOOP (4, 3, dest[0]=p[0];dest[1]=p[1];dest[2]=p[2]);
}
else
{
- for (j=0; j < (render_x1 - render_x0); j++)
- {
- const guchar *p = src + (x >> SCALE_SHIFT) * 4;
- guint32 *p32;
-
- p32 = (guint32 *) dest;
- *p32 = *((guint32 *) p);
-
- dest += 4;
- x += x_step;
- }
+ guint32 *p32;
+ INNER_LOOP(4, 4, p32=(guint32*)dest;*p32=*((guint32*)p));
}
}
}
-#undef INNER_LOOP
}
static void
double scale_y,
int overall_alpha)
{
- int i, j;
+ int i;
int x;
int x_step = (1 << SCALE_SHIFT) / scale_x;
int y_step = (1 << SCALE_SHIFT) / scale_y;
+ int xmax, xstart, xstop, x_pos, y_pos;
+ const guchar *p;
+ unsigned int a0;
for (i = 0; i < (render_y1 - render_y0); i++)
{
- const guchar *src = src_buf + (((i + render_y0) * y_step + y_step / 2) >> SCALE_SHIFT) * src_rowstride;
- guchar *dest = dest_buf + i * dest_rowstride + render_x0 * dest_channels;
+ const guchar *src;
+ guchar *dest;
+ y_pos = ((i + render_y0) * y_step + y_step / 2) >> SCALE_SHIFT;
+ y_pos = CLAMP (y_pos, 0, src_height - 1);
+ src = src_buf + y_pos * src_rowstride;
+ dest = dest_buf + i * dest_rowstride;
x = render_x0 * x_step + x_step / 2;
- for (j=0; j < (render_x1 - render_x0); j++)
- {
- const guchar *p = src + (x >> SCALE_SHIFT) * src_channels;
- unsigned int a0;
-
+ INNER_LOOP(src_channels, dest_channels,
if (src_has_alpha)
- a0 = (p[3] * overall_alpha + 0xff) >> 8;
+ a0 = (p[3] * overall_alpha) / 0xff;
else
a0 = overall_alpha;
- if (dest_has_alpha)
- {
- unsigned int a1 = dest[3];
- unsigned int total = a0 + a1;
-
- if (total)
- {
- dest[0] = (a0 * src[0] + a1 * dest[0]) / (total);
- dest[1] = (a0 * src[1] + a1 * dest[1]) / (total);
- dest[2] = (a0 * src[2] + a1 * dest[2]) / (total);
- dest[3] = total - ((a0 * a1 + 0xff) >> 8);
- }
- else
- {
- dest[0] = 0;
- dest[1] = 0;
- dest[2] = 0;
- dest[3] = 0;
- }
- }
- else
- {
- dest[0] = dest[0] + ((a0 * (p[0] - dest[0]) + 0xff) >> 8);
- dest[1] = dest[1] + ((a0 * (p[1] - dest[1]) + 0xff) >> 8);
- dest[2] = dest[2] + ((a0 * (p[2] - dest[2]) + 0xff) >> 8);
-
- if (dest_channels == 4)
- dest[3] = 0xff;
- }
-
- dest += dest_channels;
- x += x_step;
- }
+ switch (a0)
+ {
+ case 0:
+ break;
+ case 255:
+ dest[0] = p[0];
+ dest[1] = p[1];
+ dest[2] = p[2];
+ if (dest_has_alpha)
+ dest[3] = 0xff;
+ break;
+ default:
+ if (dest_has_alpha)
+ {
+ unsigned int w0 = 0xff * a0;
+ unsigned int w1 = (0xff - a0) * dest[3];
+ unsigned int w = w0 + w1;
+
+ dest[0] = (w0 * p[0] + w1 * dest[0]) / w;
+ dest[1] = (w0 * p[1] + w1 * dest[1]) / w;
+ dest[2] = (w0 * p[2] + w1 * dest[2]) / w;
+ dest[3] = w / 0xff;
+ }
+ else
+ {
+ unsigned int a1 = 0xff - a0;
+ unsigned int tmp;
+
+ tmp = a0 * p[0] + a1 * dest[0] + 0x80;
+ dest[0] = (tmp + (tmp >> 8)) >> 8;
+ tmp = a0 * p[1] + a1 * dest[1] + 0x80;
+ dest[1] = (tmp + (tmp >> 8)) >> 8;
+ tmp = a0 * p[2] + a1 * dest[2] + 0x80;
+ dest[2] = (tmp + (tmp >> 8)) >> 8;
+ }
+ break;
+ }
+ );
}
}
int y_step = (1 << SCALE_SHIFT) / scale_y;
int r1, g1, b1, r2, g2, b2;
int check_shift = get_check_shift (check_size);
+ int xmax, xstart, xstop, x_pos, y_pos;
+ const guchar *p;
+ unsigned int a0;
for (i = 0; i < (render_y1 - render_y0); i++)
{
- const guchar *src = src_buf + (((i + render_y0) * y_step + y_step/2) >> SCALE_SHIFT) * src_rowstride;
- guchar *dest = dest_buf + i * dest_rowstride;
+ const guchar *src;
+ guchar *dest;
+ y_pos = ((i + render_y0) * y_step + y_step / 2) >> SCALE_SHIFT;
+ y_pos = CLAMP (y_pos, 0, src_height - 1);
+ src = src_buf + y_pos * src_rowstride;
+ dest = dest_buf + i * dest_rowstride;
x = render_x0 * x_step + x_step / 2;
+
if (((i + check_y) >> check_shift) & 1)
{
r1 = (color2 & 0xff0000) >> 16;
b2 = color2 & 0xff;
}
- for (j=0 ; j < (render_x1 - render_x0); j++)
- {
- const guchar *p = src + (x >> SCALE_SHIFT) * src_channels;
- unsigned int a0;
-
+ j = 0;
+ INNER_LOOP(src_channels, dest_channels,
if (src_has_alpha)
a0 = (p[3] * overall_alpha + 0xff) >> 8;
else
a0 = overall_alpha;
- if (((j + check_x) >> check_shift) & 1)
- {
- dest[0] = r2 + ((a0 * ((int)p[0] - r2) + 0xff) >> 8);
- dest[1] = g2 + ((a0 * ((int)p[1] - g2) + 0xff) >> 8);
- dest[2] = b2 + ((a0 * ((int)p[2] - b2) + 0xff) >> 8);
- }
- else
- {
- dest[0] = r1 + ((a0 * ((int)p[0] - r1) + 0xff) >> 8);
- dest[1] = g1 + ((a0 * ((int)p[1] - g1) + 0xff) >> 8);
- dest[2] = b1 + ((a0 * ((int)p[2] - b1) + 0xff) >> 8);
- }
+ switch (a0)
+ {
+ case 0:
+ if (((j + check_x) >> check_shift) & 1)
+ {
+ dest[0] = r2;
+ dest[1] = g2;
+ dest[2] = b2;
+ }
+ else
+ {
+ dest[0] = r1;
+ dest[1] = g1;
+ dest[2] = b1;
+ }
+ break;
+ case 255:
+ dest[0] = p[0];
+ dest[1] = p[1];
+ dest[2] = p[2];
+ break;
+ default:
+ {
+ unsigned int tmp;
+ if (((j + check_x) >> check_shift) & 1)
+ {
+ tmp = ((int) p[0] - r2) * a0;
+ dest[0] = r2 + ((tmp + (tmp >> 8) + 0x80) >> 8);
+ tmp = ((int) p[1] - g2) * a0;
+ dest[1] = g2 + ((tmp + (tmp >> 8) + 0x80) >> 8);
+ tmp = ((int) p[2] - b2) * a0;
+ dest[2] = b2 + ((tmp + (tmp >> 8) + 0x80) >> 8);
+ }
+ else
+ {
+ tmp = ((int) p[0] - r1) * a0;
+ dest[0] = r1 + ((tmp + (tmp >> 8) + 0x80) >> 8);
+ tmp = ((int) p[1] - g1) * a0;
+ dest[1] = g1 + ((tmp + (tmp >> 8) + 0x80) >> 8);
+ tmp = ((int) p[2] - b1) * a0;
+ dest[2] = b1 + ((tmp + (tmp >> 8) + 0x80) >> 8);
+ }
+ }
+ break;
+ }
if (dest_channels == 4)
dest[3] = 0xff;
- dest += dest_channels;
- x += x_step;
- }
+ j++;
+ );
}
}
+#undef INNER_LOOP
static void
composite_pixel (guchar *dest, int dest_x, int dest_channels, int dest_has_alpha,
int src_has_alpha, int check_size, guint32 color1, guint32 color2,
- int r, int g, int b, int a)
+ guint r, guint g, guint b, guint a)
{
if (dest_has_alpha)
{
- unsigned int w = (((1 << 16) - a) * dest[3]) >> 8;
- unsigned int total = a + w;
-
- if (total)
+ unsigned int w0 = a - (a >> 8);
+ unsigned int w1 = ((0xff0000 - a) >> 8) * dest[3];
+ unsigned int w = w0 + w1;
+
+ if (w != 0)
{
- dest[0] = (r + w * dest[0]) / total;
- dest[1] = (g + w * dest[1]) / total;
- dest[2] = (b + w * dest[2]) / total;
- dest[3] = (r * w) >> 16;
+ dest[0] = (r - (r >> 8) + w1 * dest[0]) / w;
+ dest[1] = (g - (g >> 8) + w1 * dest[1]) / w;
+ dest[2] = (b - (b >> 8) + w1 * dest[2]) / w;
+ dest[3] = w / 0xff00;
}
else
{
}
else
{
- dest[0] = ((0xff0000 - a) * dest[0] + r) >> 24;
- dest[1] = ((0xff0000 - a) * dest[1] + g) >> 24;
- dest[2] = ((0xff0000 - a) * dest[2] + b) >> 24;
+ dest[0] = (r + (0xff0000 - a) * dest[0]) / 0xff0000;
+ dest[1] = (g + (0xff0000 - a) * dest[1]) / 0xff0000;
+ dest[2] = (b + (0xff0000 - a) * dest[2]) / 0xff0000;
}
}
int *pixel_weights;
pixel_weights = weights + ((x >> (SCALE_SHIFT - SUBSAMPLE_BITS)) & SUBSAMPLE_MASK) * n_x * n_y;
-
+
for (i=0; i<n_y; i++)
{
guchar *q = src[i] + x_scaled * src_channels;
ta = q[3] * line_weights[j];
else
ta = 0xff * line_weights[j];
-
+
r += ta * q[0];
g += ta * q[1];
b += ta * q[2];
if (dest_has_alpha)
{
- unsigned int w = (((1 << 16) - a) * dest[3]) >> 8;
- unsigned int total = a + w;
+ unsigned int w0 = a - (a >> 8);
+ unsigned int w1 = ((0xff0000 - a) >> 8) * dest[3];
+ unsigned int w = w0 + w1;
- if (total)
+ if (w != 0)
{
- dest[0] = (r + w * dest[0]) / total;
- dest[1] = (r + w * dest[1]) / total;
- dest[2] = (r + w * dest[2]) / total;
- dest[3] = (r * w) >> 16;
+ dest[0] = (r - (r >> 8) + w1 * dest[0]) / w;
+ dest[1] = (g - (g >> 8) + w1 * dest[1]) / w;
+ dest[2] = (b - (b >> 8) + w1 * dest[2]) / w;
+ dest[3] = w / 0xff00;
}
else
{
}
else
{
- dest[0] = ((0xff0000 - a) * dest[0] + r) >> 24;
- dest[1] = ((0xff0000 - a) * dest[1] + g) >> 24;
- dest[2] = ((0xff0000 - a) * dest[2] + b) >> 24;
+ dest[0] = (r + (0xff0000 - a) * dest[0]) / 0xff0000;
+ dest[1] = (g + (0xff0000 - a) * dest[1]) / 0xff0000;
+ dest[2] = (b + (0xff0000 - a) * dest[2]) / 0xff0000;
}
dest += dest_channels;
q0 = src0 + x_scaled * 4;
q1 = src1 + x_scaled * 4;
- pixel_weights = (int *)((char *)weights + ((x >> (SCALE_SHIFT - SUBSAMPLE_BITS - 4)) & (SUBSAMPLE_MASK << 4)));
+ pixel_weights = (int *)((char *)weights +
+ ((x >> (SCALE_SHIFT - SUBSAMPLE_BITS - 4)) & (SUBSAMPLE_MASK << 4)));
w1 = pixel_weights[0];
w2 = pixel_weights[1];
b += ta * q0[6];
a += ta;
- ta = w3 * q0[3];
- r += ta * q0[0];
- g += ta * q0[1];
- b += ta * q0[2];
+ ta = w3 * q1[3];
+ r += ta * q1[0];
+ g += ta * q1[1];
+ b += ta * q1[2];
a += ta;
- ta += w4 * q1[7];
+ ta = w4 * q1[7];
r += ta * q1[4];
g += ta * q1[5];
b += ta * q1[6];
#ifdef USE_MMX
static guchar *
-composite_line_22_4a4_mmx_stub (int *weights, int n_x, int n_y,
- guchar *dest, int dest_x, guchar *dest_end, int dest_channels, int dest_has_alpha,
- guchar **src, int src_channels, gboolean src_has_alpha,
- int x_init, int x_step, int src_width,
- int check_size, guint32 color1, guint32 color2)
+composite_line_22_4a4_mmx_stub (int *weights, int n_x, int n_y, guchar *dest,
+ int dest_x, guchar *dest_end,
+ int dest_channels, int dest_has_alpha,
+ guchar **src, int src_channels,
+ gboolean src_has_alpha, int x_init,
+ int x_step, int src_width, int check_size,
+ guint32 color1, guint32 color2)
{
guint32 mmx_weights[16][8];
int j;
mmx_weights[j][7] = 0x00010001 * (weights[4*j + 3] >> 8);
}
- return pixops_composite_line_22_4a4_mmx (mmx_weights, dest, src[0], src[1], x_step, dest_end, x_init);
+ return _pixops_composite_line_22_4a4_mmx (mmx_weights, dest, src[0], src[1],
+ x_step, dest_end, x_init);
}
#endif /* USE_MMX */
static void
-composite_pixel_color (guchar *dest, int dest_x, int dest_channels, int dest_has_alpha,
- int src_has_alpha, int check_size, guint32 color1, guint32 color2,
- int r, int g, int b, int a)
+composite_pixel_color (guchar *dest, int dest_x, int dest_channels,
+ int dest_has_alpha, int src_has_alpha, int check_size,
+ guint32 color1, guint32 color2, guint r, guint g,
+ guint b, guint a)
{
int dest_r, dest_g, dest_b;
int check_shift = get_check_shift (check_size);
}
static guchar *
-composite_line_color (int *weights, int n_x, int n_y,
- guchar *dest, int dest_x, guchar *dest_end, int dest_channels, int dest_has_alpha,
- guchar **src, int src_channels, gboolean src_has_alpha,
- int x_init, int x_step, int src_width,
- int check_size, guint32 color1, guint32 color2)
+composite_line_color (int *weights, int n_x, int n_y, guchar *dest,
+ int dest_x, guchar *dest_end, int dest_channels,
+ int dest_has_alpha, guchar **src, int src_channels,
+ gboolean src_has_alpha, int x_init, int x_step,
+ int src_width, int check_size, guint32 color1,
+ guint32 color2)
{
int x = x_init;
int i, j;
#ifdef USE_MMX
static guchar *
composite_line_color_22_4a4_mmx_stub (int *weights, int n_x, int n_y,
- guchar *dest, int dest_x, guchar *dest_end, int dest_channels, int dest_has_alpha,
- guchar **src, int src_channels, gboolean src_has_alpha,
+ guchar *dest, int dest_x,
+ guchar *dest_end, int dest_channels,
+ int dest_has_alpha, guchar **src,
+ int src_channels, gboolean src_has_alpha,
int x_init, int x_step, int src_width,
- int check_size, guint32 color1, guint32 color2)
+ int check_size, guint32 color1,
+ guint32 color2)
{
guint32 mmx_weights[16][8];
int check_shift = get_check_shift (check_size);
colors[2] = (color2 & 0xff00) << 8 | (color2 & 0xff);
colors[3] = (color2 & 0xff0000) >> 16;
- return pixops_composite_line_color_22_4a4_mmx (mmx_weights, dest, src[0], src[1], x_step, dest_end, x_init,
- dest_x, check_shift, colors);
+ return _pixops_composite_line_color_22_4a4_mmx (mmx_weights, dest, src[0],
+ src[1], x_step, dest_end, x_init, dest_x, check_shift, colors);
}
#endif /* USE_MMX */
static void
scale_pixel (guchar *dest, int dest_x, int dest_channels, int dest_has_alpha,
int src_has_alpha, int check_size, guint32 color1, guint32 color2,
- int r, int g, int b, int a)
+ guint r, guint g, guint b, guint a)
{
if (src_has_alpha)
{
}
static guchar *
-scale_line (int *weights, int n_x, int n_y,
- guchar *dest, int dest_x, guchar *dest_end, int dest_channels, int dest_has_alpha,
- guchar **src, int src_channels, gboolean src_has_alpha,
- int x_init, int x_step, int src_width,
- int check_size, guint32 color1, guint32 color2)
+scale_line (int *weights, int n_x, int n_y, guchar *dest, int dest_x,
+ guchar *dest_end, int dest_channels, int dest_has_alpha,
+ guchar **src, int src_channels, gboolean src_has_alpha, int x_init,
+ int x_step, int src_width, int check_size, guint32 color1,
+ guint32 color2)
{
int x = x_init;
int i, j;
int x_scaled = x >> SCALE_SHIFT;
int *pixel_weights;
- pixel_weights = weights + ((x >> (SCALE_SHIFT - SUBSAMPLE_BITS)) & SUBSAMPLE_MASK) * n_x * n_y;
+ pixel_weights = weights +
+ ((x >> (SCALE_SHIFT - SUBSAMPLE_BITS)) & SUBSAMPLE_MASK) * n_x * n_y;
if (src_has_alpha)
{
#ifdef USE_MMX
static guchar *
-scale_line_22_33_mmx_stub (int *weights, int n_x, int n_y,
- guchar *dest, int dest_x, guchar *dest_end, int dest_channels, int dest_has_alpha,
- guchar **src, int src_channels, gboolean src_has_alpha,
- int x_init, int x_step, int src_width,
- int check_size, guint32 color1, guint32 color2)
+scale_line_22_33_mmx_stub (int *weights, int n_x, int n_y, guchar *dest,
+ int dest_x, guchar *dest_end, int dest_channels,
+ int dest_has_alpha, guchar **src, int src_channels,
+ gboolean src_has_alpha, int x_init, int x_step,
+ int src_width, int check_size, guint32 color1,
+ guint32 color2)
{
guint32 mmx_weights[16][8];
int j;
mmx_weights[j][7] = 0x00010001 * (weights[4*j + 3] >> 8);
}
- return pixops_scale_line_22_33_mmx (mmx_weights, dest, src[0], src[1], x_step, dest_end, x_init);
+ return _pixops_scale_line_22_33_mmx (mmx_weights, dest, src[0], src[1],
+ x_step, dest_end, x_init);
}
#endif /* USE_MMX */
static guchar *
-scale_line_22_33 (int *weights, int n_x, int n_y,
- guchar *dest, guchar *dest_end, int dest_channels, int dest_has_alpha,
+scale_line_22_33 (int *weights, int n_x, int n_y, guchar *dest, int dest_x,
+ guchar *dest_end, int dest_channels, int dest_has_alpha,
guchar **src, int src_channels, gboolean src_has_alpha,
int x_init, int x_step, int src_width,
int check_size, guint32 color1, guint32 color2)
q0 = src0 + x_scaled * 3;
q1 = src1 + x_scaled * 3;
- pixel_weights = (int *)((char *)weights + ((x >> (SCALE_SHIFT - SUBSAMPLE_BITS - 4)) & (SUBSAMPLE_MASK << 4)));
-
+ pixel_weights = weights +
+ ((x >> (SCALE_SHIFT - SUBSAMPLE_BITS)) & SUBSAMPLE_MASK) * 4;
+
w1 = pixel_weights[0];
w2 = pixel_weights[1];
w3 = pixel_weights[2];
g += w3 * q1[1];
b += w3 * q1[2];
- r += w4 * q1[4];
- g += w4 * q1[5];
- b += w4 * q1[6];
+ r += w4 * q1[3];
+ g += w4 * q1[4];
+ b += w4 * q1[5];
- dest[0] = r >> 16;
- dest[1] = g >> 16;
- dest[2] = b >> 16;
+ dest[0] = (r + 0x8000) >> 16;
+ dest[1] = (g + 0x8000) >> 16;
+ dest[2] = (b + 0x8000) >> 16;
dest += 3;
-
x += x_step;
}
}
static void
-process_pixel (int *weights, int n_x, int n_y,
- guchar *dest, int dest_x, int dest_channels, int dest_has_alpha,
- guchar **src, int src_channels, gboolean src_has_alpha,
- int x_start, int src_width,
- int check_size, guint32 color1, guint32 color2,
+process_pixel (int *weights, int n_x, int n_y, guchar *dest, int dest_x,
+ int dest_channels, int dest_has_alpha, guchar **src,
+ int src_channels, gboolean src_has_alpha, int x_start,
+ int src_width, int check_size, guint32 color1, guint32 color2,
PixopsPixelFunc pixel_func)
{
unsigned int r = 0, g = 0, b = 0, a = 0;
}
}
- (*pixel_func) (dest, dest_x, dest_channels, dest_has_alpha, src_has_alpha, check_size, color1, color2, r, g, b, a);
+ (*pixel_func) (dest, dest_x, dest_channels, dest_has_alpha, src_has_alpha,
+ check_size, color1, color2, r, g, b, a);
+}
+
+static void
+correct_total (int *weights,
+ int n_x,
+ int n_y,
+ int total,
+ double overall_alpha)
+{
+ int correction = (int)(0.5 + 65536 * overall_alpha) - total;
+ int remaining, c, d, i;
+
+ if (correction != 0)
+ {
+ remaining = correction;
+ for (d = 1, c = correction; c != 0 && remaining != 0; d++, c = correction / d)
+ for (i = n_x * n_y - 1; i >= 0 && c != 0 && remaining != 0; i--)
+ if (*(weights + i) + c >= 0)
+ {
+ *(weights + i) += c;
+ remaining -= c;
+ if ((0 < remaining && remaining < c) ||
+ (0 > remaining && remaining > c))
+ c = remaining;
+ }
+ }
+}
+
+static int *
+make_filter_table (PixopsFilter *filter)
+{
+ int i_offset, j_offset;
+ int n_x = filter->x.n;
+ int n_y = filter->y.n;
+ int *weights = g_new (int, SUBSAMPLE * SUBSAMPLE * n_x * n_y);
+
+ for (i_offset=0; i_offset < SUBSAMPLE; i_offset++)
+ for (j_offset=0; j_offset < SUBSAMPLE; j_offset++)
+ {
+ double weight;
+ int *pixel_weights = weights + ((i_offset*SUBSAMPLE) + j_offset) * n_x * n_y;
+ int total = 0;
+ int i, j;
+
+ for (i=0; i < n_y; i++)
+ for (j=0; j < n_x; j++)
+ {
+ weight = filter->x.weights[(j_offset * n_x) + j] *
+ filter->y.weights[(i_offset * n_y) + i] *
+ filter->overall_alpha * 65536 + 0.5;
+
+ total += (int)weight;
+
+ *(pixel_weights + n_x * i + j) = weight;
+ }
+
+ correct_total (pixel_weights, n_x, n_y, total, filter->overall_alpha);
+ }
+
+ return weights;
}
static void
PixopsPixelFunc pixel_func)
{
int i, j;
- int x, y;
- guchar **line_bufs = g_new (guchar *, filter->n_y);
-
- int x_step = (1 << SCALE_SHIFT) / scale_x;
- int y_step = (1 << SCALE_SHIFT) / scale_y;
+ int x, y; /* X and Y position in source (fixed_point) */
+
+ guchar **line_bufs = g_new (guchar *, filter->y.n);
+ int *filter_weights = make_filter_table (filter);
- int dest_x;
- int scaled_x_offset = floor (filter->x_offset * (1 << SCALE_SHIFT));
+ int x_step = (1 << SCALE_SHIFT) / scale_x; /* X step in source (fixed point) */
+ int y_step = (1 << SCALE_SHIFT) / scale_y; /* Y step in source (fixed point) */
- int run_end_index = (((src_width - filter->n_x + 1) << SCALE_SHIFT) - scaled_x_offset - 1) / x_step + 1 - render_x0;
int check_shift = check_size ? get_check_shift (check_size) : 0;
- y = render_y0 * y_step + floor (filter->y_offset * (1 << SCALE_SHIFT));
+ int scaled_x_offset = floor (filter->x.offset * (1 << SCALE_SHIFT));
+
+ /* Compute the index where we run off the end of the source buffer. The
+ * furthest source pixel we access at index i is:
+ *
+ * ((render_x0 + i) * x_step + scaled_x_offset) >> SCALE_SHIFT + filter->x.n - 1
+ *
+ * So, run_end_index is the smallest i for which this pixel is src_width,
+ * i.e, for which:
+ *
+ * (i + render_x0) * x_step >= ((src_width - filter->x.n + 1) << SCALE_SHIFT) - scaled_x_offset
+ *
+ */
+#define MYDIV(a,b) ((a) > 0 ? (a) / (b) : ((a) - (b) + 1) / (b)) /* Division so that -1/5 = -1 */
+
+ int run_end_x = (((src_width - filter->x.n + 1) << SCALE_SHIFT) - scaled_x_offset);
+ int run_end_index = MYDIV (run_end_x + x_step - 1, x_step) - render_x0;
+ run_end_index = MIN (run_end_index, render_x1 - render_x0);
+
+ y = render_y0 * y_step + floor (filter->y.offset * (1 << SCALE_SHIFT));
for (i = 0; i < (render_y1 - render_y0); i++)
{
+ int dest_x;
int y_start = y >> SCALE_SHIFT;
int x_start;
- int *run_weights = filter->weights + ((y >> (SCALE_SHIFT - SUBSAMPLE_BITS)) & SUBSAMPLE_MASK) * filter->n_x * filter->n_y * SUBSAMPLE;
+ int *run_weights = filter_weights +
+ ((y >> (SCALE_SHIFT - SUBSAMPLE_BITS)) & SUBSAMPLE_MASK) *
+ filter->x.n * filter->y.n * SUBSAMPLE;
guchar *new_outbuf;
guint32 tcolor1, tcolor2;
tcolor2 = color2;
}
- for (j=0; j<filter->n_y; j++)
+ for (j=0; j<filter->y.n; j++)
{
if (y_start < 0)
line_bufs[j] = (guchar *)src_buf;
while (x_start < 0 && outbuf < outbuf_end)
{
- process_pixel (run_weights + ((x >> (SCALE_SHIFT - SUBSAMPLE_BITS)) & SUBSAMPLE_MASK) * (filter->n_x * filter->n_y), filter->n_x, filter->n_y,
+ process_pixel (run_weights + ((x >> (SCALE_SHIFT - SUBSAMPLE_BITS)) & SUBSAMPLE_MASK) * (filter->x.n * filter->y.n), filter->x.n, filter->y.n,
outbuf, dest_x, dest_channels, dest_has_alpha,
line_bufs, src_channels, src_has_alpha,
x >> SCALE_SHIFT, src_width,
outbuf += dest_channels;
}
- new_outbuf = (*line_func)(run_weights, filter->n_x, filter->n_y,
- outbuf, dest_x,
- MIN (outbuf_end, dest_buf + dest_rowstride * i + run_end_index * dest_channels),
- dest_channels, dest_has_alpha,
- line_bufs, src_channels, src_has_alpha,
- x, x_step, src_width, check_size, tcolor1, tcolor2);
+ new_outbuf = (*line_func) (run_weights, filter->x.n, filter->y.n,
+ outbuf, dest_x, dest_buf + dest_rowstride *
+ i + run_end_index * dest_channels,
+ dest_channels, dest_has_alpha,
+ line_bufs, src_channels, src_has_alpha,
+ x, x_step, src_width, check_size, tcolor1,
+ tcolor2);
dest_x += (new_outbuf - outbuf) / dest_channels;
- x = dest_x * x_step + scaled_x_offset;
+ x = (dest_x - check_x + render_x0) * x_step + scaled_x_offset;
outbuf = new_outbuf;
while (outbuf < outbuf_end)
{
- process_pixel (run_weights + ((x >> (SCALE_SHIFT - SUBSAMPLE_BITS)) & SUBSAMPLE_MASK) * (filter->n_x * filter->n_y), filter->n_x, filter->n_y,
+ process_pixel (run_weights + ((x >> (SCALE_SHIFT - SUBSAMPLE_BITS)) & SUBSAMPLE_MASK) * (filter->x.n * filter->y.n), filter->x.n, filter->y.n,
outbuf, dest_x, dest_channels, dest_has_alpha,
line_bufs, src_channels, src_has_alpha,
x >> SCALE_SHIFT, src_width,
}
g_free (line_bufs);
+ g_free (filter_weights);
}
+/* Compute weights for reconstruction by replication followed by
+ * sampling with a box filter
+ */
static void
-tile_make_weights (PixopsFilter *filter, double x_scale, double y_scale, double overall_alpha)
+tile_make_weights (PixopsFilterDimension *dim,
+ double scale)
{
- int i_offset, j_offset;
-
- int n_x = ceil(1/x_scale + 1);
- int n_y = ceil(1/y_scale + 1);
+ int n = ceil (1 / scale + 1);
+ double *pixel_weights = g_new (double, SUBSAMPLE * n);
+ int offset;
+ int i;
- filter->x_offset = 0;
- filter->y_offset = 0;
- filter->n_x = n_x;
- filter->n_y = n_y;
- filter->weights = g_new (int, SUBSAMPLE * SUBSAMPLE * n_x * n_y);
+ dim->n = n;
+ dim->offset = 0;
+ dim->weights = pixel_weights;
- for (i_offset=0; i_offset<SUBSAMPLE; i_offset++)
- for (j_offset=0; j_offset<SUBSAMPLE; j_offset++)
- {
- int *pixel_weights = filter->weights + ((i_offset*SUBSAMPLE) + j_offset) * n_x * n_y;
- double x = (double)j_offset / 16;
- double y = (double)i_offset / 16;
- int i,j;
-
- for (i = 0; i < n_y; i++)
- {
- double tw, th;
-
- if (i < y)
- {
- if (i + 1 > y)
- th = MIN(i+1, y + 1/y_scale) - y;
- else
- th = 0;
- }
- else
- {
- if (y + 1/y_scale > i)
- th = MIN(i+1, y + 1/y_scale) - i;
- else
- th = 0;
- }
-
- for (j = 0; j < n_x; j++)
- {
- if (j < x)
- {
- if (j + 1 > x)
- tw = MIN(j+1, x + 1/x_scale) - x;
- else
- tw = 0;
- }
- else
- {
- if (x + 1/x_scale > j)
- tw = MIN(j+1, x + 1/x_scale) - j;
- else
- tw = 0;
- }
-
- *(pixel_weights + n_x * i + j) = 65536 * tw * x_scale * th * y_scale * overall_alpha;
- }
- }
- }
+ for (offset = 0; offset < SUBSAMPLE; offset++)
+ {
+ double x = (double)offset / SUBSAMPLE;
+ double a = x + 1 / scale;
+
+ for (i = 0; i < n; i++)
+ {
+ if (i < x)
+ {
+ if (i + 1 > x)
+ *(pixel_weights++) = (MIN (i + 1, a) - x) * scale;
+ else
+ *(pixel_weights++) = 0;
+ }
+ else
+ {
+ if (a > i)
+ *(pixel_weights++) = (MIN (i + 1, a) - i) * scale;
+ else
+ *(pixel_weights++) = 0;
+ }
+ }
+ }
}
+/* Compute weights for a filter that, for minification
+ * is the same as 'tiles', and for magnification, is bilinear
+ * reconstruction followed by a sampling with a delta function.
+ */
static void
-bilinear_make_fast_weights (PixopsFilter *filter, double x_scale, double y_scale, double overall_alpha)
+bilinear_magnify_make_weights (PixopsFilterDimension *dim,
+ double scale)
{
- int i_offset, j_offset;
- double *x_weights, *y_weights;
- int n_x, n_y;
-
- if (x_scale > 1.0) /* Bilinear */
- {
- n_x = 2;
- filter->x_offset = 0.5 * (1/x_scale - 1);
- }
- else /* Tile */
- {
- n_x = ceil(1.0 + 1.0/x_scale);
- filter->x_offset = 0.0;
- }
+ double *pixel_weights;
+ int n;
+ int offset;
+ int i;
- if (y_scale > 1.0) /* Bilinear */
+ if (scale > 1.0) /* Linear */
{
- n_y = 2;
- filter->y_offset = 0.5 * (1/y_scale - 1);
+ n = 2;
+ dim->offset = 0.5 * (1 / scale - 1);
}
- else /* Tile */
+ else /* Tile */
{
- n_y = ceil(1.0 + 1.0/y_scale);
- filter->y_offset = 0.0;
+ n = ceil (1.0 + 1.0 / scale);
+ dim->offset = 0.0;
}
- filter->n_y = n_y;
- filter->n_x = n_x;
- filter->weights = g_new (int, SUBSAMPLE * SUBSAMPLE * n_x * n_y);
-
- x_weights = g_new (double, n_x);
- y_weights = g_new (double, n_y);
+ dim->n = n;
+ dim->weights = g_new (double, SUBSAMPLE * n);
- for (i_offset=0; i_offset<SUBSAMPLE; i_offset++)
- for (j_offset=0; j_offset<SUBSAMPLE; j_offset++)
- {
- int *pixel_weights = filter->weights + ((i_offset*SUBSAMPLE) + j_offset) * n_x * n_y;
- double x = (double)j_offset / 16;
- double y = (double)i_offset / 16;
- int i,j;
-
- if (x_scale > 1.0) /* Bilinear */
- {
- for (i = 0; i < n_x; i++)
- {
- x_weights[i] = ((i == 0) ? (1 - x) : x) / x_scale;
- }
- }
- else /* Tile */
- {
- for (i = 0; i < n_x; i++)
- {
- if (i < x)
- {
- if (i + 1 > x)
- x_weights[i] = MIN(i+1, x + 1/x_scale) - x;
- else
- x_weights[i] = 0;
- }
- else
- {
- if (x + 1/x_scale > i)
- x_weights[i] = MIN(i+1, x + 1/x_scale) - i;
- else
- x_weights[i] = 0;
- }
- }
- }
-
- if (y_scale > 1.0) /* Bilinear */
- {
- for (i = 0; i < n_y; i++)
- {
- y_weights[i] = ((i == 0) ? (1 - y) : y) / y_scale;
- }
- }
- else /* Tile */
- {
- for (i = 0; i < n_y; i++)
- {
- if (i < y)
- {
- if (i + 1 > y)
- y_weights[i] = MIN(i+1, y + 1/y_scale) - y;
- else
- y_weights[i] = 0;
- }
- else
- {
- if (y + 1/y_scale > i)
- y_weights[i] = MIN(i+1, y + 1/y_scale) - i;
- else
- y_weights[i] = 0;
- }
- }
- }
-
- for (i = 0; i < n_y; i++)
- for (j = 0; j < n_x; j++)
- *(pixel_weights + n_x * i + j) = 65536 * x_weights[j] * x_scale * y_weights[i] * y_scale * overall_alpha;
- }
+ pixel_weights = dim->weights;
- g_free (x_weights);
- g_free (y_weights);
+ for (offset=0; offset < SUBSAMPLE; offset++)
+ {
+ double x = (double)offset / SUBSAMPLE;
+
+ if (scale > 1.0) /* Linear */
+ {
+ for (i = 0; i < n; i++)
+ *(pixel_weights++) = (((i == 0) ? (1 - x) : x) / scale) * scale;
+ }
+ else /* Tile */
+ {
+ double a = x + 1 / scale;
+
+ /* x
+ * ---------|--.-|----|--.-|------- SRC
+ * ------------|---------|--------- DEST
+ */
+ for (i = 0; i < n; i++)
+ {
+ if (i < x)
+ {
+ if (i + 1 > x)
+ *(pixel_weights++) = (MIN (i + 1, a) - x) * scale;
+ else
+ *(pixel_weights++) = 0;
+ }
+ else
+ {
+ if (a > i)
+ *(pixel_weights++) = (MIN (i + 1, a) - i) * scale;
+ else
+ *(pixel_weights++) = 0;
+ }
+ }
+ }
+ }
}
+/* Computes the integral from b0 to b1 of
+ *
+ * f(x) = x; 0 <= x < 1
+ * f(x) = 0; otherwise
+ *
+ * We combine two of these to compute the convolution of
+ * a box filter with a triangular spike.
+ */
static double
-bilinear_quadrant (double bx0, double bx1, double by0, double by1)
+linear_box_half (double b0, double b1)
{
- double ax0, ax1, ay0, ay1;
- double x0, x1, y0, y1;
+ double a0, a1;
+ double x0, x1;
- ax0 = 0.;
- ax1 = 1.;
- ay0 = 0.;
- ay1 = 1.;
+ a0 = 0.;
+ a1 = 1.;
- if (ax0 < bx0)
+ if (a0 < b0)
{
- if (ax1 > bx0)
- {
- x0 = bx0;
- x1 = MIN (ax1, bx1);
- }
+ if (a1 > b0)
+ {
+ x0 = b0;
+ x1 = MIN (a1, b1);
+ }
else
- return 0;
+ return 0;
}
else
{
- if (bx1 > ax0)
- {
- x0 = ax0;
- x1 = MIN (ax1, bx1);
- }
+ if (b1 > a0)
+ {
+ x0 = a0;
+ x1 = MIN (a1, b1);
+ }
else
- return 0;
+ return 0;
}
- if (ay0 < by0)
- {
- if (ay1 > by0)
- {
- y0 = by0;
- y1 = MIN (ay1, by1);
- }
- else
- return 0;
- }
- else
- {
- if (by1 > ay0)
- {
- y0 = ay0;
- y1 = MIN (ay1, by1);
- }
- else
- return 0;
- }
-
- return 0.25 * (x1*x1 - x0*x0) * (y1*y1 - y0*y0);
+ return 0.5 * (x1*x1 - x0*x0);
}
+/* Compute weights for reconstructing with bilinear
+ * interpolation, then sampling with a box filter
+ */
static void
-bilinear_make_weights (PixopsFilter *filter, double x_scale, double y_scale, double overall_alpha)
+bilinear_box_make_weights (PixopsFilterDimension *dim,
+ double scale)
{
- int i_offset, j_offset;
+ int n = ceil (1/scale + 3.0);
+ double *pixel_weights = g_new (double, SUBSAMPLE * n);
+ double w;
+ int offset, i;
- int n_x = ceil(1/x_scale + 2.0);
- int n_y = ceil(1/y_scale + 2.0);
+ dim->offset = -1.0;
+ dim->n = n;
+ dim->weights = pixel_weights;
- filter->x_offset = -1.0;
- filter->y_offset = -1.0;
- filter->n_x = n_x;
- filter->n_y = n_y;
-
- filter->weights = g_new (int, SUBSAMPLE * SUBSAMPLE * n_x * n_y);
+ for (offset = 0; offset < SUBSAMPLE; offset++)
+ {
+ double x = (double)offset / SUBSAMPLE;
+ double a = x + 1 / scale;
- for (i_offset=0; i_offset<SUBSAMPLE; i_offset++)
- for (j_offset=0; j_offset<SUBSAMPLE; j_offset++)
- {
- int *pixel_weights = filter->weights + ((i_offset*SUBSAMPLE) + j_offset) * n_x * n_y;
- double x = (double)j_offset / 16;
- double y = (double)i_offset / 16;
- int i,j;
-
- for (i = 0; i < n_y; i++)
- for (j = 0; j < n_x; j++)
- {
- double w;
+ for (i = 0; i < n; i++)
+ {
+ w = linear_box_half (0.5 + i - a, 0.5 + i - x);
+ w += linear_box_half (1.5 + x - i, 1.5 + a - i);
+
+ *(pixel_weights++) = w * scale;
+ }
+ }
+}
- w = bilinear_quadrant (0.5 + j - (x + 1 / x_scale), 0.5 + j - x, 0.5 + i - (y + 1 / y_scale), 0.5 + i - y);
- w += bilinear_quadrant (1.5 + x - j, 1.5 + (x + 1 / x_scale) - j, 0.5 + i - (y + 1 / y_scale), 0.5 + i - y);
- w += bilinear_quadrant (0.5 + j - (x + 1 / x_scale), 0.5 + j - x, 1.5 + y - i, 1.5 + (y + 1 / y_scale) - i);
- w += bilinear_quadrant (1.5 + x - j, 1.5 + (x + 1 / x_scale) - j, 1.5 + y - i, 1.5 + (y + 1 / y_scale) - i);
-
- *(pixel_weights + n_x * i + j) = 65536 * w * x_scale * y_scale * overall_alpha;
- }
- }
+static void
+make_weights (PixopsFilter *filter,
+ PixopsInterpType interp_type,
+ double scale_x,
+ double scale_y)
+{
+ switch (interp_type)
+ {
+ case PIXOPS_INTERP_NEAREST:
+ g_assert_not_reached ();
+ break;
+
+ case PIXOPS_INTERP_TILES:
+ tile_make_weights (&filter->x, scale_x);
+ tile_make_weights (&filter->y, scale_y);
+ break;
+
+ case PIXOPS_INTERP_BILINEAR:
+ bilinear_magnify_make_weights (&filter->x, scale_x);
+ bilinear_magnify_make_weights (&filter->y, scale_y);
+ break;
+
+ case PIXOPS_INTERP_HYPER:
+ bilinear_box_make_weights (&filter->x, scale_x);
+ bilinear_box_make_weights (&filter->y, scale_y);
+ break;
+ }
}
-void
-pixops_composite_color (guchar *dest_buf,
- int render_x0,
- int render_y0,
- int render_x1,
- int render_y1,
- int dest_rowstride,
- int dest_channels,
- gboolean dest_has_alpha,
- const guchar *src_buf,
- int src_width,
- int src_height,
- int src_rowstride,
- int src_channels,
- gboolean src_has_alpha,
- double scale_x,
- double scale_y,
- GdkInterpType interp_type,
- int overall_alpha,
- int check_x,
- int check_y,
- int check_size,
- guint32 color1,
- guint32 color2)
+static void
+_pixops_composite_color_real (guchar *dest_buf,
+ int render_x0,
+ int render_y0,
+ int render_x1,
+ int render_y1,
+ int dest_rowstride,
+ int dest_channels,
+ gboolean dest_has_alpha,
+ const guchar *src_buf,
+ int src_width,
+ int src_height,
+ int src_rowstride,
+ int src_channels,
+ gboolean src_has_alpha,
+ double scale_x,
+ double scale_y,
+ PixopsInterpType interp_type,
+ int overall_alpha,
+ int check_x,
+ int check_y,
+ int check_size,
+ guint32 color1,
+ guint32 color2)
{
PixopsFilter filter;
PixopsLineFunc line_func;
#ifdef USE_MMX
- gboolean found_mmx = pixops_have_mmx();
+ gboolean found_mmx = _pixops_have_mmx ();
#endif
g_return_if_fail (!(dest_channels == 3 && dest_has_alpha));
if (scale_x == 0 || scale_y == 0)
return;
- if (!src_has_alpha && overall_alpha == 255)
- pixops_scale (dest_buf, render_x0, render_y0, render_x1, render_y1,
- dest_rowstride, dest_channels, dest_has_alpha,
- src_buf, src_width, src_height, src_rowstride, src_channels,
- src_has_alpha, scale_x, scale_y, interp_type);
-
- switch (interp_type)
+ if (interp_type == PIXOPS_INTERP_NEAREST)
{
- case GDK_INTERP_NEAREST:
- pixops_composite_color_nearest (dest_buf, render_x0, render_y0, render_x1, render_y1,
- dest_rowstride, dest_channels, dest_has_alpha,
- src_buf, src_width, src_height, src_rowstride, src_channels, src_has_alpha,
- scale_x, scale_y, overall_alpha,
- check_x, check_y, check_size, color1, color2);
+ pixops_composite_color_nearest (dest_buf, render_x0, render_y0,
+ render_x1, render_y1, dest_rowstride,
+ dest_channels, dest_has_alpha, src_buf,
+ src_width, src_height, src_rowstride,
+ src_channels, src_has_alpha, scale_x,
+ scale_y, overall_alpha, check_x, check_y,
+ check_size, color1, color2);
return;
-
- case GDK_INTERP_TILES:
- tile_make_weights (&filter, scale_x, scale_y, overall_alpha / 255.);
- break;
-
- case GDK_INTERP_BILINEAR:
- bilinear_make_fast_weights (&filter, scale_x, scale_y, overall_alpha / 255.);
- break;
-
- case GDK_INTERP_HYPER:
- bilinear_make_weights (&filter, scale_x, scale_y, overall_alpha / 255.);
- break;
}
+
+ filter.overall_alpha = overall_alpha / 255.;
+ make_weights (&filter, interp_type, scale_x, scale_y);
#ifdef USE_MMX
- if (filter.n_x == 2 && filter.n_y == 2 &&
- dest_channels == 4 && src_channels == 4 && src_has_alpha && !dest_has_alpha && found_mmx)
+ if (filter.x.n == 2 && filter.y.n == 2 &&
+ dest_channels == 4 && src_channels == 4 &&
+ src_has_alpha && !dest_has_alpha && found_mmx)
line_func = composite_line_color_22_4a4_mmx_stub;
else
-#endif
+#endif
line_func = composite_line_color;
pixops_process (dest_buf, render_x0, render_y0, render_x1, render_y1,
src_has_alpha, scale_x, scale_y, check_x, check_y, check_size, color1, color2,
&filter, line_func, composite_pixel_color);
- g_free (filter.weights);
+ g_free (filter.x.weights);
+ g_free (filter.y.weights);
}
void
-pixops_composite (guchar *dest_buf,
- int render_x0,
- int render_y0,
- int render_x1,
- int render_y1,
- int dest_rowstride,
- int dest_channels,
- gboolean dest_has_alpha,
- const guchar *src_buf,
- int src_width,
- int src_height,
- int src_rowstride,
- int src_channels,
- gboolean src_has_alpha,
- double scale_x,
- double scale_y,
- GdkInterpType interp_type,
- int overall_alpha)
+_pixops_composite_color (guchar *dest_buf,
+ int dest_width,
+ int dest_height,
+ int dest_rowstride,
+ int dest_channels,
+ gboolean dest_has_alpha,
+ const guchar *src_buf,
+ int src_width,
+ int src_height,
+ int src_rowstride,
+ int src_channels,
+ gboolean src_has_alpha,
+ int dest_x,
+ int dest_y,
+ int dest_region_width,
+ int dest_region_height,
+ double offset_x,
+ double offset_y,
+ double scale_x,
+ double scale_y,
+ PixopsInterpType interp_type,
+ int overall_alpha,
+ int check_x,
+ int check_y,
+ int check_size,
+ guint32 color1,
+ guint32 color2)
+{
+ guchar *new_dest_buf;
+ int render_x0;
+ int render_y0;
+ int render_x1;
+ int render_y1;
+
+ if (!src_has_alpha && overall_alpha == 255)
+ {
+ _pixops_scale (dest_buf, dest_width, dest_height, dest_rowstride,
+ dest_channels, dest_has_alpha, src_buf, src_width,
+ src_height, src_rowstride, src_channels, src_has_alpha,
+ dest_x, dest_y, dest_region_width, dest_region_height,
+ offset_x, offset_y, scale_x, scale_y, interp_type);
+ return;
+ }
+
+ new_dest_buf = dest_buf + dest_y * dest_rowstride + dest_x *
+ dest_channels;
+ render_x0 = dest_x - offset_x;
+ render_y0 = dest_y - offset_y;
+ render_x1 = dest_x + dest_region_width - offset_x;
+ render_y1 = dest_y + dest_region_height - offset_y;
+
+ _pixops_composite_color_real (new_dest_buf, render_x0, render_y0, render_x1,
+ render_y1, dest_rowstride, dest_channels,
+ dest_has_alpha, src_buf, src_width,
+ src_height, src_rowstride, src_channels,
+ src_has_alpha, scale_x, scale_y,
+ (PixopsInterpType)interp_type, overall_alpha,
+ check_x, check_y, check_size, color1, color2);
+}
+
+/**
+ * _pixops_composite_real:
+ * @dest_buf: pointer to location to store result
+ * @render_x0: x0 of region of scaled source to store into @dest_buf
+ * @render_y0: y0 of region of scaled source to store into @dest_buf
+ * @render_x1: x1 of region of scaled source to store into @dest_buf
+ * @render_y1: y1 of region of scaled source to store into @dest_buf
+ * @dest_rowstride: rowstride of @dest_buf
+ * @dest_channels: number of channels in @dest_buf
+ * @dest_has_alpha: whether @dest_buf has alpha
+ * @src_buf: pointer to source pixels
+ * @src_width: width of source (used for clipping)
+ * @src_height: height of source (used for clipping)
+ * @src_rowstride: rowstride of source
+ * @src_channels: number of channels in @src_buf
+ * @src_has_alpha: whether @src_buf has alpha
+ * @scale_x: amount to scale source by in X direction
+ * @scale_y: amount to scale source by in Y direction
+ * @interp_type: type of enumeration
+ * @overall_alpha: overall alpha factor to multiply source by
+ *
+ * Scale source buffer by scale_x / scale_y, then composite a given rectangle
+ * of the result into the destination buffer.
+ **/
+static void
+_pixops_composite_real (guchar *dest_buf,
+ int render_x0,
+ int render_y0,
+ int render_x1,
+ int render_y1,
+ int dest_rowstride,
+ int dest_channels,
+ gboolean dest_has_alpha,
+ const guchar *src_buf,
+ int src_width,
+ int src_height,
+ int src_rowstride,
+ int src_channels,
+ gboolean src_has_alpha,
+ double scale_x,
+ double scale_y,
+ PixopsInterpType interp_type,
+ int overall_alpha)
{
PixopsFilter filter;
PixopsLineFunc line_func;
#ifdef USE_MMX
- gboolean found_mmx = pixops_have_mmx();
+ gboolean found_mmx = _pixops_have_mmx ();
#endif
g_return_if_fail (!(dest_channels == 3 && dest_has_alpha));
if (scale_x == 0 || scale_y == 0)
return;
- if (!src_has_alpha && overall_alpha == 255)
- pixops_scale (dest_buf, render_x0, render_y0, render_x1, render_y1,
- dest_rowstride, dest_channels, dest_has_alpha,
- src_buf, src_width, src_height, src_rowstride, src_channels,
- src_has_alpha, scale_x, scale_y, interp_type);
-
- switch (interp_type)
+ if (interp_type == PIXOPS_INTERP_NEAREST)
{
- case GDK_INTERP_NEAREST:
- pixops_composite_nearest (dest_buf, render_x0, render_y0, render_x1, render_y1,
- dest_rowstride, dest_channels, dest_has_alpha,
- src_buf, src_width, src_height, src_rowstride, src_channels,
- src_has_alpha, scale_x, scale_y, overall_alpha);
+ pixops_composite_nearest (dest_buf, render_x0, render_y0, render_x1,
+ render_y1, dest_rowstride, dest_channels,
+ dest_has_alpha, src_buf, src_width, src_height,
+ src_rowstride, src_channels, src_has_alpha,
+ scale_x, scale_y, overall_alpha);
return;
-
- case GDK_INTERP_TILES:
- tile_make_weights (&filter, scale_x, scale_y, overall_alpha / 255.);
- break;
-
- case GDK_INTERP_BILINEAR:
- bilinear_make_fast_weights (&filter, scale_x, scale_y, overall_alpha / 255.);
- break;
-
- case GDK_INTERP_HYPER:
- bilinear_make_weights (&filter, scale_x, scale_y, overall_alpha / 255.);
- break;
}
+
+ filter.overall_alpha = overall_alpha / 255.;
+ make_weights (&filter, interp_type, scale_x, scale_y);
- if (filter.n_x == 2 && filter.n_y == 2 &&
- dest_channels == 4 && src_channels == 4 && src_has_alpha && !dest_has_alpha)
+ if (filter.x.n == 2 && filter.y.n == 2 && dest_channels == 4 &&
+ src_channels == 4 && src_has_alpha && !dest_has_alpha)
{
#ifdef USE_MMX
if (found_mmx)
src_has_alpha, scale_x, scale_y, 0, 0, 0, 0, 0,
&filter, line_func, composite_pixel);
- g_free (filter.weights);
+ g_free (filter.x.weights);
+ g_free (filter.y.weights);
}
void
-pixops_scale (guchar *dest_buf,
- int render_x0,
- int render_y0,
- int render_x1,
- int render_y1,
- int dest_rowstride,
- int dest_channels,
- gboolean dest_has_alpha,
- const guchar *src_buf,
- int src_width,
- int src_height,
- int src_rowstride,
- int src_channels,
- gboolean src_has_alpha,
- double scale_x,
- double scale_y,
- GdkInterpType interp_type)
+_pixops_composite (guchar *dest_buf,
+ int dest_width,
+ int dest_height,
+ int dest_rowstride,
+ int dest_channels,
+ int dest_has_alpha,
+ const guchar *src_buf,
+ int src_width,
+ int src_height,
+ int src_rowstride,
+ int src_channels,
+ int src_has_alpha,
+ int dest_x,
+ int dest_y,
+ int dest_region_width,
+ int dest_region_height,
+ double offset_x,
+ double offset_y,
+ double scale_x,
+ double scale_y,
+ PixopsInterpType interp_type,
+ int overall_alpha)
+{
+ guchar *new_dest_buf;
+ int render_x0;
+ int render_y0;
+ int render_x1;
+ int render_y1;
+
+ if (!src_has_alpha && overall_alpha == 255)
+ {
+ _pixops_scale (dest_buf, dest_width, dest_height, dest_rowstride,
+ dest_channels, dest_has_alpha, src_buf, src_width,
+ src_height, src_rowstride, src_channels, src_has_alpha,
+ dest_x, dest_y, dest_region_width, dest_region_height,
+ offset_x, offset_y, scale_x, scale_y, interp_type);
+ return;
+ }
+
+#ifdef USE_MEDIALIB
+ pixops_medialib_composite (dest_buf, dest_width, dest_height, dest_rowstride,
+ dest_channels, dest_has_alpha, src_buf,
+ src_width, src_height, src_rowstride,
+ src_channels, src_has_alpha, dest_x, dest_y,
+ dest_region_width, dest_region_height, offset_x,
+ offset_y, scale_x, scale_y,
+ (PixopsInterpType)interp_type, overall_alpha);
+ return;
+#endif
+
+ new_dest_buf = dest_buf + dest_y * dest_rowstride + dest_x * dest_channels;
+ render_x0 = dest_x - offset_x;
+ render_y0 = dest_y - offset_y;
+ render_x1 = dest_x + dest_region_width - offset_x;
+ render_y1 = dest_y + dest_region_height - offset_y;
+
+ _pixops_composite_real (new_dest_buf, render_x0, render_y0, render_x1,
+ render_y1, dest_rowstride, dest_channels,
+ dest_has_alpha, src_buf, src_width, src_height,
+ src_rowstride, src_channels, src_has_alpha, scale_x,
+ scale_y, (PixopsInterpType)interp_type,
+ overall_alpha);
+}
+
+#ifdef USE_MEDIALIB
+static void
+medialib_get_interpolation (mlInterp * ml_interp,
+ PixopsInterpType interp_type,
+ double scale_x,
+ double scale_y,
+ double overall_alpha)
+{
+ mlib_s32 leftPadding, topPadding;
+ ml_interp->interp_table = NULL;
+
+ /*
+ * medialib 2.1 and later supports scaling with user-defined interpolation
+ * tables, so this logic is used.
+ *
+ * bilinear_magnify_make_weights builds an interpolation table of size 2x2 if
+ * the scale factor >= 1.0 and "ceil (1.0 + 1.0/scale)" otherwise. These map
+ * most closely to MLIB_BILINEAR, which uses an interpolation table of size
+ * 2x2.
+ *
+ * tile_make_weights builds an interpolation table of size 2x2 if the scale
+ * factor >= 1.0 and "ceil (1.0 + 1.0/scale)" otherwise. These map most
+ * closely to MLIB_BILINEAR, which uses an interpolation table of size 2x2.
+ *
+ * bilinear_box_make_weights builds an interpolation table of size 4x4 if the
+ * scale factor >= 1.0 and "ceil (1.0 + 1.0/scale)" otherwise. These map most
+ * closely to MLIB_BICUBIC, which uses an interpolation table of size 4x4.
+ *
+ * PIXOPS_INTERP_NEAREST calls pixops_scale_nearest which does not use an
+ * interpolation table. This maps to MLIB_NEAREST.
+ */
+ switch (interp_type)
+ {
+ case PIXOPS_INTERP_BILINEAR:
+ bilinear_magnify_make_weights (&(ml_interp->po_filter.x), scale_x);
+ bilinear_magnify_make_weights (&(ml_interp->po_filter.y), scale_y);
+ leftPadding = 0;
+ topPadding = 0;
+
+ if (scale_x <= 1.0)
+ ml_interp->tx = 0.5 * (1 - scale_x);
+ else
+ ml_interp->tx = 0.0;
+
+ if (scale_y <= 1.0)
+ ml_interp->ty = 0.5 * (1 - scale_y);
+ else
+ ml_interp->ty = 0.0;
+
+ break;
+
+ case PIXOPS_INTERP_TILES:
+ tile_make_weights (&(ml_interp->po_filter.x), scale_x);
+ tile_make_weights (&(ml_interp->po_filter.y), scale_y);
+ leftPadding = 0;
+ topPadding = 0;
+ ml_interp->tx = 0.5 * (1 - scale_x);
+ ml_interp->ty = 0.5 * (1 - scale_y);
+ break;
+
+ case PIXOPS_INTERP_HYPER:
+ bilinear_box_make_weights (&(ml_interp->po_filter.x), scale_x);
+ bilinear_box_make_weights (&(ml_interp->po_filter.y), scale_y);
+ leftPadding = 1;
+ topPadding = 1;
+ ml_interp->tx = 0.5 * (1 - scale_x);
+ ml_interp->ty = 0.5 * (1 - scale_y);
+ break;
+
+ case PIXOPS_INTERP_NEAREST:
+ default:
+ /*
+ * Note that this function should not be called in the
+ * PIXOPS_INTERP_NEAREST case since it does not use an interpolation
+ * table.
+ */
+ g_assert_not_reached ();
+ break;
+ }
+
+ /*
+ * If overall_alpha is not 1.0, then multiply the vectors built by the
+ * sqrt (overall_alpha). This will cause overall_alpha to get evenly
+ * blended across both axis.
+ *
+ * Note there is no need to multiply the vectors built by the various
+ * make-weight functions by sqrt (overall_alpha) since the make-weight
+ * functions are called with overall_alpha hardcoded to 1.0.
+ */
+ if (overall_alpha != 1.0)
+ {
+ double sqrt_alpha = sqrt (overall_alpha);
+ int i;
+
+ for (i=0; i < SUBSAMPLE * ml_interp->po_filter.x.n; i++)
+ ml_interp->po_filter.x.weights[i] *= sqrt_alpha;
+ for (i=0; i < SUBSAMPLE * ml_interp->po_filter.y.n; i++)
+ ml_interp->po_filter.y.weights[i] *= sqrt_alpha;
+ }
+
+ ml_interp->interp_table = (void *) mlib_ImageInterpTableCreate (MLIB_DOUBLE,
+ ml_interp->po_filter.x.n, ml_interp->po_filter.y.n, leftPadding,
+ topPadding, SUBSAMPLE_BITS, SUBSAMPLE_BITS, 8,
+ ml_interp->po_filter.x.weights, ml_interp->po_filter.y.weights);
+
+ g_free (ml_interp->po_filter.x.weights);
+ g_free (ml_interp->po_filter.y.weights);
+}
+
+static void
+pixops_medialib_composite (guchar *dest_buf,
+ int dest_width,
+ int dest_height,
+ int dest_rowstride,
+ int dest_channels,
+ int dest_has_alpha,
+ const guchar *src_buf,
+ int src_width,
+ int src_height,
+ int src_rowstride,
+ int src_channels,
+ int src_has_alpha,
+ int dest_x,
+ int dest_y,
+ int dest_region_width,
+ int dest_region_height,
+ double offset_x,
+ double offset_y,
+ double scale_x,
+ double scale_y,
+ PixopsInterpType interp_type,
+ int overall_alpha)
+{
+ mlib_blend blend;
+ g_return_if_fail (!(dest_channels == 3 && dest_has_alpha));
+ g_return_if_fail (!(src_channels == 3 && src_has_alpha));
+
+ if (scale_x == 0 || scale_y == 0)
+ return;
+
+ if (!medialib_initialized)
+ _pixops_use_medialib ();
+
+ if (!use_medialib)
+ {
+ /* Use non-mediaLib version */
+ _pixops_composite_real (dest_buf + dest_y * dest_rowstride + dest_x *
+ dest_channels, dest_x - offset_x, dest_y -
+ offset_y, dest_x + dest_region_width - offset_x,
+ dest_y + dest_region_height - offset_y,
+ dest_rowstride, dest_channels, dest_has_alpha,
+ src_buf, src_width, src_height, src_rowstride,
+ src_channels, src_has_alpha, scale_x, scale_y,
+ interp_type, overall_alpha);
+ }
+ else
+ {
+ mlInterp ml_interp;
+ mlib_image img_src, img_dest;
+ double ml_offset_x, ml_offset_y;
+
+ if (!src_has_alpha && overall_alpha == 255 &&
+ dest_channels <= src_channels)
+ {
+ pixops_medialib_scale (dest_buf, dest_region_width,
+ dest_region_height, dest_rowstride,
+ dest_channels, dest_has_alpha, src_buf,
+ src_width, src_height, src_rowstride,
+ src_channels, src_has_alpha, dest_x, dest_y,
+ dest_region_width, dest_region_height,
+ offset_x, offset_y, scale_x, scale_y,
+ interp_type);
+ return;
+ }
+
+ mlib_ImageSetStruct (&img_src, MLIB_BYTE, src_channels,
+ src_width, src_height, src_rowstride, src_buf);
+
+ if (dest_x == 0 && dest_y == 0 &&
+ dest_width == dest_region_width &&
+ dest_height == dest_region_height)
+ {
+ mlib_ImageSetStruct (&img_dest, MLIB_BYTE, dest_channels,
+ dest_width, dest_height, dest_rowstride,
+ dest_buf);
+ }
+ else
+ {
+ mlib_u8 *data = dest_buf + (dest_y * dest_rowstride) +
+ (dest_x * dest_channels);
+
+ mlib_ImageSetStruct (&img_dest, MLIB_BYTE, dest_channels,
+ dest_region_width, dest_region_height,
+ dest_rowstride, data);
+ }
+
+ ml_offset_x = floor (offset_x) - dest_x;
+ ml_offset_y = floor (offset_y) - dest_y;
+
+ if (interp_type == PIXOPS_INTERP_NEAREST)
+ {
+ blend = src_has_alpha ? MLIB_BLEND_GTK_SRC_OVER2 : MLIB_BLEND_GTK_SRC;
+
+ mlib_ImageZoomTranslateBlend (&img_dest,
+ &img_src,
+ scale_x,
+ scale_y,
+ ml_offset_x,
+ ml_offset_y,
+ MLIB_NEAREST,
+ MLIB_EDGE_SRC_EXTEND_INDEF,
+ blend,
+ overall_alpha,
+ 1);
+ }
+ else
+ {
+ blend = src_has_alpha ? MLIB_BLEND_GTK_SRC_OVER : MLIB_BLEND_GTK_SRC;
+
+ if (interp_type == PIXOPS_INTERP_BILINEAR &&
+ scale_x > 1.0 && scale_y > 1.0)
+ {
+ mlib_ImageZoomTranslateBlend (&img_dest,
+ &img_src,
+ scale_x,
+ scale_y,
+ ml_offset_x,
+ ml_offset_y,
+ MLIB_BILINEAR,
+ MLIB_EDGE_SRC_EXTEND_INDEF,
+ blend,
+ overall_alpha,
+ 1);
+ }
+ else
+ {
+ medialib_get_interpolation (&ml_interp, interp_type, scale_x,
+ scale_y, overall_alpha/255.0);
+
+ if (ml_interp.interp_table != NULL)
+ {
+ mlib_ImageZoomTranslateTableBlend (&img_dest,
+ &img_src,
+ scale_x,
+ scale_y,
+ ml_offset_x + ml_interp.tx,
+ ml_offset_y + ml_interp.ty,
+ ml_interp.interp_table,
+ MLIB_EDGE_SRC_EXTEND_INDEF,
+ blend,
+ 1);
+ mlib_ImageInterpTableDelete (ml_interp.interp_table);
+ }
+ else
+ {
+ /* Should not happen - Use non-mediaLib version */
+ _pixops_composite_real (dest_buf + dest_y * dest_rowstride +
+ dest_x * dest_channels,
+ dest_x - offset_x, dest_y - offset_y,
+ dest_x + dest_region_width - offset_x,
+ dest_y + dest_region_height - offset_y,
+ dest_rowstride, dest_channels,
+ dest_has_alpha, src_buf, src_width,
+ src_height, src_rowstride,
+ src_channels, src_has_alpha, scale_x,
+ scale_y, interp_type, overall_alpha);
+ }
+ }
+ }
+ }
+}
+#endif
+
+static void
+_pixops_scale_real (guchar *dest_buf,
+ int render_x0,
+ int render_y0,
+ int render_x1,
+ int render_y1,
+ int dest_rowstride,
+ int dest_channels,
+ gboolean dest_has_alpha,
+ const guchar *src_buf,
+ int src_width,
+ int src_height,
+ int src_rowstride,
+ int src_channels,
+ gboolean src_has_alpha,
+ double scale_x,
+ double scale_y,
+ PixopsInterpType interp_type)
{
PixopsFilter filter;
PixopsLineFunc line_func;
#ifdef USE_MMX
- gboolean found_mmx = pixops_have_mmx();
+ gboolean found_mmx = _pixops_have_mmx ();
#endif
g_return_if_fail (!(dest_channels == 3 && dest_has_alpha));
if (scale_x == 0 || scale_y == 0)
return;
- switch (interp_type)
+ if (interp_type == PIXOPS_INTERP_NEAREST)
{
- case GDK_INTERP_NEAREST:
- pixops_scale_nearest (dest_buf, render_x0, render_y0, render_x1, render_y1,
- dest_rowstride, dest_channels, dest_has_alpha,
- src_buf, src_width, src_height, src_rowstride, src_channels, src_has_alpha,
+ pixops_scale_nearest (dest_buf, render_x0, render_y0, render_x1,
+ render_y1, dest_rowstride, dest_channels,
+ dest_has_alpha, src_buf, src_width, src_height,
+ src_rowstride, src_channels, src_has_alpha,
scale_x, scale_y);
return;
-
- case GDK_INTERP_TILES:
- tile_make_weights (&filter, scale_x, scale_y, 1.0);
- break;
-
- case GDK_INTERP_BILINEAR:
- bilinear_make_fast_weights (&filter, scale_x, scale_y, 1.0);
- break;
-
- case GDK_INTERP_HYPER:
- bilinear_make_weights (&filter, scale_x, scale_y, 1.0);
- break;
}
+
+ filter.overall_alpha = 1.0;
+ make_weights (&filter, interp_type, scale_x, scale_y);
+ if (filter.x.n == 2 && filter.y.n == 2 && dest_channels == 3 && src_channels == 3)
+ {
#ifdef USE_MMX
- if (filter.n_x == 2 && filter.n_y == 2 &&
- found_mmx && dest_channels == 3 && src_channels == 3)
- line_func = scale_line_22_33_mmx_stub;
+ if (found_mmx)
+ line_func = scale_line_22_33_mmx_stub;
+ else
+#endif
+ line_func = scale_line_22_33;
+ }
else
-#endif
line_func = scale_line;
pixops_process (dest_buf, render_x0, render_y0, render_x1, render_y1,
src_has_alpha, scale_x, scale_y, 0, 0, 0, 0, 0,
&filter, line_func, scale_pixel);
- g_free (filter.weights);
+ g_free (filter.x.weights);
+ g_free (filter.y.weights);
}
+void
+_pixops_scale (guchar *dest_buf,
+ int dest_width,
+ int dest_height,
+ int dest_rowstride,
+ int dest_channels,
+ int dest_has_alpha,
+ const guchar *src_buf,
+ int src_width,
+ int src_height,
+ int src_rowstride,
+ int src_channels,
+ int src_has_alpha,
+ int dest_x,
+ int dest_y,
+ int dest_region_width,
+ int dest_region_height,
+ double offset_x,
+ double offset_y,
+ double scale_x,
+ double scale_y,
+ PixopsInterpType interp_type)
+{
+ guchar *new_dest_buf;
+ int render_x0;
+ int render_y0;
+ int render_x1;
+ int render_y1;
+
+#ifdef USE_MEDIALIB
+ pixops_medialib_scale (dest_buf, dest_width, dest_height, dest_rowstride,
+ dest_channels, dest_has_alpha, src_buf, src_width,
+ src_height, src_rowstride, src_channels,
+ src_has_alpha, dest_x, dest_y, dest_region_width,
+ dest_region_height, offset_x, offset_y, scale_x,
+ scale_y, (PixopsInterpType)interp_type);
+ return;
+#endif
+
+ new_dest_buf = dest_buf + dest_y * dest_rowstride + dest_x * dest_channels;
+ render_x0 = dest_x - offset_x;
+ render_y0 = dest_y - offset_y;
+ render_x1 = dest_x + dest_region_width - offset_x;
+ render_y1 = dest_y + dest_region_height - offset_y;
+
+ _pixops_scale_real (new_dest_buf, render_x0, render_y0, render_x1,
+ render_y1, dest_rowstride, dest_channels,
+ dest_has_alpha, src_buf, src_width, src_height,
+ src_rowstride, src_channels, src_has_alpha,
+ scale_x, scale_y, (PixopsInterpType)interp_type);
+}
+
+#ifdef USE_MEDIALIB
+static void
+pixops_medialib_scale (guchar *dest_buf,
+ int dest_width,
+ int dest_height,
+ int dest_rowstride,
+ int dest_channels,
+ int dest_has_alpha,
+ const guchar *src_buf,
+ int src_width,
+ int src_height,
+ int src_rowstride,
+ int src_channels,
+ int src_has_alpha,
+ int dest_x,
+ int dest_y,
+ int dest_region_width,
+ int dest_region_height,
+ double offset_x,
+ double offset_y,
+ double scale_x,
+ double scale_y,
+ PixopsInterpType interp_type)
+{
+ if (scale_x == 0 || scale_y == 0)
+ return;
+
+ if (!medialib_initialized)
+ _pixops_use_medialib ();
+
+ /*
+ * We no longer support mediaLib 2.1 because it has a core dumping problem
+ * in the mlib_ImageZoomTranslateTable function that has been corrected in
+ * 2.2. Although the mediaLib_zoom function could be used, it does not
+ * work properly if the source and destination images have different
+ * values for "has_alpha" or "num_channels". The complicated if-logic
+ * required to support both versions is not worth supporting
+ * mediaLib 2.1 moving forward.
+ */
+ if (!use_medialib)
+ {
+ _pixops_scale_real (dest_buf + dest_y * dest_rowstride + dest_x *
+ dest_channels, dest_x - offset_x, dest_y - offset_y,
+ dest_x + dest_region_width - offset_x,
+ dest_y + dest_region_height - offset_y,
+ dest_rowstride, dest_channels, dest_has_alpha,
+ src_buf, src_width, src_height, src_rowstride,
+ src_channels, src_has_alpha, scale_x, scale_y,
+ interp_type);
+ }
+ else
+ {
+ mlInterp ml_interp;
+ mlib_image img_orig_src, img_src, img_dest;
+ double ml_offset_x, ml_offset_y;
+ guchar *tmp_buf = NULL;
+
+ mlib_ImageSetStruct (&img_orig_src, MLIB_BYTE, src_channels, src_width,
+ src_height, src_rowstride, src_buf);
+
+ if (dest_x == 0 && dest_y == 0 &&
+ dest_width == dest_region_width &&
+ dest_height == dest_region_height)
+ {
+ mlib_ImageSetStruct (&img_dest, MLIB_BYTE, dest_channels,
+ dest_width, dest_height, dest_rowstride,
+ dest_buf);
+ }
+ else
+ {
+ mlib_u8 *data = dest_buf + (dest_y * dest_rowstride) +
+ (dest_x * dest_channels);
+
+ mlib_ImageSetStruct (&img_dest, MLIB_BYTE, dest_channels,
+ dest_region_width, dest_region_height,
+ dest_rowstride, data);
+ }
+
+ ml_offset_x = floor (offset_x) - dest_x;
+ ml_offset_y = floor (offset_y) - dest_y;
+
+ /*
+ * Note that zoomTranslate and zoomTranslateTable are faster
+ * than zoomTranslateBlend and zoomTranslateTableBlend. However
+ * the faster functions only work in the following case:
+ *
+ * if (src_channels == dest_channels &&
+ * (!src_alpha && interp_table != PIXOPS_INTERP_NEAREST))
+ *
+ * We use the faster versions if we can.
+ *
+ * Note when the interp_type is BILINEAR and the interpolation
+ * table will be size 2x2 (when both x/y scale factors > 1.0),
+ * then we do not bother building the interpolation table. In
+ * this case we can just use MLIB_BILINEAR, which is faster than
+ * using a specified interpolation table.
+ */
+ img_src = img_orig_src;
+
+ if (!src_has_alpha)
+ {
+ if (src_channels > dest_channels)
+ {
+ int channels = 3;
+ int rowstride = (channels * src_width + 3) & ~3;
+
+ tmp_buf = g_malloc (src_rowstride * src_height);
+
+ if (src_buf != NULL)
+ {
+ src_channels = channels;
+ src_rowstride = rowstride;
+
+ mlib_ImageSetStruct (&img_src, MLIB_BYTE, src_channels,
+ src_width, src_height, src_rowstride,
+ tmp_buf);
+ mlib_ImageChannelExtract (&img_src, &img_orig_src, 0xE);
+ }
+ }
+ }
+
+ if (interp_type == PIXOPS_INTERP_NEAREST)
+ {
+ if (src_channels == dest_channels)
+ {
+ mlib_ImageZoomTranslate (&img_dest,
+ &img_src,
+ scale_x,
+ scale_y,
+ ml_offset_x,
+ ml_offset_y,
+ MLIB_NEAREST,
+ MLIB_EDGE_SRC_EXTEND_INDEF);
+ }
+ else
+ {
+ mlib_ImageZoomTranslateBlend (&img_dest,
+ &img_src,
+ scale_x,
+ scale_y,
+ ml_offset_x,
+ ml_offset_y,
+ MLIB_NEAREST,
+ MLIB_EDGE_SRC_EXTEND_INDEF,
+ MLIB_BLEND_GTK_SRC,
+ 1.0,
+ 1);
+ }
+ }
+ else if (src_channels == dest_channels && !src_has_alpha)
+ {
+ if (interp_type == PIXOPS_INTERP_BILINEAR &&
+ scale_x > 1.0 && scale_y > 1.0)
+ {
+ mlib_ImageZoomTranslate (&img_dest,
+ &img_src,
+ scale_x,
+ scale_y,
+ ml_offset_x,
+ ml_offset_y,
+ MLIB_BILINEAR,
+ MLIB_EDGE_SRC_EXTEND_INDEF);
+ }
+ else
+ {
+ medialib_get_interpolation (&ml_interp, interp_type,
+ scale_x, scale_y, 1.0);
+
+ if (ml_interp.interp_table != NULL)
+ {
+ mlib_ImageZoomTranslateTable (&img_dest,
+ &img_src,
+ scale_x,
+ scale_y,
+ ml_offset_x + ml_interp.tx,
+ ml_offset_y + ml_interp.ty,
+ ml_interp.interp_table,
+ MLIB_EDGE_SRC_EXTEND_INDEF);
+
+ mlib_ImageInterpTableDelete (ml_interp.interp_table);
+ }
+ else
+ {
+ /* Should not happen. */
+ mlib_filter ml_filter;
+
+ switch (interp_type)
+ {
+ case PIXOPS_INTERP_BILINEAR:
+ ml_filter = MLIB_BILINEAR;
+ break;
+
+ case PIXOPS_INTERP_TILES:
+ ml_filter = MLIB_BILINEAR;
+ break;
+
+ case PIXOPS_INTERP_HYPER:
+ ml_filter = MLIB_BICUBIC;
+ break;
+ }
+
+ mlib_ImageZoomTranslate (&img_dest,
+ &img_src,
+ scale_x,
+ scale_y,
+ ml_offset_x,
+ ml_offset_y,
+ ml_filter,
+ MLIB_EDGE_SRC_EXTEND_INDEF);
+ }
+ }
+ }
+
+ /* Deal with case where src_channels != dest_channels || src_has_alpha */
+ else if (interp_type == PIXOPS_INTERP_BILINEAR &&
+ scale_x > 1.0 && scale_y > 1.0)
+ {
+ mlib_ImageZoomTranslateBlend (&img_dest,
+ &img_src,
+ scale_x,
+ scale_y,
+ ml_offset_x,
+ ml_offset_y,
+ MLIB_BILINEAR,
+ MLIB_EDGE_SRC_EXTEND_INDEF,
+ MLIB_BLEND_GTK_SRC,
+ 1.0,
+ 1);
+ }
+ else
+ {
+ medialib_get_interpolation (&ml_interp, interp_type,
+ scale_x, scale_y, 1.0);
+
+ if (ml_interp.interp_table != NULL)
+ {
+ mlib_ImageZoomTranslateTableBlend (&img_dest,
+ &img_src,
+ scale_x,
+ scale_y,
+ ml_offset_x + ml_interp.tx,
+ ml_offset_y + ml_interp.ty,
+ ml_interp.interp_table,
+ MLIB_EDGE_SRC_EXTEND_INDEF,
+ MLIB_BLEND_GTK_SRC,
+ 1);
+ mlib_ImageInterpTableDelete (ml_interp.interp_table);
+ }
+ else
+ {
+ mlib_filter ml_filter;
+
+ switch (interp_type)
+ {
+ case PIXOPS_INTERP_BILINEAR:
+ ml_filter = MLIB_BILINEAR;
+ break;
+
+ case PIXOPS_INTERP_TILES:
+ ml_filter = MLIB_BILINEAR;
+ break;
+
+ case PIXOPS_INTERP_HYPER:
+ ml_filter = MLIB_BICUBIC;
+ break;
+ }
+
+ mlib_ImageZoomTranslate (&img_dest,
+ &img_src,
+ scale_x,
+ scale_y,
+ ml_offset_x,
+ ml_offset_y,
+ ml_filter,
+ MLIB_EDGE_SRC_EXTEND_INDEF);
+ }
+ }
+
+ if (tmp_buf != NULL)
+ g_free (tmp_buf);
+ }
+}
+#endif