projects / mesa.git / blob
? search:
summary | shortlog | log | commit | commitdiff | tree
history | raw | HEAD
Merge remote branch 'origin/master' into nv50-compiler
[mesa.git] / src / mesa / main / mipmap.c
1 /*
2 * Mesa 3-D graphics library
3 * Version: 7.1
4 *
5 * Copyright (C) 1999-2007 Brian Paul All Rights Reserved.
6 *
7 * Permission is hereby granted, free of charge, to any person obtaining a
8 * copy of this software and associated documentation files (the "Software"),
9 * to deal in the Software without restriction, including without limitation
10 * the rights to use, copy, modify, merge, publish, distribute, sublicense,
11 * and/or sell copies of the Software, and to permit persons to whom the
12 * Software is furnished to do so, subject to the following conditions:
13 *
14 * The above copyright notice and this permission notice shall be included
15 * in all copies or substantial portions of the Software.
16 *
17 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
18 * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
19 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
20 * BRIAN PAUL BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN
21 * AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN
22 * CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
23 */
24
25
26 /**
27 * \file mipmap.c mipmap generation and teximage resizing functions.
28 */
29
30 #include "imports.h"
31 #include "formats.h"
32 #include "mipmap.h"
33 #include "teximage.h"
34 #include "texstore.h"
35 #include "image.h"
36
37
38
39 static GLint
40 bytes_per_pixel(GLenum datatype, GLuint comps)
41 {
42 GLint b = _mesa_sizeof_packed_type(datatype);
43 assert(b >= 0);
44
45 if (_mesa_type_is_packed(datatype))
46 return b;
47 else
48 return b * comps;
49 }
50
51
52 /**
53 * \name Support macros for do_row and do_row_3d
54 *
55 * The macro madness is here for two reasons. First, it compacts the code
56 * slightly. Second, it makes it much easier to adjust the specifics of the
57 * filter to tune the rounding characteristics.
58 */
59 /*@{*/
60 #define DECLARE_ROW_POINTERS(t, e) \
61 const t(*rowA)[e] = (const t(*)[e]) srcRowA; \
62 const t(*rowB)[e] = (const t(*)[e]) srcRowB; \
63 const t(*rowC)[e] = (const t(*)[e]) srcRowC; \
64 const t(*rowD)[e] = (const t(*)[e]) srcRowD; \
65 t(*dst)[e] = (t(*)[e]) dstRow
66
67 #define DECLARE_ROW_POINTERS0(t) \
68 const t *rowA = (const t *) srcRowA; \
69 const t *rowB = (const t *) srcRowB; \
70 const t *rowC = (const t *) srcRowC; \
71 const t *rowD = (const t *) srcRowD; \
72 t *dst = (t *) dstRow
73
74 #define FILTER_SUM_3D(Aj, Ak, Bj, Bk, Cj, Ck, Dj, Dk) \
75 ((unsigned) Aj + (unsigned) Ak \
76 + (unsigned) Bj + (unsigned) Bk \
77 + (unsigned) Cj + (unsigned) Ck \
78 + (unsigned) Dj + (unsigned) Dk \
79 + 4) >> 3
80
81 #define FILTER_3D(e) \
82 do { \
83 dst[i][e] = FILTER_SUM_3D(rowA[j][e], rowA[k][e], \
84 rowB[j][e], rowB[k][e], \
85 rowC[j][e], rowC[k][e], \
86 rowD[j][e], rowD[k][e]); \
87 } while(0)
88
89 #define FILTER_SUM_3D_SIGNED(Aj, Ak, Bj, Bk, Cj, Ck, Dj, Dk) \
90 (Aj + Ak \
91 + Bj + Bk \
92 + Cj + Ck \
93 + Dj + Dk \
94 + 4) / 8
95
96 #define FILTER_3D_SIGNED(e) \
97 do { \
98 dst[i][e] = FILTER_SUM_3D_SIGNED(rowA[j][e], rowA[k][e], \
99 rowB[j][e], rowB[k][e], \
100 rowC[j][e], rowC[k][e], \
101 rowD[j][e], rowD[k][e]); \
102 } while(0)
103
104 #define FILTER_F_3D(e) \
105 do { \
106 dst[i][e] = (rowA[j][e] + rowA[k][e] \
107 + rowB[j][e] + rowB[k][e] \
108 + rowC[j][e] + rowC[k][e] \
109 + rowD[j][e] + rowD[k][e]) * 0.125F; \
110 } while(0)
111
112 #define FILTER_HF_3D(e) \
113 do { \
114 const GLfloat aj = _mesa_half_to_float(rowA[j][e]); \
115 const GLfloat ak = _mesa_half_to_float(rowA[k][e]); \
116 const GLfloat bj = _mesa_half_to_float(rowB[j][e]); \
117 const GLfloat bk = _mesa_half_to_float(rowB[k][e]); \
118 const GLfloat cj = _mesa_half_to_float(rowC[j][e]); \
119 const GLfloat ck = _mesa_half_to_float(rowC[k][e]); \
120 const GLfloat dj = _mesa_half_to_float(rowD[j][e]); \
121 const GLfloat dk = _mesa_half_to_float(rowD[k][e]); \
122 dst[i][e] = _mesa_float_to_half((aj + ak + bj + bk + cj + ck + dj + dk) \
123 * 0.125F); \
124 } while(0)
125 /*@}*/
126
127
128 /**
129 * Average together two rows of a source image to produce a single new
130 * row in the dest image. It's legal for the two source rows to point
131 * to the same data. The source width must be equal to either the
132 * dest width or two times the dest width.
133 * \param datatype GL_UNSIGNED_BYTE, GL_UNSIGNED_SHORT, GL_FLOAT, etc.
134 * \param comps number of components per pixel (1..4)
135 */
136 static void
137 do_row(GLenum datatype, GLuint comps, GLint srcWidth,
138 const GLvoid *srcRowA, const GLvoid *srcRowB,
139 GLint dstWidth, GLvoid *dstRow)
140 {
141 const GLuint k0 = (srcWidth == dstWidth) ? 0 : 1;
142 const GLuint colStride = (srcWidth == dstWidth) ? 1 : 2;
143
144 ASSERT(comps >= 1);
145 ASSERT(comps <= 4);
146
147 /* This assertion is no longer valid with non-power-of-2 textures
148 assert(srcWidth == dstWidth || srcWidth == 2 * dstWidth);
149 */
150
151 if (datatype == GL_UNSIGNED_BYTE && comps == 4) {
152 GLuint i, j, k;
153 const GLubyte(*rowA)[4] = (const GLubyte(*)[4]) srcRowA;
154 const GLubyte(*rowB)[4] = (const GLubyte(*)[4]) srcRowB;
155 GLubyte(*dst)[4] = (GLubyte(*)[4]) dstRow;
156 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
157 i++, j += colStride, k += colStride) {
158 dst[i][0] = (rowA[j][0] + rowA[k][0] + rowB[j][0] + rowB[k][0]) / 4;
159 dst[i][1] = (rowA[j][1] + rowA[k][1] + rowB[j][1] + rowB[k][1]) / 4;
160 dst[i][2] = (rowA[j][2] + rowA[k][2] + rowB[j][2] + rowB[k][2]) / 4;
161 dst[i][3] = (rowA[j][3] + rowA[k][3] + rowB[j][3] + rowB[k][3]) / 4;
162 }
163 }
164 else if (datatype == GL_UNSIGNED_BYTE && comps == 3) {
165 GLuint i, j, k;
166 const GLubyte(*rowA)[3] = (const GLubyte(*)[3]) srcRowA;
167 const GLubyte(*rowB)[3] = (const GLubyte(*)[3]) srcRowB;
168 GLubyte(*dst)[3] = (GLubyte(*)[3]) dstRow;
169 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
170 i++, j += colStride, k += colStride) {
171 dst[i][0] = (rowA[j][0] + rowA[k][0] + rowB[j][0] + rowB[k][0]) / 4;
172 dst[i][1] = (rowA[j][1] + rowA[k][1] + rowB[j][1] + rowB[k][1]) / 4;
173 dst[i][2] = (rowA[j][2] + rowA[k][2] + rowB[j][2] + rowB[k][2]) / 4;
174 }
175 }
176 else if (datatype == GL_UNSIGNED_BYTE && comps == 2) {
177 GLuint i, j, k;
178 const GLubyte(*rowA)[2] = (const GLubyte(*)[2]) srcRowA;
179 const GLubyte(*rowB)[2] = (const GLubyte(*)[2]) srcRowB;
180 GLubyte(*dst)[2] = (GLubyte(*)[2]) dstRow;
181 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
182 i++, j += colStride, k += colStride) {
183 dst[i][0] = (rowA[j][0] + rowA[k][0] + rowB[j][0] + rowB[k][0]) >> 2;
184 dst[i][1] = (rowA[j][1] + rowA[k][1] + rowB[j][1] + rowB[k][1]) >> 2;
185 }
186 }
187 else if (datatype == GL_UNSIGNED_BYTE && comps == 1) {
188 GLuint i, j, k;
189 const GLubyte *rowA = (const GLubyte *) srcRowA;
190 const GLubyte *rowB = (const GLubyte *) srcRowB;
191 GLubyte *dst = (GLubyte *) dstRow;
192 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
193 i++, j += colStride, k += colStride) {
194 dst[i] = (rowA[j] + rowA[k] + rowB[j] + rowB[k]) >> 2;
195 }
196 }
197
198 else if (datatype == GL_BYTE && comps == 4) {
199 GLuint i, j, k;
200 const GLbyte(*rowA)[4] = (const GLbyte(*)[4]) srcRowA;
201 const GLbyte(*rowB)[4] = (const GLbyte(*)[4]) srcRowB;
202 GLbyte(*dst)[4] = (GLbyte(*)[4]) dstRow;
203 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
204 i++, j += colStride, k += colStride) {
205 dst[i][0] = (rowA[j][0] + rowA[k][0] + rowB[j][0] + rowB[k][0]) / 4;
206 dst[i][1] = (rowA[j][1] + rowA[k][1] + rowB[j][1] + rowB[k][1]) / 4;
207 dst[i][2] = (rowA[j][2] + rowA[k][2] + rowB[j][2] + rowB[k][2]) / 4;
208 dst[i][3] = (rowA[j][3] + rowA[k][3] + rowB[j][3] + rowB[k][3]) / 4;
209 }
210 }
211 else if (datatype == GL_BYTE && comps == 3) {
212 GLuint i, j, k;
213 const GLbyte(*rowA)[3] = (const GLbyte(*)[3]) srcRowA;
214 const GLbyte(*rowB)[3] = (const GLbyte(*)[3]) srcRowB;
215 GLbyte(*dst)[3] = (GLbyte(*)[3]) dstRow;
216 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
217 i++, j += colStride, k += colStride) {
218 dst[i][0] = (rowA[j][0] + rowA[k][0] + rowB[j][0] + rowB[k][0]) / 4;
219 dst[i][1] = (rowA[j][1] + rowA[k][1] + rowB[j][1] + rowB[k][1]) / 4;
220 dst[i][2] = (rowA[j][2] + rowA[k][2] + rowB[j][2] + rowB[k][2]) / 4;
221 }
222 }
223 else if (datatype == GL_BYTE && comps == 2) {
224 GLuint i, j, k;
225 const GLbyte(*rowA)[2] = (const GLbyte(*)[2]) srcRowA;
226 const GLbyte(*rowB)[2] = (const GLbyte(*)[2]) srcRowB;
227 GLbyte(*dst)[2] = (GLbyte(*)[2]) dstRow;
228 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
229 i++, j += colStride, k += colStride) {
230 dst[i][0] = (rowA[j][0] + rowA[k][0] + rowB[j][0] + rowB[k][0]) / 4;
231 dst[i][1] = (rowA[j][1] + rowA[k][1] + rowB[j][1] + rowB[k][1]) / 4;
232 }
233 }
234 else if (datatype == GL_BYTE && comps == 1) {
235 GLuint i, j, k;
236 const GLbyte *rowA = (const GLbyte *) srcRowA;
237 const GLbyte *rowB = (const GLbyte *) srcRowB;
238 GLbyte *dst = (GLbyte *) dstRow;
239 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
240 i++, j += colStride, k += colStride) {
241 dst[i] = (rowA[j] + rowA[k] + rowB[j] + rowB[k]) / 4;
242 }
243 }
244
245 else if (datatype == GL_UNSIGNED_SHORT && comps == 4) {
246 GLuint i, j, k;
247 const GLushort(*rowA)[4] = (const GLushort(*)[4]) srcRowA;
248 const GLushort(*rowB)[4] = (const GLushort(*)[4]) srcRowB;
249 GLushort(*dst)[4] = (GLushort(*)[4]) dstRow;
250 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
251 i++, j += colStride, k += colStride) {
252 dst[i][0] = (rowA[j][0] + rowA[k][0] + rowB[j][0] + rowB[k][0]) / 4;
253 dst[i][1] = (rowA[j][1] + rowA[k][1] + rowB[j][1] + rowB[k][1]) / 4;
254 dst[i][2] = (rowA[j][2] + rowA[k][2] + rowB[j][2] + rowB[k][2]) / 4;
255 dst[i][3] = (rowA[j][3] + rowA[k][3] + rowB[j][3] + rowB[k][3]) / 4;
256 }
257 }
258 else if (datatype == GL_UNSIGNED_SHORT && comps == 3) {
259 GLuint i, j, k;
260 const GLushort(*rowA)[3] = (const GLushort(*)[3]) srcRowA;
261 const GLushort(*rowB)[3] = (const GLushort(*)[3]) srcRowB;
262 GLushort(*dst)[3] = (GLushort(*)[3]) dstRow;
263 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
264 i++, j += colStride, k += colStride) {
265 dst[i][0] = (rowA[j][0] + rowA[k][0] + rowB[j][0] + rowB[k][0]) / 4;
266 dst[i][1] = (rowA[j][1] + rowA[k][1] + rowB[j][1] + rowB[k][1]) / 4;
267 dst[i][2] = (rowA[j][2] + rowA[k][2] + rowB[j][2] + rowB[k][2]) / 4;
268 }
269 }
270 else if (datatype == GL_UNSIGNED_SHORT && comps == 2) {
271 GLuint i, j, k;
272 const GLushort(*rowA)[2] = (const GLushort(*)[2]) srcRowA;
273 const GLushort(*rowB)[2] = (const GLushort(*)[2]) srcRowB;
274 GLushort(*dst)[2] = (GLushort(*)[2]) dstRow;
275 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
276 i++, j += colStride, k += colStride) {
277 dst[i][0] = (rowA[j][0] + rowA[k][0] + rowB[j][0] + rowB[k][0]) / 4;
278 dst[i][1] = (rowA[j][1] + rowA[k][1] + rowB[j][1] + rowB[k][1]) / 4;
279 }
280 }
281 else if (datatype == GL_UNSIGNED_SHORT && comps == 1) {
282 GLuint i, j, k;
283 const GLushort *rowA = (const GLushort *) srcRowA;
284 const GLushort *rowB = (const GLushort *) srcRowB;
285 GLushort *dst = (GLushort *) dstRow;
286 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
287 i++, j += colStride, k += colStride) {
288 dst[i] = (rowA[j] + rowA[k] + rowB[j] + rowB[k]) / 4;
289 }
290 }
291 else if (datatype == GL_FLOAT && comps == 4) {
292 GLuint i, j, k;
293 const GLfloat(*rowA)[4] = (const GLfloat(*)[4]) srcRowA;
294 const GLfloat(*rowB)[4] = (const GLfloat(*)[4]) srcRowB;
295 GLfloat(*dst)[4] = (GLfloat(*)[4]) dstRow;
296 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
297 i++, j += colStride, k += colStride) {
298 dst[i][0] = (rowA[j][0] + rowA[k][0] +
299 rowB[j][0] + rowB[k][0]) * 0.25F;
300 dst[i][1] = (rowA[j][1] + rowA[k][1] +
301 rowB[j][1] + rowB[k][1]) * 0.25F;
302 dst[i][2] = (rowA[j][2] + rowA[k][2] +
303 rowB[j][2] + rowB[k][2]) * 0.25F;
304 dst[i][3] = (rowA[j][3] + rowA[k][3] +
305 rowB[j][3] + rowB[k][3]) * 0.25F;
306 }
307 }
308 else if (datatype == GL_FLOAT && comps == 3) {
309 GLuint i, j, k;
310 const GLfloat(*rowA)[3] = (const GLfloat(*)[3]) srcRowA;
311 const GLfloat(*rowB)[3] = (const GLfloat(*)[3]) srcRowB;
312 GLfloat(*dst)[3] = (GLfloat(*)[3]) dstRow;
313 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
314 i++, j += colStride, k += colStride) {
315 dst[i][0] = (rowA[j][0] + rowA[k][0] +
316 rowB[j][0] + rowB[k][0]) * 0.25F;
317 dst[i][1] = (rowA[j][1] + rowA[k][1] +
318 rowB[j][1] + rowB[k][1]) * 0.25F;
319 dst[i][2] = (rowA[j][2] + rowA[k][2] +
320 rowB[j][2] + rowB[k][2]) * 0.25F;
321 }
322 }
323 else if (datatype == GL_FLOAT && comps == 2) {
324 GLuint i, j, k;
325 const GLfloat(*rowA)[2] = (const GLfloat(*)[2]) srcRowA;
326 const GLfloat(*rowB)[2] = (const GLfloat(*)[2]) srcRowB;
327 GLfloat(*dst)[2] = (GLfloat(*)[2]) dstRow;
328 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
329 i++, j += colStride, k += colStride) {
330 dst[i][0] = (rowA[j][0] + rowA[k][0] +
331 rowB[j][0] + rowB[k][0]) * 0.25F;
332 dst[i][1] = (rowA[j][1] + rowA[k][1] +
333 rowB[j][1] + rowB[k][1]) * 0.25F;
334 }
335 }
336 else if (datatype == GL_FLOAT && comps == 1) {
337 GLuint i, j, k;
338 const GLfloat *rowA = (const GLfloat *) srcRowA;
339 const GLfloat *rowB = (const GLfloat *) srcRowB;
340 GLfloat *dst = (GLfloat *) dstRow;
341 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
342 i++, j += colStride, k += colStride) {
343 dst[i] = (rowA[j] + rowA[k] + rowB[j] + rowB[k]) * 0.25F;
344 }
345 }
346
347 else if (datatype == GL_HALF_FLOAT_ARB && comps == 4) {
348 GLuint i, j, k, comp;
349 const GLhalfARB(*rowA)[4] = (const GLhalfARB(*)[4]) srcRowA;
350 const GLhalfARB(*rowB)[4] = (const GLhalfARB(*)[4]) srcRowB;
351 GLhalfARB(*dst)[4] = (GLhalfARB(*)[4]) dstRow;
352 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
353 i++, j += colStride, k += colStride) {
354 for (comp = 0; comp < 4; comp++) {
355 GLfloat aj, ak, bj, bk;
356 aj = _mesa_half_to_float(rowA[j][comp]);
357 ak = _mesa_half_to_float(rowA[k][comp]);
358 bj = _mesa_half_to_float(rowB[j][comp]);
359 bk = _mesa_half_to_float(rowB[k][comp]);
360 dst[i][comp] = _mesa_float_to_half((aj + ak + bj + bk) * 0.25F);
361 }
362 }
363 }
364 else if (datatype == GL_HALF_FLOAT_ARB && comps == 3) {
365 GLuint i, j, k, comp;
366 const GLhalfARB(*rowA)[3] = (const GLhalfARB(*)[3]) srcRowA;
367 const GLhalfARB(*rowB)[3] = (const GLhalfARB(*)[3]) srcRowB;
368 GLhalfARB(*dst)[3] = (GLhalfARB(*)[3]) dstRow;
369 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
370 i++, j += colStride, k += colStride) {
371 for (comp = 0; comp < 3; comp++) {
372 GLfloat aj, ak, bj, bk;
373 aj = _mesa_half_to_float(rowA[j][comp]);
374 ak = _mesa_half_to_float(rowA[k][comp]);
375 bj = _mesa_half_to_float(rowB[j][comp]);
376 bk = _mesa_half_to_float(rowB[k][comp]);
377 dst[i][comp] = _mesa_float_to_half((aj + ak + bj + bk) * 0.25F);
378 }
379 }
380 }
381 else if (datatype == GL_HALF_FLOAT_ARB && comps == 2) {
382 GLuint i, j, k, comp;
383 const GLhalfARB(*rowA)[2] = (const GLhalfARB(*)[2]) srcRowA;
384 const GLhalfARB(*rowB)[2] = (const GLhalfARB(*)[2]) srcRowB;
385 GLhalfARB(*dst)[2] = (GLhalfARB(*)[2]) dstRow;
386 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
387 i++, j += colStride, k += colStride) {
388 for (comp = 0; comp < 2; comp++) {
389 GLfloat aj, ak, bj, bk;
390 aj = _mesa_half_to_float(rowA[j][comp]);
391 ak = _mesa_half_to_float(rowA[k][comp]);
392 bj = _mesa_half_to_float(rowB[j][comp]);
393 bk = _mesa_half_to_float(rowB[k][comp]);
394 dst[i][comp] = _mesa_float_to_half((aj + ak + bj + bk) * 0.25F);
395 }
396 }
397 }
398 else if (datatype == GL_HALF_FLOAT_ARB && comps == 1) {
399 GLuint i, j, k;
400 const GLhalfARB *rowA = (const GLhalfARB *) srcRowA;
401 const GLhalfARB *rowB = (const GLhalfARB *) srcRowB;
402 GLhalfARB *dst = (GLhalfARB *) dstRow;
403 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
404 i++, j += colStride, k += colStride) {
405 GLfloat aj, ak, bj, bk;
406 aj = _mesa_half_to_float(rowA[j]);
407 ak = _mesa_half_to_float(rowA[k]);
408 bj = _mesa_half_to_float(rowB[j]);
409 bk = _mesa_half_to_float(rowB[k]);
410 dst[i] = _mesa_float_to_half((aj + ak + bj + bk) * 0.25F);
411 }
412 }
413
414 else if (datatype == GL_UNSIGNED_INT && comps == 1) {
415 GLuint i, j, k;
416 const GLuint *rowA = (const GLuint *) srcRowA;
417 const GLuint *rowB = (const GLuint *) srcRowB;
418 GLfloat *dst = (GLfloat *) dstRow;
419 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
420 i++, j += colStride, k += colStride) {
421 dst[i] = (GLfloat)(rowA[j] / 4 + rowA[k] / 4 + rowB[j] / 4 + rowB[k] / 4);
422 }
423 }
424
425 else if (datatype == GL_UNSIGNED_SHORT_5_6_5 && comps == 3) {
426 GLuint i, j, k;
427 const GLushort *rowA = (const GLushort *) srcRowA;
428 const GLushort *rowB = (const GLushort *) srcRowB;
429 GLushort *dst = (GLushort *) dstRow;
430 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
431 i++, j += colStride, k += colStride) {
432 const GLint rowAr0 = rowA[j] & 0x1f;
433 const GLint rowAr1 = rowA[k] & 0x1f;
434 const GLint rowBr0 = rowB[j] & 0x1f;
435 const GLint rowBr1 = rowB[k] & 0x1f;
436 const GLint rowAg0 = (rowA[j] >> 5) & 0x3f;
437 const GLint rowAg1 = (rowA[k] >> 5) & 0x3f;
438 const GLint rowBg0 = (rowB[j] >> 5) & 0x3f;
439 const GLint rowBg1 = (rowB[k] >> 5) & 0x3f;
440 const GLint rowAb0 = (rowA[j] >> 11) & 0x1f;
441 const GLint rowAb1 = (rowA[k] >> 11) & 0x1f;
442 const GLint rowBb0 = (rowB[j] >> 11) & 0x1f;
443 const GLint rowBb1 = (rowB[k] >> 11) & 0x1f;
444 const GLint red = (rowAr0 + rowAr1 + rowBr0 + rowBr1) >> 2;
445 const GLint green = (rowAg0 + rowAg1 + rowBg0 + rowBg1) >> 2;
446 const GLint blue = (rowAb0 + rowAb1 + rowBb0 + rowBb1) >> 2;
447 dst[i] = (blue << 11) | (green << 5) | red;
448 }
449 }
450 else if (datatype == GL_UNSIGNED_SHORT_4_4_4_4 && comps == 4) {
451 GLuint i, j, k;
452 const GLushort *rowA = (const GLushort *) srcRowA;
453 const GLushort *rowB = (const GLushort *) srcRowB;
454 GLushort *dst = (GLushort *) dstRow;
455 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
456 i++, j += colStride, k += colStride) {
457 const GLint rowAr0 = rowA[j] & 0xf;
458 const GLint rowAr1 = rowA[k] & 0xf;
459 const GLint rowBr0 = rowB[j] & 0xf;
460 const GLint rowBr1 = rowB[k] & 0xf;
461 const GLint rowAg0 = (rowA[j] >> 4) & 0xf;
462 const GLint rowAg1 = (rowA[k] >> 4) & 0xf;
463 const GLint rowBg0 = (rowB[j] >> 4) & 0xf;
464 const GLint rowBg1 = (rowB[k] >> 4) & 0xf;
465 const GLint rowAb0 = (rowA[j] >> 8) & 0xf;
466 const GLint rowAb1 = (rowA[k] >> 8) & 0xf;
467 const GLint rowBb0 = (rowB[j] >> 8) & 0xf;
468 const GLint rowBb1 = (rowB[k] >> 8) & 0xf;
469 const GLint rowAa0 = (rowA[j] >> 12) & 0xf;
470 const GLint rowAa1 = (rowA[k] >> 12) & 0xf;
471 const GLint rowBa0 = (rowB[j] >> 12) & 0xf;
472 const GLint rowBa1 = (rowB[k] >> 12) & 0xf;
473 const GLint red = (rowAr0 + rowAr1 + rowBr0 + rowBr1) >> 2;
474 const GLint green = (rowAg0 + rowAg1 + rowBg0 + rowBg1) >> 2;
475 const GLint blue = (rowAb0 + rowAb1 + rowBb0 + rowBb1) >> 2;
476 const GLint alpha = (rowAa0 + rowAa1 + rowBa0 + rowBa1) >> 2;
477 dst[i] = (alpha << 12) | (blue << 8) | (green << 4) | red;
478 }
479 }
480 else if (datatype == GL_UNSIGNED_SHORT_1_5_5_5_REV && comps == 4) {
481 GLuint i, j, k;
482 const GLushort *rowA = (const GLushort *) srcRowA;
483 const GLushort *rowB = (const GLushort *) srcRowB;
484 GLushort *dst = (GLushort *) dstRow;
485 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
486 i++, j += colStride, k += colStride) {
487 const GLint rowAr0 = rowA[j] & 0x1f;
488 const GLint rowAr1 = rowA[k] & 0x1f;
489 const GLint rowBr0 = rowB[j] & 0x1f;
490 const GLint rowBr1 = rowB[k] & 0x1f;
491 const GLint rowAg0 = (rowA[j] >> 5) & 0x1f;
492 const GLint rowAg1 = (rowA[k] >> 5) & 0x1f;
493 const GLint rowBg0 = (rowB[j] >> 5) & 0x1f;
494 const GLint rowBg1 = (rowB[k] >> 5) & 0x1f;
495 const GLint rowAb0 = (rowA[j] >> 10) & 0x1f;
496 const GLint rowAb1 = (rowA[k] >> 10) & 0x1f;
497 const GLint rowBb0 = (rowB[j] >> 10) & 0x1f;
498 const GLint rowBb1 = (rowB[k] >> 10) & 0x1f;
499 const GLint rowAa0 = (rowA[j] >> 15) & 0x1;
500 const GLint rowAa1 = (rowA[k] >> 15) & 0x1;
501 const GLint rowBa0 = (rowB[j] >> 15) & 0x1;
502 const GLint rowBa1 = (rowB[k] >> 15) & 0x1;
503 const GLint red = (rowAr0 + rowAr1 + rowBr0 + rowBr1) >> 2;
504 const GLint green = (rowAg0 + rowAg1 + rowBg0 + rowBg1) >> 2;
505 const GLint blue = (rowAb0 + rowAb1 + rowBb0 + rowBb1) >> 2;
506 const GLint alpha = (rowAa0 + rowAa1 + rowBa0 + rowBa1) >> 2;
507 dst[i] = (alpha << 15) | (blue << 10) | (green << 5) | red;
508 }
509 }
510 else if (datatype == GL_UNSIGNED_BYTE_3_3_2 && comps == 3) {
511 GLuint i, j, k;
512 const GLubyte *rowA = (const GLubyte *) srcRowA;
513 const GLubyte *rowB = (const GLubyte *) srcRowB;
514 GLubyte *dst = (GLubyte *) dstRow;
515 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
516 i++, j += colStride, k += colStride) {
517 const GLint rowAr0 = rowA[j] & 0x3;
518 const GLint rowAr1 = rowA[k] & 0x3;
519 const GLint rowBr0 = rowB[j] & 0x3;
520 const GLint rowBr1 = rowB[k] & 0x3;
521 const GLint rowAg0 = (rowA[j] >> 2) & 0x7;
522 const GLint rowAg1 = (rowA[k] >> 2) & 0x7;
523 const GLint rowBg0 = (rowB[j] >> 2) & 0x7;
524 const GLint rowBg1 = (rowB[k] >> 2) & 0x7;
525 const GLint rowAb0 = (rowA[j] >> 5) & 0x7;
526 const GLint rowAb1 = (rowA[k] >> 5) & 0x7;
527 const GLint rowBb0 = (rowB[j] >> 5) & 0x7;
528 const GLint rowBb1 = (rowB[k] >> 5) & 0x7;
529 const GLint red = (rowAr0 + rowAr1 + rowBr0 + rowBr1) >> 2;
530 const GLint green = (rowAg0 + rowAg1 + rowBg0 + rowBg1) >> 2;
531 const GLint blue = (rowAb0 + rowAb1 + rowBb0 + rowBb1) >> 2;
532 dst[i] = (blue << 5) | (green << 2) | red;
533 }
534 }
535 else {
536 _mesa_problem(NULL, "bad format in do_row()");
537 }
538 }
539
540
541 /**
542 * Average together four rows of a source image to produce a single new
543 * row in the dest image. It's legal for the two source rows to point
544 * to the same data. The source width must be equal to either the
545 * dest width or two times the dest width.
546 *
547 * \param datatype GL pixel type \c GL_UNSIGNED_BYTE, \c GL_UNSIGNED_SHORT,
548 * \c GL_FLOAT, etc.
549 * \param comps number of components per pixel (1..4)
550 * \param srcWidth Width of a row in the source data
551 * \param srcRowA Pointer to one of the rows of source data
552 * \param srcRowB Pointer to one of the rows of source data
553 * \param srcRowC Pointer to one of the rows of source data
554 * \param srcRowD Pointer to one of the rows of source data
555 * \param dstWidth Width of a row in the destination data
556 * \param srcRowA Pointer to the row of destination data
557 */
558 static void
559 do_row_3D(GLenum datatype, GLuint comps, GLint srcWidth,
560 const GLvoid *srcRowA, const GLvoid *srcRowB,
561 const GLvoid *srcRowC, const GLvoid *srcRowD,
562 GLint dstWidth, GLvoid *dstRow)
563 {
564 const GLuint k0 = (srcWidth == dstWidth) ? 0 : 1;
565 const GLuint colStride = (srcWidth == dstWidth) ? 1 : 2;
566 GLuint i, j, k;
567
568 ASSERT(comps >= 1);
569 ASSERT(comps <= 4);
570
571 if ((datatype == GL_UNSIGNED_BYTE) && (comps == 4)) {
572 DECLARE_ROW_POINTERS(GLubyte, 4);
573
574 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
575 i++, j += colStride, k += colStride) {
576 FILTER_3D(0);
577 FILTER_3D(1);
578 FILTER_3D(2);
579 FILTER_3D(3);
580 }
581 }
582 else if ((datatype == GL_UNSIGNED_BYTE) && (comps == 3)) {
583 DECLARE_ROW_POINTERS(GLubyte, 3);
584
585 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
586 i++, j += colStride, k += colStride) {
587 FILTER_3D(0);
588 FILTER_3D(1);
589 FILTER_3D(2);
590 }
591 }
592 else if ((datatype == GL_UNSIGNED_BYTE) && (comps == 2)) {
593 DECLARE_ROW_POINTERS(GLubyte, 2);
594
595 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
596 i++, j += colStride, k += colStride) {
597 FILTER_3D(0);
598 FILTER_3D(1);
599 }
600 }
601 else if ((datatype == GL_UNSIGNED_BYTE) && (comps == 1)) {
602 DECLARE_ROW_POINTERS(GLubyte, 1);
603
604 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
605 i++, j += colStride, k += colStride) {
606 FILTER_3D(0);
607 }
608 }
609 if ((datatype == GL_BYTE) && (comps == 4)) {
610 DECLARE_ROW_POINTERS(GLbyte, 4);
611
612 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
613 i++, j += colStride, k += colStride) {
614 FILTER_3D_SIGNED(0);
615 FILTER_3D_SIGNED(1);
616 FILTER_3D_SIGNED(2);
617 FILTER_3D_SIGNED(3);
618 }
619 }
620 else if ((datatype == GL_BYTE) && (comps == 3)) {
621 DECLARE_ROW_POINTERS(GLbyte, 3);
622
623 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
624 i++, j += colStride, k += colStride) {
625 FILTER_3D_SIGNED(0);
626 FILTER_3D_SIGNED(1);
627 FILTER_3D_SIGNED(2);
628 }
629 }
630 else if ((datatype == GL_BYTE) && (comps == 2)) {
631 DECLARE_ROW_POINTERS(GLbyte, 2);
632
633 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
634 i++, j += colStride, k += colStride) {
635 FILTER_3D_SIGNED(0);
636 FILTER_3D_SIGNED(1);
637 }
638 }
639 else if ((datatype == GL_BYTE) && (comps == 1)) {
640 DECLARE_ROW_POINTERS(GLbyte, 1);
641
642 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
643 i++, j += colStride, k += colStride) {
644 FILTER_3D_SIGNED(0);
645 }
646 }
647 else if ((datatype == GL_UNSIGNED_SHORT) && (comps == 4)) {
648 DECLARE_ROW_POINTERS(GLushort, 4);
649
650 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
651 i++, j += colStride, k += colStride) {
652 FILTER_3D(0);
653 FILTER_3D(1);
654 FILTER_3D(2);
655 FILTER_3D(3);
656 }
657 }
658 else if ((datatype == GL_UNSIGNED_SHORT) && (comps == 3)) {
659 DECLARE_ROW_POINTERS(GLushort, 3);
660
661 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
662 i++, j += colStride, k += colStride) {
663 FILTER_3D(0);
664 FILTER_3D(1);
665 FILTER_3D(2);
666 }
667 }
668 else if ((datatype == GL_UNSIGNED_SHORT) && (comps == 2)) {
669 DECLARE_ROW_POINTERS(GLushort, 2);
670
671 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
672 i++, j += colStride, k += colStride) {
673 FILTER_3D(0);
674 FILTER_3D(1);
675 }
676 }
677 else if ((datatype == GL_UNSIGNED_SHORT) && (comps == 1)) {
678 DECLARE_ROW_POINTERS(GLushort, 1);
679
680 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
681 i++, j += colStride, k += colStride) {
682 FILTER_3D(0);
683 }
684 }
685 else if ((datatype == GL_FLOAT) && (comps == 4)) {
686 DECLARE_ROW_POINTERS(GLfloat, 4);
687
688 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
689 i++, j += colStride, k += colStride) {
690 FILTER_F_3D(0);
691 FILTER_F_3D(1);
692 FILTER_F_3D(2);
693 FILTER_F_3D(3);
694 }
695 }
696 else if ((datatype == GL_FLOAT) && (comps == 3)) {
697 DECLARE_ROW_POINTERS(GLfloat, 3);
698
699 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
700 i++, j += colStride, k += colStride) {
701 FILTER_F_3D(0);
702 FILTER_F_3D(1);
703 FILTER_F_3D(2);
704 }
705 }
706 else if ((datatype == GL_FLOAT) && (comps == 2)) {
707 DECLARE_ROW_POINTERS(GLfloat, 2);
708
709 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
710 i++, j += colStride, k += colStride) {
711 FILTER_F_3D(0);
712 FILTER_F_3D(1);
713 }
714 }
715 else if ((datatype == GL_FLOAT) && (comps == 1)) {
716 DECLARE_ROW_POINTERS(GLfloat, 1);
717
718 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
719 i++, j += colStride, k += colStride) {
720 FILTER_F_3D(0);
721 }
722 }
723 else if ((datatype == GL_HALF_FLOAT_ARB) && (comps == 4)) {
724 DECLARE_ROW_POINTERS(GLhalfARB, 4);
725
726 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
727 i++, j += colStride, k += colStride) {
728 FILTER_HF_3D(0);
729 FILTER_HF_3D(1);
730 FILTER_HF_3D(2);
731 FILTER_HF_3D(3);
732 }
733 }
734 else if ((datatype == GL_HALF_FLOAT_ARB) && (comps == 3)) {
735 DECLARE_ROW_POINTERS(GLhalfARB, 4);
736
737 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
738 i++, j += colStride, k += colStride) {
739 FILTER_HF_3D(0);
740 FILTER_HF_3D(1);
741 FILTER_HF_3D(2);
742 }
743 }
744 else if ((datatype == GL_HALF_FLOAT_ARB) && (comps == 2)) {
745 DECLARE_ROW_POINTERS(GLhalfARB, 4);
746
747 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
748 i++, j += colStride, k += colStride) {
749 FILTER_HF_3D(0);
750 FILTER_HF_3D(1);
751 }
752 }
753 else if ((datatype == GL_HALF_FLOAT_ARB) && (comps == 1)) {
754 DECLARE_ROW_POINTERS(GLhalfARB, 4);
755
756 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
757 i++, j += colStride, k += colStride) {
758 FILTER_HF_3D(0);
759 }
760 }
761 else if ((datatype == GL_UNSIGNED_INT) && (comps == 1)) {
762 const GLuint *rowA = (const GLuint *) srcRowA;
763 const GLuint *rowB = (const GLuint *) srcRowB;
764 const GLuint *rowC = (const GLuint *) srcRowC;
765 const GLuint *rowD = (const GLuint *) srcRowD;
766 GLfloat *dst = (GLfloat *) dstRow;
767
768 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
769 i++, j += colStride, k += colStride) {
770 const uint64_t tmp = (((uint64_t) rowA[j] + (uint64_t) rowA[k])
771 + ((uint64_t) rowB[j] + (uint64_t) rowB[k])
772 + ((uint64_t) rowC[j] + (uint64_t) rowC[k])
773 + ((uint64_t) rowD[j] + (uint64_t) rowD[k]));
774 dst[i] = (GLfloat)((double) tmp * 0.125);
775 }
776 }
777 else if ((datatype == GL_UNSIGNED_SHORT_5_6_5) && (comps == 3)) {
778 DECLARE_ROW_POINTERS0(GLushort);
779
780 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
781 i++, j += colStride, k += colStride) {
782 const GLint rowAr0 = rowA[j] & 0x1f;
783 const GLint rowAr1 = rowA[k] & 0x1f;
784 const GLint rowBr0 = rowB[j] & 0x1f;
785 const GLint rowBr1 = rowB[k] & 0x1f;
786 const GLint rowCr0 = rowC[j] & 0x1f;
787 const GLint rowCr1 = rowC[k] & 0x1f;
788 const GLint rowDr0 = rowD[j] & 0x1f;
789 const GLint rowDr1 = rowD[k] & 0x1f;
790 const GLint rowAg0 = (rowA[j] >> 5) & 0x3f;
791 const GLint rowAg1 = (rowA[k] >> 5) & 0x3f;
792 const GLint rowBg0 = (rowB[j] >> 5) & 0x3f;
793 const GLint rowBg1 = (rowB[k] >> 5) & 0x3f;
794 const GLint rowCg0 = (rowC[j] >> 5) & 0x3f;
795 const GLint rowCg1 = (rowC[k] >> 5) & 0x3f;
796 const GLint rowDg0 = (rowD[j] >> 5) & 0x3f;
797 const GLint rowDg1 = (rowD[k] >> 5) & 0x3f;
798 const GLint rowAb0 = (rowA[j] >> 11) & 0x1f;
799 const GLint rowAb1 = (rowA[k] >> 11) & 0x1f;
800 const GLint rowBb0 = (rowB[j] >> 11) & 0x1f;
801 const GLint rowBb1 = (rowB[k] >> 11) & 0x1f;
802 const GLint rowCb0 = (rowC[j] >> 11) & 0x1f;
803 const GLint rowCb1 = (rowC[k] >> 11) & 0x1f;
804 const GLint rowDb0 = (rowD[j] >> 11) & 0x1f;
805 const GLint rowDb1 = (rowD[k] >> 11) & 0x1f;
806 const GLint r = FILTER_SUM_3D(rowAr0, rowAr1, rowBr0, rowBr1,
807 rowCr0, rowCr1, rowDr0, rowDr1);
808 const GLint g = FILTER_SUM_3D(rowAg0, rowAg1, rowBg0, rowBg1,
809 rowCg0, rowCg1, rowDg0, rowDg1);
810 const GLint b = FILTER_SUM_3D(rowAb0, rowAb1, rowBb0, rowBb1,
811 rowCb0, rowCb1, rowDb0, rowDb1);
812 dst[i] = (b << 11) | (g << 5) | r;
813 }
814 }
815 else if ((datatype == GL_UNSIGNED_SHORT_4_4_4_4) && (comps == 4)) {
816 DECLARE_ROW_POINTERS0(GLushort);
817
818 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
819 i++, j += colStride, k += colStride) {
820 const GLint rowAr0 = rowA[j] & 0xf;
821 const GLint rowAr1 = rowA[k] & 0xf;
822 const GLint rowBr0 = rowB[j] & 0xf;
823 const GLint rowBr1 = rowB[k] & 0xf;
824 const GLint rowCr0 = rowC[j] & 0xf;
825 const GLint rowCr1 = rowC[k] & 0xf;
826 const GLint rowDr0 = rowD[j] & 0xf;
827 const GLint rowDr1 = rowD[k] & 0xf;
828 const GLint rowAg0 = (rowA[j] >> 4) & 0xf;
829 const GLint rowAg1 = (rowA[k] >> 4) & 0xf;
830 const GLint rowBg0 = (rowB[j] >> 4) & 0xf;
831 const GLint rowBg1 = (rowB[k] >> 4) & 0xf;
832 const GLint rowCg0 = (rowC[j] >> 4) & 0xf;
833 const GLint rowCg1 = (rowC[k] >> 4) & 0xf;
834 const GLint rowDg0 = (rowD[j] >> 4) & 0xf;
835 const GLint rowDg1 = (rowD[k] >> 4) & 0xf;
836 const GLint rowAb0 = (rowA[j] >> 8) & 0xf;
837 const GLint rowAb1 = (rowA[k] >> 8) & 0xf;
838 const GLint rowBb0 = (rowB[j] >> 8) & 0xf;
839 const GLint rowBb1 = (rowB[k] >> 8) & 0xf;
840 const GLint rowCb0 = (rowC[j] >> 8) & 0xf;
841 const GLint rowCb1 = (rowC[k] >> 8) & 0xf;
842 const GLint rowDb0 = (rowD[j] >> 8) & 0xf;
843 const GLint rowDb1 = (rowD[k] >> 8) & 0xf;
844 const GLint rowAa0 = (rowA[j] >> 12) & 0xf;
845 const GLint rowAa1 = (rowA[k] >> 12) & 0xf;
846 const GLint rowBa0 = (rowB[j] >> 12) & 0xf;
847 const GLint rowBa1 = (rowB[k] >> 12) & 0xf;
848 const GLint rowCa0 = (rowC[j] >> 12) & 0xf;
849 const GLint rowCa1 = (rowC[k] >> 12) & 0xf;
850 const GLint rowDa0 = (rowD[j] >> 12) & 0xf;
851 const GLint rowDa1 = (rowD[k] >> 12) & 0xf;
852 const GLint r = FILTER_SUM_3D(rowAr0, rowAr1, rowBr0, rowBr1,
853 rowCr0, rowCr1, rowDr0, rowDr1);
854 const GLint g = FILTER_SUM_3D(rowAg0, rowAg1, rowBg0, rowBg1,
855 rowCg0, rowCg1, rowDg0, rowDg1);
856 const GLint b = FILTER_SUM_3D(rowAb0, rowAb1, rowBb0, rowBb1,
857 rowCb0, rowCb1, rowDb0, rowDb1);
858 const GLint a = FILTER_SUM_3D(rowAa0, rowAa1, rowBa0, rowBa1,
859 rowCa0, rowCa1, rowDa0, rowDa1);
860
861 dst[i] = (a << 12) | (b << 8) | (g << 4) | r;
862 }
863 }
864 else if ((datatype == GL_UNSIGNED_SHORT_1_5_5_5_REV) && (comps == 4)) {
865 DECLARE_ROW_POINTERS0(GLushort);
866
867 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
868 i++, j += colStride, k += colStride) {
869 const GLint rowAr0 = rowA[j] & 0x1f;
870 const GLint rowAr1 = rowA[k] & 0x1f;
871 const GLint rowBr0 = rowB[j] & 0x1f;
872 const GLint rowBr1 = rowB[k] & 0x1f;
873 const GLint rowCr0 = rowC[j] & 0x1f;
874 const GLint rowCr1 = rowC[k] & 0x1f;
875 const GLint rowDr0 = rowD[j] & 0x1f;
876 const GLint rowDr1 = rowD[k] & 0x1f;
877 const GLint rowAg0 = (rowA[j] >> 5) & 0x1f;
878 const GLint rowAg1 = (rowA[k] >> 5) & 0x1f;
879 const GLint rowBg0 = (rowB[j] >> 5) & 0x1f;
880 const GLint rowBg1 = (rowB[k] >> 5) & 0x1f;
881 const GLint rowCg0 = (rowC[j] >> 5) & 0x1f;
882 const GLint rowCg1 = (rowC[k] >> 5) & 0x1f;
883 const GLint rowDg0 = (rowD[j] >> 5) & 0x1f;
884 const GLint rowDg1 = (rowD[k] >> 5) & 0x1f;
885 const GLint rowAb0 = (rowA[j] >> 10) & 0x1f;
886 const GLint rowAb1 = (rowA[k] >> 10) & 0x1f;
887 const GLint rowBb0 = (rowB[j] >> 10) & 0x1f;
888 const GLint rowBb1 = (rowB[k] >> 10) & 0x1f;
889 const GLint rowCb0 = (rowC[j] >> 10) & 0x1f;
890 const GLint rowCb1 = (rowC[k] >> 10) & 0x1f;
891 const GLint rowDb0 = (rowD[j] >> 10) & 0x1f;
892 const GLint rowDb1 = (rowD[k] >> 10) & 0x1f;
893 const GLint rowAa0 = (rowA[j] >> 15) & 0x1;
894 const GLint rowAa1 = (rowA[k] >> 15) & 0x1;
895 const GLint rowBa0 = (rowB[j] >> 15) & 0x1;
896 const GLint rowBa1 = (rowB[k] >> 15) & 0x1;
897 const GLint rowCa0 = (rowC[j] >> 15) & 0x1;
898 const GLint rowCa1 = (rowC[k] >> 15) & 0x1;
899 const GLint rowDa0 = (rowD[j] >> 15) & 0x1;
900 const GLint rowDa1 = (rowD[k] >> 15) & 0x1;
901 const GLint r = FILTER_SUM_3D(rowAr0, rowAr1, rowBr0, rowBr1,
902 rowCr0, rowCr1, rowDr0, rowDr1);
903 const GLint g = FILTER_SUM_3D(rowAg0, rowAg1, rowBg0, rowBg1,
904 rowCg0, rowCg1, rowDg0, rowDg1);
905 const GLint b = FILTER_SUM_3D(rowAb0, rowAb1, rowBb0, rowBb1,
906 rowCb0, rowCb1, rowDb0, rowDb1);
907 const GLint a = FILTER_SUM_3D(rowAa0, rowAa1, rowBa0, rowBa1,
908 rowCa0, rowCa1, rowDa0, rowDa1);
909
910 dst[i] = (a << 15) | (b << 10) | (g << 5) | r;
911 }
912 }
913 else if ((datatype == GL_UNSIGNED_BYTE_3_3_2) && (comps == 3)) {
914 DECLARE_ROW_POINTERS0(GLushort);
915
916 for (i = j = 0, k = k0; i < (GLuint) dstWidth;
917 i++, j += colStride, k += colStride) {
918 const GLint rowAr0 = rowA[j] & 0x3;
919 const GLint rowAr1 = rowA[k] & 0x3;
920 const GLint rowBr0 = rowB[j] & 0x3;
921 const GLint rowBr1 = rowB[k] & 0x3;
922 const GLint rowCr0 = rowC[j] & 0x3;
923 const GLint rowCr1 = rowC[k] & 0x3;
924 const GLint rowDr0 = rowD[j] & 0x3;
925 const GLint rowDr1 = rowD[k] & 0x3;
926 const GLint rowAg0 = (rowA[j] >> 2) & 0x7;
927 const GLint rowAg1 = (rowA[k] >> 2) & 0x7;
928 const GLint rowBg0 = (rowB[j] >> 2) & 0x7;
929 const GLint rowBg1 = (rowB[k] >> 2) & 0x7;
930 const GLint rowCg0 = (rowC[j] >> 2) & 0x7;
931 const GLint rowCg1 = (rowC[k] >> 2) & 0x7;
932 const GLint rowDg0 = (rowD[j] >> 2) & 0x7;
933 const GLint rowDg1 = (rowD[k] >> 2) & 0x7;
934 const GLint rowAb0 = (rowA[j] >> 5) & 0x7;
935 const GLint rowAb1 = (rowA[k] >> 5) & 0x7;
936 const GLint rowBb0 = (rowB[j] >> 5) & 0x7;
937 const GLint rowBb1 = (rowB[k] >> 5) & 0x7;
938 const GLint rowCb0 = (rowC[j] >> 5) & 0x7;
939 const GLint rowCb1 = (rowC[k] >> 5) & 0x7;
940 const GLint rowDb0 = (rowD[j] >> 5) & 0x7;
941 const GLint rowDb1 = (rowD[k] >> 5) & 0x7;
942 const GLint r = FILTER_SUM_3D(rowAr0, rowAr1, rowBr0, rowBr1,
943 rowCr0, rowCr1, rowDr0, rowDr1);
944 const GLint g = FILTER_SUM_3D(rowAg0, rowAg1, rowBg0, rowBg1,
945 rowCg0, rowCg1, rowDg0, rowDg1);
946 const GLint b = FILTER_SUM_3D(rowAb0, rowAb1, rowBb0, rowBb1,
947 rowCb0, rowCb1, rowDb0, rowDb1);
948 dst[i] = (b << 5) | (g << 2) | r;
949 }
950 }
951 else {
952 _mesa_problem(NULL, "bad format in do_row()");
953 }
954 }
955
956
957 /*
958 * These functions generate a 1/2-size mipmap image from a source image.
959 * Texture borders are handled by copying or averaging the source image's
960 * border texels, depending on the scale-down factor.
961 */
962
963 static void
964 make_1d_mipmap(GLenum datatype, GLuint comps, GLint border,
965 GLint srcWidth, const GLubyte *srcPtr,
966 GLint dstWidth, GLubyte *dstPtr)
967 {
968 const GLint bpt = bytes_per_pixel(datatype, comps);
969 const GLubyte *src;
970 GLubyte *dst;
971
972 /* skip the border pixel, if any */
973 src = srcPtr + border * bpt;
974 dst = dstPtr + border * bpt;
975
976 /* we just duplicate the input row, kind of hack, saves code */
977 do_row(datatype, comps, srcWidth - 2 * border, src, src,
978 dstWidth - 2 * border, dst);
979
980 if (border) {
981 /* copy left-most pixel from source */
982 assert(dstPtr);
983 assert(srcPtr);
984 memcpy(dstPtr, srcPtr, bpt);
985 /* copy right-most pixel from source */
986 memcpy(dstPtr + (dstWidth - 1) * bpt,
987 srcPtr + (srcWidth - 1) * bpt,
988 bpt);
989 }
990 }
991
992
993 static void
994 make_2d_mipmap(GLenum datatype, GLuint comps, GLint border,
995 GLint srcWidth, GLint srcHeight,
996 const GLubyte *srcPtr, GLint srcRowStride,
997 GLint dstWidth, GLint dstHeight,
998 GLubyte *dstPtr, GLint dstRowStride)
999 {
1000 const GLint bpt = bytes_per_pixel(datatype, comps);
1001 const GLint srcWidthNB = srcWidth - 2 * border; /* sizes w/out border */
1002 const GLint dstWidthNB = dstWidth - 2 * border;
1003 const GLint dstHeightNB = dstHeight - 2 * border;
1004 const GLint srcRowBytes = bpt * srcRowStride;
1005 const GLint dstRowBytes = bpt * dstRowStride;
1006 const GLubyte *srcA, *srcB;
1007 GLubyte *dst;
1008 GLint row, srcRowStep;
1009
1010 /* Compute src and dst pointers, skipping any border */
1011 srcA = srcPtr + border * ((srcWidth + 1) * bpt);
1012 if (srcHeight > 1 && srcHeight > dstHeight) {
1013 /* sample from two source rows */
1014 srcB = srcA + srcRowBytes;
1015 srcRowStep = 2;
1016 }
1017 else {
1018 /* sample from one source row */
1019 srcB = srcA;
1020 srcRowStep = 1;
1021 }
1022
1023 dst = dstPtr + border * ((dstWidth + 1) * bpt);
1024
1025 for (row = 0; row < dstHeightNB; row++) {
1026 do_row(datatype, comps, srcWidthNB, srcA, srcB,
1027 dstWidthNB, dst);
1028 srcA += srcRowStep * srcRowBytes;
1029 srcB += srcRowStep * srcRowBytes;
1030 dst += dstRowBytes;
1031 }
1032
1033 /* This is ugly but probably won't be used much */
1034 if (border > 0) {
1035 /* fill in dest border */
1036 /* lower-left border pixel */
1037 assert(dstPtr);
1038 assert(srcPtr);
1039 memcpy(dstPtr, srcPtr, bpt);
1040 /* lower-right border pixel */
1041 memcpy(dstPtr + (dstWidth - 1) * bpt,
1042 srcPtr + (srcWidth - 1) * bpt, bpt);
1043 /* upper-left border pixel */
1044 memcpy(dstPtr + dstWidth * (dstHeight - 1) * bpt,
1045 srcPtr + srcWidth * (srcHeight - 1) * bpt, bpt);
1046 /* upper-right border pixel */
1047 memcpy(dstPtr + (dstWidth * dstHeight - 1) * bpt,
1048 srcPtr + (srcWidth * srcHeight - 1) * bpt, bpt);
1049 /* lower border */
1050 do_row(datatype, comps, srcWidthNB,
1051 srcPtr + bpt,
1052 srcPtr + bpt,
1053 dstWidthNB, dstPtr + bpt);
1054 /* upper border */
1055 do_row(datatype, comps, srcWidthNB,
1056 srcPtr + (srcWidth * (srcHeight - 1) + 1) * bpt,
1057 srcPtr + (srcWidth * (srcHeight - 1) + 1) * bpt,
1058 dstWidthNB,
1059 dstPtr + (dstWidth * (dstHeight - 1) + 1) * bpt);
1060 /* left and right borders */
1061 if (srcHeight == dstHeight) {
1062 /* copy border pixel from src to dst */
1063 for (row = 1; row < srcHeight; row++) {
1064 memcpy(dstPtr + dstWidth * row * bpt,
1065 srcPtr + srcWidth * row * bpt, bpt);
1066 memcpy(dstPtr + (dstWidth * row + dstWidth - 1) * bpt,
1067 srcPtr + (srcWidth * row + srcWidth - 1) * bpt, bpt);
1068 }
1069 }
1070 else {
1071 /* average two src pixels each dest pixel */
1072 for (row = 0; row < dstHeightNB; row += 2) {
1073 do_row(datatype, comps, 1,
1074 srcPtr + (srcWidth * (row * 2 + 1)) * bpt,
1075 srcPtr + (srcWidth * (row * 2 + 2)) * bpt,
1076 1, dstPtr + (dstWidth * row + 1) * bpt);
1077 do_row(datatype, comps, 1,
1078 srcPtr + (srcWidth * (row * 2 + 1) + srcWidth - 1) * bpt,
1079 srcPtr + (srcWidth * (row * 2 + 2) + srcWidth - 1) * bpt,
1080 1, dstPtr + (dstWidth * row + 1 + dstWidth - 1) * bpt);
1081 }
1082 }
1083 }
1084 }
1085
1086
1087 static void
1088 make_3d_mipmap(GLenum datatype, GLuint comps, GLint border,
1089 GLint srcWidth, GLint srcHeight, GLint srcDepth,
1090 const GLubyte *srcPtr, GLint srcRowStride,
1091 GLint dstWidth, GLint dstHeight, GLint dstDepth,
1092 GLubyte *dstPtr, GLint dstRowStride)
1093 {
1094 const GLint bpt = bytes_per_pixel(datatype, comps);
1095 const GLint srcWidthNB = srcWidth - 2 * border; /* sizes w/out border */
1096 const GLint srcDepthNB = srcDepth - 2 * border;
1097 const GLint dstWidthNB = dstWidth - 2 * border;
1098 const GLint dstHeightNB = dstHeight - 2 * border;
1099 const GLint dstDepthNB = dstDepth - 2 * border;
1100 GLint img, row;
1101 GLint bytesPerSrcImage, bytesPerDstImage;
1102 GLint bytesPerSrcRow, bytesPerDstRow;
1103 GLint srcImageOffset, srcRowOffset;
1104
1105 (void) srcDepthNB; /* silence warnings */
1106
1107
1108 bytesPerSrcImage = srcWidth * srcHeight * bpt;
1109 bytesPerDstImage = dstWidth * dstHeight * bpt;
1110
1111 bytesPerSrcRow = srcWidth * bpt;
1112 bytesPerDstRow = dstWidth * bpt;
1113
1114 /* Offset between adjacent src images to be averaged together */
1115 srcImageOffset = (srcDepth == dstDepth) ? 0 : bytesPerSrcImage;
1116
1117 /* Offset between adjacent src rows to be averaged together */
1118 srcRowOffset = (srcHeight == dstHeight) ? 0 : srcWidth * bpt;
1119
1120 /*
1121 * Need to average together up to 8 src pixels for each dest pixel.
1122 * Break that down into 3 operations:
1123 * 1. take two rows from source image and average them together.
1124 * 2. take two rows from next source image and average them together.
1125 * 3. take the two averaged rows and average them for the final dst row.
1126 */
1127
1128 /*
1129 printf("mip3d %d x %d x %d -> %d x %d x %d\n",
1130 srcWidth, srcHeight, srcDepth, dstWidth, dstHeight, dstDepth);
1131 */
1132
1133 for (img = 0; img < dstDepthNB; img++) {
1134 /* first source image pointer, skipping border */
1135 const GLubyte *imgSrcA = srcPtr
1136 + (bytesPerSrcImage + bytesPerSrcRow + border) * bpt * border
1137 + img * (bytesPerSrcImage + srcImageOffset);
1138 /* second source image pointer, skipping border */
1139 const GLubyte *imgSrcB = imgSrcA + srcImageOffset;
1140 /* address of the dest image, skipping border */
1141 GLubyte *imgDst = dstPtr
1142 + (bytesPerDstImage + bytesPerDstRow + border) * bpt * border
1143 + img * bytesPerDstImage;
1144
1145 /* setup the four source row pointers and the dest row pointer */
1146 const GLubyte *srcImgARowA = imgSrcA;
1147 const GLubyte *srcImgARowB = imgSrcA + srcRowOffset;
1148 const GLubyte *srcImgBRowA = imgSrcB;
1149 const GLubyte *srcImgBRowB = imgSrcB + srcRowOffset;
1150 GLubyte *dstImgRow = imgDst;
1151
1152 for (row = 0; row < dstHeightNB; row++) {
1153 do_row_3D(datatype, comps, srcWidthNB,
1154 srcImgARowA, srcImgARowB,
1155 srcImgBRowA, srcImgBRowB,
1156 dstWidthNB, dstImgRow);
1157
1158 /* advance to next rows */
1159 srcImgARowA += bytesPerSrcRow + srcRowOffset;
1160 srcImgARowB += bytesPerSrcRow + srcRowOffset;
1161 srcImgBRowA += bytesPerSrcRow + srcRowOffset;
1162 srcImgBRowB += bytesPerSrcRow + srcRowOffset;
1163 dstImgRow += bytesPerDstRow;
1164 }
1165 }
1166
1167
1168 /* Luckily we can leverage the make_2d_mipmap() function here! */
1169 if (border > 0) {
1170 /* do front border image */
1171 make_2d_mipmap(datatype, comps, 1, srcWidth, srcHeight, srcPtr, srcRowStride,
1172 dstWidth, dstHeight, dstPtr, dstRowStride);
1173 /* do back border image */
1174 make_2d_mipmap(datatype, comps, 1, srcWidth, srcHeight,
1175 srcPtr + bytesPerSrcImage * (srcDepth - 1), srcRowStride,
1176 dstWidth, dstHeight,
1177 dstPtr + bytesPerDstImage * (dstDepth - 1), dstRowStride);
1178 /* do four remaining border edges that span the image slices */
1179 if (srcDepth == dstDepth) {
1180 /* just copy border pixels from src to dst */
1181 for (img = 0; img < dstDepthNB; img++) {
1182 const GLubyte *src;
1183 GLubyte *dst;
1184
1185 /* do border along [img][row=0][col=0] */
1186 src = srcPtr + (img + 1) * bytesPerSrcImage;
1187 dst = dstPtr + (img + 1) * bytesPerDstImage;
1188 memcpy(dst, src, bpt);
1189
1190 /* do border along [img][row=dstHeight-1][col=0] */
1191 src = srcPtr + (img * 2 + 1) * bytesPerSrcImage
1192 + (srcHeight - 1) * bytesPerSrcRow;
1193 dst = dstPtr + (img + 1) * bytesPerDstImage
1194 + (dstHeight - 1) * bytesPerDstRow;
1195 memcpy(dst, src, bpt);
1196
1197 /* do border along [img][row=0][col=dstWidth-1] */
1198 src = srcPtr + (img * 2 + 1) * bytesPerSrcImage
1199 + (srcWidth - 1) * bpt;
1200 dst = dstPtr + (img + 1) * bytesPerDstImage
1201 + (dstWidth - 1) * bpt;
1202 memcpy(dst, src, bpt);
1203
1204 /* do border along [img][row=dstHeight-1][col=dstWidth-1] */
1205 src = srcPtr + (img * 2 + 1) * bytesPerSrcImage
1206 + (bytesPerSrcImage - bpt);
1207 dst = dstPtr + (img + 1) * bytesPerDstImage
1208 + (bytesPerDstImage - bpt);
1209 memcpy(dst, src, bpt);
1210 }
1211 }
1212 else {
1213 /* average border pixels from adjacent src image pairs */
1214 ASSERT(srcDepthNB == 2 * dstDepthNB);
1215 for (img = 0; img < dstDepthNB; img++) {
1216 const GLubyte *src;
1217 GLubyte *dst;
1218
1219 /* do border along [img][row=0][col=0] */
1220 src = srcPtr + (img * 2 + 1) * bytesPerSrcImage;
1221 dst = dstPtr + (img + 1) * bytesPerDstImage;
1222 do_row(datatype, comps, 1, src, src + srcImageOffset, 1, dst);
1223
1224 /* do border along [img][row=dstHeight-1][col=0] */
1225 src = srcPtr + (img * 2 + 1) * bytesPerSrcImage
1226 + (srcHeight - 1) * bytesPerSrcRow;
1227 dst = dstPtr + (img + 1) * bytesPerDstImage
1228 + (dstHeight - 1) * bytesPerDstRow;
1229 do_row(datatype, comps, 1, src, src + srcImageOffset, 1, dst);
1230
1231 /* do border along [img][row=0][col=dstWidth-1] */
1232 src = srcPtr + (img * 2 + 1) * bytesPerSrcImage
1233 + (srcWidth - 1) * bpt;
1234 dst = dstPtr + (img + 1) * bytesPerDstImage
1235 + (dstWidth - 1) * bpt;
1236 do_row(datatype, comps, 1, src, src + srcImageOffset, 1, dst);
1237
1238 /* do border along [img][row=dstHeight-1][col=dstWidth-1] */
1239 src = srcPtr + (img * 2 + 1) * bytesPerSrcImage
1240 + (bytesPerSrcImage - bpt);
1241 dst = dstPtr + (img + 1) * bytesPerDstImage
1242 + (bytesPerDstImage - bpt);
1243 do_row(datatype, comps, 1, src, src + srcImageOffset, 1, dst);
1244 }
1245 }
1246 }
1247 }
1248
1249
1250 static void
1251 make_1d_stack_mipmap(GLenum datatype, GLuint comps, GLint border,
1252 GLint srcWidth, const GLubyte *srcPtr, GLuint srcRowStride,
1253 GLint dstWidth, GLint dstHeight,
1254 GLubyte *dstPtr, GLuint dstRowStride )
1255 {
1256 const GLint bpt = bytes_per_pixel(datatype, comps);
1257 const GLint srcWidthNB = srcWidth - 2 * border; /* sizes w/out border */
1258 const GLint dstWidthNB = dstWidth - 2 * border;
1259 const GLint dstHeightNB = dstHeight - 2 * border;
1260 const GLint srcRowBytes = bpt * srcRowStride;
1261 const GLint dstRowBytes = bpt * dstRowStride;
1262 const GLubyte *src;
1263 GLubyte *dst;
1264 GLint row;
1265
1266 /* Compute src and dst pointers, skipping any border */
1267 src = srcPtr + border * ((srcWidth + 1) * bpt);
1268 dst = dstPtr + border * ((dstWidth + 1) * bpt);
1269
1270 for (row = 0; row < dstHeightNB; row++) {
1271 do_row(datatype, comps, srcWidthNB, src, src,
1272 dstWidthNB, dst);
1273 src += srcRowBytes;
1274 dst += dstRowBytes;
1275 }
1276
1277 if (border) {
1278 /* copy left-most pixel from source */
1279 assert(dstPtr);
1280 assert(srcPtr);
1281 memcpy(dstPtr, srcPtr, bpt);
1282 /* copy right-most pixel from source */
1283 memcpy(dstPtr + (dstWidth - 1) * bpt,
1284 srcPtr + (srcWidth - 1) * bpt,
1285 bpt);
1286 }
1287 }
1288
1289
1290 /**
1291 * \bug
1292 * There is quite a bit of refactoring that could be done with this function
1293 * and \c make_2d_mipmap.
1294 */
1295 static void
1296 make_2d_stack_mipmap(GLenum datatype, GLuint comps, GLint border,
1297 GLint srcWidth, GLint srcHeight,
1298 const GLubyte *srcPtr, GLint srcRowStride,
1299 GLint dstWidth, GLint dstHeight, GLint dstDepth,
1300 GLubyte *dstPtr, GLint dstRowStride)
1301 {
1302 const GLint bpt = bytes_per_pixel(datatype, comps);
1303 const GLint srcWidthNB = srcWidth - 2 * border; /* sizes w/out border */
1304 const GLint dstWidthNB = dstWidth - 2 * border;
1305 const GLint dstHeightNB = dstHeight - 2 * border;
1306 const GLint dstDepthNB = dstDepth - 2 * border;
1307 const GLint srcRowBytes = bpt * srcRowStride;
1308 const GLint dstRowBytes = bpt * dstRowStride;
1309 const GLubyte *srcA, *srcB;
1310 GLubyte *dst;
1311 GLint layer;
1312 GLint row;
1313
1314 /* Compute src and dst pointers, skipping any border */
1315 srcA = srcPtr + border * ((srcWidth + 1) * bpt);
1316 if (srcHeight > 1)
1317 srcB = srcA + srcRowBytes;
1318 else
1319 srcB = srcA;
1320 dst = dstPtr + border * ((dstWidth + 1) * bpt);
1321
1322 for (layer = 0; layer < dstDepthNB; layer++) {
1323 for (row = 0; row < dstHeightNB; row++) {
1324 do_row(datatype, comps, srcWidthNB, srcA, srcB,
1325 dstWidthNB, dst);
1326 srcA += 2 * srcRowBytes;
1327 srcB += 2 * srcRowBytes;
1328 dst += dstRowBytes;
1329 }
1330
1331 /* This is ugly but probably won't be used much */
1332 if (border > 0) {
1333 /* fill in dest border */
1334 /* lower-left border pixel */
1335 assert(dstPtr);
1336 assert(srcPtr);
1337 memcpy(dstPtr, srcPtr, bpt);
1338 /* lower-right border pixel */
1339 memcpy(dstPtr + (dstWidth - 1) * bpt,
1340 srcPtr + (srcWidth - 1) * bpt, bpt);
1341 /* upper-left border pixel */
1342 memcpy(dstPtr + dstWidth * (dstHeight - 1) * bpt,
1343 srcPtr + srcWidth * (srcHeight - 1) * bpt, bpt);
1344 /* upper-right border pixel */
1345 memcpy(dstPtr + (dstWidth * dstHeight - 1) * bpt,
1346 srcPtr + (srcWidth * srcHeight - 1) * bpt, bpt);
1347 /* lower border */
1348 do_row(datatype, comps, srcWidthNB,
1349 srcPtr + bpt,
1350 srcPtr + bpt,
1351 dstWidthNB, dstPtr + bpt);
1352 /* upper border */
1353 do_row(datatype, comps, srcWidthNB,
1354 srcPtr + (srcWidth * (srcHeight - 1) + 1) * bpt,
1355 srcPtr + (srcWidth * (srcHeight - 1) + 1) * bpt,
1356 dstWidthNB,
1357 dstPtr + (dstWidth * (dstHeight - 1) + 1) * bpt);
1358 /* left and right borders */
1359 if (srcHeight == dstHeight) {
1360 /* copy border pixel from src to dst */
1361 for (row = 1; row < srcHeight; row++) {
1362 memcpy(dstPtr + dstWidth * row * bpt,
1363 srcPtr + srcWidth * row * bpt, bpt);
1364 memcpy(dstPtr + (dstWidth * row + dstWidth - 1) * bpt,
1365 srcPtr + (srcWidth * row + srcWidth - 1) * bpt, bpt);
1366 }
1367 }
1368 else {
1369 /* average two src pixels each dest pixel */
1370 for (row = 0; row < dstHeightNB; row += 2) {
1371 do_row(datatype, comps, 1,
1372 srcPtr + (srcWidth * (row * 2 + 1)) * bpt,
1373 srcPtr + (srcWidth * (row * 2 + 2)) * bpt,
1374 1, dstPtr + (dstWidth * row + 1) * bpt);
1375 do_row(datatype, comps, 1,
1376 srcPtr + (srcWidth * (row * 2 + 1) + srcWidth - 1) * bpt,
1377 srcPtr + (srcWidth * (row * 2 + 2) + srcWidth - 1) * bpt,
1378 1, dstPtr + (dstWidth * row + 1 + dstWidth - 1) * bpt);
1379 }
1380 }
1381 }
1382 }
1383 }
1384
1385
1386 /**
1387 * Down-sample a texture image to produce the next lower mipmap level.
1388 * \param comps components per texel (1, 2, 3 or 4)
1389 * \param srcRowStride stride between source rows, in texels
1390 * \param dstRowStride stride between destination rows, in texels
1391 */
1392 void
1393 _mesa_generate_mipmap_level(GLenum target,
1394 GLenum datatype, GLuint comps,
1395 GLint border,
1396 GLint srcWidth, GLint srcHeight, GLint srcDepth,
1397 const GLubyte *srcData,
1398 GLint srcRowStride,
1399 GLint dstWidth, GLint dstHeight, GLint dstDepth,
1400 GLubyte *dstData,
1401 GLint dstRowStride)
1402 {
1403 /*
1404 * We use simple 2x2 averaging to compute the next mipmap level.
1405 */
1406 switch (target) {
1407 case GL_TEXTURE_1D:
1408 make_1d_mipmap(datatype, comps, border,
1409 srcWidth, srcData,
1410 dstWidth, dstData);
1411 break;
1412 case GL_TEXTURE_2D:
1413 case GL_TEXTURE_CUBE_MAP_POSITIVE_X_ARB:
1414 case GL_TEXTURE_CUBE_MAP_NEGATIVE_X_ARB:
1415 case GL_TEXTURE_CUBE_MAP_POSITIVE_Y_ARB:
1416 case GL_TEXTURE_CUBE_MAP_NEGATIVE_Y_ARB:
1417 case GL_TEXTURE_CUBE_MAP_POSITIVE_Z_ARB:
1418 case GL_TEXTURE_CUBE_MAP_NEGATIVE_Z_ARB:
1419 make_2d_mipmap(datatype, comps, border,
1420 srcWidth, srcHeight, srcData, srcRowStride,
1421 dstWidth, dstHeight, dstData, dstRowStride);
1422 break;
1423 case GL_TEXTURE_3D:
1424 make_3d_mipmap(datatype, comps, border,
1425 srcWidth, srcHeight, srcDepth,
1426 srcData, srcRowStride,
1427 dstWidth, dstHeight, dstDepth,
1428 dstData, dstRowStride);
1429 break;
1430 case GL_TEXTURE_1D_ARRAY_EXT:
1431 make_1d_stack_mipmap(datatype, comps, border,
1432 srcWidth, srcData, srcRowStride,
1433 dstWidth, dstHeight,
1434 dstData, dstRowStride);
1435 break;
1436 case GL_TEXTURE_2D_ARRAY_EXT:
1437 make_2d_stack_mipmap(datatype, comps, border,
1438 srcWidth, srcHeight,
1439 srcData, srcRowStride,
1440 dstWidth, dstHeight,
1441 dstDepth, dstData, dstRowStride);
1442 break;
1443 case GL_TEXTURE_RECTANGLE_NV:
1444 /* no mipmaps, do nothing */
1445 break;
1446 default:
1447 _mesa_problem(NULL, "bad dimensions in _mesa_generate_mipmaps");
1448 return;
1449 }
1450 }
1451
1452
1453 /**
1454 * compute next (level+1) image size
1455 * \return GL_FALSE if no smaller size can be generated (eg. src is 1x1x1 size)
1456 */
1457 static GLboolean
1458 next_mipmap_level_size(GLenum target, GLint border,
1459 GLint srcWidth, GLint srcHeight, GLint srcDepth,
1460 GLint *dstWidth, GLint *dstHeight, GLint *dstDepth)
1461 {
1462 if (srcWidth - 2 * border > 1) {
1463 *dstWidth = (srcWidth - 2 * border) / 2 + 2 * border;
1464 }
1465 else {
1466 *dstWidth = srcWidth; /* can't go smaller */
1467 }
1468
1469 if ((srcHeight - 2 * border > 1) &&
1470 (target != GL_TEXTURE_1D_ARRAY_EXT)) {
1471 *dstHeight = (srcHeight - 2 * border) / 2 + 2 * border;
1472 }
1473 else {
1474 *dstHeight = srcHeight; /* can't go smaller */
1475 }
1476
1477 if ((srcDepth - 2 * border > 1) &&
1478 (target != GL_TEXTURE_2D_ARRAY_EXT)) {
1479 *dstDepth = (srcDepth - 2 * border) / 2 + 2 * border;
1480 }
1481 else {
1482 *dstDepth = srcDepth; /* can't go smaller */
1483 }
1484
1485 if (*dstWidth == srcWidth &&
1486 *dstHeight == srcHeight &&
1487 *dstDepth == srcDepth) {
1488 return GL_FALSE;
1489 }
1490 else {
1491 return GL_TRUE;
1492 }
1493 }
1494
1495
1496
1497
1498 /**
1499 * Automatic mipmap generation.
1500 * This is the fallback/default function for ctx->Driver.GenerateMipmap().
1501 * Generate a complete set of mipmaps from texObj's BaseLevel image.
1502 * Stop at texObj's MaxLevel or when we get to the 1x1 texture.
1503 * For cube maps, target will be one of
1504 * GL_TEXTURE_CUBE_MAP_POSITIVE/NEGATIVE_X/Y/Z; never GL_TEXTURE_CUBE_MAP.
1505 */
1506 void
1507 _mesa_generate_mipmap(GLcontext *ctx, GLenum target,
1508 struct gl_texture_object *texObj)
1509 {
1510 const struct gl_texture_image *srcImage;
1511 gl_format convertFormat;
1512 const GLubyte *srcData = NULL;
1513 GLubyte *dstData = NULL;
1514 GLint level, maxLevels;
1515 GLenum datatype;
1516 GLuint comps;
1517
1518 ASSERT(texObj);
1519 srcImage = _mesa_select_tex_image(ctx, texObj, target, texObj->BaseLevel);
1520 ASSERT(srcImage);
1521
1522 maxLevels = _mesa_max_texture_levels(ctx, texObj->Target);
1523 ASSERT(maxLevels > 0); /* bad target */
1524
1525 /* Find convertFormat - the format that do_row() will process */
1526
1527 if (_mesa_is_format_compressed(srcImage->TexFormat)) {
1528 /* setup for compressed textures - need to allocate temporary
1529 * image buffers to hold uncompressed images.
1530 */
1531 GLuint row;
1532 GLint components, size;
1533 GLchan *dst;
1534
1535 assert(texObj->Target == GL_TEXTURE_2D ||
1536 texObj->Target == GL_TEXTURE_CUBE_MAP_ARB);
1537
1538 if (srcImage->_BaseFormat == GL_RGB) {
1539 convertFormat = MESA_FORMAT_RGB888;
1540 components = 3;
1541 }
1542 else if (srcImage->_BaseFormat == GL_RGBA) {
1543 convertFormat = MESA_FORMAT_RGBA8888;
1544 components = 4;
1545 }
1546 else {
1547 _mesa_problem(ctx, "bad srcImage->_BaseFormat in _mesa_generate_mipmaps");
1548 return;
1549 }
1550
1551 /* allocate storage for uncompressed GL_RGB or GL_RGBA images */
1552 size = _mesa_bytes_per_pixel(srcImage->_BaseFormat, CHAN_TYPE)
1553 * srcImage->Width * srcImage->Height * srcImage->Depth + 20;
1554 /* 20 extra bytes, just be safe when calling last FetchTexel */
1555 srcData = (GLubyte *) malloc(size);
1556 if (!srcData) {
1557 _mesa_error(ctx, GL_OUT_OF_MEMORY, "generate mipmaps");
1558 return;
1559 }
1560 dstData = (GLubyte *) malloc(size / 2); /* 1/4 would probably be OK */
1561 if (!dstData) {
1562 _mesa_error(ctx, GL_OUT_OF_MEMORY, "generate mipmaps");
1563 free((void *) srcData);
1564 return;
1565 }
1566
1567 /* decompress base image here */
1568 dst = (GLchan *) srcData;
1569 for (row = 0; row < srcImage->Height; row++) {
1570 GLuint col;
1571 for (col = 0; col < srcImage->Width; col++) {
1572 srcImage->FetchTexelc(srcImage, col, row, 0, dst);
1573 dst += components;
1574 }
1575 }
1576 }
1577 else {
1578 /* uncompressed */
1579 convertFormat = srcImage->TexFormat;
1580 }
1581
1582 _mesa_format_to_type_and_comps(convertFormat, &datatype, &comps);
1583
1584 for (level = texObj->BaseLevel; level < texObj->MaxLevel
1585 && level < maxLevels - 1; level++) {
1586 /* generate image[level+1] from image[level] */
1587 const struct gl_texture_image *srcImage;
1588 struct gl_texture_image *dstImage;
1589 GLint srcWidth, srcHeight, srcDepth;
1590 GLint dstWidth, dstHeight, dstDepth;
1591 GLint border;
1592 GLboolean nextLevel;
1593
1594 /* get src image parameters */
1595 srcImage = _mesa_select_tex_image(ctx, texObj, target, level);
1596 ASSERT(srcImage);
1597 srcWidth = srcImage->Width;
1598 srcHeight = srcImage->Height;
1599 srcDepth = srcImage->Depth;
1600 border = srcImage->Border;
1601
1602 nextLevel = next_mipmap_level_size(target, border,
1603 srcWidth, srcHeight, srcDepth,
1604 &dstWidth, &dstHeight, &dstDepth);
1605 if (!nextLevel) {
1606 /* all done */
1607 if (_mesa_is_format_compressed(srcImage->TexFormat)) {
1608 free((void *) srcData);
1609 free(dstData);
1610 }
1611 return;
1612 }
1613
1614 /* get dest gl_texture_image */
1615 dstImage = _mesa_get_tex_image(ctx, texObj, target, level + 1);
1616 if (!dstImage) {
1617 _mesa_error(ctx, GL_OUT_OF_MEMORY, "generating mipmaps");
1618 return;
1619 }
1620
1621 /* Free old image data */
1622 if (dstImage->Data)
1623 ctx->Driver.FreeTexImageData(ctx, dstImage);
1624
1625 /* initialize new image */
1626 _mesa_init_teximage_fields(ctx, target, dstImage, dstWidth, dstHeight,
1627 dstDepth, border, srcImage->InternalFormat);
1628 dstImage->DriverData = NULL;
1629 dstImage->TexFormat = srcImage->TexFormat;
1630 dstImage->FetchTexelc = srcImage->FetchTexelc;
1631 dstImage->FetchTexelf = srcImage->FetchTexelf;
1632
1633 /* Alloc new teximage data buffer */
1634 {
1635 GLuint size = _mesa_format_image_size(dstImage->TexFormat,
1636 dstWidth, dstHeight, dstDepth);
1637 dstImage->Data = _mesa_alloc_texmemory(size);
1638 if (!dstImage->Data) {
1639 _mesa_error(ctx, GL_OUT_OF_MEMORY, "generating mipmaps");
1640 return;
1641 }
1642 }
1643
1644 /* Setup src and dest data pointers */
1645 if (_mesa_is_format_compressed(dstImage->TexFormat)) {
1646 /* srcData and dstData are already set */
1647 ASSERT(srcData);
1648 ASSERT(dstData);
1649 }
1650 else {
1651 srcData = (const GLubyte *) srcImage->Data;
1652 dstData = (GLubyte *) dstImage->Data;
1653 }
1654
1655 ASSERT(dstImage->TexFormat);
1656 ASSERT(dstImage->FetchTexelc);
1657 ASSERT(dstImage->FetchTexelf);
1658
1659 _mesa_generate_mipmap_level(target, datatype, comps, border,
1660 srcWidth, srcHeight, srcDepth,
1661 srcData, srcImage->RowStride,
1662 dstWidth, dstHeight, dstDepth,
1663 dstData, dstImage->RowStride);
1664
1665
1666 if (_mesa_is_format_compressed(dstImage->TexFormat)) {
1667 GLubyte *temp;
1668 /* compress image from dstData into dstImage->Data */
1669 const GLenum srcFormat = _mesa_get_format_base_format(convertFormat);
1670 GLint dstRowStride
1671 = _mesa_format_row_stride(dstImage->TexFormat, dstWidth);
1672 ASSERT(srcFormat == GL_RGB || srcFormat == GL_RGBA);
1673
1674 _mesa_texstore(ctx, 2, dstImage->_BaseFormat,
1675 dstImage->TexFormat,
1676 dstImage->Data,
1677 0, 0, 0, /* dstX/Y/Zoffset */
1678 dstRowStride, 0, /* strides */
1679 dstWidth, dstHeight, 1, /* size */
1680 srcFormat, CHAN_TYPE,
1681 dstData, /* src data, actually */
1682 &ctx->DefaultPacking);
1683
1684 /* swap src and dest pointers */
1685 temp = (GLubyte *) srcData;
1686 srcData = dstData;
1687 dstData = temp;
1688 }
1689
1690 } /* loop over mipmap levels */
1691 }
1692
1693
1694 /**
1695 * Helper function for drivers which need to rescale texture images to
1696 * certain aspect ratios.
1697 * Nearest filtering only (for broken hardware that can't support
1698 * all aspect ratios). This can be made a lot faster, but I don't
1699 * really care enough...
1700 */
1701 void
1702 _mesa_rescale_teximage2d(GLuint bytesPerPixel,
1703 GLuint srcStrideInPixels,
1704 GLuint dstRowStride,
1705 GLint srcWidth, GLint srcHeight,
1706 GLint dstWidth, GLint dstHeight,
1707 const GLvoid *srcImage, GLvoid *dstImage)
1708 {
1709 GLint row, col;
1710
1711 #define INNER_LOOP( TYPE, HOP, WOP ) \
1712 for ( row = 0 ; row < dstHeight ; row++ ) { \
1713 GLint srcRow = row HOP hScale; \
1714 for ( col = 0 ; col < dstWidth ; col++ ) { \
1715 GLint srcCol = col WOP wScale; \
1716 dst[col] = src[srcRow * srcStrideInPixels + srcCol]; \
1717 } \
1718 dst = (TYPE *) ((GLubyte *) dst + dstRowStride); \
1719 } \
1720
1721 #define RESCALE_IMAGE( TYPE ) \
1722 do { \
1723 const TYPE *src = (const TYPE *)srcImage; \
1724 TYPE *dst = (TYPE *)dstImage; \
1725 \
1726 if ( srcHeight < dstHeight ) { \
1727 const GLint hScale = dstHeight / srcHeight; \
1728 if ( srcWidth < dstWidth ) { \
1729 const GLint wScale = dstWidth / srcWidth; \
1730 INNER_LOOP( TYPE, /, / ); \
1731 } \
1732 else { \
1733 const GLint wScale = srcWidth / dstWidth; \
1734 INNER_LOOP( TYPE, /, * ); \
1735 } \
1736 } \
1737 else { \
1738 const GLint hScale = srcHeight / dstHeight; \
1739 if ( srcWidth < dstWidth ) { \
1740 const GLint wScale = dstWidth / srcWidth; \
1741 INNER_LOOP( TYPE, *, / ); \
1742 } \
1743 else { \
1744 const GLint wScale = srcWidth / dstWidth; \
1745 INNER_LOOP( TYPE, *, * ); \
1746 } \
1747 } \
1748 } while (0)
1749
1750 switch ( bytesPerPixel ) {
1751 case 4:
1752 RESCALE_IMAGE( GLuint );
1753 break;
1754
1755 case 2:
1756 RESCALE_IMAGE( GLushort );
1757 break;
1758
1759 case 1:
1760 RESCALE_IMAGE( GLubyte );
1761 break;
1762 default:
1763 _mesa_problem(NULL,"unexpected bytes/pixel in _mesa_rescale_teximage2d");
1764 }
1765 }
1766
1767
1768 /**
1769 * Upscale an image by replication, not (typical) stretching.
1770 * We use this when the image width or height is less than a
1771 * certain size (4, 8) and we need to upscale an image.
1772 */
1773 void
1774 _mesa_upscale_teximage2d(GLsizei inWidth, GLsizei inHeight,
1775 GLsizei outWidth, GLsizei outHeight,
1776 GLint comps, const GLchan *src, GLint srcRowStride,
1777 GLchan *dest )
1778 {
1779 GLint i, j, k;
1780
1781 ASSERT(outWidth >= inWidth);
1782 ASSERT(outHeight >= inHeight);
1783 #if 0
1784 ASSERT(inWidth == 1 || inWidth == 2 || inHeight == 1 || inHeight == 2);
1785 ASSERT((outWidth & 3) == 0);
1786 ASSERT((outHeight & 3) == 0);
1787 #endif
1788
1789 for (i = 0; i < outHeight; i++) {
1790 const GLint ii = i % inHeight;
1791 for (j = 0; j < outWidth; j++) {
1792 const GLint jj = j % inWidth;
1793 for (k = 0; k < comps; k++) {
1794 dest[(i * outWidth + j) * comps + k]
1795 = src[ii * srcRowStride + jj * comps + k];
1796 }
1797 }
1798 }
1799 }
1800