summaryrefslogtreecommitdiff
path: root/gegl/gegl-algorithms-2x2-downscale.inc
blob: 51d4e1e990e0321a0d0339027e883db884098695 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
#define S(a)   ((DOWNSCALE_SUM)(a))

void
DOWNSCALE_FUNCNAME (gint    bpp,
                    gint    src_width,
                    gint    src_height,
                    guchar *src_data,
                    gint    src_rowstride,
                    guchar *dst_data,
                    gint    dst_rowstride)
{
  gint y;
  gint diag = src_rowstride + bpp;
  const gint components = bpp / sizeof(DOWNSCALE_TYPE);

  if (!src_data || !dst_data)
    return;

  switch (components)
  {
    case 1:
      for (y = 0; y < src_height / 2; y++)
        {
          gint    x;
          guchar *src = src_data + src_rowstride * y * 2;
          guchar *dst = dst_data + dst_rowstride * y;

          for (x = 0; x < src_width / 2; x++)
            {
              DOWNSCALE_TYPE * aa = ((DOWNSCALE_TYPE *)(src));
              DOWNSCALE_TYPE * ab = ((DOWNSCALE_TYPE *)(src + bpp));
              DOWNSCALE_TYPE * ba = ((DOWNSCALE_TYPE *)(src + src_rowstride));
              DOWNSCALE_TYPE * bb = ((DOWNSCALE_TYPE *)(src + diag));

              ((DOWNSCALE_TYPE *)dst)[0] = (S(aa[0]) + S(ab[0]) + S(ba[0]) + S(bb[0])) / DOWNSCALE_DIVISOR;

              dst += bpp;
              src += bpp * 2;
            }
        }
    break;
    case 2:
      for (y = 0; y < src_height / 2; y++)
        {
          gint    x;
          guchar *src = src_data + src_rowstride * y * 2;
          guchar *dst = dst_data + dst_rowstride * y;

          for (x = 0; x < src_width / 2; x++)
            {
              DOWNSCALE_TYPE * aa = ((DOWNSCALE_TYPE *)(src));
              DOWNSCALE_TYPE * ab = ((DOWNSCALE_TYPE *)(src + bpp));
              DOWNSCALE_TYPE * ba = ((DOWNSCALE_TYPE *)(src + src_rowstride));
              DOWNSCALE_TYPE * bb = ((DOWNSCALE_TYPE *)(src + diag));

              ((DOWNSCALE_TYPE *)dst)[0] = (S(aa[0]) + S(ab[0]) + S(ba[0]) +
                S(bb[0])) / DOWNSCALE_DIVISOR;
              ((DOWNSCALE_TYPE *)dst)[1] = (S(aa[1]) + S(ab[1]) + S(ba[1]) +
                  S(bb[1])) / DOWNSCALE_DIVISOR;

              dst += bpp;
              src += bpp * 2;
            }
        }
    break;
    case 3:
      for (y = 0; y < src_height / 2; y++)
        {
          gint    x;
          guchar *src = src_data + src_rowstride * y * 2;
          guchar *dst = dst_data + dst_rowstride * y;

          for (x = 0; x < src_width / 2; x++)
            {
              DOWNSCALE_TYPE * aa = ((DOWNSCALE_TYPE *)(src));
              DOWNSCALE_TYPE * ab = ((DOWNSCALE_TYPE *)(src + bpp));
              DOWNSCALE_TYPE * ba = ((DOWNSCALE_TYPE *)(src + src_rowstride));
              DOWNSCALE_TYPE * bb = ((DOWNSCALE_TYPE *)(src + diag));

              ((DOWNSCALE_TYPE *)dst)[0] =
               (S(aa[0]) + S(ab[0]) + S(ba[0]) + S(bb[0])) / DOWNSCALE_DIVISOR;
              ((DOWNSCALE_TYPE *)dst)[1] =
               (S(aa[1]) + S(ab[1]) + S(ba[1]) + S(bb[1])) / DOWNSCALE_DIVISOR;
              ((DOWNSCALE_TYPE *)dst)[2] =
               (S(aa[2]) + S(ab[2]) + S(ba[2]) + S(bb[2])) / DOWNSCALE_DIVISOR;

              dst += bpp;
              src += bpp * 2;
            }
        }
    break;
    case 4:
      for (y = 0; y < src_height / 2; y++)
        {
          gint    x;
          guchar *src = src_data + src_rowstride * y * 2;
          guchar *dst = dst_data + dst_rowstride * y;

          for (x = 0; x < src_width / 2; x++)
            {
              DOWNSCALE_TYPE * aa = ((DOWNSCALE_TYPE *)(src));
              DOWNSCALE_TYPE * ab = ((DOWNSCALE_TYPE *)(src + bpp));
              DOWNSCALE_TYPE * ba = ((DOWNSCALE_TYPE *)(src + src_rowstride));
              DOWNSCALE_TYPE * bb = ((DOWNSCALE_TYPE *)(src + diag));

              ((DOWNSCALE_TYPE *)dst)[0] =
               (S(aa[0]) + S(ab[0]) + S(ba[0]) + S(bb[0])) / DOWNSCALE_DIVISOR;
              ((DOWNSCALE_TYPE *)dst)[1] =
               (S(aa[1]) + S(ab[1]) + S(ba[1]) + S(bb[1])) / DOWNSCALE_DIVISOR;
              ((DOWNSCALE_TYPE *)dst)[2] =
               (S(aa[2]) + S(ab[2]) + S(ba[2]) + S(bb[2])) / DOWNSCALE_DIVISOR;
              ((DOWNSCALE_TYPE *)dst)[3] =
               (S(aa[3]) + S(ab[3]) + S(ba[3]) + S(bb[3])) / DOWNSCALE_DIVISOR;

              dst += bpp;
              src += bpp * 2;
            }
        }
    break;
    default:
      for (y = 0; y < src_height / 2; y++)
        {
          gint    x;
          guchar *src = src_data + src_rowstride * y * 2;
          guchar *dst = dst_data + dst_rowstride * y;

          for (x = 0; x < src_width / 2; x++)
            {
              gint i;
              DOWNSCALE_TYPE * aa = ((DOWNSCALE_TYPE *)(src));
              DOWNSCALE_TYPE * ab = ((DOWNSCALE_TYPE *)(src + bpp));
              DOWNSCALE_TYPE * ba = ((DOWNSCALE_TYPE *)(src + src_rowstride));
              DOWNSCALE_TYPE * bb = ((DOWNSCALE_TYPE *)(src + diag));

              for (i = 0; i < components; i++)
                ((DOWNSCALE_TYPE *)dst)[i] =
                  (S(aa[i]) + S(ab[i]) + S(ba[i]) + S(bb[i]))/DOWNSCALE_DIVISOR;

              dst += bpp;
              src += bpp * 2;
            }
        }
    break;
  }
}

#undef S