Subversion Repositories Kolibri OS

/* -*- Mode: c; c-basic-offset: 4; tab-width: 8; indent-tabs-mode: t; -*- */

/*

 * Copyright © 2010, 2012 Soren Sandmann Pedersen

 * Copyright © 2010, 2012 Red Hat, Inc.

 *

 * Permission is hereby granted, free of charge, to any person obtaining a

 * copy of this software and associated documentation files (the "Software"),

 * to deal in the Software without restriction, including without limitation

 * the rights to use, copy, modify, merge, publish, distribute, sublicense,

 * and/or sell copies of the Software, and to permit persons to whom the

 * Software is furnished to do so, subject to the following conditions:

 *

 * The above copyright notice and this permission notice (including the next

 * paragraph) shall be included in all copies or substantial portions of the

 * Software.

 *

 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR

 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,

 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.  IN NO EVENT SHALL

 * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER

 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING

 * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER

 * DEALINGS IN THE SOFTWARE.

 *

 * Author: Soren Sandmann Pedersen (sandmann@cs.au.dk)

 */

#ifdef HAVE_CONFIG_H

#include 

#endif

#include 

#include 

#include 

#include "pixman-private.h"

/* Workaround for http://gcc.gnu.org/PR54965 */

/* GCC 4.6 has problems with force_inline, so just use normal inline instead */

#if defined(__GNUC__) && (__GNUC__ == 4) && (__GNUC_MINOR__ == 6)

#undef force_inline

#define force_inline __inline__

#endif

typedef float (* combine_channel_t) (float sa, float s, float da, float d);

static force_inline void

combine_inner (pixman_bool_t component,

	       float *dest, const float *src, const float *mask, int n_pixels,

	       combine_channel_t combine_a, combine_channel_t combine_c)

{

    int i;

    if (!mask)

    {

	for (i = 0; i < 4 * n_pixels; i += 4)

	{

	    float sa = src[i + 0];

	    float sr = src[i + 1];

	    float sg = src[i + 2];

	    float sb = src[i + 3];

	    float da = dest[i + 0];

	    float dr = dest[i + 1];

	    float dg = dest[i + 2];

	    float db = dest[i + 3];

	    dest[i + 0] = combine_a (sa, sa, da, da);

	    dest[i + 1] = combine_c (sa, sr, da, dr);

	    dest[i + 2] = combine_c (sa, sg, da, dg);

	    dest[i + 3] = combine_c (sa, sb, da, db);

	}

    }

    else

    {

	for (i = 0; i < 4 * n_pixels; i += 4)

	{

	    float sa, sr, sg, sb;

	    float ma, mr, mg, mb;

	    float da, dr, dg, db;

	    sa = src[i + 0];

	    sr = src[i + 1];

	    sg = src[i + 2];

	    sb = src[i + 3];

	    if (component)

	    {

		ma = mask[i + 0];

		mr = mask[i + 1];

		mg = mask[i + 2];

		mb = mask[i + 3];

		sr *= mr;

		sg *= mg;

		sb *= mb;

		ma *= sa;

		mr *= sa;

		mg *= sa;

		mb *= sa;

		sa = ma;

	    }

	    else

	    {

		ma = mask[i + 0];

		sa *= ma;

		sr *= ma;

		sg *= ma;

		sb *= ma;

		ma = mr = mg = mb = sa;

	    }

	    da = dest[i + 0];

	    dr = dest[i + 1];

	    dg = dest[i + 2];

	    db = dest[i + 3];

	    dest[i + 0] = combine_a (ma, sa, da, da);

	    dest[i + 1] = combine_c (mr, sr, da, dr);

	    dest[i + 2] = combine_c (mg, sg, da, dg);

	    dest[i + 3] = combine_c (mb, sb, da, db);

	}

    }

}

#define MAKE_COMBINER(name, component, combine_a, combine_c)		\

    static void								\

    combine_ ## name ## _float (pixman_implementation_t *imp,		\

				pixman_op_t              op,		\

				float                   *dest,		\

				const float             *src,		\

				const float             *mask,		\

				int		         n_pixels)	\

    {									\

	combine_inner (component, dest, src, mask, n_pixels,		\

		       combine_a, combine_c);				\

    }

#define MAKE_COMBINERS(name, combine_a, combine_c)			\

    MAKE_COMBINER(name ## _ca, TRUE, combine_a, combine_c)		\

    MAKE_COMBINER(name ## _u, FALSE, combine_a, combine_c)

/*

 * Porter/Duff operators

 */

typedef enum

{

    ZERO,

    ONE,

    SRC_ALPHA,

    DEST_ALPHA,

    INV_SA,

    INV_DA,

    SA_OVER_DA,

    DA_OVER_SA,

    INV_SA_OVER_DA,

    INV_DA_OVER_SA,

    ONE_MINUS_SA_OVER_DA,

    ONE_MINUS_DA_OVER_SA,

    ONE_MINUS_INV_DA_OVER_SA,

    ONE_MINUS_INV_SA_OVER_DA

} combine_factor_t;

#define CLAMP(f)					\

    (((f) < 0)? 0 : (((f) > 1.0) ? 1.0 : (f)))

static force_inline float

get_factor (combine_factor_t factor, float sa, float da)

{

    float f = -1;

    switch (factor)

    {

    case ZERO:

	f = 0.0f;

	break;

    case ONE:

	f = 1.0f;

	break;

    case SRC_ALPHA:

	f = sa;

	break;

    case DEST_ALPHA:

	f = da;

	break;

    case INV_SA:

	f = 1 - sa;

	break;

    case INV_DA:

	f = 1 - da;

	break;

    case SA_OVER_DA:

	if (FLOAT_IS_ZERO (da))

	    f = 1.0f;

	else

	    f = CLAMP (sa / da);

	break;

    case DA_OVER_SA:

	if (FLOAT_IS_ZERO (sa))

	    f = 1.0f;

	else

	    f = CLAMP (da / sa);

	break;

    case INV_SA_OVER_DA:

	if (FLOAT_IS_ZERO (da))

	    f = 1.0f;

	else

	    f = CLAMP ((1.0f - sa) / da);

	break;

    case INV_DA_OVER_SA:

	if (FLOAT_IS_ZERO (sa))

	    f = 1.0f;

	else

	    f = CLAMP ((1.0f - da) / sa);

	break;

    case ONE_MINUS_SA_OVER_DA:

	if (FLOAT_IS_ZERO (da))

	    f = 0.0f;

	else

	    f = CLAMP (1.0f - sa / da);

	break;

    case ONE_MINUS_DA_OVER_SA:

	if (FLOAT_IS_ZERO (sa))

	    f = 0.0f;

	else

	    f = CLAMP (1.0f - da / sa);

	break;

    case ONE_MINUS_INV_DA_OVER_SA:

	if (FLOAT_IS_ZERO (sa))

	    f = 0.0f;

	else

	    f = CLAMP (1.0f - (1.0f - da) / sa);

	break;

    case ONE_MINUS_INV_SA_OVER_DA:

	if (FLOAT_IS_ZERO (da))

	    f = 0.0f;

	else

	    f = CLAMP (1.0f - (1.0f - sa) / da);

	break;

    }

    return f;

}

#define MAKE_PD_COMBINERS(name, a, b)					\

    static float force_inline						\

    pd_combine_ ## name (float sa, float s, float da, float d)		\

    {									\

	const float fa = get_factor (a, sa, da);			\

	const float fb = get_factor (b, sa, da);			\

									\

	return MIN (1.0f, s * fa + d * fb);				\

    }									\

    									\

    MAKE_COMBINERS(name, pd_combine_ ## name, pd_combine_ ## name)

MAKE_PD_COMBINERS (clear,			ZERO,				ZERO)

MAKE_PD_COMBINERS (src,				ONE,				ZERO)

MAKE_PD_COMBINERS (dst,				ZERO,				ONE)

MAKE_PD_COMBINERS (over,			ONE,				INV_SA)

MAKE_PD_COMBINERS (over_reverse,		INV_DA,				ONE)

MAKE_PD_COMBINERS (in,				DEST_ALPHA,			ZERO)

MAKE_PD_COMBINERS (in_reverse,			ZERO,				SRC_ALPHA)

MAKE_PD_COMBINERS (out,				INV_DA,				ZERO)

MAKE_PD_COMBINERS (out_reverse,			ZERO,				INV_SA)

MAKE_PD_COMBINERS (atop,			DEST_ALPHA,			INV_SA)

MAKE_PD_COMBINERS (atop_reverse,		INV_DA,				SRC_ALPHA)

MAKE_PD_COMBINERS (xor,				INV_DA,				INV_SA)

MAKE_PD_COMBINERS (add,				ONE,				ONE)

MAKE_PD_COMBINERS (saturate,			INV_DA_OVER_SA,			ONE)

MAKE_PD_COMBINERS (disjoint_clear,		ZERO,				ZERO)

MAKE_PD_COMBINERS (disjoint_src,		ONE,				ZERO)

MAKE_PD_COMBINERS (disjoint_dst,		ZERO,				ONE)

MAKE_PD_COMBINERS (disjoint_over,		ONE,				INV_SA_OVER_DA)

MAKE_PD_COMBINERS (disjoint_over_reverse,	INV_DA_OVER_SA,			ONE)

MAKE_PD_COMBINERS (disjoint_in,			ONE_MINUS_INV_DA_OVER_SA,	ZERO)

MAKE_PD_COMBINERS (disjoint_in_reverse,		ZERO,				ONE_MINUS_INV_SA_OVER_DA)

MAKE_PD_COMBINERS (disjoint_out,		INV_DA_OVER_SA,			ZERO)

MAKE_PD_COMBINERS (disjoint_out_reverse,	ZERO,				INV_SA_OVER_DA)

MAKE_PD_COMBINERS (disjoint_atop,		ONE_MINUS_INV_DA_OVER_SA,	INV_SA_OVER_DA)

MAKE_PD_COMBINERS (disjoint_atop_reverse,	INV_DA_OVER_SA,			ONE_MINUS_INV_SA_OVER_DA)

MAKE_PD_COMBINERS (disjoint_xor,		INV_DA_OVER_SA,			INV_SA_OVER_DA)

MAKE_PD_COMBINERS (conjoint_clear,		ZERO,				ZERO)

MAKE_PD_COMBINERS (conjoint_src,		ONE,				ZERO)

MAKE_PD_COMBINERS (conjoint_dst,		ZERO,				ONE)

MAKE_PD_COMBINERS (conjoint_over,		ONE,				ONE_MINUS_SA_OVER_DA)

MAKE_PD_COMBINERS (conjoint_over_reverse,	ONE_MINUS_DA_OVER_SA,		ONE)

MAKE_PD_COMBINERS (conjoint_in,			DA_OVER_SA,			ZERO)

MAKE_PD_COMBINERS (conjoint_in_reverse,		ZERO,				SA_OVER_DA)

MAKE_PD_COMBINERS (conjoint_out,		ONE_MINUS_DA_OVER_SA,		ZERO)

MAKE_PD_COMBINERS (conjoint_out_reverse,	ZERO,				ONE_MINUS_SA_OVER_DA)

MAKE_PD_COMBINERS (conjoint_atop,		DA_OVER_SA,			ONE_MINUS_SA_OVER_DA)

MAKE_PD_COMBINERS (conjoint_atop_reverse,	ONE_MINUS_DA_OVER_SA,		SA_OVER_DA)

MAKE_PD_COMBINERS (conjoint_xor,		ONE_MINUS_DA_OVER_SA,		ONE_MINUS_SA_OVER_DA)

/*

 * PDF blend modes:

 *

 * The following blend modes have been taken from the PDF ISO 32000

 * specification, which at this point in time is available from

 * http://www.adobe.com/devnet/acrobat/pdfs/PDF32000_2008.pdf

 * The relevant chapters are 11.3.5 and 11.3.6.

 * The formula for computing the final pixel color given in 11.3.6 is:

 * αr × Cr = (1 – αs) × αb × Cb + (1 – αb) × αs × Cs + αb × αs × B(Cb, Cs)

 * with B() being the blend function.

 * Note that OVER is a special case of this operation, using B(Cb, Cs) = Cs

 *

 * These blend modes should match the SVG filter draft specification, as

 * it has been designed to mirror ISO 32000. Note that at the current point

 * no released draft exists that shows this, as the formulas have not been

 * updated yet after the release of ISO 32000.

 *

 * The default implementation here uses the PDF_SEPARABLE_BLEND_MODE and

 * PDF_NON_SEPARABLE_BLEND_MODE macros, which take the blend function as an

 * argument. Note that this implementation operates on premultiplied colors,

 * while the PDF specification does not. Therefore the code uses the formula

 * ar.Cra = (1 – as) . Dca + (1 – ad) . Sca + B(Dca, ad, Sca, as)

 */

#define MAKE_SEPARABLE_PDF_COMBINERS(name)				\

    static force_inline float						\

    combine_ ## name ## _a (float sa, float s, float da, float d)	\

    {									\

	return da + sa - da * sa;					\

    }									\

    									\

    static force_inline float						\

    combine_ ## name ## _c (float sa, float s, float da, float d)	\

    {									\

	float f = (1 - sa) * d + (1 - da) * s;				\

									\

	return f + blend_ ## name (sa, s, da, d);			\

    }									\

    									\

    MAKE_COMBINERS (name, combine_ ## name ## _a, combine_ ## name ## _c)

static force_inline float

blend_multiply (float sa, float s, float da, float d)

{

    return d * s;

}

static force_inline float

blend_screen (float sa, float s, float da, float d)

{

    return d * sa + s * da - s * d;

}

static force_inline float

blend_overlay (float sa, float s, float da, float d)

{

    if (2 * d < da)

	return 2 * s * d;

    else

	return sa * da - 2 * (da - d) * (sa - s);

}

static force_inline float

blend_darken (float sa, float s, float da, float d)

{

    s = s * da;

    d = d * sa;

    if (s > d)

	return d;

    else

	return s;

}

static force_inline float

blend_lighten (float sa, float s, float da, float d)

{

    s = s * da;

    d = d * sa;

    if (s > d)

	return s;

    else

	return d;

}

static force_inline float

blend_color_dodge (float sa, float s, float da, float d)

{

    if (FLOAT_IS_ZERO (d))

	return 0.0f;

    else if (d * sa >= sa * da - s * da)

	return sa * da;

    else if (FLOAT_IS_ZERO (sa - s))

	return sa * da;

    else

	return sa * sa * d / (sa - s);

}

static force_inline float

blend_color_burn (float sa, float s, float da, float d)

{

    if (d >= da)

	return sa * da;

    else if (sa * (da - d) >= s * da)

	return 0.0f;

    else if (FLOAT_IS_ZERO (s))

	return 0.0f;

    else

	return sa * (da - sa * (da - d) / s);

}

static force_inline float

blend_hard_light (float sa, float s, float da, float d)

{

    if (2 * s < sa)

	return 2 * s * d;

    else

	return sa * da - 2 * (da - d) * (sa - s);

}

static force_inline float

blend_soft_light (float sa, float s, float da, float d)

{

    if (2 * s < sa)

    {

	if (FLOAT_IS_ZERO (da))

	    return d * sa;

	else

	    return d * sa - d * (da - d) * (sa - 2 * s) / da;

    }

    else

    {

	if (FLOAT_IS_ZERO (da))

	{

	    return 0.0f;

	}

	else

	{

	    if (4 * d <= da)

		return d * sa + (2 * s - sa) * d * ((16 * d / da - 12) * d / da + 3);

	    else

		return d * sa + (sqrtf (d * da) - d) * (2 * s - sa);

	}

    }

}

static force_inline float

blend_difference (float sa, float s, float da, float d)

{

    float dsa = d * sa;

    float sda = s * da;

    if (sda < dsa)

	return dsa - sda;

    else

	return sda - dsa;

}

static force_inline float

blend_exclusion (float sa, float s, float da, float d)

{

    return s * da + d * sa - 2 * d * s;

}

MAKE_SEPARABLE_PDF_COMBINERS (multiply)

MAKE_SEPARABLE_PDF_COMBINERS (screen)

MAKE_SEPARABLE_PDF_COMBINERS (overlay)

MAKE_SEPARABLE_PDF_COMBINERS (darken)

MAKE_SEPARABLE_PDF_COMBINERS (lighten)

MAKE_SEPARABLE_PDF_COMBINERS (color_dodge)

MAKE_SEPARABLE_PDF_COMBINERS (color_burn)

MAKE_SEPARABLE_PDF_COMBINERS (hard_light)

MAKE_SEPARABLE_PDF_COMBINERS (soft_light)

MAKE_SEPARABLE_PDF_COMBINERS (difference)

MAKE_SEPARABLE_PDF_COMBINERS (exclusion)

/*

 * PDF nonseperable blend modes.

 *

 * These are implemented using the following functions to operate in Hsl

 * space, with Cmax, Cmid, Cmin referring to the max, mid and min value

 * of the red, green and blue components.

 *

 * LUM (C) = 0.3 × Cred + 0.59 × Cgreen + 0.11 × Cblue

 *

 * clip_color (C):

 *   l = LUM (C)

 *   min = Cmin

 *   max = Cmax

 *   if n < 0.0

 *     C = l + (((C – l) × l) ⁄     (l – min))

 *   if x > 1.0

 *     C = l + (((C – l) × (1 – l)) (max – l))

 *   return C

 *

 * set_lum (C, l):

 *   d = l – LUM (C)

 *   C += d

 *   return clip_color (C)

 *

 * SAT (C) = CH_MAX (C) - CH_MIN (C)

 *

 * set_sat (C, s):

 *  if Cmax > Cmin

 *    Cmid = ( ( ( Cmid – Cmin ) × s ) ⁄ ( Cmax – Cmin ) )

 *    Cmax = s

 *  else

 *    Cmid = Cmax = 0.0

 *  Cmin = 0.0

 *  return C

 */

/* For premultiplied colors, we need to know what happens when C is

 * multiplied by a real number. LUM and SAT are linear:

 *

 *    LUM (r × C) = r × LUM (C)		SAT (r × C) = r × SAT (C)

 *

 * If we extend clip_color with an extra argument a and change

 *

 *        if x >= 1.0

 *

 * into

 *

 *        if x >= a

 *

 * then clip_color is also linear:

 *

 *     r * clip_color (C, a) = clip_color (r_c, ra);

 *

 * for positive r.

 *

 * Similarly, we can extend set_lum with an extra argument that is just passed

 * on to clip_color:

 *

 *     r × set_lum ( C, l, a)

 *

 *   = r × clip_color ( C + l - LUM (C), a)

 *

 *   = clip_color ( r * C + r × l - LUM (r × C), r * a)

 *

 *   = set_lum ( r * C, r * l, r * a)

 *

 * Finally, set_sat:

 *

 *     r * set_sat (C, s) = set_sat (x * C, r * s)

 *

 * The above holds for all non-zero x because they x'es in the fraction for

 * C_mid cancel out. Specifically, it holds for x = r:

 *

 *     r * set_sat (C, s) = set_sat (r_c, rs)

 *

 *

 *

 *

 * So, for the non-separable PDF blend modes, we have (using s, d for

 * non-premultiplied colors, and S, D for premultiplied:

 *

 *   Color:

 *

 *     a_s * a_d * B(s, d)

 *   = a_s * a_d * set_lum (S/a_s, LUM (D/a_d), 1)

 *   = set_lum (S * a_d, a_s * LUM (D), a_s * a_d)

 *

 *

 *   Luminosity:

 *

 *     a_s * a_d * B(s, d)

 *   = a_s * a_d * set_lum (D/a_d, LUM(S/a_s), 1)

 *   = set_lum (a_s * D, a_d * LUM(S), a_s * a_d)

 *

 *

 *   Saturation:

 *

 *     a_s * a_d * B(s, d)

 *   = a_s * a_d * set_lum (set_sat (D/a_d, SAT (S/a_s)), LUM (D/a_d), 1)

 *   = set_lum (a_s * a_d * set_sat (D/a_d, SAT (S/a_s)),

 *                                        a_s * LUM (D), a_s * a_d)

 *   = set_lum (set_sat (a_s * D, a_d * SAT (S), a_s * LUM (D), a_s * a_d))

 *

 *   Hue:

 *

 *     a_s * a_d * B(s, d)

 *   = a_s * a_d * set_lum (set_sat (S/a_s, SAT (D/a_d)), LUM (D/a_d), 1)

 *   = set_lum (set_sat (a_d * S, a_s * SAT (D)), a_s * LUM (D), a_s * a_d)

 *

 */

typedef struct

{

    float	r;

    float	g;

    float	b;

} rgb_t;

static force_inline float

minf (float a, float b)

{

    return a < b? a : b;

}

static force_inline float

maxf (float a, float b)

{

    return a > b? a : b;

}

static force_inline float

channel_min (const rgb_t *c)

{

    return minf (minf (c->r, c->g), c->b);

}

static force_inline float

channel_max (const rgb_t *c)

{

    return maxf (maxf (c->r, c->g), c->b);

}

static force_inline float

get_lum (const rgb_t *c)

{

    return c->r * 0.3f + c->g * 0.59f + c->b * 0.11f;

}

static force_inline float

get_sat (const rgb_t *c)

{

    return channel_max (c) - channel_min (c);

}

static void

clip_color (rgb_t *color, float a)

{

    float l = get_lum (color);

    float n = channel_min (color);

    float x = channel_max (color);

    float t;

    if (n < 0.0f)

    {

	t = l - n;

	if (FLOAT_IS_ZERO (t))

	{

	    color->r = 0.0f;

	    color->g = 0.0f;

	    color->b = 0.0f;

	}

	else

	{

	    color->r = l + (((color->r - l) * l) / t);

	    color->g = l + (((color->g - l) * l) / t);

	    color->b = l + (((color->b - l) * l) / t);

	}

    }

    if (x > a)

    {

	t = x - l;

	if (FLOAT_IS_ZERO (t))

	{

	    color->r = a;

	    color->g = a;

	    color->b = a;

	}

	else

	{

	    color->r = l + (((color->r - l) * (a - l) / t));

	    color->g = l + (((color->g - l) * (a - l) / t));

	    color->b = l + (((color->b - l) * (a - l) / t));

	}

    }

}

static void

set_lum (rgb_t *color, float sa, float l)

{

    float d = l - get_lum (color);

    color->r = color->r + d;

    color->g = color->g + d;

    color->b = color->b + d;

    clip_color (color, sa);

}

static void

set_sat (rgb_t *src, float sat)

{

    float *max, *mid, *min;

    float t;

    if (src->r > src->g)

    {

	if (src->r > src->b)

	{

	    max = &(src->r);

	    if (src->g > src->b)

	    {

		mid = &(src->g);

		min = &(src->b);

	    }

	    else

	    {

		mid = &(src->b);

		min = &(src->g);

	    }

	}

	else

	{

	    max = &(src->b);

	    mid = &(src->r);

	    min = &(src->g);

	}

    }

    else

    {

	if (src->r > src->b)

	{

	    max = &(src->g);

	    mid = &(src->r);

	    min = &(src->b);

	}

	else

	{

	    min = &(src->r);

	    if (src->g > src->b)

	    {

		max = &(src->g);

		mid = &(src->b);

	    }

	    else

	    {

		max = &(src->b);

		mid = &(src->g);

	    }

	}

    }

    t = *max - *min;

    if (FLOAT_IS_ZERO (t))

    {

	*mid = *max = 0.0f;

    }

    else

    {

	*mid = ((*mid - *min) * sat) / t;

	*max = sat;

    }

    *min = 0.0f;

}

/*

 * Hue:

 * B(Cb, Cs) = set_lum (set_sat (Cs, SAT (Cb)), LUM (Cb))

 */

static force_inline void

blend_hsl_hue (rgb_t *res,

	       const rgb_t *dest, float da,

	       const rgb_t *src, float sa)

{

    res->r = src->r * da;

    res->g = src->g * da;

    res->b = src->b * da;

    set_sat (res, get_sat (dest) * sa);

    set_lum (res, sa * da, get_lum (dest) * sa);

}

/*

 * Saturation:

 * B(Cb, Cs) = set_lum (set_sat (Cb, SAT (Cs)), LUM (Cb))

 */

static force_inline void

blend_hsl_saturation (rgb_t *res,

		      const rgb_t *dest, float da,

		      const rgb_t *src, float sa)

{

    res->r = dest->r * sa;

    res->g = dest->g * sa;

    res->b = dest->b * sa;

    set_sat (res, get_sat (src) * da);

    set_lum (res, sa * da, get_lum (dest) * sa);

}

/*

 * Color:

 * B(Cb, Cs) = set_lum (Cs, LUM (Cb))

 */

static force_inline void

blend_hsl_color (rgb_t *res,

		 const rgb_t *dest, float da,

		 const rgb_t *src, float sa)

{

    res->r = src->r * da;

    res->g = src->g * da;

    res->b = src->b * da;

    set_lum (res, sa * da, get_lum (dest) * sa);

}

/*

 * Luminosity:

 * B(Cb, Cs) = set_lum (Cb, LUM (Cs))

 */

static force_inline void

blend_hsl_luminosity (rgb_t *res,

		      const rgb_t *dest, float da,

		      const rgb_t *src, float sa)

{

    res->r = dest->r * sa;

    res->g = dest->g * sa;

    res->b = dest->b * sa;

    set_lum (res, sa * da, get_lum (src) * da);

}

#define MAKE_NON_SEPARABLE_PDF_COMBINERS(name)				\

    static void								\

    combine_ ## name ## _u_float (pixman_implementation_t *imp,		\

				  pixman_op_t              op,		\

				  float                   *dest,	\

				  const float             *src,		\

				  const float             *mask,	\

				  int		           n_pixels)	\

    {									\

    	int i;								\

									\

	for (i = 0; i < 4 * n_pixels; i += 4)				\

	{								\

	    float sa, da;						\

	    rgb_t sc, dc, rc;						\

									\

	    sa = src[i + 0];						\

	    sc.r = src[i + 1];						\

	    sc.g = src[i + 2];						\

	    sc.b = src[i + 3];						\

									\

	    da = dest[i + 0];						\

	    dc.r = dest[i + 1];						\

	    dc.g = dest[i + 2];						\

	    dc.b = dest[i + 3];						\

									\

	    if (mask)							\

	    {								\

		float ma = mask[i + 0];					\

									\

		/* Component alpha is not supported for HSL modes */	\

		sa *= ma;						\

		sc.r *= ma;						\

		sc.g *= ma;						\

		sc.g *= ma;						\

	    }								\

									\

	    blend_ ## name (&rc, &dc, da, &sc, sa);			\

									\

	    dest[i + 0] = sa + da - sa * da;				\

	    dest[i + 1] = (1 - sa) * dc.r + (1 - da) * sc.r + rc.r;	\

	    dest[i + 2] = (1 - sa) * dc.g + (1 - da) * sc.g + rc.g;	\

	    dest[i + 3] = (1 - sa) * dc.b + (1 - da) * sc.b + rc.b;	\

	}								\

    }

MAKE_NON_SEPARABLE_PDF_COMBINERS(hsl_hue)

MAKE_NON_SEPARABLE_PDF_COMBINERS(hsl_saturation)

MAKE_NON_SEPARABLE_PDF_COMBINERS(hsl_color)

MAKE_NON_SEPARABLE_PDF_COMBINERS(hsl_luminosity)

void

_pixman_setup_combiner_functions_float (pixman_implementation_t *imp)

{

    /* Unified alpha */

    imp->combine_float[PIXMAN_OP_CLEAR] = combine_clear_u_float;

    imp->combine_float[PIXMAN_OP_SRC] = combine_src_u_float;

    imp->combine_float[PIXMAN_OP_DST] = combine_dst_u_float;

    imp->combine_float[PIXMAN_OP_OVER] = combine_over_u_float;

    imp->combine_float[PIXMAN_OP_OVER_REVERSE] = combine_over_reverse_u_float;

    imp->combine_float[PIXMAN_OP_IN] = combine_in_u_float;

    imp->combine_float[PIXMAN_OP_IN_REVERSE] = combine_in_reverse_u_float;

    imp->combine_float[PIXMAN_OP_OUT] = combine_out_u_float;

    imp->combine_float[PIXMAN_OP_OUT_REVERSE] = combine_out_reverse_u_float;

    imp->combine_float[PIXMAN_OP_ATOP] = combine_atop_u_float;

    imp->combine_float[PIXMAN_OP_ATOP_REVERSE] = combine_atop_reverse_u_float;

    imp->combine_float[PIXMAN_OP_XOR] = combine_xor_u_float;

    imp->combine_float[PIXMAN_OP_ADD] = combine_add_u_float;

    imp->combine_float[PIXMAN_OP_SATURATE] = combine_saturate_u_float;

    /* Disjoint, unified */

    imp->combine_float[PIXMAN_OP_DISJOINT_CLEAR] = combine_disjoint_clear_u_float;

    imp->combine_float[PIXMAN_OP_DISJOINT_SRC] = combine_disjoint_src_u_float;

    imp->combine_float[PIXMAN_OP_DISJOINT_DST] = combine_disjoint_dst_u_float;

    imp->combine_float[PIXMAN_OP_DISJOINT_OVER] = combine_disjoint_over_u_float;

    imp->combine_float[PIXMAN_OP_DISJOINT_OVER_REVERSE] = combine_disjoint_over_reverse_u_float;

    imp->combine_float[PIXMAN_OP_DISJOINT_IN] = combine_disjoint_in_u_float;

    imp->combine_float[PIXMAN_OP_DISJOINT_IN_REVERSE] = combine_disjoint_in_reverse_u_float;

    imp->combine_float[PIXMAN_OP_DISJOINT_OUT] = combine_disjoint_out_u_float;

    imp->combine_float[PIXMAN_OP_DISJOINT_OUT_REVERSE] = combine_disjoint_out_reverse_u_float;

    imp->combine_float[PIXMAN_OP_DISJOINT_ATOP] = combine_disjoint_atop_u_float;

    imp->combine_float[PIXMAN_OP_DISJOINT_ATOP_REVERSE] = combine_disjoint_atop_reverse_u_float;

    imp->combine_float[PIXMAN_OP_DISJOINT_XOR] = combine_disjoint_xor_u_float;

    /* Conjoint, unified */

    imp->combine_float[PIXMAN_OP_CONJOINT_CLEAR] = combine_conjoint_clear_u_float;

    imp->combine_float[PIXMAN_OP_CONJOINT_SRC] = combine_conjoint_src_u_float;

    imp->combine_float[PIXMAN_OP_CONJOINT_DST] = combine_conjoint_dst_u_float;

    imp->combine_float[PIXMAN_OP_CONJOINT_OVER] = combine_conjoint_over_u_float;

    imp->combine_float[PIXMAN_OP_CONJOINT_OVER_REVERSE] = combine_conjoint_over_reverse_u_float;

    imp->combine_float[PIXMAN_OP_CONJOINT_IN] = combine_conjoint_in_u_float;

    imp->combine_float[PIXMAN_OP_CONJOINT_IN_REVERSE] = combine_conjoint_in_reverse_u_float;

    imp->combine_float[PIXMAN_OP_CONJOINT_OUT] = combine_conjoint_out_u_float;

    imp->combine_float[PIXMAN_OP_CONJOINT_OUT_REVERSE] = combine_conjoint_out_reverse_u_float;

    imp->combine_float[PIXMAN_OP_CONJOINT_ATOP] = combine_conjoint_atop_u_float;

    imp->combine_float[PIXMAN_OP_CONJOINT_ATOP_REVERSE] = combine_conjoint_atop_reverse_u_float;

    imp->combine_float[PIXMAN_OP_CONJOINT_XOR] = combine_conjoint_xor_u_float;

    /* PDF operators, unified */

    imp->combine_float[PIXMAN_OP_MULTIPLY] = combine_multiply_u_float;

    imp->combine_float[PIXMAN_OP_SCREEN] = combine_screen_u_float;

    imp->combine_float[PIXMAN_OP_OVERLAY] = combine_overlay_u_float;

    imp->combine_float[PIXMAN_OP_DARKEN] = combine_darken_u_float;

    imp->combine_float[PIXMAN_OP_LIGHTEN] = combine_lighten_u_float;

    imp->combine_float[PIXMAN_OP_COLOR_DODGE] = combine_color_dodge_u_float;

    imp->combine_float[PIXMAN_OP_COLOR_BURN] = combine_color_burn_u_float;

    imp->combine_float[PIXMAN_OP_HARD_LIGHT] = combine_hard_light_u_float;

    imp->combine_float[PIXMAN_OP_SOFT_LIGHT] = combine_soft_light_u_float;

    imp->combine_float[PIXMAN_OP_DIFFERENCE] = combine_difference_u_float;

    imp->combine_float[PIXMAN_OP_EXCLUSION] = combine_exclusion_u_float;

    imp->combine_float[PIXMAN_OP_HSL_HUE] = combine_hsl_hue_u_float;

    imp->combine_float[PIXMAN_OP_HSL_SATURATION] = combine_hsl_saturation_u_float;

    imp->combine_float[PIXMAN_OP_HSL_COLOR] = combine_hsl_color_u_float;

    imp->combine_float[PIXMAN_OP_HSL_LUMINOSITY] = combine_hsl_luminosity_u_float;

    /* Component alpha combiners */

    imp->combine_float_ca[PIXMAN_OP_CLEAR] = combine_clear_ca_float;

    imp->combine_float_ca[PIXMAN_OP_SRC] = combine_src_ca_float;

    imp->combine_float_ca[PIXMAN_OP_DST] = combine_dst_ca_float;

    imp->combine_float_ca[PIXMAN_OP_OVER] = combine_over_ca_float;

    imp->combine_float_ca[PIXMAN_OP_OVER_REVERSE] = combine_over_reverse_ca_float;

    imp->combine_float_ca[PIXMAN_OP_IN] = combine_in_ca_float;

    imp->combine_float_ca[PIXMAN_OP_IN_REVERSE] = combine_in_reverse_ca_float;

    imp->combine_float_ca[PIXMAN_OP_OUT] = combine_out_ca_float;

    imp->combine_float_ca[PIXMAN_OP_OUT_REVERSE] = combine_out_reverse_ca_float;

    imp->combine_float_ca[PIXMAN_OP_ATOP] = combine_atop_ca_float;

    imp->combine_float_ca[PIXMAN_OP_ATOP_REVERSE] = combine_atop_reverse_ca_float;

    imp->combine_float_ca[PIXMAN_OP_XOR] = combine_xor_ca_float;

    imp->combine_float_ca[PIXMAN_OP_ADD] = combine_add_ca_float;

    imp->combine_float_ca[PIXMAN_OP_SATURATE] = combine_saturate_ca_float;

    /* Disjoint CA */

    imp->combine_float_ca[PIXMAN_OP_DISJOINT_CLEAR] = combine_disjoint_clear_ca_float;

    imp->combine_float_ca[PIXMAN_OP_DISJOINT_SRC] = combine_disjoint_src_ca_float;

    imp->combine_float_ca[PIXMAN_OP_DISJOINT_DST] = combine_disjoint_dst_ca_float;

    imp->combine_float_ca[PIXMAN_OP_DISJOINT_OVER] = combine_disjoint_over_ca_float;

    imp->combine_float_ca[PIXMAN_OP_DISJOINT_OVER_REVERSE] = combine_disjoint_over_reverse_ca_float;

    imp->combine_float_ca[PIXMAN_OP_DISJOINT_IN] = combine_disjoint_in_ca_float;

    imp->combine_float_ca[PIXMAN_OP_DISJOINT_IN_REVERSE] = combine_disjoint_in_reverse_ca_float;

    imp->combine_float_ca[PIXMAN_OP_DISJOINT_OUT] = combine_disjoint_out_ca_float;

    imp->combine_float_ca[PIXMAN_OP_DISJOINT_OUT_REVERSE] = combine_disjoint_out_reverse_ca_float;

    imp->combine_float_ca[PIXMAN_OP_DISJOINT_ATOP] = combine_disjoint_atop_ca_float;

    imp->combine_float_ca[PIXMAN_OP_DISJOINT_ATOP_REVERSE] = combine_disjoint_atop_reverse_ca_float;

    imp->combine_float_ca[PIXMAN_OP_DISJOINT_XOR] = combine_disjoint_xor_ca_float;

    /* Conjoint CA */

    imp->combine_float_ca[PIXMAN_OP_CONJOINT_CLEAR] = combine_conjoint_clear_ca_float;

    imp->combine_float_ca[PIXMAN_OP_CONJOINT_SRC] = combine_conjoint_src_ca_float;

    imp->combine_float_ca[PIXMAN_OP_CONJOINT_DST] = combine_conjoint_dst_ca_float;

    imp->combine_float_ca[PIXMAN_OP_CONJOINT_OVER] = combine_conjoint_over_ca_float;

    imp->combine_float_ca[PIXMAN_OP_CONJOINT_OVER_REVERSE] = combine_conjoint_over_reverse_ca_float;

    imp->combine_float_ca[PIXMAN_OP_CONJOINT_IN] = combine_conjoint_in_ca_float;

    imp->combine_float_ca[PIXMAN_OP_CONJOINT_IN_REVERSE] = combine_conjoint_in_reverse_ca_float;

    imp->combine_float_ca[PIXMAN_OP_CONJOINT_OUT] = combine_conjoint_out_ca_float;

    imp->combine_float_ca[PIXMAN_OP_CONJOINT_OUT_REVERSE] = combine_conjoint_out_reverse_ca_float;

    imp->combine_float_ca[PIXMAN_OP_CONJOINT_ATOP] = combine_conjoint_atop_ca_float;

    imp->combine_float_ca[PIXMAN_OP_CONJOINT_ATOP_REVERSE] = combine_conjoint_atop_reverse_ca_float;

    imp->combine_float_ca[PIXMAN_OP_CONJOINT_XOR] = combine_conjoint_xor_ca_float;

    /* PDF operators CA */

    imp->combine_float_ca[PIXMAN_OP_MULTIPLY] = combine_multiply_ca_float;

    imp->combine_float_ca[PIXMAN_OP_SCREEN] = combine_screen_ca_float;

    imp->combine_float_ca[PIXMAN_OP_OVERLAY] = combine_overlay_ca_float;

    imp->combine_float_ca[PIXMAN_OP_DARKEN] = combine_darken_ca_float;

    imp->combine_float_ca[PIXMAN_OP_LIGHTEN] = combine_lighten_ca_float;

    imp->combine_float_ca[PIXMAN_OP_COLOR_DODGE] = combine_color_dodge_ca_float;

    imp->combine_float_ca[PIXMAN_OP_COLOR_BURN] = combine_color_burn_ca_float;

    imp->combine_float_ca[PIXMAN_OP_HARD_LIGHT] = combine_hard_light_ca_float;

    imp->combine_float_ca[PIXMAN_OP_SOFT_LIGHT] = combine_soft_light_ca_float;

    imp->combine_float_ca[PIXMAN_OP_DIFFERENCE] = combine_difference_ca_float;

    imp->combine_float_ca[PIXMAN_OP_EXCLUSION] = combine_exclusion_ca_float;

    /* It is not clear that these make sense, so make them noops for now */

    imp->combine_float_ca[PIXMAN_OP_HSL_HUE] = combine_dst_u_float;

    imp->combine_float_ca[PIXMAN_OP_HSL_SATURATION] = combine_dst_u_float;

    imp->combine_float_ca[PIXMAN_OP_HSL_COLOR] = combine_dst_u_float;

    imp->combine_float_ca[PIXMAN_OP_HSL_LUMINOSITY] = combine_dst_u_float;

}