Logo Search packages:      
Sourcecode: inkscape version File versions  Download package

nr-compose.cpp

#define __NR_COMPOSE_C__

/*
 * Pixel buffer rendering library
 *
 * Authors:
 *   Lauris Kaplinski <lauris@kaplinski.com>
 *
 * This code is in public domain
 */

#ifdef HAVE_CONFIG_H
# include "config.h"
#endif

#include <string.h>
#include "nr-pixops.h"

#ifdef WITH_MMX
/* fixme: */
#ifdef __cplusplus
extern "C" {
#endif /* __cplusplus */
int nr_have_mmx (void);
void nr_mmx_R8G8B8A8_P_EMPTY_A8_RGBAP (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, unsigned char *c);
void nr_mmx_R8G8B8A8_P_R8G8B8A8_P_A8_RGBAP (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, unsigned char *c);
void nr_mmx_R8G8B8_R8G8B8_R8G8B8A8_P (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, unsigned int alpha);
#define NR_PIXOPS_MMX nr_have_mmx ()
#ifdef __cplusplus
}
#endif /* __cplusplus */
#endif

// Naming: nr_RESULT_BACKGROUND_FOREGROUND_extra

void
nr_R8G8B8A8_N_EMPTY_R8G8B8A8_N (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, unsigned int alpha)
{
      unsigned int r, c;

      for (r = h; r > 0; r--) {
            if (alpha == 0) {
                  memset(px, 0x0, 4 * w);
            } else if (alpha == 255) {
                  memcpy(px, spx, 4 * w);
            } else {
                  unsigned char *d = px;
                  const unsigned char *s = spx;
                  for (c = w; c > 0; c--) {
                        *d++ = *s++;
                        *d++ = *s++;
                        *d++ = *s++;
                        *d++ = NR_PREMUL_111(*s, alpha);
                        s++;
                  }
            }
            px += rs;
            spx += srs;
      }
}

void
nr_R8G8B8A8_N_EMPTY_R8G8B8A8_P (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, unsigned int alpha)
{
      unsigned int r, c;

      for (r = h; r > 0; r--) {
            if (alpha == 0) {
                  memset(px, 0x0, 4 * w);
            } else {
                  unsigned char *d = px;
                  const unsigned char *s = spx;
                  for (c = w; c > 0; c--) {
                        if (s[3] == 0) {
                              d[3] = 0;
                        } else if (s[3] == 255) {
                              memcpy(d, s, 4);
                        } else {
                              d[0] = NR_DEMUL_111(s[0], s[3]);
                              d[1] = NR_DEMUL_111(s[1], s[3]);
                              d[2] = NR_DEMUL_111(s[2], s[3]);
                              d[3] = NR_PREMUL_111(s[3], alpha);
                        }
                        d += 4;
                        s += 4;
                  }
            }
            px += rs;
            spx += srs;
      }
}

void
nr_R8G8B8A8_P_EMPTY_R8G8B8A8_N (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, unsigned int alpha)
{
      unsigned int r, c;

      for (r = h; r > 0; r--) {
            unsigned char *d = px;
            const unsigned char *s = spx;
            if (alpha == 0) {
                  memset(px, 0x0, 4 * w);
            } else if (alpha == 255) {
                  for (c = w; c > 0; c--) {
                        d[0] = NR_PREMUL_111(s[0], s[3]);
                        d[1] = NR_PREMUL_111(s[1], s[3]);
                        d[2] = NR_PREMUL_111(s[2], s[3]);
                        d[3] = s[3];
                        d += 4;
                        s += 4;
                  }
            } else {
                  for (c = w; c > 0; c--) {
                        if (s[3] == 0) {
                              memset(d, 0, 4);
                        } else {
                              unsigned int a;
                              a = NR_PREMUL_112(s[3], alpha);
                              d[0] = NR_PREMUL_121(s[0], a);
                              d[1] = NR_PREMUL_121(s[1], a);
                              d[2] = NR_PREMUL_121(s[2], a);
                              d[3] = NR_NORMALIZE_21(a);
                        }
                        d += 4;
                        s += 4;
                  }
            }
            px += rs;
            spx += srs;
      }
}

void
nr_R8G8B8A8_P_EMPTY_R8G8B8A8_P (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, unsigned int alpha)
{
      unsigned int r, c;

      for (r = h; r > 0; r--) {
            if (alpha == 0) {
                  memset(px, 0x0, 4 * w);
            } else if (alpha == 255) {
                  memcpy(px, spx, 4 * w);
            } else {
                  unsigned char *d = px;
                  const unsigned char *s = spx;
                  for (c = w; c > 0; c--) {
                        d[0] = NR_PREMUL_111(s[0], alpha);
                        d[1] = NR_PREMUL_111(s[1], alpha);
                        d[2] = NR_PREMUL_111(s[2], alpha);
                        d[3] = NR_PREMUL_111(s[3], alpha);
                        d += 4;
                        s += 4;
                  }
            }
            px += rs;
            spx += srs;
      }
}

void
nr_R8G8B8A8_N_R8G8B8A8_N_R8G8B8A8_N (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, unsigned int alpha)
{
      unsigned int r, c;

      if (alpha == 0) {
            /* NOP */
      } else if (alpha == 255) {
            for (r = h; r > 0; r--) {
                  unsigned char *d = px;
                  const unsigned char *s = spx;
                  for (c = w; c > 0; c--) {
                        if (s[3] == 0) {
                              /* Transparent FG, NOP */
                        } else if ((s[3] == 255) || (d[3] == 0)) {
                              /* Full coverage, COPY */
                              memcpy(d, s, 4);
                        } else {
                              /* Full composition */
                              unsigned int ca;
                              ca = NR_COMPOSEA_112(s[3], d[3]);
                              d[0] = NR_COMPOSENNN_111121(s[0], s[3], d[0], d[3], ca);
                              d[1] = NR_COMPOSENNN_111121(s[1], s[3], d[1], d[3], ca);
                              d[2] = NR_COMPOSENNN_111121(s[2], s[3], d[2], d[3], ca);
                              d[3] = NR_NORMALIZE_21(ca);
                        }
                        d += 4;
                        s += 4;
                  }
                  px += rs;
                  spx += srs;
            }
      } else {
            for (r = h; r > 0; r--) {
                  unsigned char *d = px;
                  const unsigned char *s = spx;
                  for (c = w; c > 0; c--) {
                        unsigned int a;
                        a = NR_PREMUL_112(s[3], alpha);
                        if (a == 0) {
                              /* Transparent FG, NOP */
                        } else if ((a == 255*255) || (d[3] == 0)) {
                              /* Full coverage, COPY */
                              d[0] = s[0];
                              d[1] = s[1];
                              d[2] = s[2];
                              d[3] = NR_NORMALIZE_21(a);
                        } else {
                              /* Full composition */
                              unsigned int ca;
                              ca = NR_COMPOSEA_213(a, d[3]);
                              d[0] = NR_COMPOSENNN_121131(s[0], a, d[0], d[3], ca);
                              d[1] = NR_COMPOSENNN_121131(s[1], a, d[1], d[3], ca);
                              d[2] = NR_COMPOSENNN_121131(s[2], a, d[2], d[3], ca);
                              d[3] = NR_NORMALIZE_31(ca);
                        }
                        d += 4;
                        s += 4;
                  }
                  px += rs;
                  spx += srs;
            }
      }
}

void
nr_R8G8B8A8_N_R8G8B8A8_N_R8G8B8A8_P (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, unsigned int alpha)
{
      unsigned int r, c;

      if (alpha == 0) {
            /* NOP */
      } else if (alpha == 255) {
            for (r = h; r > 0; r--) {
                  unsigned char *d = px;
                  const unsigned char *s = spx;
                  for (c = w; c > 0; c--) {
                        if (s[3] == 0) {
                              /* Transparent FG, NOP */
                        } else if (s[3] == 255) {
                              /* Full coverage, demul src */
                              //   dc' = ((1 - sa) * da*dc + sc)/da' = sc/da' = sc
                              //   da' = 1 - (1 - sa) * (1 - da) = 1 - 0 * (1 - da) = 1
                              memcpy(d, s, 4);
                        } else if (d[3] == 0) {
                              /* Full coverage, demul src */
                              //   dc' = ((1 - sa) * da*dc + sc)/da' = sc/da' = sc/sa = sc/sa
                              //   da' = 1 - (1 - sa) * (1 - da) = 1 - (1 - sa) = sa
                              d[0] = NR_DEMUL_111(s[0], s[3]);
                              d[1] = NR_DEMUL_111(s[1], s[3]);
                              d[2] = NR_DEMUL_111(s[2], s[3]);
                              d[3] = s[3];
                        } else {
                              /* Full composition */
                              //   dc' = ((1 - sa) * da*dc + sc)/da' = ((1 - sa) * da*dc + sc)/da'
                              //   da' = 1 - (1 - sa) * (1 - da) = 1 - (1 - sa) * (1 - da)
                              unsigned int da = NR_COMPOSEA_112(s[3], d[3]);
                              d[0] = NR_COMPOSEPNN_111121(s[0], s[3], d[0], d[3], da);
                              d[1] = NR_COMPOSEPNN_111121(s[1], s[3], d[1], d[3], da);
                              d[2] = NR_COMPOSEPNN_111121(s[2], s[3], d[2], d[3], da);
                              d[3] = NR_NORMALIZE_21(da);
                        }
                        d += 4;
                        s += 4;
                  }
                  px += rs;
                  spx += srs;
            }
      } else {
            for (r = h; r > 0; r--) {
                  unsigned char *d = px;
                  const unsigned char *s = spx;
                  for (c = w; c > 0; c--) {
                        unsigned int a;
                        a = NR_PREMUL_112(s[3], alpha);
                        if (a == 0) {
                              /* Transparent FG, NOP */
                        } else if (d[3] == 0) {
                              /* Full coverage, demul src */
                              //   dc' = ((1 - alpha*sa) * da*dc + alpha*sc)/da' = alpha*sc/da' = alpha*sc/(alpha*sa) = sc/sa
                              //   da' = 1 - (1 - alpha*sa) * (1 - da) = 1 - (1 - alpha*sa) = alpha*sa
                              d[0] = NR_DEMUL_111(s[0], s[3]);
                              d[1] = NR_DEMUL_111(s[1], s[3]);
                              d[2] = NR_DEMUL_111(s[2], s[3]);
                              d[3] = NR_NORMALIZE_21(a);
                        } else {
                              //   dc' = ((1 - alpha*sa) * da*dc + alpha*sc)/da'
                              //   da' = 1 - (1 - alpha*sa) * (1 - da)
                              unsigned int da = NR_COMPOSEA_213(a, d[3]);
                              d[0] = NR_COMPOSEPNN_221131(NR_PREMUL_112(s[0], alpha), a, d[0], d[3], da);
                              d[1] = NR_COMPOSEPNN_221131(NR_PREMUL_112(s[1], alpha), a, d[1], d[3], da);
                              d[2] = NR_COMPOSEPNN_221131(NR_PREMUL_112(s[2], alpha), a, d[2], d[3], da);
                              d[3] = NR_NORMALIZE_31(da);
                        }
                        d += 4;
                        s += 4;
                  }
                  px += rs;
                  spx += srs;
            }
      }
}

void
nr_R8G8B8A8_P_R8G8B8A8_P_R8G8B8A8_N (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, unsigned int alpha)
{
      unsigned int r, c;

      if (alpha == 0) {
            /* NOP */
      } else if (alpha == 255) {
            for (r = h; r > 0; r--) {
                  unsigned char *d = px;
                  const unsigned char *s = spx;
                  for (c = w; c > 0; c--) {
                        if (s[3] == 0) {
                              /* Transparent FG, NOP */
                        } else if (s[3] == 255) {
                              /* Opaque FG, COPY */
                              //   dc' = (1 - sa) * dc + sa*sc = sa*sc = sc
                              //   da' = 1 - (1 - sa) * (1 - da) = 1 - 0 * (1 - da) = 1 (= sa)
                              memcpy(d, s, 4);
                        } else if (d[3] == 0) {
                              /* Transparent BG, premul src */
                              //   dc' = (1 - sa) * dc + sa*sc = sa*sc
                              //   da' = 1 - (1 - sa) * (1 - da) = 1 - (1 - sa) = sa
                              d[0] = NR_PREMUL_111(s[0], s[3]);
                              d[1] = NR_PREMUL_111(s[1], s[3]);
                              d[2] = NR_PREMUL_111(s[2], s[3]);
                              d[3] = s[3];
                        } else {
                              //   dc' = (1 - sa) * dc + sa*sc
                              //   da' = 1 - (1 - sa) * (1 - da)
                              d[0] = NR_COMPOSENPP_1111(s[0], s[3], d[0]);
                              d[1] = NR_COMPOSENPP_1111(s[1], s[3], d[1]);
                              d[2] = NR_COMPOSENPP_1111(s[2], s[3], d[2]);
                              d[3] = NR_COMPOSEA_111(s[3], d[3]);
                        }
                        d += 4;
                        s += 4;
                  }
                  px += rs;
                  spx += srs;
            }
      } else {
            for (r = h; r > 0; r--) {
                  unsigned char *d = px;
                  const unsigned char *s = spx;
                  for (c = w; c > 0; c--) {
                        unsigned int a;
                        a = NR_PREMUL_112 (s[3], alpha);
                        if (a == 0) {
                              /* Transparent FG, NOP */
                        } else if (d[3] == 0) {
                              /* Transparent BG, premul src */
                              //   dc' = (1 - alpha*sa) * dc + alpha*sa*sc = alpha*sa*sc
                              //   da' = 1 - (1 - alpha*sa) * (1 - da) = 1 - (1 - alpha*sa) = alpha*sa
                              d[0] = NR_PREMUL_121(s[0], a);
                              d[1] = NR_PREMUL_121(s[1], a);
                              d[2] = NR_PREMUL_121(s[2], a);
                              d[3] = NR_NORMALIZE_21(a);
                        } else {
                              //   dc' = (1 - alpha*sa) * dc + alpha*sa*sc
                              //   da' = 1 - (1 - alpha*sa) * (1 - da)
                              d[0] = NR_COMPOSENPP_1211(s[0], a, d[0]);
                              d[1] = NR_COMPOSENPP_1211(s[1], a, d[1]);
                              d[2] = NR_COMPOSENPP_1211(s[2], a, d[2]);
                              d[3] = NR_COMPOSEA_211(a, d[3]);
                        }
                        d += 4;
                        s += 4;
                  }
                  px += rs;
                  spx += srs;
            }
      }
}

void
nr_R8G8B8A8_P_R8G8B8A8_P_R8G8B8A8_P (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, unsigned int alpha)
{
      unsigned int r, c;

      if (alpha == 0) {
            /* Transparent FG, NOP */
      } else if (alpha == 255) {
            /* Simple */
            for (r = h; r > 0; r--) {
                  unsigned char *d = px;
                  const unsigned char *s = spx;
                  for (c = w; c > 0; c--) {
                        if (s[3] == 0) {
                              /* Transparent FG, NOP */
                        } else if ((s[3] == 255) || (d[3] == 0)) {
                              /* Transparent BG, COPY */
                              memcpy(d, s, 4);
                        } else {
                              d[0] = NR_COMPOSEPPP_1111(s[0], s[3], d[0]);
                              d[1] = NR_COMPOSEPPP_1111(s[1], s[3], d[1]);
                              d[2] = NR_COMPOSEPPP_1111(s[2], s[3], d[2]);
                              d[3] = NR_COMPOSEA_111(s[3], d[3]);
                        }
                        d += 4;
                        s += 4;
                  }
                  px += rs;
                  spx += srs;
            }
      } else {
            for (r = h; r > 0; r--) {
                  unsigned char *d = px;
                  const unsigned char *s = spx;
                  for (c = w; c > 0; c--) {
                        if (s[3] == 0) {
                              /* Transparent FG, NOP */
                        } else if (d[3] == 0) {
                              /* Transparent BG, COPY */
                              d[0] = NR_PREMUL_111(s[0], alpha);
                              d[1] = NR_PREMUL_111(s[1], alpha);
                              d[2] = NR_PREMUL_111(s[2], alpha);
                              d[3] = NR_PREMUL_111(s[3], alpha);
                        } else {
                              //   dc' = (1 - alpha*sa) * dc + alpha*sc
                              //   da' = 1 - (1 - alpha*sa) * (1 - da)
                              unsigned int a;
                              a = NR_PREMUL_112(s[3], alpha);
                              d[0] = NR_COMPOSEPPP_2211(NR_PREMUL_112(alpha, s[0]), a, d[0]);
                              d[1] = NR_COMPOSEPPP_2211(NR_PREMUL_112(alpha, s[1]), a, d[1]);
                              d[2] = NR_COMPOSEPPP_2211(NR_PREMUL_112(alpha, s[2]), a, d[2]);
                              d[3] = NR_COMPOSEA_211(a, d[3]);
                        }
                        d += 4;
                        s += 4;
                  }
                  px += rs;
                  spx += srs;
            }
      }
}

/* Masked operations */

void
nr_R8G8B8A8_N_EMPTY_R8G8B8A8_N_A8 (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, const unsigned char *mpx, int mrs)
{
      unsigned int r, c;

      for (r = h; r > 0; r--) {
            unsigned char *d = px;
            const unsigned char *s = spx;
            const unsigned char *m = mpx;
            for (c = w; c > 0; c--) {
                  d[0] = s[0];
                  d[1] = s[1];
                  d[2] = s[2];
                  d[3] = NR_PREMUL_111(s[3], m[0]);
                  d += 4;
                  s += 4;
                  m += 1;
            }
            px += rs;
            spx += srs;
            mpx += mrs;
      }
}

void
nr_R8G8B8A8_N_EMPTY_R8G8B8A8_P_A8 (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, const unsigned char *mpx, int mrs)
{
      unsigned int r, c;

      for (r = h; r > 0; r--) {
            unsigned char *d = px;
            const unsigned char *s = spx;
            const unsigned char *m = mpx;
            for (c = w; c > 0; c--) {
                  unsigned int a;
                  a = NR_PREMUL_112 (s[3], m[0]);
                  if (a == 0) {
                        d[3] = 0;
                  } else if (a == 255*255) {
                        memcpy(d, s, 4);
                  } else {
                        //   dc' = ((1 - m*sa) * da*dc + m*sc)/da' = m*sc/da' = m*sc/(m*sa) = sc/sa
                        //   da' = 1 - (1 - m*sa) * (1 - da) = 1 - (1 - m*sa) = m*sa
                        d[0] = NR_DEMUL_111(s[0], s[3]);
                        d[1] = NR_DEMUL_111(s[1], s[3]);
                        d[2] = NR_DEMUL_111(s[2], s[3]);
                        d[3] = NR_NORMALIZE_21(a);
                  }
                  d += 4;
                  s += 4;
                  m += 1;
            }
            px += rs;
            spx += srs;
            mpx += mrs;
      }
}

void
nr_R8G8B8A8_P_EMPTY_R8G8B8A8_N_A8 (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, const unsigned char *mpx, int mrs)
{
      unsigned int r, c;

      for (r = h; r > 0; r--) {
            unsigned char *d = px;
            const unsigned char *s = spx;
            const unsigned char *m = mpx;
            for (c = w; c > 0; c--) {
                  unsigned int a;
                  a = NR_PREMUL_112(s[3], m[0]);
                  if (a == 0) {
                        memset(d, 0, 4);
                  } else if (a == 255*255) {
                        memcpy(d, s, 4);
                  } else {
                        d[0] = NR_PREMUL_121(s[0], a);
                        d[1] = NR_PREMUL_121(s[1], a);
                        d[2] = NR_PREMUL_121(s[2], a);
                        d[3] = NR_NORMALIZE_21(a);
                  }
                  d += 4;
                  s += 4;
                  m += 1;
            }
            px += rs;
            spx += srs;
            mpx += mrs;
      }
}

void
nr_R8G8B8A8_P_EMPTY_R8G8B8A8_P_A8 (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, const unsigned char *mpx, int mrs)
{
      unsigned int r, c;

      for (r = h; r > 0; r--) {
            unsigned char *d = px;
            const unsigned char *s = spx;
            const unsigned char *m = mpx;
            for (c = w; c > 0; c--) {
                  d[0] = NR_PREMUL_111(s[0], m[0]);
                  d[1] = NR_PREMUL_111(s[1], m[0]);
                  d[2] = NR_PREMUL_111(s[2], m[0]);
                  d[3] = NR_PREMUL_111(s[3], m[0]);
                  d += 4;
                  s += 4;
                  m += 1;
            }
            px += rs;
            spx += srs;
            mpx += mrs;
      }
}

void
nr_R8G8B8A8_N_R8G8B8A8_N_R8G8B8A8_N_A8 (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, const unsigned char *mpx, int mrs)
{
      unsigned int r, c;

      for (r = h; r > 0; r--) {
            unsigned char *d = px;
            const unsigned char *s = spx;
            const unsigned char *m = mpx;
            for (c = w; c > 0; c--) {
                  unsigned int a;
                  a = NR_PREMUL_112(s[3], m[0]);
                  if (a == 0) {
                        /* Transparent FG, NOP */
                  } else if ((a == 255*255) || (d[3] == 0)) {
                        /* Full coverage, COPY */
                        d[0] = s[0];
                        d[1] = s[1];
                        d[2] = s[2];
                        d[3] = NR_NORMALIZE_21(a);
                  } else {
                        /* Full composition */
                        unsigned int ca;
                        ca = NR_COMPOSEA_213(a, d[3]);
                        d[0] = NR_COMPOSENNN_121131(s[0], a, d[0], d[3], ca);
                        d[1] = NR_COMPOSENNN_121131(s[1], a, d[1], d[3], ca);
                        d[2] = NR_COMPOSENNN_121131(s[2], a, d[2], d[3], ca);
                        d[3] = NR_NORMALIZE_31(ca);
                  }
                  d += 4;
                  s += 4;
                  m += 1;
            }
            px += rs;
            spx += srs;
            mpx += mrs;
      }
}

void
nr_R8G8B8A8_N_R8G8B8A8_N_R8G8B8A8_P_A8 (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, const unsigned char *mpx, int mrs)
{
      unsigned int r, c;

      for (r = h; r > 0; r--) {
            unsigned char *d = px;
            const unsigned char *s = spx;
            const unsigned char *m = mpx;
            for (c = w; c > 0; c--) {
                  unsigned int a;
                  a = NR_PREMUL_112(s[3], m[0]);
                  if (a == 0) {
                        /* Transparent FG, NOP */
                  } else if (a == 255*255) {
                        /* Opaque FG, COPY */
                        memcpy(d, s, 4);
                  } else if (d[3] == 0) {
                        /* Full coverage, demul src */
                        //   dc' = ((1 - m*sa) * da*dc + m*sc)/da' = m*sc/da' = m*sc/(m*sa) = sc/sa
                        //   da' = 1 - (1 - m*sa) * (1 - da) = 1 - (1 - m*sa) = m*sa
                        d[0] = NR_DEMUL_111(s[0], s[3]);
                        d[1] = NR_DEMUL_111(s[1], s[3]);
                        d[2] = NR_DEMUL_111(s[2], s[3]);
                        d[3] = NR_NORMALIZE_21(a);
                  } else if (m[0] == 255) {
                        /* Full composition */
                        //   dc' = ((1 - m*sa) * da*dc + m*sc)/da' = ((1 - sa) * da*dc + sc)/da'
                        //   da' = 1 - (1 - m*sa) * (1 - da) = 1 - (1 - sa) * (1 - da)
                        unsigned int da = NR_COMPOSEA_112(s[3], d[3]);
                        d[0] = NR_COMPOSEPNN_111121(s[0], s[3], d[0], d[3], da);
                        d[1] = NR_COMPOSEPNN_111121(s[1], s[3], d[1], d[3], da);
                        d[2] = NR_COMPOSEPNN_111121(s[2], s[3], d[2], d[3], da);
                        d[3] = NR_NORMALIZE_21(da);
                  } else {
                        //   dc' = ((1 - m*sa) * da*dc + m*sc)/da'
                        //   da' = 1 - (1 - m*sa) * (1 - da)
                        unsigned int da = NR_COMPOSEA_213(a, d[3]);
                        d[0] = NR_COMPOSEPNN_221131(NR_PREMUL_112(s[0], m[0]), a, d[0], d[3], da);
                        d[1] = NR_COMPOSEPNN_221131(NR_PREMUL_112(s[1], m[0]), a, d[1], d[3], da);
                        d[2] = NR_COMPOSEPNN_221131(NR_PREMUL_112(s[2], m[0]), a, d[2], d[3], da);
                        d[3] = NR_NORMALIZE_31(da);
                  }
                  d += 4;
                  s += 4;
                  m += 1;
            }
            px += rs;
            spx += srs;
            mpx += mrs;
      }
}

void
nr_R8G8B8A8_P_R8G8B8A8_P_R8G8B8A8_N_A8 (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, const unsigned char *mpx, int mrs)
{
      unsigned int r, c;

      for (r = h; r>0; r--) {
            unsigned char *d = px;
            const unsigned char *s = spx;
            const unsigned char *m = mpx;
            for (c = w; c>0; c--) {
                  unsigned int a;
                  a = NR_PREMUL_112(s[3], m[0]);
                  if (a == 0) {
                        /* Transparent FG, NOP */
                  } else if (a == 255*255) {
                        memcpy(d, s, 4);
                  } else {
                        d[0] = NR_COMPOSENPP_1211(s[0], a, d[0]);
                        d[1] = NR_COMPOSENPP_1211(s[1], a, d[1]);
                        d[2] = NR_COMPOSENPP_1211(s[2], a, d[2]);
                        d[3] = NR_COMPOSEA_211(a, d[3]);
                  }
                  d += 4;
                  s += 4;
                  m += 1;
            }
            px += rs;
            spx += srs;
            mpx += mrs;
      }
}

void
nr_R8G8B8A8_P_R8G8B8A8_P_R8G8B8A8_P_A8 (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, const unsigned char *mpx, int mrs)
{
      unsigned int r, c;

      for (r = h; r > 0; r--) {
            unsigned char *d = px;
            const unsigned char *s = spx;
            const unsigned char *m = mpx;
            for (c = w; c > 0; c--) {
                  unsigned int a;
                  a = NR_PREMUL_112 (s[3], m[0]);
                  if (a == 0) {
                        /* Transparent FG, NOP */
                  } else if (a == 255*255) {
                        /* Opaque FG, COPY */
                        memcpy(d, s, 4);
                  } else if (d[3] == 0) {
                        /* Transparent BG, COPY */
                        //   dc' = (1 - m*sa) * dc + m*sc = m*sc
                        //   da' = 1 - (1 - m*sa) * (1 - da) = 1 - (1 - m*sa)  = m*sa
                        d[0] = NR_PREMUL_111 (s[0], m[0]);
                        d[1] = NR_PREMUL_111 (s[1], m[0]);
                        d[2] = NR_PREMUL_111 (s[2], m[0]);
                        d[3] = NR_NORMALIZE_21(a);
                  } else {
                        //   dc' = (1 - m*sa) * dc + m*sc
                        //   da' = 1 - (1 - m*sa) * (1 - da)
                        d[0] = NR_COMPOSEPPP_2211 (NR_PREMUL_112 (s[0], m[0]), a, d[0]);
                        d[1] = NR_COMPOSEPPP_2211 (NR_PREMUL_112 (s[1], m[0]), a, d[1]);
                        d[2] = NR_COMPOSEPPP_2211 (NR_PREMUL_112 (s[2], m[0]), a, d[2]);
                        d[3] = NR_COMPOSEA_211(a, d[3]);
                  }
                  d += 4;
                  s += 4;
                  m += 1;
            }
            px += rs;
            spx += srs;
            mpx += mrs;
      }
}

/* FINAL DST MASK COLOR */

void
nr_R8G8B8A8_N_EMPTY_A8_RGBA32 (unsigned char *px, int w, int h, int rs, const unsigned char *mpx, int mrs, unsigned long rgba)
{
      unsigned int r, g, b, a;
      unsigned int x, y;

      r = NR_RGBA32_R (rgba);
      g = NR_RGBA32_G (rgba);
      b = NR_RGBA32_B (rgba);
      a = NR_RGBA32_A (rgba);

      for (y = h; y > 0; y--) {
            if (a == 0) {
                  memset(px, 0, w*4);
            } else {
                  unsigned char *d = px;
                  const unsigned char *m = mpx;
                  for (x = w; x > 0; x--) {
                        d[0] = r;
                        d[1] = g;
                        d[2] = b;
                        d[3] = NR_PREMUL_111 (m[0], a);
                        d += 4;
                        m += 1;
                  }
            }
            px += rs;
            mpx += mrs;
      }
}

void
nr_R8G8B8A8_P_EMPTY_A8_RGBA32 (unsigned char *px, int w, int h, int rs, const unsigned char *mpx, int mrs, unsigned long rgba)
{
      unsigned int r, g, b, a;
      unsigned int x, y;

      r = NR_RGBA32_R (rgba);
      g = NR_RGBA32_G (rgba);
      b = NR_RGBA32_B (rgba);
      a = NR_RGBA32_A (rgba);

#ifdef WITH_MMX
      if (NR_PIXOPS_MMX) {
            unsigned char c[4];
            c[0] = NR_PREMUL_111 (r, a);
            c[1] = NR_PREMUL_111 (g, a);
            c[2] = NR_PREMUL_111 (b, a);
            c[3] = a;
            /* WARNING: MMX composer REQUIRES w > 0 and h > 0 */
            nr_mmx_R8G8B8A8_P_EMPTY_A8_RGBAP (px, w, h, rs, mpx, mrs, c);
            // This mmx optimized code is approx. 2x faster than the non-optimized code below (Measured by Diederik van Lierop, 2009-12-17)
            return;
      }
#endif

      if ( a != 255 ){
            // Pre-premultiply color values
            r *= a;
            g *= a;
            b *= a;
      }

      for (y = h; y > 0; y--) {
            unsigned char *d = px;
            const unsigned char *m = mpx;
            if (a == 0) {
                  memset(px, 0, w*4);
            } else if (a == 255) {
                  for (x = w; x > 0; x--) {
                        d[0] = NR_PREMUL_111(m[0], r);
                        d[1] = NR_PREMUL_111(m[0], g);
                        d[2] = NR_PREMUL_111(m[0], b);
                        d[3] = m[0];
                        d += 4;
                        m += 1;
                  }
            } else {
                  for (x = w; x > 0; x--) {
                        // Color values are already premultiplied with a
                        d[0] = NR_PREMUL_121(m[0], r);
                        d[1] = NR_PREMUL_121(m[0], g);
                        d[2] = NR_PREMUL_121(m[0], b);
                        d[3] = NR_PREMUL_111(m[0], a);
                        d += 4;
                        m += 1;
                  }
            }
            px += rs;
            mpx += mrs;
      }
}

void
nr_R8G8B8_R8G8B8_A8_RGBA32 (unsigned char *px, int w, int h, int rs, const unsigned char *mpx, int mrs, unsigned long rgba)
{
      unsigned int r, g, b, a;
      unsigned int x, y;

      r = NR_RGBA32_R (rgba);
      g = NR_RGBA32_G (rgba);
      b = NR_RGBA32_B (rgba);
      a = NR_RGBA32_A (rgba);

      if (a == 0) {
            /* NOP */
      } else if (a == 255) {
            for (y = h; y > 0; y--) {
                  unsigned char *d = px;
                  const unsigned char *m = mpx;
                  for (x = w; x > 0; x--) {
                        d[0] = NR_COMPOSEN11_1111 (r, m[0], d[0]);
                        d[1] = NR_COMPOSEN11_1111 (g, m[0], d[1]);
                        d[2] = NR_COMPOSEN11_1111 (b, m[0], d[2]);
                        d += 3;
                        m += 1;
                  }
                  px += rs;
                  mpx += mrs;
            }
      } else {
            for (y = h; y > 0; y--) {
                  unsigned char *d = px;
                  const unsigned char *m = mpx;
                  for (x = w; x > 0; x--) {
                        //   dc' = (1 - m*sa) * dc + m*sa*sc
                        unsigned int alpha;
                        alpha = NR_PREMUL_112 (a, m[0]);
                        d[0] = NR_COMPOSEN11_1211 (r, alpha, d[0]);
                        d[1] = NR_COMPOSEN11_1211 (g, alpha, d[1]);
                        d[2] = NR_COMPOSEN11_1211 (b, alpha, d[2]);
                        d += 3;
                        m += 1;
                  }
                  px += rs;
                  mpx += mrs;
            }
      }
}

void
nr_R8G8B8A8_N_R8G8B8A8_N_A8_RGBA32 (unsigned char *px, int w, int h, int rs, const unsigned char *mpx, int mrs, unsigned long rgba)
{
      unsigned int r, g, b, a;
      unsigned int x, y;

      r = NR_RGBA32_R (rgba);
      g = NR_RGBA32_G (rgba);
      b = NR_RGBA32_B (rgba);
      a = NR_RGBA32_A (rgba);

      if (a == 0) {
            /* NOP */
      } else if (a == 255) {
            for (y = h; y > 0; y--) {
                  unsigned char *d = px;
                  const unsigned char *m = mpx;
                  for (x = w; x > 0; x--) {
                        if (m[0] == 0) {
                              /* Transparent FG, NOP */
                        } else if (m[0] == 255 || d[3] == 0) {
                              /* Full coverage, COPY */
                              d[0] = r;
                              d[1] = g;
                              d[2] = b;
                              d[3] = m[0];
                        } else {
                              /* Full composition */
                              unsigned int da = NR_COMPOSEA_112(m[0], d[3]);
                              d[0] = NR_COMPOSENNN_111121(r, m[0], d[0], d[3], da);
                              d[1] = NR_COMPOSENNN_111121(g, m[0], d[1], d[3], da);
                              d[2] = NR_COMPOSENNN_111121(b, m[0], d[2], d[3], da);
                              d[3] = NR_NORMALIZE_21(da);
                        }
                        d += 4;
                        m += 1;
                  }
                  px += rs;
                  mpx += mrs;
            }
      } else {
            for (y = h; y > 0; y--) {
                  unsigned char *d = px;
                  const unsigned char *m = mpx;
                  for (x = w; x > 0; x--) {
                        unsigned int ca;
                        ca = NR_PREMUL_112 (m[0], a);
                        if (ca == 0) {
                              /* Transparent FG, NOP */
                        } else if (d[3] == 0) {
                              /* Full coverage, COPY */
                              d[0] = r;
                              d[1] = g;
                              d[2] = b;
                              d[3] = NR_NORMALIZE_21(ca);
                        } else {
                              /* Full composition */
                              unsigned int da = NR_COMPOSEA_213(ca, d[3]);
                              d[0] = NR_COMPOSENNN_121131(r, ca, d[0], d[3], da);
                              d[1] = NR_COMPOSENNN_121131(g, ca, d[1], d[3], da);
                              d[2] = NR_COMPOSENNN_121131(b, ca, d[2], d[3], da);
                              d[3] = NR_NORMALIZE_31(da);
                        }
                        d += 4;
                        m += 1;
                  }
                  px += rs;
                  mpx += mrs;
            }
      }
}

void
nr_R8G8B8A8_P_R8G8B8A8_P_A8_RGBA32 (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, unsigned long rgba)
{
      unsigned int r, g, b, a;
      unsigned int x, y;

      r = NR_RGBA32_R (rgba);
      g = NR_RGBA32_G (rgba);
      b = NR_RGBA32_B (rgba);
      a = NR_RGBA32_A (rgba);

#ifdef WITH_MMX
      if (NR_PIXOPS_MMX && a != 0) {
            unsigned char c[4];
            c[0] = NR_PREMUL_111 (r, a);
            c[1] = NR_PREMUL_111 (g, a);
            c[2] = NR_PREMUL_111 (b, a);
            c[3] = a;
            /* WARNING: MMX composer REQUIRES w > 0 and h > 0 */
            nr_mmx_R8G8B8A8_P_R8G8B8A8_P_A8_RGBAP (px, w, h, rs, spx, srs, c);
            return;
      }
#endif

      if (a == 0) {
            /* Transparent FG, NOP */
      } else if (a == 255) {
            /* Simple */
            for (y = h; y > 0; y--) {
                  unsigned char *d, *s;
                  d = (unsigned char *) px;
                  s = (unsigned char *) spx;
                  for (x = w; x > 0; x--) {
                        if (s[0] == 0) {
                              /* Transparent FG, NOP */
                        } else {
                              /* Full composition */
                              unsigned int invca = 255-s[0]; // By swapping the arguments GCC can better optimize these calls
                              d[0] = NR_COMPOSENPP_1111(d[0], invca, r);
                              d[1] = NR_COMPOSENPP_1111(d[1], invca, g);
                              d[2] = NR_COMPOSENPP_1111(d[2], invca, b);
                              d[3] = NR_COMPOSEA_111(s[0], d[3]);
                        }
                        d += 4;
                        s += 1;
                  }
                  px += rs;
                  spx += srs;
            }
      } else {
            for (y = h; y > 0; y--) {
                  unsigned char *d, *s;
                  d = (unsigned char *) px;
                  s = (unsigned char *) spx;
                  for (x = w; x > 0; x--) {
                        unsigned int ca;
                        ca = NR_PREMUL_112 (s[0], a);
                        if (ca == 0) {
                              /* Transparent FG, NOP */
                        } else {
                              /* Full composition */
                              unsigned int invca = 255*255-ca; // By swapping the arguments GCC can better optimize these calls
                              d[0] = NR_COMPOSENPP_1211(d[0], invca, r);
                              d[1] = NR_COMPOSENPP_1211(d[1], invca, g);
                              d[2] = NR_COMPOSENPP_1211(d[2], invca, b);
                              d[3] = NR_COMPOSEA_211(ca, d[3]);
                        }
                        d += 4;
                        s += 1;
                  }
                  px += rs;
                  spx += srs;
            }
      }
}

/* RGB */

void
nr_R8G8B8_R8G8B8_R8G8B8A8_P (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, unsigned int alpha)
{
      unsigned int r, c;

#ifdef WITH_MMX
      if (NR_PIXOPS_MMX && alpha != 0) {
            /* WARNING: MMX composer REQUIRES w > 0 and h > 0 */
            nr_mmx_R8G8B8_R8G8B8_R8G8B8A8_P (px, w, h, rs, spx, srs, alpha);
            return;
      }
#endif

      if (alpha == 0) {
            /* NOP */
      } else if (alpha == 255) {
            for (r = h; r > 0; r--) {
                  unsigned char *d = px;
                  const unsigned char *s = spx;
                  for (c = w; c > 0; c--) {
                        //   dc' = (1 - alpha*sa) * dc + alpha*sc = (1 - sa) * dc + sc
                        if (s[3] == 0) {
                              /* NOP */
                        } else if (s[3] == 255) {
                              d[0] = s[0];
                              d[1] = s[1];
                              d[2] = s[2];
                        } else {
                              d[0] = NR_COMPOSEP11_1111(s[0], s[3], d[0]);
                              d[1] = NR_COMPOSEP11_1111(s[1], s[3], d[1]);
                              d[2] = NR_COMPOSEP11_1111(s[2], s[3], d[2]);
                        }
                        d += 3;
                        s += 4;
                  }
                  px += rs;
                  spx += srs;
            }
      } else {
            for (r = h; r > 0; r--) {
                  unsigned char *d = px;
                  const unsigned char *s = spx;
                  for (c = w; c > 0; c--) {
                        unsigned int a;
                        a = NR_PREMUL_112(s[3], alpha);
                        //   dc' = (1 - alpha*sa) * dc + alpha*sc
                        if (a == 0) {
                              /* NOP */
                        } else {
                              d[0] = NR_COMPOSEP11_2211(NR_PREMUL_112(s[0], alpha), a, d[0]);
                              d[1] = NR_COMPOSEP11_2211(NR_PREMUL_112(s[1], alpha), a, d[1]);
                              d[2] = NR_COMPOSEP11_2211(NR_PREMUL_112(s[2], alpha), a, d[2]);
                        }
                        /* a == 255 is impossible, because alpha < 255 */
                        d += 3;
                        s += 4;
                  }
                  px += rs;
                  spx += srs;
            }
      }
}

void
nr_R8G8B8_R8G8B8_R8G8B8A8_N (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, unsigned int alpha)
{
      unsigned int r, c;

      if (alpha == 0) {
            /* NOP */
      } else if (alpha == 255) {
            for (r = h; r > 0; r--) {
                  unsigned char *d = px;
                  const unsigned char *s = spx;
                  for (c = w; c > 0; c--) {
                        //   dc' = (1 - alpha*sa) * dc + alpha*sa*sc = (1 - sa) * dc + sa*sc
                        if (s[3] == 0) {
                              /* NOP */
                        } else if (s[3] == 255) {
                              d[0] = s[0];
                              d[1] = s[1];
                              d[2] = s[2];
                        } else {
                              d[0] = NR_COMPOSEN11_1111(s[0], s[3], d[0]);
                              d[1] = NR_COMPOSEN11_1111(s[1], s[3], d[1]);
                              d[2] = NR_COMPOSEN11_1111(s[2], s[3], d[2]);
                        }
                        d += 3;
                        s += 4;
                  }
                  px += rs;
                  spx += srs;
            }
      } else {
            for (r = h; r > 0; r--) {
                  unsigned char *d = px;
                  const unsigned char *s = spx;
                  for (c = w; c > 0; c--) {
                        unsigned int a;
                        a = NR_PREMUL_112(s[3], alpha);
                        //   dc' = (1 - alpha*sa) * dc + alpha*sa*sc
                        if (a == 0) {
                              /* NOP */
                        } else {
                              d[0] = NR_COMPOSEN11_1211(s[0], a, d[0]);
                              d[1] = NR_COMPOSEN11_1211(s[1], a, d[1]);
                              d[2] = NR_COMPOSEN11_1211(s[2], a, d[2]);
                        }
                        /* a == 255 is impossible, because alpha < 255 */
                        d += 3;
                        s += 4;
                  }
                  px += rs;
                  spx += srs;
            }
      }
}

void
nr_R8G8B8_R8G8B8_R8G8B8A8_P_A8 (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, const unsigned char *mpx, int mrs)
{
      unsigned int x, y;

      for (y = h; y > 0; y--) {
            unsigned char* d = px;
            const unsigned char* s = spx;
            const unsigned char* m = mpx;
            for (x = w; x > 0; x--) {
                  unsigned int a;
                  a = NR_PREMUL_112(s[3], m[0]);
                  if (a == 0) {
                        /* NOP */
                  } else if (a == 255*255) {
                        memcpy(d, s, 3);
                  } else {
                        //   dc' = (1 - m*sa) * dc + m*sc
                        d[0] = NR_COMPOSEP11_2211(NR_PREMUL_112(s[0], m[0]), a, d[0]);
                        d[1] = NR_COMPOSEP11_2211(NR_PREMUL_112(s[1], m[0]), a, d[1]);
                        d[2] = NR_COMPOSEP11_2211(NR_PREMUL_112(s[2], m[0]), a, d[2]);
                  }
                  d += 3;
                  s += 4;
                  m += 1;
            }
            px += rs;
            spx += srs;
            mpx += mrs;
      }
}

void
nr_R8G8B8_R8G8B8_R8G8B8A8_N_A8 (unsigned char *px, int w, int h, int rs, const unsigned char *spx, int srs, const unsigned char *mpx, int mrs)
{
      unsigned int x, y;

      for (y = h; y > 0; y--) {
            unsigned char* d = px;
            const unsigned char* s = spx;
            const unsigned char* m = mpx;
            for (x = w; x > 0; x--) {
                  unsigned int a;
                  a = NR_PREMUL_112(s[3], m[0]);
                  if (a == 0) {
                        /* NOP */
                  } else if (a == 255*255) {
                        memcpy(d, s, 3);
                  } else {
                        //   dc' = (1 - m*sa) * dc + m*sa*sc
                        d[0] = NR_COMPOSEN11_1211(s[0], a, d[0]);
                        d[1] = NR_COMPOSEN11_1211(s[1], a, d[1]);
                        d[2] = NR_COMPOSEN11_1211(s[2], a, d[2]);
                  }
                  d += 3;
                  s += 4;
                  m += 1;
            }
            px += rs;
            spx += srs;
            mpx += mrs;
      }
}



Generated by  Doxygen 1.6.0   Back to index