1 324891a5 2006-06-25 devnull #include <u.h>
2 324891a5 2006-06-25 devnull #include <libc.h>
3 324891a5 2006-06-25 devnull #include <bio.h>
4 324891a5 2006-06-25 devnull #include <draw.h>
5 324891a5 2006-06-25 devnull #include <memdraw.h>
7 324891a5 2006-06-25 devnull #define DBG if(0)
8 324891a5 2006-06-25 devnull #define RGB2K(r,g,b) ((299*((u32int)(r))+587*((u32int)(g))+114*((u32int)(b)))/1000)
11 324891a5 2006-06-25 devnull * This program tests the 'memimagedraw' primitive stochastically.
12 324891a5 2006-06-25 devnull * It tests the combination aspects of it thoroughly, but since the
13 324891a5 2006-06-25 devnull * three images it uses are disjoint, it makes no check of the
14 324891a5 2006-06-25 devnull * correct behavior when images overlap. That is, however, much
15 324891a5 2006-06-25 devnull * easier to get right and to test.
18 324891a5 2006-06-25 devnull void drawonepixel(Memimage*, Point, Memimage*, Point, Memimage*, Point);
19 324891a5 2006-06-25 devnull void verifyone(void);
20 324891a5 2006-06-25 devnull void verifyline(void);
21 324891a5 2006-06-25 devnull void verifyrect(void);
22 324891a5 2006-06-25 devnull void verifyrectrepl(int, int);
23 324891a5 2006-06-25 devnull void putpixel(Memimage *img, Point pt, u32int nv);
24 324891a5 2006-06-25 devnull u32int rgbatopix(uchar, uchar, uchar, uchar);
26 324891a5 2006-06-25 devnull char *dchan, *schan, *mchan;
27 324891a5 2006-06-25 devnull int dbpp, sbpp, mbpp;
29 324891a5 2006-06-25 devnull int drawdebug=0;
30 324891a5 2006-06-25 devnull int seed;
31 324891a5 2006-06-25 devnull int niters = 100;
32 324891a5 2006-06-25 devnull int dbpp; /* bits per pixel in destination */
33 324891a5 2006-06-25 devnull int sbpp; /* bits per pixel in src */
34 324891a5 2006-06-25 devnull int mbpp; /* bits per pixel in mask */
35 324891a5 2006-06-25 devnull int dpm; /* pixel mask at high part of byte, in destination */
36 324891a5 2006-06-25 devnull int nbytes; /* in destination */
38 324891a5 2006-06-25 devnull int Xrange = 64;
39 324891a5 2006-06-25 devnull int Yrange = 8;
41 324891a5 2006-06-25 devnull Memimage *dst;
42 324891a5 2006-06-25 devnull Memimage *src;
43 324891a5 2006-06-25 devnull Memimage *mask;
44 324891a5 2006-06-25 devnull Memimage *stmp;
45 324891a5 2006-06-25 devnull Memimage *mtmp;
46 324891a5 2006-06-25 devnull Memimage *ones;
47 324891a5 2006-06-25 devnull uchar *dstbits;
48 324891a5 2006-06-25 devnull uchar *srcbits;
49 324891a5 2006-06-25 devnull uchar *maskbits;
50 324891a5 2006-06-25 devnull u32int *savedstbits;
53 324891a5 2006-06-25 devnull rdb(void)
58 324891a5 2006-06-25 devnull iprint(char *fmt, ...)
61 324891a5 2006-06-25 devnull va_list va;
62 324891a5 2006-06-25 devnull char buf[1024];
64 324891a5 2006-06-25 devnull va_start(va, fmt);
65 324891a5 2006-06-25 devnull n = doprint(buf, buf+sizeof buf, fmt, va) - buf;
66 324891a5 2006-06-25 devnull va_end(va);
68 324891a5 2006-06-25 devnull write(1,buf,n);
69 324891a5 2006-06-25 devnull return 1;
73 324891a5 2006-06-25 devnull main(int argc, char *argv[])
75 324891a5 2006-06-25 devnull memimageinit();
76 324891a5 2006-06-25 devnull seed = time(0);
78 324891a5 2006-06-25 devnull ARGBEGIN{
79 324891a5 2006-06-25 devnull case 'x':
80 324891a5 2006-06-25 devnull Xrange = atoi(ARGF());
82 324891a5 2006-06-25 devnull case 'y':
83 324891a5 2006-06-25 devnull Yrange = atoi(ARGF());
85 324891a5 2006-06-25 devnull case 'n':
86 324891a5 2006-06-25 devnull niters = atoi(ARGF());
88 324891a5 2006-06-25 devnull case 's':
89 324891a5 2006-06-25 devnull seed = atoi(ARGF());
93 324891a5 2006-06-25 devnull dchan = "r8g8b8";
94 324891a5 2006-06-25 devnull schan = "r8g8b8";
95 324891a5 2006-06-25 devnull mchan = "r8g8b8";
96 324891a5 2006-06-25 devnull switch(argc){
97 324891a5 2006-06-25 devnull case 3: mchan = argv[2];
98 324891a5 2006-06-25 devnull case 2: schan = argv[1];
99 324891a5 2006-06-25 devnull case 1: dchan = argv[0];
100 324891a5 2006-06-25 devnull case 0: break;
101 324891a5 2006-06-25 devnull default: goto Usage;
103 324891a5 2006-06-25 devnull fprint(2, "usage: dtest [dchan [schan [mchan]]]\n");
104 324891a5 2006-06-25 devnull exits("usage");
107 324891a5 2006-06-25 devnull fmtinstall('b', numbconv); /* binary! */
109 324891a5 2006-06-25 devnull fprint(2, "%s -x %d -y %d -s 0x%x %s %s %s\n", argv0, Xrange, Yrange, seed, dchan, schan, mchan);
110 324891a5 2006-06-25 devnull srand(seed);
112 324891a5 2006-06-25 devnull dst = allocmemimage(Rect(0, 0, Xrange, Yrange), strtochan(dchan));
113 324891a5 2006-06-25 devnull src = allocmemimage(Rect(0, 0, Xrange, Yrange), strtochan(schan));
114 324891a5 2006-06-25 devnull mask = allocmemimage(Rect(0, 0, Xrange, Yrange), strtochan(mchan));
115 324891a5 2006-06-25 devnull stmp = allocmemimage(Rect(0, 0, Xrange, Yrange), strtochan(schan));
116 324891a5 2006-06-25 devnull mtmp = allocmemimage(Rect(0, 0, Xrange, Yrange), strtochan(mchan));
117 324891a5 2006-06-25 devnull ones = allocmemimage(Rect(0, 0, Xrange, Yrange), strtochan(mchan));
118 324891a5 2006-06-25 devnull /* print("chan %lux %lux %lux %lux %lux %lux\n", dst->chan, src->chan, mask->chan, stmp->chan, mtmp->chan, ones->chan); */
119 324891a5 2006-06-25 devnull if(dst==0 || src==0 || mask==0 || mtmp==0 || ones==0) {
121 324891a5 2006-06-25 devnull fprint(2, "dtest: allocation failed: %r\n");
122 324891a5 2006-06-25 devnull exits("alloc");
124 324891a5 2006-06-25 devnull nbytes = (4*Xrange+4)*Yrange;
125 324891a5 2006-06-25 devnull srcbits = malloc(nbytes);
126 324891a5 2006-06-25 devnull dstbits = malloc(nbytes);
127 324891a5 2006-06-25 devnull maskbits = malloc(nbytes);
128 324891a5 2006-06-25 devnull savedstbits = malloc(nbytes);
129 324891a5 2006-06-25 devnull if(dstbits==0 || srcbits==0 || maskbits==0 || savedstbits==0)
130 324891a5 2006-06-25 devnull goto Alloc;
131 324891a5 2006-06-25 devnull dbpp = dst->depth;
132 324891a5 2006-06-25 devnull sbpp = src->depth;
133 324891a5 2006-06-25 devnull mbpp = mask->depth;
134 324891a5 2006-06-25 devnull dpm = 0xFF ^ (0xFF>>dbpp);
135 324891a5 2006-06-25 devnull memset(ones->data->bdata, 0xFF, ones->width*sizeof(u32int)*Yrange);
138 324891a5 2006-06-25 devnull fprint(2, "dtest: verify single pixel operation\n");
139 324891a5 2006-06-25 devnull verifyone();
141 324891a5 2006-06-25 devnull fprint(2, "dtest: verify full line non-replicated\n");
142 324891a5 2006-06-25 devnull verifyline();
144 324891a5 2006-06-25 devnull fprint(2, "dtest: verify full rectangle non-replicated\n");
145 324891a5 2006-06-25 devnull verifyrect();
147 324891a5 2006-06-25 devnull fprint(2, "dtest: verify full rectangle source replicated\n");
148 324891a5 2006-06-25 devnull verifyrectrepl(1, 0);
150 324891a5 2006-06-25 devnull fprint(2, "dtest: verify full rectangle mask replicated\n");
151 324891a5 2006-06-25 devnull verifyrectrepl(0, 1);
153 324891a5 2006-06-25 devnull fprint(2, "dtest: verify full rectangle source and mask replicated\n");
154 324891a5 2006-06-25 devnull verifyrectrepl(1, 1);
156 324891a5 2006-06-25 devnull exits(0);
160 324891a5 2006-06-25 devnull * Dump out an ASCII representation of an image. The label specifies
161 324891a5 2006-06-25 devnull * a list of characters to put at various points in the picture.
163 324891a5 2006-06-25 devnull static void
164 324891a5 2006-06-25 devnull Bprintr5g6b5(Biobuf *bio, char*, u32int v)
166 324891a5 2006-06-25 devnull int r,g,b;
167 324891a5 2006-06-25 devnull r = (v>>11)&31;
168 324891a5 2006-06-25 devnull g = (v>>5)&63;
169 324891a5 2006-06-25 devnull b = v&31;
170 324891a5 2006-06-25 devnull Bprint(bio, "%.2x%.2x%.2x", r,g,b);
173 324891a5 2006-06-25 devnull static void
174 324891a5 2006-06-25 devnull Bprintr5g5b5a1(Biobuf *bio, char*, u32int v)
176 324891a5 2006-06-25 devnull int r,g,b,a;
177 324891a5 2006-06-25 devnull r = (v>>11)&31;
178 324891a5 2006-06-25 devnull g = (v>>6)&31;
179 324891a5 2006-06-25 devnull b = (v>>1)&31;
180 324891a5 2006-06-25 devnull a = v&1;
181 324891a5 2006-06-25 devnull Bprint(bio, "%.2x%.2x%.2x%.2x", r,g,b,a);
185 324891a5 2006-06-25 devnull dumpimage(char *name, Memimage *img, void *vdata, Point labelpt)
187 324891a5 2006-06-25 devnull Biobuf b;
188 324891a5 2006-06-25 devnull uchar *data;
189 324891a5 2006-06-25 devnull uchar *p;
190 324891a5 2006-06-25 devnull char *arg;
191 324891a5 2006-06-25 devnull void (*fmt)(Biobuf*, char*, u32int);
192 324891a5 2006-06-25 devnull int npr, x, y, nb, bpp;
193 324891a5 2006-06-25 devnull u32int v, mask;
194 324891a5 2006-06-25 devnull Rectangle r;
196 324891a5 2006-06-25 devnull fmt = nil;
197 324891a5 2006-06-25 devnull arg = nil;
198 324891a5 2006-06-25 devnull switch(img->depth){
202 324891a5 2006-06-25 devnull fmt = (void(*)(Biobuf*,char*,u32int))Bprint;
203 324891a5 2006-06-25 devnull arg = "%.1ux";
206 324891a5 2006-06-25 devnull fmt = (void(*)(Biobuf*,char*,u32int))Bprint;
207 324891a5 2006-06-25 devnull arg = "%.2ux";
209 324891a5 2006-06-25 devnull case 16:
210 324891a5 2006-06-25 devnull arg = nil;
211 324891a5 2006-06-25 devnull if(img->chan == RGB16)
212 324891a5 2006-06-25 devnull fmt = Bprintr5g6b5;
214 324891a5 2006-06-25 devnull fmt = (void(*)(Biobuf*,char*,u32int))Bprint;
215 324891a5 2006-06-25 devnull arg = "%.4ux";
218 324891a5 2006-06-25 devnull case 24:
219 324891a5 2006-06-25 devnull fmt = (void(*)(Biobuf*,char*,u32int))Bprint;
220 324891a5 2006-06-25 devnull arg = "%.6lux";
222 324891a5 2006-06-25 devnull case 32:
223 324891a5 2006-06-25 devnull fmt = (void(*)(Biobuf*,char*,u32int))Bprint;
224 324891a5 2006-06-25 devnull arg = "%.8lux";
227 324891a5 2006-06-25 devnull if(fmt == nil){
228 324891a5 2006-06-25 devnull fprint(2, "bad format\n");
229 324891a5 2006-06-25 devnull abort();
232 324891a5 2006-06-25 devnull r = img->r;
233 324891a5 2006-06-25 devnull Binit(&b, 2, OWRITE);
234 324891a5 2006-06-25 devnull data = vdata;
235 324891a5 2006-06-25 devnull bpp = img->depth;
236 324891a5 2006-06-25 devnull Bprint(&b, "%s\t%d\tr %R clipr %R repl %d data %p *%P\n", name, r.min.x, r, img->clipr, (img->flags&Frepl) ? 1 : 0, vdata, labelpt);
237 324891a5 2006-06-25 devnull mask = (1ULL<<bpp)-1;
238 324891a5 2006-06-25 devnull /* for(y=r.min.y; y<r.max.y; y++){ */
239 324891a5 2006-06-25 devnull for(y=0; y<Yrange; y++){
242 324891a5 2006-06-25 devnull p = data+(byteaddr(img, Pt(0,y))-(uchar*)img->data->bdata);
243 324891a5 2006-06-25 devnull Bprint(&b, "%-4d\t", y);
244 324891a5 2006-06-25 devnull /* for(x=r.min.x; x<r.max.x; x++){ */
245 324891a5 2006-06-25 devnull for(x=0; x<Xrange; x++){
246 324891a5 2006-06-25 devnull if(x==0)
247 324891a5 2006-06-25 devnull Bprint(&b, "\t");
249 324891a5 2006-06-25 devnull if(x != 0 && (x%8)==0)
250 324891a5 2006-06-25 devnull Bprint(&b, " ");
252 324891a5 2006-06-25 devnull npr = 0;
253 324891a5 2006-06-25 devnull if(x==labelpt.x && y==labelpt.y){
254 324891a5 2006-06-25 devnull Bprint(&b, "*");
257 324891a5 2006-06-25 devnull if(npr == 0)
258 324891a5 2006-06-25 devnull Bprint(&b, " ");
260 324891a5 2006-06-25 devnull while(nb < bpp){
261 324891a5 2006-06-25 devnull v &= (1<<nb)-1;
262 324891a5 2006-06-25 devnull v |= (u32int)(*p++) << nb;
263 324891a5 2006-06-25 devnull nb += 8;
265 324891a5 2006-06-25 devnull nb -= bpp;
266 324891a5 2006-06-25 devnull /* print("bpp %d v %.8lux mask %.8lux nb %d\n", bpp, v, mask, nb); */
267 324891a5 2006-06-25 devnull fmt(&b, arg, (v>>nb)&mask);
269 324891a5 2006-06-25 devnull Bprint(&b, "\n");
271 324891a5 2006-06-25 devnull Bterm(&b);
275 324891a5 2006-06-25 devnull * Verify that the destination pixel has the specified value.
276 324891a5 2006-06-25 devnull * The value is in the high bits of v, suitably masked, but must
277 324891a5 2006-06-25 devnull * be extracted from the destination Memimage.
280 324891a5 2006-06-25 devnull checkone(Point p, Point sp, Point mp)
282 324891a5 2006-06-25 devnull int delta;
283 324891a5 2006-06-25 devnull uchar *dp, *sdp;
285 324891a5 2006-06-25 devnull delta = (uchar*)byteaddr(dst, p)-(uchar*)dst->data->bdata;
286 324891a5 2006-06-25 devnull dp = (uchar*)dst->data->bdata+delta;
287 324891a5 2006-06-25 devnull sdp = (uchar*)savedstbits+delta;
289 324891a5 2006-06-25 devnull if(memcmp(dp, sdp, (dst->depth+7)/8) != 0) {
290 324891a5 2006-06-25 devnull fprint(2, "dtest: one bad pixel drawing at dst %P from source %P mask %P\n", p, sp, mp);
291 324891a5 2006-06-25 devnull fprint(2, " %.2ux %.2ux %.2ux %.2ux should be %.2ux %.2ux %.2ux %.2ux\n",
292 324891a5 2006-06-25 devnull dp[0], dp[1], dp[2], dp[3], sdp[0], sdp[1], sdp[2], sdp[3]);
293 324891a5 2006-06-25 devnull fprint(2, "addresses dst %p src %p mask %p\n", dp, byteaddr(src, sp), byteaddr(mask, mp));
294 324891a5 2006-06-25 devnull dumpimage("src", src, src->data->bdata, sp);
295 324891a5 2006-06-25 devnull dumpimage("mask", mask, mask->data->bdata, mp);
296 324891a5 2006-06-25 devnull dumpimage("origdst", dst, dstbits, p);
297 324891a5 2006-06-25 devnull dumpimage("dst", dst, dst->data->bdata, p);
298 324891a5 2006-06-25 devnull dumpimage("gooddst", dst, savedstbits, p);
299 324891a5 2006-06-25 devnull abort();
304 324891a5 2006-06-25 devnull * Verify that the destination line has the same value as the saved line.
306 324891a5 2006-06-25 devnull #define RECTPTS(r) (r).min.x, (r).min.y, (r).max.x, (r).max.y
308 324891a5 2006-06-25 devnull checkline(Rectangle r, Point sp, Point mp, int y, Memimage *stmp, Memimage *mtmp)
310 324891a5 2006-06-25 devnull u32int *dp;
312 324891a5 2006-06-25 devnull u32int *saved;
314 324891a5 2006-06-25 devnull dp = wordaddr(dst, Pt(0, y));
315 324891a5 2006-06-25 devnull saved = savedstbits + y*dst->width;
316 324891a5 2006-06-25 devnull if(dst->depth < 8)
317 324891a5 2006-06-25 devnull nb = Xrange/(8/dst->depth);
319 324891a5 2006-06-25 devnull nb = Xrange*(dst->depth/8);
320 324891a5 2006-06-25 devnull if(memcmp(dp, saved, nb) != 0){
321 324891a5 2006-06-25 devnull fprint(2, "dtest: bad line at y=%d; saved %p dp %p\n", y, saved, dp);
322 324891a5 2006-06-25 devnull fprint(2, "draw dst %R src %P mask %P\n", r, sp, mp);
323 324891a5 2006-06-25 devnull dumpimage("src", src, src->data->bdata, sp);
324 324891a5 2006-06-25 devnull if(stmp) dumpimage("stmp", stmp, stmp->data->bdata, sp);
325 324891a5 2006-06-25 devnull dumpimage("mask", mask, mask->data->bdata, mp);
326 324891a5 2006-06-25 devnull if(mtmp) dumpimage("mtmp", mtmp, mtmp->data->bdata, mp);
327 324891a5 2006-06-25 devnull dumpimage("origdst", dst, dstbits, r.min);
328 324891a5 2006-06-25 devnull dumpimage("dst", dst, dst->data->bdata, r.min);
329 324891a5 2006-06-25 devnull dumpimage("gooddst", dst, savedstbits, r.min);
330 324891a5 2006-06-25 devnull abort();
335 324891a5 2006-06-25 devnull * Fill the bits of an image with random data.
336 324891a5 2006-06-25 devnull * The Memimage parameter is used only to make sure
337 324891a5 2006-06-25 devnull * the data is well formatted: only ucbits is written.
340 324891a5 2006-06-25 devnull fill(Memimage *img, uchar *ucbits)
342 324891a5 2006-06-25 devnull int i, x, y;
343 324891a5 2006-06-25 devnull ushort *up;
344 324891a5 2006-06-25 devnull uchar alpha, r, g, b;
345 324891a5 2006-06-25 devnull void *data;
347 324891a5 2006-06-25 devnull if((img->flags&Falpha) == 0){
348 324891a5 2006-06-25 devnull up = (ushort*)ucbits;
349 324891a5 2006-06-25 devnull for(i=0; i<nbytes/2; i++)
350 324891a5 2006-06-25 devnull *up++ = lrand() >> 7;
351 324891a5 2006-06-25 devnull if(i+i != nbytes)
352 324891a5 2006-06-25 devnull *(uchar*)up = lrand() >> 7;
354 324891a5 2006-06-25 devnull data = img->data->bdata;
355 324891a5 2006-06-25 devnull img->data->bdata = ucbits;
357 324891a5 2006-06-25 devnull for(x=img->r.min.x; x<img->r.max.x; x++)
358 324891a5 2006-06-25 devnull for(y=img->r.min.y; y<img->r.max.y; y++){
359 324891a5 2006-06-25 devnull alpha = rand() >> 4;
360 324891a5 2006-06-25 devnull r = rand()%(alpha+1);
361 324891a5 2006-06-25 devnull g = rand()%(alpha+1);
362 324891a5 2006-06-25 devnull b = rand()%(alpha+1);
363 324891a5 2006-06-25 devnull putpixel(img, Pt(x,y), rgbatopix(r,g,b,alpha));
365 324891a5 2006-06-25 devnull img->data->bdata = data;
371 324891a5 2006-06-25 devnull * Mask is preset; do the rest
374 324891a5 2006-06-25 devnull verifyonemask(void)
376 324891a5 2006-06-25 devnull Point dp, sp, mp;
378 324891a5 2006-06-25 devnull fill(dst, dstbits);
379 324891a5 2006-06-25 devnull fill(src, srcbits);
380 324891a5 2006-06-25 devnull memmove(dst->data->bdata, dstbits, dst->width*sizeof(u32int)*Yrange);
381 324891a5 2006-06-25 devnull memmove(src->data->bdata, srcbits, src->width*sizeof(u32int)*Yrange);
382 324891a5 2006-06-25 devnull memmove(mask->data->bdata, maskbits, mask->width*sizeof(u32int)*Yrange);
384 324891a5 2006-06-25 devnull dp.x = nrand(Xrange);
385 324891a5 2006-06-25 devnull dp.y = nrand(Yrange);
387 324891a5 2006-06-25 devnull sp.x = nrand(Xrange);
388 324891a5 2006-06-25 devnull sp.y = nrand(Yrange);
390 324891a5 2006-06-25 devnull mp.x = nrand(Xrange);
391 324891a5 2006-06-25 devnull mp.y = nrand(Yrange);
393 324891a5 2006-06-25 devnull drawonepixel(dst, dp, src, sp, mask, mp);
394 324891a5 2006-06-25 devnull memmove(mask->data->bdata, maskbits, mask->width*sizeof(u32int)*Yrange);
395 324891a5 2006-06-25 devnull memmove(savedstbits, dst->data->bdata, dst->width*sizeof(u32int)*Yrange);
397 324891a5 2006-06-25 devnull memmove(dst->data->bdata, dstbits, dst->width*sizeof(u32int)*Yrange);
398 324891a5 2006-06-25 devnull memimagedraw(dst, Rect(dp.x, dp.y, dp.x+1, dp.y+1), src, sp, mask, mp, SoverD);
399 324891a5 2006-06-25 devnull memmove(mask->data->bdata, maskbits, mask->width*sizeof(u32int)*Yrange);
401 324891a5 2006-06-25 devnull checkone(dp, sp, mp);
405 324891a5 2006-06-25 devnull verifyone(void)
409 324891a5 2006-06-25 devnull /* mask all zeros */
410 324891a5 2006-06-25 devnull memset(maskbits, 0, nbytes);
411 324891a5 2006-06-25 devnull for(i=0; i<niters; i++)
412 324891a5 2006-06-25 devnull verifyonemask();
414 324891a5 2006-06-25 devnull /* mask all ones */
415 324891a5 2006-06-25 devnull memset(maskbits, 0xFF, nbytes);
416 324891a5 2006-06-25 devnull for(i=0; i<niters; i++)
417 324891a5 2006-06-25 devnull verifyonemask();
419 324891a5 2006-06-25 devnull /* random mask */
420 324891a5 2006-06-25 devnull for(i=0; i<niters; i++){
421 324891a5 2006-06-25 devnull fill(mask, maskbits);
422 324891a5 2006-06-25 devnull verifyonemask();
427 324891a5 2006-06-25 devnull * Mask is preset; do the rest
430 324891a5 2006-06-25 devnull verifylinemask(void)
432 324891a5 2006-06-25 devnull Point sp, mp, tp, up;
433 324891a5 2006-06-25 devnull Rectangle dr;
436 324891a5 2006-06-25 devnull fill(dst, dstbits);
437 324891a5 2006-06-25 devnull fill(src, srcbits);
438 324891a5 2006-06-25 devnull memmove(dst->data->bdata, dstbits, dst->width*sizeof(u32int)*Yrange);
439 324891a5 2006-06-25 devnull memmove(src->data->bdata, srcbits, src->width*sizeof(u32int)*Yrange);
440 324891a5 2006-06-25 devnull memmove(mask->data->bdata, maskbits, mask->width*sizeof(u32int)*Yrange);
442 324891a5 2006-06-25 devnull dr.min.x = nrand(Xrange-1);
443 324891a5 2006-06-25 devnull dr.min.y = nrand(Yrange-1);
444 324891a5 2006-06-25 devnull dr.max.x = dr.min.x + 1 + nrand(Xrange-1-dr.min.x);
445 324891a5 2006-06-25 devnull dr.max.y = dr.min.y + 1;
447 324891a5 2006-06-25 devnull sp.x = nrand(Xrange);
448 324891a5 2006-06-25 devnull sp.y = nrand(Yrange);
450 324891a5 2006-06-25 devnull mp.x = nrand(Xrange);
451 324891a5 2006-06-25 devnull mp.y = nrand(Yrange);
453 324891a5 2006-06-25 devnull tp = sp;
454 324891a5 2006-06-25 devnull up = mp;
455 324891a5 2006-06-25 devnull for(x=dr.min.x; x<dr.max.x && tp.x<Xrange && up.x<Xrange; x++,tp.x++,up.x++)
456 324891a5 2006-06-25 devnull memimagedraw(dst, Rect(x, dr.min.y, x+1, dr.min.y+1), src, tp, mask, up, SoverD);
457 324891a5 2006-06-25 devnull memmove(savedstbits, dst->data->bdata, dst->width*sizeof(u32int)*Yrange);
459 324891a5 2006-06-25 devnull memmove(dst->data->bdata, dstbits, dst->width*sizeof(u32int)*Yrange);
461 324891a5 2006-06-25 devnull memimagedraw(dst, dr, src, sp, mask, mp, SoverD);
462 324891a5 2006-06-25 devnull checkline(dr, drawrepl(src->r, sp), drawrepl(mask->r, mp), dr.min.y, nil, nil);
466 324891a5 2006-06-25 devnull verifyline(void)
470 324891a5 2006-06-25 devnull /* mask all ones */
471 324891a5 2006-06-25 devnull memset(maskbits, 0xFF, nbytes);
472 324891a5 2006-06-25 devnull for(i=0; i<niters; i++)
473 324891a5 2006-06-25 devnull verifylinemask();
475 324891a5 2006-06-25 devnull /* mask all zeros */
476 324891a5 2006-06-25 devnull memset(maskbits, 0, nbytes);
477 324891a5 2006-06-25 devnull for(i=0; i<niters; i++)
478 324891a5 2006-06-25 devnull verifylinemask();
480 324891a5 2006-06-25 devnull /* random mask */
481 324891a5 2006-06-25 devnull for(i=0; i<niters; i++){
482 324891a5 2006-06-25 devnull fill(mask, maskbits);
483 324891a5 2006-06-25 devnull verifylinemask();
488 324891a5 2006-06-25 devnull * Mask is preset; do the rest
491 324891a5 2006-06-25 devnull verifyrectmask(void)
493 324891a5 2006-06-25 devnull Point sp, mp, tp, up;
494 324891a5 2006-06-25 devnull Rectangle dr;
495 324891a5 2006-06-25 devnull int x, y;
497 324891a5 2006-06-25 devnull fill(dst, dstbits);
498 324891a5 2006-06-25 devnull fill(src, srcbits);
499 324891a5 2006-06-25 devnull memmove(dst->data->bdata, dstbits, dst->width*sizeof(u32int)*Yrange);
500 324891a5 2006-06-25 devnull memmove(src->data->bdata, srcbits, src->width*sizeof(u32int)*Yrange);
501 324891a5 2006-06-25 devnull memmove(mask->data->bdata, maskbits, mask->width*sizeof(u32int)*Yrange);
503 324891a5 2006-06-25 devnull dr.min.x = nrand(Xrange-1);
504 324891a5 2006-06-25 devnull dr.min.y = nrand(Yrange-1);
505 324891a5 2006-06-25 devnull dr.max.x = dr.min.x + 1 + nrand(Xrange-1-dr.min.x);
506 324891a5 2006-06-25 devnull dr.max.y = dr.min.y + 1 + nrand(Yrange-1-dr.min.y);
508 324891a5 2006-06-25 devnull sp.x = nrand(Xrange);
509 324891a5 2006-06-25 devnull sp.y = nrand(Yrange);
511 324891a5 2006-06-25 devnull mp.x = nrand(Xrange);
512 324891a5 2006-06-25 devnull mp.y = nrand(Yrange);
514 324891a5 2006-06-25 devnull tp = sp;
515 324891a5 2006-06-25 devnull up = mp;
516 324891a5 2006-06-25 devnull for(y=dr.min.y; y<dr.max.y && tp.y<Yrange && up.y<Yrange; y++,tp.y++,up.y++){
517 324891a5 2006-06-25 devnull for(x=dr.min.x; x<dr.max.x && tp.x<Xrange && up.x<Xrange; x++,tp.x++,up.x++)
518 324891a5 2006-06-25 devnull memimagedraw(dst, Rect(x, y, x+1, y+1), src, tp, mask, up, SoverD);
519 324891a5 2006-06-25 devnull tp.x = sp.x;
520 324891a5 2006-06-25 devnull up.x = mp.x;
522 324891a5 2006-06-25 devnull memmove(savedstbits, dst->data->bdata, dst->width*sizeof(u32int)*Yrange);
524 324891a5 2006-06-25 devnull memmove(dst->data->bdata, dstbits, dst->width*sizeof(u32int)*Yrange);
526 324891a5 2006-06-25 devnull memimagedraw(dst, dr, src, sp, mask, mp, SoverD);
527 324891a5 2006-06-25 devnull for(y=0; y<Yrange; y++)
528 324891a5 2006-06-25 devnull checkline(dr, drawrepl(src->r, sp), drawrepl(mask->r, mp), y, nil, nil);
532 324891a5 2006-06-25 devnull verifyrect(void)
536 324891a5 2006-06-25 devnull /* mask all zeros */
537 324891a5 2006-06-25 devnull memset(maskbits, 0, nbytes);
538 324891a5 2006-06-25 devnull for(i=0; i<niters; i++)
539 324891a5 2006-06-25 devnull verifyrectmask();
541 324891a5 2006-06-25 devnull /* mask all ones */
542 324891a5 2006-06-25 devnull memset(maskbits, 0xFF, nbytes);
543 324891a5 2006-06-25 devnull for(i=0; i<niters; i++)
544 324891a5 2006-06-25 devnull verifyrectmask();
546 324891a5 2006-06-25 devnull /* random mask */
547 324891a5 2006-06-25 devnull for(i=0; i<niters; i++){
548 324891a5 2006-06-25 devnull fill(mask, maskbits);
549 324891a5 2006-06-25 devnull verifyrectmask();
553 324891a5 2006-06-25 devnull Rectangle
554 324891a5 2006-06-25 devnull randrect(void)
556 324891a5 2006-06-25 devnull Rectangle r;
558 324891a5 2006-06-25 devnull r.min.x = nrand(Xrange-1);
559 324891a5 2006-06-25 devnull r.min.y = nrand(Yrange-1);
560 324891a5 2006-06-25 devnull r.max.x = r.min.x + 1 + nrand(Xrange-1-r.min.x);
561 324891a5 2006-06-25 devnull r.max.y = r.min.y + 1 + nrand(Yrange-1-r.min.y);
562 324891a5 2006-06-25 devnull return r;
566 324891a5 2006-06-25 devnull * Return coordinate corresponding to x withing range [minx, maxx)
569 324891a5 2006-06-25 devnull tilexy(int minx, int maxx, int x)
573 324891a5 2006-06-25 devnull sx = (x-minx) % (maxx-minx);
574 324891a5 2006-06-25 devnull if(sx < 0)
575 324891a5 2006-06-25 devnull sx += maxx-minx;
576 324891a5 2006-06-25 devnull return sx+minx;
580 324891a5 2006-06-25 devnull replicate(Memimage *i, Memimage *tmp)
582 324891a5 2006-06-25 devnull Rectangle r, r1;
583 324891a5 2006-06-25 devnull int x, y, nb;
585 324891a5 2006-06-25 devnull /* choose the replication window (i->r) */
586 324891a5 2006-06-25 devnull r.min.x = nrand(Xrange-1);
587 324891a5 2006-06-25 devnull r.min.y = nrand(Yrange-1);
588 324891a5 2006-06-25 devnull /* make it trivial more often than pure chance allows */
589 324891a5 2006-06-25 devnull switch(lrand()&0){
591 324891a5 2006-06-25 devnull r.max.x = r.min.x + 2;
592 324891a5 2006-06-25 devnull r.max.y = r.min.y + 2;
593 324891a5 2006-06-25 devnull if(r.max.x < Xrange && r.max.y < Yrange)
595 324891a5 2006-06-25 devnull /* fall through */
597 324891a5 2006-06-25 devnull r.max.x = r.min.x + 1;
598 324891a5 2006-06-25 devnull r.max.y = r.min.y + 1;
600 324891a5 2006-06-25 devnull default:
601 324891a5 2006-06-25 devnull if(r.min.x+3 >= Xrange)
602 324891a5 2006-06-25 devnull r.max.x = Xrange;
604 324891a5 2006-06-25 devnull r.max.x = r.min.x+3 + nrand(Xrange-(r.min.x+3));
606 324891a5 2006-06-25 devnull if(r.min.y+3 >= Yrange)
607 324891a5 2006-06-25 devnull r.max.y = Yrange;
609 324891a5 2006-06-25 devnull r.max.y = r.min.y+3 + nrand(Yrange-(r.min.y+3));
611 324891a5 2006-06-25 devnull assert(r.min.x >= 0);
612 324891a5 2006-06-25 devnull assert(r.max.x <= Xrange);
613 324891a5 2006-06-25 devnull assert(r.min.y >= 0);
614 324891a5 2006-06-25 devnull assert(r.max.y <= Yrange);
615 324891a5 2006-06-25 devnull /* copy from i to tmp so we have just the replicated bits */
616 324891a5 2006-06-25 devnull nb = tmp->width*sizeof(u32int)*Yrange;
617 324891a5 2006-06-25 devnull memset(tmp->data->bdata, 0, nb);
618 324891a5 2006-06-25 devnull memimagedraw(tmp, r, i, r.min, ones, r.min, SoverD);
619 324891a5 2006-06-25 devnull memmove(i->data->bdata, tmp->data->bdata, nb);
620 324891a5 2006-06-25 devnull /* i is now a non-replicated instance of the replication */
621 324891a5 2006-06-25 devnull /* replicate it by hand through tmp */
622 324891a5 2006-06-25 devnull memset(tmp->data->bdata, 0, nb);
623 324891a5 2006-06-25 devnull x = -(tilexy(r.min.x, r.max.x, 0)-r.min.x);
624 324891a5 2006-06-25 devnull for(; x<Xrange; x+=Dx(r)){
625 324891a5 2006-06-25 devnull y = -(tilexy(r.min.y, r.max.y, 0)-r.min.y);
626 324891a5 2006-06-25 devnull for(; y<Yrange; y+=Dy(r)){
627 324891a5 2006-06-25 devnull /* set r1 to instance of tile by translation */
628 324891a5 2006-06-25 devnull r1.min.x = x;
629 324891a5 2006-06-25 devnull r1.min.y = y;
630 324891a5 2006-06-25 devnull r1.max.x = r1.min.x+Dx(r);
631 324891a5 2006-06-25 devnull r1.max.y = r1.min.y+Dy(r);
632 324891a5 2006-06-25 devnull memimagedraw(tmp, r1, i, r.min, ones, r.min, SoverD);
635 324891a5 2006-06-25 devnull i->flags |= Frepl;
636 324891a5 2006-06-25 devnull i->r = r;
637 324891a5 2006-06-25 devnull i->clipr = randrect();
638 324891a5 2006-06-25 devnull /* fprint(2, "replicate [[%d %d] [%d %d]] [[%d %d][%d %d]]\n", r.min.x, r.min.y, r.max.x, r.max.y, */
639 324891a5 2006-06-25 devnull /* i->clipr.min.x, i->clipr.min.y, i->clipr.max.x, i->clipr.max.y); */
640 324891a5 2006-06-25 devnull tmp->clipr = i->clipr;
644 324891a5 2006-06-25 devnull * Mask is preset; do the rest
647 324891a5 2006-06-25 devnull verifyrectmaskrepl(int srcrepl, int maskrepl)
649 324891a5 2006-06-25 devnull Point sp, mp, tp, up;
650 324891a5 2006-06-25 devnull Rectangle dr;
651 324891a5 2006-06-25 devnull int x, y;
652 324891a5 2006-06-25 devnull Memimage *s, *m;
654 324891a5 2006-06-25 devnull /* print("verfrect %d %d\n", srcrepl, maskrepl); */
655 324891a5 2006-06-25 devnull src->flags &= ~Frepl;
656 324891a5 2006-06-25 devnull src->r = Rect(0, 0, Xrange, Yrange);
657 324891a5 2006-06-25 devnull src->clipr = src->r;
658 324891a5 2006-06-25 devnull stmp->flags &= ~Frepl;
659 324891a5 2006-06-25 devnull stmp->r = Rect(0, 0, Xrange, Yrange);
660 324891a5 2006-06-25 devnull stmp->clipr = src->r;
661 324891a5 2006-06-25 devnull mask->flags &= ~Frepl;
662 324891a5 2006-06-25 devnull mask->r = Rect(0, 0, Xrange, Yrange);
663 324891a5 2006-06-25 devnull mask->clipr = mask->r;
664 324891a5 2006-06-25 devnull mtmp->flags &= ~Frepl;
665 324891a5 2006-06-25 devnull mtmp->r = Rect(0, 0, Xrange, Yrange);
666 324891a5 2006-06-25 devnull mtmp->clipr = mask->r;
668 324891a5 2006-06-25 devnull fill(dst, dstbits);
669 324891a5 2006-06-25 devnull fill(src, srcbits);
671 324891a5 2006-06-25 devnull memmove(dst->data->bdata, dstbits, dst->width*sizeof(u32int)*Yrange);
672 324891a5 2006-06-25 devnull memmove(src->data->bdata, srcbits, src->width*sizeof(u32int)*Yrange);
673 324891a5 2006-06-25 devnull memmove(mask->data->bdata, maskbits, mask->width*sizeof(u32int)*Yrange);
675 324891a5 2006-06-25 devnull if(srcrepl){
676 324891a5 2006-06-25 devnull replicate(src, stmp);
677 324891a5 2006-06-25 devnull s = stmp;
679 324891a5 2006-06-25 devnull s = src;
680 324891a5 2006-06-25 devnull if(maskrepl){
681 324891a5 2006-06-25 devnull replicate(mask, mtmp);
682 324891a5 2006-06-25 devnull m = mtmp;
684 324891a5 2006-06-25 devnull m = mask;
686 324891a5 2006-06-25 devnull dr = randrect();
688 324891a5 2006-06-25 devnull sp.x = nrand(Xrange);
689 324891a5 2006-06-25 devnull sp.y = nrand(Yrange);
691 324891a5 2006-06-25 devnull mp.x = nrand(Xrange);
692 324891a5 2006-06-25 devnull mp.y = nrand(Yrange);
694 324891a5 2006-06-25 devnull DBG print("smalldraws\n");
695 324891a5 2006-06-25 devnull for(tp.y=sp.y,up.y=mp.y,y=dr.min.y; y<dr.max.y && tp.y<Yrange && up.y<Yrange; y++,tp.y++,up.y++)
696 324891a5 2006-06-25 devnull for(tp.x=sp.x,up.x=mp.x,x=dr.min.x; x<dr.max.x && tp.x<Xrange && up.x<Xrange; x++,tp.x++,up.x++)
697 324891a5 2006-06-25 devnull memimagedraw(dst, Rect(x, y, x+1, y+1), s, tp, m, up, SoverD);
698 324891a5 2006-06-25 devnull memmove(savedstbits, dst->data->bdata, dst->width*sizeof(u32int)*Yrange);
700 324891a5 2006-06-25 devnull memmove(dst->data->bdata, dstbits, dst->width*sizeof(u32int)*Yrange);
702 324891a5 2006-06-25 devnull DBG print("bigdraw\n");
703 324891a5 2006-06-25 devnull memimagedraw(dst, dr, src, sp, mask, mp, SoverD);
704 324891a5 2006-06-25 devnull for(y=0; y<Yrange; y++)
705 324891a5 2006-06-25 devnull checkline(dr, drawrepl(src->r, sp), drawrepl(mask->r, mp), y, srcrepl?stmp:nil, maskrepl?mtmp:nil);
709 324891a5 2006-06-25 devnull verifyrectrepl(int srcrepl, int maskrepl)
713 324891a5 2006-06-25 devnull /* mask all ones */
714 324891a5 2006-06-25 devnull memset(maskbits, 0xFF, nbytes);
715 324891a5 2006-06-25 devnull for(i=0; i<niters; i++)
716 324891a5 2006-06-25 devnull verifyrectmaskrepl(srcrepl, maskrepl);
718 324891a5 2006-06-25 devnull /* mask all zeros */
719 324891a5 2006-06-25 devnull memset(maskbits, 0, nbytes);
720 324891a5 2006-06-25 devnull for(i=0; i<niters; i++)
721 324891a5 2006-06-25 devnull verifyrectmaskrepl(srcrepl, maskrepl);
723 324891a5 2006-06-25 devnull /* random mask */
724 324891a5 2006-06-25 devnull for(i=0; i<niters; i++){
725 324891a5 2006-06-25 devnull fill(mask, maskbits);
726 324891a5 2006-06-25 devnull verifyrectmaskrepl(srcrepl, maskrepl);
731 324891a5 2006-06-25 devnull * Trivial draw implementation.
732 324891a5 2006-06-25 devnull * Color values are passed around as u32ints containing ααRRGGBB
736 324891a5 2006-06-25 devnull * Convert v, which is nhave bits wide, into its nwant bits wide equivalent.
737 324891a5 2006-06-25 devnull * Replicates to widen the value, truncates to narrow it.
740 324891a5 2006-06-25 devnull replbits(u32int v, int nhave, int nwant)
742 324891a5 2006-06-25 devnull v &= (1<<nhave)-1;
743 324891a5 2006-06-25 devnull for(; nhave<nwant; nhave*=2)
744 324891a5 2006-06-25 devnull v |= v<<nhave;
745 324891a5 2006-06-25 devnull v >>= (nhave-nwant);
746 324891a5 2006-06-25 devnull return v & ((1<<nwant)-1);
750 324891a5 2006-06-25 devnull * Decode a pixel into the uchar* values.
753 324891a5 2006-06-25 devnull pixtorgba(u32int v, uchar *r, uchar *g, uchar *b, uchar *a)
755 324891a5 2006-06-25 devnull *a = v>>24;
756 324891a5 2006-06-25 devnull *r = v>>16;
757 324891a5 2006-06-25 devnull *g = v>>8;
762 324891a5 2006-06-25 devnull * Convert uchar channels into u32int pixel.
765 324891a5 2006-06-25 devnull rgbatopix(uchar r, uchar g, uchar b, uchar a)
767 324891a5 2006-06-25 devnull return (a<<24)|(r<<16)|(g<<8)|b;
771 324891a5 2006-06-25 devnull * Retrieve the pixel value at pt in the image.
774 324891a5 2006-06-25 devnull getpixel(Memimage *img, Point pt)
776 324891a5 2006-06-25 devnull uchar r, g, b, a, *p;
777 324891a5 2006-06-25 devnull int nbits, npack, bpp;
778 324891a5 2006-06-25 devnull u32int v, c, rbits, bits;
780 324891a5 2006-06-25 devnull r = g = b = 0;
781 324891a5 2006-06-25 devnull a = ~0; /* default alpha is full */
783 324891a5 2006-06-25 devnull p = byteaddr(img, pt);
784 324891a5 2006-06-25 devnull v = p[0]|(p[1]<<8)|(p[2]<<16)|(p[3]<<24);
785 324891a5 2006-06-25 devnull bpp = img->depth;
786 324891a5 2006-06-25 devnull if(bpp<8){
788 324891a5 2006-06-25 devnull * Sub-byte greyscale pixels.
790 324891a5 2006-06-25 devnull * We want to throw away the top pt.x%npack pixels and then use the next bpp bits
791 324891a5 2006-06-25 devnull * in the bottom byte of v. This madness is due to having big endian bits
792 324891a5 2006-06-25 devnull * but little endian bytes.
794 324891a5 2006-06-25 devnull npack = 8/bpp;
795 324891a5 2006-06-25 devnull v >>= 8 - bpp*(pt.x%npack+1);
796 324891a5 2006-06-25 devnull v &= (1<<bpp)-1;
797 324891a5 2006-06-25 devnull r = g = b = replbits(v, bpp, 8);
800 324891a5 2006-06-25 devnull * General case. We need to parse the channel descriptor and do what it says.
801 324891a5 2006-06-25 devnull * In all channels but the color map, we replicate to 8 bits because that's the
802 324891a5 2006-06-25 devnull * precision that all calculations are done at.
804 324891a5 2006-06-25 devnull * In the case of the color map, we leave the bits alone, in case a color map
805 324891a5 2006-06-25 devnull * with less than 8 bits of index is used. This is currently disallowed, so it's
806 324891a5 2006-06-25 devnull * sort of silly.
809 324891a5 2006-06-25 devnull for(c=img->chan; c; c>>=8){
810 324891a5 2006-06-25 devnull nbits = NBITS(c);
811 324891a5 2006-06-25 devnull bits = v & ((1<<nbits)-1);
812 324891a5 2006-06-25 devnull rbits = replbits(bits, nbits, 8);
813 324891a5 2006-06-25 devnull v >>= nbits;
814 324891a5 2006-06-25 devnull switch(TYPE(c)){
815 324891a5 2006-06-25 devnull case CRed:
816 324891a5 2006-06-25 devnull r = rbits;
818 324891a5 2006-06-25 devnull case CGreen:
819 324891a5 2006-06-25 devnull g = rbits;
821 324891a5 2006-06-25 devnull case CBlue:
822 324891a5 2006-06-25 devnull b = rbits;
824 324891a5 2006-06-25 devnull case CGrey:
825 324891a5 2006-06-25 devnull r = g = b = rbits;
827 324891a5 2006-06-25 devnull case CAlpha:
828 324891a5 2006-06-25 devnull a = rbits;
830 324891a5 2006-06-25 devnull case CMap:
831 324891a5 2006-06-25 devnull p = img->cmap->cmap2rgb + 3*bits;
832 324891a5 2006-06-25 devnull r = p[0];
833 324891a5 2006-06-25 devnull g = p[1];
834 324891a5 2006-06-25 devnull b = p[2];
836 324891a5 2006-06-25 devnull case CIgnore:
838 324891a5 2006-06-25 devnull default:
839 324891a5 2006-06-25 devnull fprint(2, "unknown channel type %lud\n", TYPE(c));
840 324891a5 2006-06-25 devnull abort();
844 324891a5 2006-06-25 devnull return rgbatopix(r, g, b, a);
848 324891a5 2006-06-25 devnull * Return the greyscale equivalent of a pixel.
851 324891a5 2006-06-25 devnull getgrey(Memimage *img, Point pt)
853 324891a5 2006-06-25 devnull uchar r, g, b, a;
854 324891a5 2006-06-25 devnull pixtorgba(getpixel(img, pt), &r, &g, &b, &a);
855 324891a5 2006-06-25 devnull return RGB2K(r, g, b);
859 324891a5 2006-06-25 devnull * Return the value at pt in image, if image is interpreted
860 324891a5 2006-06-25 devnull * as a mask. This means the alpha channel if present, else
861 324891a5 2006-06-25 devnull * the greyscale or its computed equivalent.
864 324891a5 2006-06-25 devnull getmask(Memimage *img, Point pt)
866 324891a5 2006-06-25 devnull if(img->flags&Falpha)
867 324891a5 2006-06-25 devnull return getpixel(img, pt)>>24;
869 324891a5 2006-06-25 devnull return getgrey(img, pt);
871 324891a5 2006-06-25 devnull #undef DBG
873 324891a5 2006-06-25 devnull #define DBG if(0)
875 324891a5 2006-06-25 devnull * Write a pixel to img at point pt.
877 324891a5 2006-06-25 devnull * We do this by reading a 32-bit little endian
878 324891a5 2006-06-25 devnull * value from p and then writing it back
879 324891a5 2006-06-25 devnull * after tweaking the appropriate bits. Because
880 324891a5 2006-06-25 devnull * the data is little endian, we don't have to worry
881 324891a5 2006-06-25 devnull * about what the actual depth is, as long as it is
882 324891a5 2006-06-25 devnull * less than 32 bits.
885 324891a5 2006-06-25 devnull putpixel(Memimage *img, Point pt, u32int nv)
887 324891a5 2006-06-25 devnull uchar r, g, b, a, *p, *q;
888 324891a5 2006-06-25 devnull u32int c, mask, bits, v;
889 324891a5 2006-06-25 devnull int bpp, sh, npack, nbits;
891 324891a5 2006-06-25 devnull pixtorgba(nv, &r, &g, &b, &a);
893 324891a5 2006-06-25 devnull p = byteaddr(img, pt);
894 324891a5 2006-06-25 devnull v = p[0]|(p[1]<<8)|(p[2]<<16)|(p[3]<<24);
895 324891a5 2006-06-25 devnull bpp = img->depth;
896 324891a5 2006-06-25 devnull DBG print("v %.8lux...", v);
897 324891a5 2006-06-25 devnull if(bpp < 8){
899 324891a5 2006-06-25 devnull * Sub-byte greyscale pixels. We need to skip the leftmost pt.x%npack pixels,
900 324891a5 2006-06-25 devnull * which is equivalent to skipping the rightmost npack - pt.x%npack - 1 pixels.
902 324891a5 2006-06-25 devnull npack = 8/bpp;
903 324891a5 2006-06-25 devnull sh = bpp*(npack - pt.x%npack - 1);
904 324891a5 2006-06-25 devnull bits = RGB2K(r,g,b);
905 324891a5 2006-06-25 devnull DBG print("repl %lux 8 %d = %lux...", bits, bpp, replbits(bits, 8, bpp));
906 324891a5 2006-06-25 devnull bits = replbits(bits, 8, bpp);
907 324891a5 2006-06-25 devnull mask = (1<<bpp)-1;
908 324891a5 2006-06-25 devnull DBG print("bits %lux mask %lux sh %d...", bits, mask, sh);
909 324891a5 2006-06-25 devnull mask <<= sh;
910 324891a5 2006-06-25 devnull bits <<= sh;
911 324891a5 2006-06-25 devnull DBG print("(%lux & %lux) | (%lux & %lux)", v, ~mask, bits, mask);
912 324891a5 2006-06-25 devnull v = (v & ~mask) | (bits & mask);
913 324891a5 2006-06-25 devnull } else {
915 324891a5 2006-06-25 devnull * General case. We need to parse the channel descriptor again.
918 324891a5 2006-06-25 devnull for(c=img->chan; c; c>>=8){
919 324891a5 2006-06-25 devnull nbits = NBITS(c);
920 324891a5 2006-06-25 devnull switch(TYPE(c)){
921 324891a5 2006-06-25 devnull case CRed:
922 324891a5 2006-06-25 devnull bits = r;
924 324891a5 2006-06-25 devnull case CGreen:
925 324891a5 2006-06-25 devnull bits = g;
927 324891a5 2006-06-25 devnull case CBlue:
928 324891a5 2006-06-25 devnull bits = b;
930 324891a5 2006-06-25 devnull case CGrey:
931 324891a5 2006-06-25 devnull bits = RGB2K(r, g, b);
933 324891a5 2006-06-25 devnull case CAlpha:
934 324891a5 2006-06-25 devnull bits = a;
936 324891a5 2006-06-25 devnull case CIgnore:
937 324891a5 2006-06-25 devnull bits = 0;
939 324891a5 2006-06-25 devnull case CMap:
940 324891a5 2006-06-25 devnull q = img->cmap->rgb2cmap;
941 324891a5 2006-06-25 devnull bits = q[(r>>4)*16*16+(g>>4)*16+(b>>4)];
943 324891a5 2006-06-25 devnull default:
944 324891a5 2006-06-25 devnull SET(bits);
945 324891a5 2006-06-25 devnull fprint(2, "unknown channel type %lud\n", TYPE(c));
946 324891a5 2006-06-25 devnull abort();
949 324891a5 2006-06-25 devnull DBG print("repl %lux 8 %d = %lux...", bits, nbits, replbits(bits, 8, nbits));
950 324891a5 2006-06-25 devnull if(TYPE(c) != CMap)
951 324891a5 2006-06-25 devnull bits = replbits(bits, 8, nbits);
952 324891a5 2006-06-25 devnull mask = (1<<nbits)-1;
953 324891a5 2006-06-25 devnull DBG print("bits %lux mask %lux sh %d...", bits, mask, sh);
954 324891a5 2006-06-25 devnull bits <<= sh;
955 324891a5 2006-06-25 devnull mask <<= sh;
956 324891a5 2006-06-25 devnull v = (v & ~mask) | (bits & mask);
957 324891a5 2006-06-25 devnull sh += nbits;
960 324891a5 2006-06-25 devnull DBG print("v %.8lux\n", v);
961 324891a5 2006-06-25 devnull p[0] = v;
962 324891a5 2006-06-25 devnull p[1] = v>>8;
963 324891a5 2006-06-25 devnull p[2] = v>>16;
964 324891a5 2006-06-25 devnull p[3] = v>>24;
966 324891a5 2006-06-25 devnull #undef DBG
968 324891a5 2006-06-25 devnull #define DBG if(0)
970 324891a5 2006-06-25 devnull drawonepixel(Memimage *dst, Point dp, Memimage *src, Point sp, Memimage *mask, Point mp)
972 324891a5 2006-06-25 devnull uchar m, M, sr, sg, sb, sa, sk, dr, dg, db, da, dk;
974 324891a5 2006-06-25 devnull pixtorgba(getpixel(dst, dp), &dr, &dg, &db, &da);
975 324891a5 2006-06-25 devnull pixtorgba(getpixel(src, sp), &sr, &sg, &sb, &sa);
976 324891a5 2006-06-25 devnull m = getmask(mask, mp);
977 324891a5 2006-06-25 devnull M = 255-(sa*m)/255;
979 324891a5 2006-06-25 devnull DBG print("dst %x %x %x %x src %x %x %x %x m %x = ", dr,dg,db,da, sr,sg,sb,sa, m);
980 324891a5 2006-06-25 devnull if(dst->flags&Fgrey){
982 324891a5 2006-06-25 devnull * We need to do the conversion to grey before the alpha calculation
983 324891a5 2006-06-25 devnull * because the draw operator does this, and we need to be operating
984 324891a5 2006-06-25 devnull * at the same precision so we get exactly the same answers.
986 324891a5 2006-06-25 devnull sk = RGB2K(sr, sg, sb);
987 324891a5 2006-06-25 devnull dk = RGB2K(dr, dg, db);
988 324891a5 2006-06-25 devnull dk = (sk*m + dk*M)/255;
989 324891a5 2006-06-25 devnull dr = dg = db = dk;
990 324891a5 2006-06-25 devnull da = (sa*m + da*M)/255;
993 324891a5 2006-06-25 devnull * True color alpha calculation treats all channels (including alpha)
994 324891a5 2006-06-25 devnull * the same. It might have been nice to use an array, but oh well.
996 324891a5 2006-06-25 devnull dr = (sr*m + dr*M)/255;
997 324891a5 2006-06-25 devnull dg = (sg*m + dg*M)/255;
998 324891a5 2006-06-25 devnull db = (sb*m + db*M)/255;
999 324891a5 2006-06-25 devnull da = (sa*m + da*M)/255;
1002 324891a5 2006-06-25 devnull DBG print("%x %x %x %x\n", dr,dg,db,da);
1003 324891a5 2006-06-25 devnull putpixel(dst, dp, rgbatopix(dr, dg, db, da));
