From http://www.jwz.org/xscreensaver/xscreensaver-5.30.tar.gz
[xscreensaver] / hacks / abstractile.c
1 /*
2  * Copyright (c) 2004-2009 Steve Sundstrom
3  *
4  * Permission to use, copy, modify, distribute, and sell this software and its
5  * documentation for any purpose is hereby granted without fee, provided that
6  * the above copyright notice appear in all copies and that both that
7  * copyright notice and this permission notice appear in supporting
8  * documentation.  No representations are made about the suitability of this
9  * software for any purpose.  It is provided "as is" without express or
10  * implied warranty.
11  */
12
13 #include "screenhack.h"
14 #include "colors.h"
15 #include "hsv.h"
16 #include <stdio.h>
17 #include <math.h>
18 #include <sys/time.h>
19 /*#include <sys/utsname.h>*/
20
21 #define MODE_CREATE             0   /* init, create, then finish sleep */
22 #define MODE_ERASE              1   /* erase, then reset colors */
23 #define MODE_DRAW               2
24
25 #define DIR_NONE                0
26 #define DIR_UP                  1
27 #define DIR_DOWN                2
28 #define DIR_LEFT                3
29 #define DIR_RIGHT               4
30
31 #define LINE_FORCE              1
32 #define LINE_NEW                2
33 #define LINE_BRIN               3
34 #define LINE_BROUT              4
35
36 #define PT_UL                   0
37 #define PT_MP                   1
38 #define PT_LR                   2
39 #define PT_NL                   3
40
41 #define D3D_NONE                0
42 #define D3D_BLOCK               1
43 #define D3D_NEON                2
44 #define D3D_TILED               3
45
46 #define TILE_RANDOM             0
47 #define TILE_FLAT               1
48 #define TILE_THIN               2
49 #define TILE_OUTLINE            3
50 #define TILE_BLOCK              4
51 #define TILE_NEON               5
52 #define TILE_TILED              6
53
54 #define BASECOLORS              30
55 #define SHADES                  12
56 #define MAXCOLORS               40
57 #define LAYERS                  4
58 #define PATTERNS                40
59 #define SHAPES                  18
60 #define DRAWORDERS              40
61 #define COLORMAPS               20
62 #define WAVES                   6
63 #define STRETCHES               8
64
65 struct lineStruct {
66     unsigned int x, y, len, obj, color, ndol;
67     int deo;
68     Bool hv;
69 };
70
71 struct gridStruct {
72     unsigned int line, hl, hr, vu, vd, dhl, dhr, dvu, dvd;
73 };
74
75 /* basically the same as global variables, but used to keep them in a bucket
76    and pass them around easier like the original C++ implementation */
77 struct state {
78   /* window values */
79   Display *display;
80   Window window;
81   XWindowAttributes xgwa;
82   GC fgc, bgc;
83   XColor colors[255];
84
85   /* memory values */
86   struct lineStruct *dline, *eline;
87   struct gridStruct *grid;
88   unsigned int *zlist, *fdol;
89   Bool *odi;
90     /* draw, erase, fill, init, line, object, z indexes */
91   unsigned int di, ei, fi, ii, bi, li, eli, oi, zi;
92   /* size variables */
93   unsigned int gridx, gridy, gridn;                     /* grid size */
94   int lwid, bwid, swid;/* line width, background width, shadow width */
95   int narray, max_wxh;
96   int elwid, elpu, egridx, egridy; /* for now */
97   /* fill variables */
98   int bnratio;                 /* ratio of branch lines to new lines */
99   int maxlen;                              /* maximum length of line */
100   int forcemax;                  /* make line be max possible length */
101   int olen;                           /* open length set by findopen */
102   int bln;          /* blocking line number set by findopen, -1=edge */
103   /* color variables */
104   int ncolors;                        /* number of colors for screen */
105   int shades;
106   int rco[MAXCOLORS];           /* random ordering of colors for deo */
107   int cmap;
108   int layers;
109   Bool newcols;      /* can we create new colormaps with each screen */
110   /* draw variables */
111   int dmap, emap;  /* pattern by which line draw order is determined */
112   int dvar, evar;             /* random number added to .deo to vary */
113   int ddir, edir;      /* draw/erase in forward direction or reverse */
114   int lpu;            /* lines drawn per update used to adjust speed */
115   int d3d;
116   int round;
117   int outline;
118   /* layered draw variables */
119   int pattern[LAYERS], shape[LAYERS], mix[LAYERS];
120   int csw[LAYERS], wsx[LAYERS], wsy[LAYERS], sec[LAYERS];
121   int cs1[LAYERS], cs2[LAYERS], cs3[LAYERS]; int cs4[LAYERS];
122   int wave[LAYERS], waveh[LAYERS], wavel[LAYERS];
123   int rx1[LAYERS], rx2[LAYERS], rx3[LAYERS];
124   int ry1[LAYERS], ry2[LAYERS], ry3[LAYERS];
125   /* misc variables */
126   int mode, sleep, speed, tile, dialog;
127   Bool grid_full, resized;
128   struct timeval time;
129 };
130
131 static int
132 _min(int a, int b)
133 {
134   if (a<=b)
135     return(a);
136   return(b);
137 }
138
139 static int
140 _max(int a, int b)
141 {
142   if (a>=b)
143     return(a);
144   return(b);
145 }
146
147 static int
148 _dist(struct state *st, int x1, int x2, int y1, int y2, int s)
149 {
150   double xd=x1-x2;
151   double yd=y1-y2;
152   switch(s) {
153     case 0:
154       return((int)sqrt(xd*xd+yd*yd));
155     case 1:
156       return((int)sqrt(xd*xd*st->cs1[0]*2+yd*yd));
157     case 2:
158       return((int)sqrt(xd*xd+yd*yd*st->cs2[0]*2));
159     default:
160       return((int)sqrt(xd*xd*st->cs1[0]/st->cs2[0]+yd*yd*st->cs3[0]/st->cs4[0]));
161   }
162 }
163
164 static int
165 _wave(struct state *st, int x, int h, int l, int wave)
166 {
167   l+=1;
168   switch(wave) {
169     case 0:                                         /* cos wave*/
170       return((int)(cos((double)x*M_PI/l)*h));
171     case 1:                                      /* double wave*/
172     case 2:                                      /* double wave*/
173       return((int)(cos((double)x*M_PI/l)*h)+(int)(sin((double)x*M_PI/l/st->cs1[1])*h));
174     case 3:                                         /* zig zag */
175       return(abs((x%(l*2)-l))*h/l);
176     case 4:                                   /* giant zig zag */
177       return(abs((x%(l*4)-l*2))*h*3/l);
178     case 5:                                        /* sawtooth */
179       return((x%(l))*h/l);
180     default:                                        /* no wave */
181       return(0);
182   }
183 }
184
185 static int
186 _triangle(struct state *st, int x, int y, int rx, int ry, int t)
187 {
188   switch(t) {
189     case 1:
190       return(_min(_min(x+y+rx-(st->gridx/2),st->gridx-x+y),(st->gridy-y+(ry/2))*3/2));
191     case 2:
192       return(_min(_min(x-rx,y-ry),(rx+ry-x-y)*2/3));
193     case 3:
194       return(_min(_min(st->gridx-x-rx,y-ry),(rx+ry-st->gridx+x-y)*2/3));
195     case 4:
196       return(_min(_min(x-rx,st->gridy-y-ry),(rx+ry-x-st->gridy+y)*2/3));
197   }
198   return(_min(_min(st->gridx-x-rx,st->gridy-y-ry),(rx+ry-st->gridx+x-st->gridy+y)*2/3));
199 }
200
201 static void
202 _init_zlist(struct state *st)
203 {
204   unsigned int tmp, y, z;
205
206   st->gridx=st->xgwa.width/st->lwid;
207   st->gridy=st->xgwa.height/st->lwid;
208   st->gridn=st->gridx*st->gridy;
209   /* clear grid */
210   for (z=0; z<st->gridn; z++) {
211     st->grid[z].line=st->grid[z].hl=st->grid[z].hr=st->grid[z].vu=st->grid[z].vd=st->grid[z].dhl=st->grid[z].dhr=st->grid[z].dvu=st->grid[z].dvd=0;
212     st->zlist[z]=z;
213   }
214   /* rather than pull x,y points randomly and wait to hit final empy cells a
215      list of all points is created and mixed so empty cells do get hit last */
216   for (z=0; z<st->gridn; z++) {
217     y=random()%st->gridn;
218     tmp=st->zlist[y];
219     st->zlist[y]=st->zlist[z];
220     st->zlist[z]=tmp;
221   }
222 }
223
224 static void
225 make_color_ramp_rgb (Screen *screen, Visual *visual, Colormap cmap,
226     int r1, int g1, int b1,  int r2, int g2, int b2,
227     XColor *colors, int *ncolorsP, Bool closed_p)
228 {
229     int h1, h2;
230     double s1, s2, v1, v2;
231     rgb_to_hsv(r1, g1, b1, &h1, &s1, &v1);
232     rgb_to_hsv(r2, g2, b2, &h2, &s2, &v2);
233     make_color_ramp(screen, visual, cmap, h1, s1, v1, h2, s2, v2,
234         colors, ncolorsP, False, True, 0);
235 }
236
237
238 static void
239 _init_colors(struct state *st)
240 {
241   int col[BASECOLORS];
242   int c1, c2, c3, h1, h2, h3;
243   int r1, g1, b1, r2, g2, b2, r3, g3, b3;
244   double s1, s2, s3, v1, v2, v3;
245   XColor tmp_col1[16], tmp_col2[16], tmp_col3[16];
246
247   unsigned short basecol[BASECOLORS][3]={
248     /* 0  dgray */    {0x3333,0x3333,0x3333},
249     /* 1  dbrown */   {0x6666,0x3333,0x0000},
250     /* 2  dred */     {0x9999,0x0000,0x0000},
251     /* 3  orange */   {0xFFFF,0x6666,0x0000},
252     /* 4  gold */     {0xFFFF,0xCCCC,0x0000},
253     /* 5  olive */    {0x6666,0x6666,0x0000},
254     /* 6  ivy */      {0x0000,0x6666,0x0000},
255     /* 7  dgreen */   {0x0000,0x9999,0x0000},
256     /* 8  bluegray */ {0x3333,0x6666,0x6666},
257     /* 9  dblue */    {0x0000,0x0000,0x9999},
258     /* 10 blue */     {0x3333,0x3333,0xFFFF},
259     /* 11 dpurple */  {0x6666,0x0000,0xCCCC},
260     /* 12 purple */   {0x6666,0x3333,0xFFFF},
261     /* 13 violet */   {0x9999,0x3333,0x9999},
262     /* 14 magenta */  {0xCCCC,0x3333,0xCCCC},
263     /* lights */
264     /* 15 gray */     {0x3333,0x3333,0x3333},
265     /* 16 brown */    {0x9999,0x6666,0x3333},
266     /* 17 tan */      {0xCCCC,0x9999,0x3333},
267     /* 18 red */      {0xFFFF,0x0000,0x0000},
268     /* 19 lorange */  {0xFFFF,0x9999,0x0000},
269     /* 20 yellow */   {0xFFFF,0xFFFF,0x0000},
270     /* 21 lolive */   {0x9999,0x9999,0x0000},
271     /* 22 green */    {0x3333,0xCCCC,0x0000},
272     /* 23 lgreen */   {0x3333,0xFFFF,0x3333},
273     /* 24 cyan */     {0x0000,0xCCCC,0xCCCC},
274     /* 25 sky */      {0x3333,0xFFFF,0xFFFF},
275     /* 26 marine */   {0x3333,0x6666,0xFFFF},
276     /* 27 lblue */    {0x3333,0xCCCC,0xFFFF},
277     /* 28 lpurple */  {0x9999,0x9999,0xFFFF},
278     /* 29 pink */     {0xFFFF,0x9999,0xFFFF}};
279
280   if (st->d3d) {
281     st->shades = (st->d3d==D3D_TILED) ? 5 : st->lwid/2+1;
282     st->ncolors=4+random()%4;
283     if (st->cmap>0) {                      /* tint the basecolors a bit */
284       for (c1=0; c1<BASECOLORS; c1++)
285         for (c2=0; c2<2; c2++)
286           if (!basecol[c1][c2]) {
287             basecol[c1][c2]+=random()%16000;
288           } else if (basecol[c1][c2]==0xFFFF) {
289             basecol[c1][c2]-=random()%16000;
290           } else {
291             basecol[c1][c2]-=8000;
292             basecol[c1][c2]+=random()%16000;
293           }
294     }
295     switch(st->cmap%4) {
296       case 0:                                            /* all */
297         for (c1=0; c1<st->ncolors; c1++)
298           col[c1]=random()%BASECOLORS;
299         break;
300       case 1:                                          /* darks */
301         for (c1=0; c1<st->ncolors; c1++)
302           col[c1]=random()%15;
303         break;
304       case 2:                                   /* semi consecutive darks */
305         col[0]=random()%15;
306         for (c1=1; c1<st->ncolors; c1++)
307           col[c1]=(col[c1-1]+1+random()%2)%15;
308         break;
309       case 3:                                   /* consecutive darks */
310         col[0]=random()%(15-st->ncolors);
311         for (c1=1; c1<st->ncolors; c1++)
312           col[c1]=col[c1-1]+1;
313         break;
314     }
315     for (c1=0; c1<st->ncolors; c1++) {
316       /* adjust colors already set */
317       for (h1=c1*st->shades-1; h1>=0; h1--)
318         st->colors[h1+st->shades]=st->colors[h1];
319       make_color_ramp_rgb(st->xgwa.screen, st->xgwa.visual, st->xgwa.colormap,
320         basecol[col[c1]][0], basecol[col[c1]][1], basecol[col[c1]][2],
321         0xFFFF, 0xFFFF, 0xFFFF, st->colors, &st->shades,
322         False);
323     }
324     return;
325   }
326   /* not 3d */
327   st->shades=1;
328   if (st->cmap%2) {                                  /* basecolors */
329     if (random()%3) {
330       c1=random()%15;
331       c2=(c1+3+(random()%5))%15;
332       c3=(c2+3+(random()%5))%15;
333     } else {
334       c1=random()%BASECOLORS;
335       c2=(c1+5+(random()%10))%BASECOLORS;
336       c3=(c2+5+(random()%10))%BASECOLORS;
337     }
338     r1=basecol[c1][0];
339     g1=basecol[c1][1];
340     b1=basecol[c1][2];
341     r2=basecol[c2][0];
342     g2=basecol[c2][1];
343     b2=basecol[c2][2];
344     r3=basecol[c3][0];
345     g3=basecol[c3][1];
346     b3=basecol[c3][2];
347   } else {                                             /* random rgb's */
348     r1=random()%65535;
349     g1=random()%65535;
350     b1=random()%65535;
351     r2=(r1+16384+random()%32768)%65535;
352     g2=(g1+16384+random()%32768)%65535;
353     b2=(b1+16384+random()%32768)%65535;
354     r3=(r2+16384+random()%32768)%65535;
355     g3=(g2+16384+random()%32768)%65535;
356     b3=(b2+16384+random()%32768)%65535;
357  }
358  switch(st->cmap) {
359     case 0:                           /* make_color_ramp color->color */
360     case 1:
361     case 2:                           /* make_color_ramp color->white */
362     case 3:
363       st->ncolors=5+random()%5;
364       if (st->cmap>1)
365         r2=g2=b2=0xFFFF;
366       make_color_ramp_rgb(st->xgwa.screen, st->xgwa.visual, st->xgwa.colormap,
367         r1, g1, b1, r2, g2, b2,
368         st->colors, &st->ncolors, random()%2);
369       break;
370     case 4:                                /* 3 color make_color_loop */
371     case 5:
372     case 6:
373     case 7:
374       st->ncolors=8+random()%12;
375       rgb_to_hsv(r1, g1, b1, &h1, &s1, &v1);
376       rgb_to_hsv(r2, g2, b2, &h2, &s2, &v2);
377       rgb_to_hsv(r3, g3, b3, &h3, &s3, &v3);
378
379       make_color_loop(st->xgwa.screen, st->xgwa.visual, st->xgwa.colormap,
380         h1, s1, v1, h2, s2, v2, h3, s3, v3,
381         st->colors, &st->ncolors, True, False);
382       break;
383     case 8:                                            /* random smooth */
384     case 9:
385       st->ncolors=(random()%4)*6+12;
386       make_smooth_colormap (st->xgwa.screen, st->xgwa.visual,
387         st->xgwa.colormap, st->colors, &st->ncolors,
388         True, False, True);
389       break;
390     case 10:                                                 /* rainbow */
391       st->ncolors=(random()%4)*6+12;
392       make_uniform_colormap (st->xgwa.screen, st->xgwa.visual,
393         st->xgwa.colormap, st->colors, &st->ncolors,
394         True, False, True);
395       break;
396     case 11:                                     /* dark to light blend */
397     case 12:
398     case 13:
399     case 14:
400       st->ncolors=7;
401       make_color_ramp_rgb(st->xgwa.screen, st->xgwa.visual,  st->xgwa.colormap,
402         r1, g1, b1, 0xFFFF, 0xFFFF, 0xFFFF,
403         tmp_col1, &st->ncolors, False);
404       make_color_ramp_rgb(st->xgwa.screen, st->xgwa.visual,  st->xgwa.colormap,
405         r2, g2, b2, 0xFFFF, 0xFFFF, 0xFFFF,
406         tmp_col2, &st->ncolors, False);
407       if (st->cmap<13) {
408         for(c1=0; c1<=4; c1++) {
409            st->colors[c1*2]=tmp_col1[c1];
410            st->colors[c1*2+1]=tmp_col2[c1];
411         }
412         st->ncolors=10;
413       } else {
414         make_color_ramp_rgb(st->xgwa.screen, st->xgwa.visual,  st->xgwa.colormap,
415           r3, g3, b3, 0xFFFF, 0xFFFF, 0xFFFF,
416           tmp_col3, &st->ncolors, False);
417         for(c1=0; c1<=4; c1++) {
418            st->colors[c1*3]=tmp_col1[c1];
419            st->colors[c1*3+1]=tmp_col2[c1];
420            st->colors[c1*3+2]=tmp_col3[c1];
421         }
422         st->ncolors=15;
423       }
424       break;
425     default:                                                  /* random */
426       st->ncolors=(random()%4)*6+12;
427       make_random_colormap (st->xgwa.screen, st->xgwa.visual,
428         st->xgwa.colormap, st->colors, &st->ncolors,
429         False, True, False, True);
430       break;
431   }
432
433   /* set random color order for drawing and erasing */
434   for (c1=0; c1<MAXCOLORS; c1++)
435     st->rco[c1]=c1;
436   for (c1=0; c1<MAXCOLORS; c1++) {
437     c3=random()%MAXCOLORS;
438     c2=st->rco[c1];
439     st->rco[c1]=st->rco[c3];
440     st->rco[c3]=c2;
441   }
442 }
443
444 static int _comparedeo(const void *i, const void *j)
445 {
446         struct lineStruct *h1, *h2;
447
448         h1=(struct lineStruct *)i;
449         h2=(struct lineStruct *)j;
450         if (h1->deo > h2->deo)
451                 return(1);
452         if (h1->deo < h2->deo)
453                 return(-1);
454         return(0);
455 }
456
457 static int
458 _hv(struct state *st, int x, int y, int d1, int d2, int pn, Bool de)
459 {
460   int v1, v2, r;
461
462   switch (d1) {
463     case 0:
464       v1 = (de) ? st->egridx-x : st->gridx-x;
465       break;
466     case 1:
467       v1 = y;
468       break;
469     case 2:
470       v1 = x;
471       break;
472     default:
473       v1 = (de) ? st->egridy-y : st->gridy-y;
474       break;
475   }
476   switch (d2) {
477     case 0:
478       v2 = (de) ? st->egridx-x : st->gridx-x;
479       break;
480     case 1:
481       v2 = y;
482       break;
483     case 2:
484       v2 = x;
485       break;
486     default:
487       v2 = (de) ? st->egridy-y : st->gridy-y;
488       break;
489   }
490   r = (de) ? (st->dline[st->li].hv) ? (v1+10000)*pn : (v2+10000)*-pn :
491     (st->eline[st->li].hv) ? (v1+10000)*pn : (v2+10000)*-pn;
492   return(r);
493 }
494         
495 static int
496 _getdeo(struct state *st, int x, int y, int map, int de)
497 {
498   int cr;
499   switch(map) {
500     case 0:                                        /* horizontal one side */
501       return(x);
502     case 1:                                          /* vertical one side */
503       return(y);
504     case 2:                                        /* horizontal two side */
505       return(_min(x,st->gridx-x)+1);
506     case 3:                                          /* vertical two side */
507       return(_min(y,st->gridy-y)+1);
508     case 4:                                                     /* square */
509       return(_max(abs(x-st->rx3[de]),abs(y-st->ry3[de]))+1);
510     case 5:                                                /* two squares */
511       return(_min(_max(abs(x-(st->rx3[de]/2)),abs(y-st->ry3[de])),_max(abs(x-(st->gridx-(st->rx2[de]/2))),abs(y-st->ry2[de])))+1);
512     case 6:                                       /* horizontal rectangle */
513       return(_max(abs(x-st->rx3[de]),abs(y-(st->ry3[de]))*st->cs1[de])+1);
514     case 7:                                         /* vertical rectangle */
515       return(_max(abs(x-st->rx3[de])*st->cs1[de],abs(y-(st->ry3[de])))+1);
516     case 8:                                                    /* + cross */
517       return(_min(abs(x-st->rx3[de]),abs(y-(st->ry3[de])))+1);
518     case 9:                                                   /* diagonal */
519       return((x*3/4+y)+1);
520     case 10:                                         /* opposite diagonal */
521       return((x*3/4+st->gridy-y)+1);
522     case 11:                                                   /* diamond */
523       return((abs(x-st->rx3[de])+abs(y-st->ry3[de]))/2+1);
524     case 12:                                              /* two diamonds */
525       return(_min(abs(x-(st->rx3[de]/2))+abs(y-st->ry3[de]),abs(x-(st->gridx-(st->rx2[de]/2)))+abs(y-st->ry2[de]))/2+1);
526     case 13:                                                    /* circle */
527       return(_dist(st,x,st->rx3[de],y,st->ry3[de],0)+1);
528     case 14:                                        /* horizontal ellipse */
529       return(_dist(st,x,st->rx3[de],y,st->ry3[de],1)+1);
530     case 15:                                          /* vertical ellipse */
531       return(_dist(st,x,st->rx3[de],y,st->ry3[de],2)+1);
532     case 16:                                               /* two circles */
533       return(_min(_dist(st,x,st->rx3[de]/2,y,st->ry3[de],0),_dist(st,x,st->gridx-(st->rx2[de]/2),y,st->ry2[de],0))+1);
534     case 17:                                  /* horizontal straight wave */
535       return(x+_wave(st,st->gridy+y,st->csw[0]*st->cs1[0],st->csw[0]*st->cs2[0],st->wave[de]));
536     case 18:                                    /* vertical straight wave */
537       return(y+_wave(st,st->gridx+x,st->csw[0]*st->cs1[0],st->csw[0]*st->cs2[0],st->wave[de]));
538     case 19:                                     /* horizontal wavey wave */
539       return(x+_wave(st,st->gridy+y+((x/5)*st->edir),st->csw[de]*st->cs1[de],st->csw[de]*st->cs2[de],st->wave[de])+1);
540     case 20:                                       /* vertical wavey wave */
541       return(y+_wave(st,st->gridx+x+((y/5)*st->edir),st->csw[de]*st->cs1[de],st->csw[de]*st->cs2[de],st->wave[de])+1);
542 /* no d3d for 21,22 */
543     case 21:                                  /* simultaneous directional */
544       return(_hv(st,x,y,st->cs1[0]%2,st->cs2[0]%2,1,de));
545     case 22:                                       /* reverse directional */
546       return(_hv(st,x,y,st->cs1[0]%2,st->cs2[0]%2,-1,de));
547     case 23:                                                    /* length */
548       if (de)
549         return(st->dline[st->li].len*1000+random()%5000);
550       else
551         return(st->eline[st->li].len*1000+random()%5000);
552     case 24:                                                    /* object */
553     case 25:
554     case 26:
555     case 27:
556       if (de)
557         return(st->dline[st->li].obj*100);
558       else
559         return(st->eline[st->li].obj*100);
560     default:                                                     /* color */
561       cr = (de) ? st->dline[st->li].color : st->eline[st->li].color;
562       if (map<34) cr=st->rco[cr];
563       if ((map%6<4) || (de))  {                               /* by color */
564         cr*=1000;
565         cr+=random()%1000;
566       } else if (map%6==4) {                      /* by color horizontaly */
567         cr*=st->gridx;
568         cr+=(x+random()%(st->gridx/2));
569       } else {                                     /* by color vertically */
570         cr*=st->gridy;
571         cr+=(y+random()%(st->gridy/2));
572       }
573       return(cr);
574   }
575   return(1);
576 }
577
578 static void
579 _init_screen(struct state *st)
580 {
581   int nstr, x;
582   struct lineStruct *tmp;
583
584   /* malloc memory in case of resize */
585   if (st->resized) {
586     st->max_wxh=st->xgwa.width*st->xgwa.height;
587     if (st->dline!=NULL)
588       free(st->dline);
589     if (st->eline!=NULL)
590       free(st->eline);
591     if (st->grid!=NULL)
592       free(st->grid);
593     if (st->zlist!=NULL)
594       free(st->zlist);
595     if (st->fdol!=NULL)
596       free(st->fdol);
597     if (st->odi!=NULL)
598       free(st->odi);
599     st->narray = (st->xgwa.width+1)*(st->xgwa.height+1)/4+1;
600     st->dline = calloc(st->narray, sizeof(struct lineStruct));
601     st->eline = calloc(st->narray, sizeof(struct lineStruct));
602     st->grid = calloc(st->narray, sizeof(struct gridStruct));
603     st->zlist = calloc(st->narray, sizeof(unsigned int));
604     st->fdol = calloc(st->narray, sizeof(unsigned int));
605     st->odi = calloc(st->narray, sizeof(Bool));
606     if ((st->dline == NULL) || (st->eline == NULL) ||
607       (st->grid == NULL) || (st->zlist == NULL) ||
608       (st->fdol == NULL) || (st->odi == NULL)) {
609       fprintf(stderr, "not enough memory\n");
610       exit(1);
611     }
612     st->dialog = (st->xgwa.width<500) ? 1 : 0;
613     st->resized=False;
614   }
615   if (st->ii) {
616     /* swap st->dline and st->eline pointers to resort and erase */
617     tmp=st->eline;
618     st->eline=st->dline;
619     st->dline=tmp;
620     st->eli=st->li;
621     st->elwid=st->lwid;
622     st->elpu=st->lpu;
623     st->egridx=st->gridx;
624     st->egridy=st->gridy;
625
626     /* create new erase order */
627     for (st->li=1; st->li<=st->eli; st->li++)
628       st->eline[st->li].deo=(_getdeo(st,st->eline[st->li].x,st->eline[st->li].y,st->emap,0) + (random()%st->evar) + (random()%st->evar))*st->edir;
629     qsort(st->eline, st->eli+1, sizeof(struct lineStruct), _comparedeo);
630   }
631   st->ii++;
632
633   /* clear arrays and other counters */
634   st->di=st->ei=st->fi=st->li=st->oi=st->zi=0;
635   st->grid_full=False;
636   /* li starts from 1 */
637   st->dline[0].x=st->dline[0].y=st->dline[0].len=0;
638   /* to keep it first after sorting so di is never null */
639   st->dline[0].deo=-999999999;
640
641   /* set random screen variables */
642   st->lwid = (st->ii==1) ? 3 : 2+((random()%6)%4);
643   st->d3d = ((st->tile==TILE_FLAT) || (st->tile==TILE_THIN) ||
644     (st->tile==TILE_OUTLINE)) ? D3D_NONE :
645     (st->tile==TILE_BLOCK) ? D3D_BLOCK :
646     (st->tile==TILE_NEON) ? D3D_NEON :
647     (st->tile==TILE_TILED) ? D3D_TILED :
648     /* force TILE_D3D on first screen to properly load all shades */
649     ((st->ii==1) && (!st->newcols)) ? D3D_TILED : (random()%5)%4;
650 /* st->d3d=D3D_BLOCK; st->lwid=2; */
651   st->outline = (st->tile==TILE_OUTLINE) ? 1 :
652      ((st->tile!=TILE_RANDOM) || (random()%5)) ? 0 : 1;
653   st->round = (st->d3d==D3D_NEON) ? 1 :
654     ((st->d3d==D3D_BLOCK) || (st->outline) || (random()%6)) ? 0 : 1;
655   if ((st->d3d) || (st->outline) || (st->round))
656     st->lwid+=2;
657   if ((!st->d3d) && (!st->round) && (!st->outline) && (st->lwid>3))
658     st->lwid-=2;
659   if (st->d3d==D3D_TILED)
660     st->lwid++;
661   if (st->tile==TILE_THIN)
662     st->lwid=2;
663
664   _init_zlist(st);
665
666   st->maxlen=(st->lwid>6) ? 2+(random()%4) :
667                 (st->lwid>4) ? 2+(random()%8)%6 :
668                 (st->lwid>2) ? 2+(random()%12)%8 : 2+(random()%15)%10;
669   st->bnratio = 4+(random()%4)+(random()%4);
670   st->forcemax = (random()%6) ? 0 : 1;
671
672   if ((st->ii==1) || (st->newcols))
673     _init_colors(st);
674
675   st->dmap = (st->emap+5+(random()%5))%DRAWORDERS;
676
677   st->dmap=20+random()%20;
678
679   st->dvar = (st->dmap>22) ? 100 : 10+(st->csw[0]*(random()%5));
680   st->ddir= (random()%2) ? 1 : -1;
681
682   st->emap = (st->dmap+10+(random()%10))%20;
683   st->evar = (st->emap>22) ? 100 : 10+(st->csw[0]*(random()%5));
684   st->edir= (random()%2) ? 1 : -1;
685
686   st->layers= (random()%2) ? 2 : (random()%2) ? 1 : (random()%2) ? 3 : 4;
687   st->cmap=(st->cmap+5+(random()%10))%COLORMAPS;
688
689   for (x=0; x<LAYERS; x++) {
690     st->pattern[x]=random()%PATTERNS;
691     st->shape[x]=random()%SHAPES;
692     st->mix[x]=random()%20;
693     nstr = (st->lwid==2) ? 20+random()%12 :
694       (st->lwid==3) ? 16+random()%8 :
695       (st->lwid==4) ? 12+random()%6 :
696       (st->lwid==5) ? 10+random()%5 :
697       (st->lwid==6) ? 8+random()%4 :
698         5+random()%5;
699     st->csw[x] = _max(5,st->gridy/nstr);
700     st->wsx[x] = (st->wsx[x]+3+(random()%3))%STRETCHES;
701     st->wsy[x] = (st->wsy[x]+3+(random()%3))%STRETCHES;
702     st->sec[x] = random()%5;
703     if ((!st->dialog) && (st->sec[x]<2)) st->csw[x]/=2;
704     st->cs1[x] = (st->dialog) ? 1+random()%3 : 2+random()%5;
705     st->cs2[x] = (st->dialog) ? 1+random()%3 : 2+random()%5;
706     st->cs3[x] = (st->dialog) ? 1+random()%3 : 2+random()%5;
707     st->cs4[x] = (st->dialog) ? 1+random()%3 : 2+random()%5;
708     st->wave[x]=random()%WAVES;
709     st->wavel[x]=st->csw[x]*(2+random()%6);
710     st->waveh[x]=st->csw[x]*(1+random()%3);
711     st->rx1[x]=(st->gridx/10+random()%(st->gridx*8/10));
712     st->ry1[x]=(st->gridy/10+random()%(st->gridy*8/10));
713     st->rx2[x]=(st->gridx*2/10+random()%(st->gridx*6/10));
714     st->ry2[x]=(st->gridy*2/10+random()%(st->gridy*6/10));
715     st->rx3[x]=(st->gridx*3/10+random()%(st->gridx*4/10));
716     st->ry3[x]=(st->gridy*3/10+random()%(st->gridy*4/10));
717   }
718 }
719
720 static int
721 _shape(struct state *st, int x, int y, int rx, int ry, int n)
722 {
723   switch(st->shape[n]) {
724     case 0:                                        /* square/rectangle */
725     case 1:
726     case 2:
727       return(1+_max(abs(x-rx)*st->cs1[n]/st->cs2[n],abs(y-ry)*st->cs3[n]/st->cs4[n]));
728     case 3:                                                 /* diamond */
729     case 4:
730       return(1+(abs(x-rx)*st->cs1[n]/st->cs2[n]+abs(y-ry)*st->cs3[n]/st->cs4[n]));
731     case 5:                                            /* 8 point star */
732       return(1+_min(_max(abs(x-rx),abs(y-ry))*3/2,abs(x-rx)+abs(y-ry)));
733     case 6:                                             /* circle/oval */
734     case 7:
735     case 8:
736       return(1+_dist(st,x,rx,y,ry,st->cs1[n]));
737     case 9:                                       /* black hole circle */
738       return(1+(st->gridx*st->gridy/(1+(_dist(st,x,rx,y,ry,st->cs2[n])))));
739     case 10:                                                   /* sun */
740       return(1+_min(abs(x-rx)*st->gridx/(abs(y-ry)+1),abs(y-ry)*st->gridx/(abs(x-rx)+1)));
741     case 11:                             /* 2 circles+inverted circle */
742       return(1+(_dist(st,x,rx,y,ry,st->cs1[n])*_dist(st,x,(rx*3)%st->gridx,y,(ry*5)%st->gridy,st->cs1[n])/(1+_dist(st,x,(rx*4)%st->gridx,y,(ry*7)%st->gridy,st->cs1[n]))));
743     case 12:                                                   /* star */
744       return(1+(int)sqrt(abs((x-rx)*(y-ry))));
745     case 13:                                       /* centered ellipse */
746       return(1+_dist(st,x,rx,y,ry,0)+_dist(st,x,st->gridx-rx,y,st->gridy-ry,0));
747     default:                                               /* triangle */
748       return(1+_triangle(st,x,y,rx,ry,st->cs4[n]));
749   }
750   return(1+_triangle(st,x,y,rx,ry,st->cs4[n]));
751 }
752
753 static int
754 _pattern(struct state *st, int x, int y, int n)
755 {
756   int v=0, ox;
757   ox=x;
758   switch(st->wsx[n]) {
759     case 0:                                             /* slants */
760       x+=y/(1+st->cs4[n]);
761       break;
762     case 1:
763       x+=(st->gridy-y)/(1+st->cs4[n]);
764       break;
765     case 2:                                             /* curves */
766       x+=_wave(st,y,st->gridx/(1+st->cs1[n]),st->gridy,0);
767       break;
768     case 3:
769       x+=_wave(st,st->gridy-y,st->gridy/(1+st->cs1[n]),st->gridy,0);
770       break;
771     case 4:                                           /* U curves */
772       x+=_wave(st,y,st->cs1[n]*st->csw[n]/2,st->gridy*2/M_PI,0);
773       break;
774     case 5:
775       x-=_wave(st,y,st->cs1[n]*st->csw[n]/2,st->gridy*2/M_PI,0);
776       break;
777   }
778   switch(st->wsy[0]) {
779     case 0:                                          /* slants */
780       y+=ox/(1+st->cs1[n]);
781       break;
782     case 1:
783       y+=(st->gridx-ox)/(1+st->cs1[n]);
784       break;
785     case 2:                                           /* curves */
786       y+=_wave(st,ox,st->gridx/(1+st->cs1[n]),st->gridx,0);
787       break;
788     case 3:
789       y+=_wave(st,st->gridx-ox,st->gridx/(1+st->cs1[n]),st->gridx,0);
790       break;
791     case 4:                                         /* U curves */
792       y+=_wave(st,ox,st->cs1[n]*st->csw[n]/2,st->gridy*2/M_PI,0);
793       break;
794     case 5:
795       y-=_wave(st,ox,st->cs1[n]*st->csw[n]/2,st->gridy*2/M_PI,0);
796       break;
797   }
798   switch(st->pattern[n]) {
799     case 0:                                          /* horizontal stripes */
800       v=y;
801       break;
802     case 1:                                            /* vertical stripes */
803       v=x;
804       break;
805     case 2:                                            /* diagonal stripes */
806       v=(x+(y*st->cs1[n]/st->cs2[n]));
807       break;
808     case 3:                                    /* reverse diagonal stripes */
809       v=(x-(y*st->cs1[n]/st->cs2[n]));
810       break;
811     case 4:                                                /* checkerboard */
812       v=(y/st->csw[n]*3+x/st->csw[n])*st->csw[n];
813       break;
814     case 5:                                       /* diagonal checkerboard */
815       v=((x+y)/2/st->csw[n]+(x+st->gridy-y)/2/st->csw[n]*3)*st->csw[n];
816       break;
817     case 6:                                                     /* + cross */
818       v=st->gridx+(_min(abs(x-st->rx3[n]),abs(y-st->ry3[n]))*2);
819       break;
820     case 7:                                              /* double + cross */
821       v=_min(_min(abs(x-st->rx2[n]),abs(y-st->ry2[n])),_min(abs(x-st->rx1[n]),abs(y-st->ry1[n])))*2;
822       break;
823     case 8:                                                     /* X cross */
824       v=st->gridx+(_min(abs(x-st->rx3[n])*st->cs1[n]/st->cs2[n]+abs(y-st->ry2[n])*st->cs3[n]/st->cs4[n],abs(x-st->rx3[n])*st->cs1[n]/st->cs2[n]-abs(y-st->ry3[n])*st->cs3[n]/st->cs4[n])*2);
825       break;
826     case 9:                                              /* double X cross */
827       v=_min(_min(abs(x-st->rx2[n])+abs(y-st->ry2[n]),abs(x-st->rx2[n])-abs(y-st->ry2[n])),_min(abs(x-st->rx1[n])+abs(y-st->ry1[n]),abs(x-st->rx1[n])-abs(y-st->ry1[n])))*2;
828       break;
829     case 10:                                   /* horizontal stripes/waves */
830       v=st->gridy+(y+_wave(st,x,st->waveh[n],st->wavel[n],st->wave[n]));
831       break;
832     case 11:                                     /* vertical stripes/waves */
833       v=st->gridx+(x+_wave(st,y,st->waveh[n],st->wavel[n],st->wave[n]));
834       break;
835     case 12:                                     /* diagonal stripes/waves */
836       v=st->gridx+(x+(y*st->cs1[n]/st->cs2[n])+_wave(st,x,st->waveh[n],st->wavel[n],st->wave[n]));
837       break;
838     case 13:                                     /* diagonal stripes/waves */
839       v=st->gridx+(x-(y*st->cs1[n]/st->cs2[n])+_wave(st,y,st->waveh[n],st->wavel[n],st->wave[n]));
840       break;
841     case 14:                                    /* horizontal spikey waves */
842       v=y+(st->csw[n]*st->cs4[n]/st->cs3[n])+_wave(st,x+((y/st->cs3[n])*st->edir),st->csw[n]/2*st->cs1[n]/st->cs2[n],st->csw[n]/2*st->cs2[n]/st->cs1[n],st->wave[n]);
843       break;
844     case 15:                                      /* vertical spikey waves */
845       v=x+(st->csw[n]*st->cs1[n]/st->cs2[n])+_wave(st,y+((x/st->cs3[n])*st->edir),st->csw[n]/2*st->cs1[n]/st->cs2[n],st->csw[n]/2*st->cs3[n]/st->cs4[n],st->wave[n]);
846       break;
847     case 16:                                    /* big slanted hwaves */
848       v=st->gridy-y-(x*st->cs1[n]/st->cs3[n])+(st->csw[n]*st->cs1[n]*st->cs2[n]) +_wave(st,x,st->csw[n]/3*st->cs1[n]*st->cs2[n],st->csw[n]/3*st->cs3[n]*st->cs2[n],st->wave[n]);
849       break;
850     case 17:                                    /* big slanted vwaves */
851       v=x-(y*st->cs1[n]/st->cs3[n])+(st->csw[n]*st->cs1[n]*st->cs2[n]) +_wave(st,y, st->csw[n]/3*st->cs1[n]*st->cs2[n], st->csw[n]/3*st->cs3[n]*st->cs2[n], st->wave[n]);
852       break;
853     case 18:                                          /* double hwave */
854       v=y+(y+st->csw[n]*st->cs3[n])+_wave(st,x,st->csw[n]/3*st->cs3[n],st->csw[n]/3*st->cs2[n],st->wave[n])+_wave(st,x,st->csw[n]/3*st->cs4[n],st->csw[n]/3*st->cs1[n]*3/2,st->wave[n]);
855       break;
856     case 19:                                          /* double vwave */
857       v=x+(x+st->csw[n]*st->cs1[n])+_wave(st,y,st->csw[n]/3*st->cs1[n],st->csw[n]/3*st->cs3[n],st->wave[n])+_wave(st,y,st->csw[n]/3*st->cs2[n],st->csw[n]/3*st->cs4[n]*3/2,st->wave[n]);
858       break;
859     case 20:                                                  /* one shape */
860     case 21:
861     case 22:
862       v=_shape(st,x, y, st->rx3[n], st->ry3[n], n);
863       break;
864     case 23:                                                 /* two shapes */
865     case 24:
866     case 25:
867       v=_min(_shape(st,x, y, st->rx1[n], st->ry1[n], n),_shape(st,x, y, st->rx2[n], st->ry2[n], n));
868       break;
869     case 26:                                       /* two shapes opposites */
870     case 27:
871       v=_min(_shape(st,x, y, st->rx2[n], st->ry2[n], n),_shape(st,x, y, st->gridx-st->rx2[n], st->gridy-st->rx2[n], n));
872       break;
873     case 28:                                     /* two shape checkerboard */
874     case 29:
875       v=((_shape(st,x, y, st->rx1[n], st->ry1[n], n)/st->csw[n])+(_shape(st,x, y, st->rx2[n], st->ry2[n], n)/st->csw[n]))*st->csw[n];
876       break;
877     case 30:                                             /* two shape blob */
878     case 31:
879       v=(_shape(st,x, y, st->rx1[n], st->ry1[n], n)+_shape(st,x, y, st->rx2[n], st->ry2[n], n))/2;
880       break;
881     case 32:                                    /* inverted two shape blob */
882     case 33:
883       v=(_shape(st,x, y, st->rx1[n], st->ry1[n], n)+_shape(st,st->gridx-x, st->gridy-y, st->rx1[n], st->ry1[n], n))/2;
884       break;
885     case 34:                                               /* three shapes */
886     case 35:
887       v=_min(_shape(st,x, y, st->rx3[n], st->ry3[n], n),_min(_shape(st,x, y, st->rx1[n], st->ry1[n], n),_shape(st,x, y, st->rx2[n], st->ry2[n], n)));
888       break;
889     case 36:                                           /* three shape blob */
890     case 37:
891       v=(_shape(st,x, y, st->rx1[n], st->ry1[n], n)+_shape(st,x, y, st->rx2[n], st->ry2[n], n)+_shape(st,x, y, st->rx3[n], st->ry3[n], n))/3;
892       break;
893     case 38:                                                  /* 4 shapes */    
894       v=(_min(_shape(st,x, y, st->rx2[n], st->ry2[n], n),_shape(st,x, y, st->gridx-st->rx2[n], st->gridy-st->ry2[n], n)),_min(_shape(st,x, y, st->gridx-st->rx2[n], st->ry2[n], n),_shape(st,x, y, st->rx2[n], st->gridy-st->ry2[n], n)));
895       break;
896     case 39:                                            /* four rainbows */
897       v=(_min(_shape(st,x, y, st->gridx-st->rx2[n]/2, st->csw[n], n),_shape(st,x, y, st->csw[n], st->ry2[n]/2, n)),_min(_shape(st,x, y, st->rx2[n]/2, st->gridy-st->csw[n], n),_shape(st,x, y, st->gridx-st->csw[n], st->gridy-(st->ry2[n]/2), n)));
898       break;
899   }
900   /* stretch or contract stripe */
901   switch(st->sec[n]) {
902     case 0:
903       v=(int)sqrt((int)sqrt(abs(v)*st->gridx)*st->gridx);
904       break;
905     case 1:
906       v=((int)pow(v,2)/st->gridx);
907       break;
908   }
909   return (abs(v));
910 }
911
912 static int
913 _getcolor(struct state *st, int x, int y)
914 {
915   int n, cv[LAYERS];
916
917   for (n=0; n<st->layers; n++) {
918     cv[n]=_pattern(st,x,y,n);
919                   /* first wave/shape */
920     cv[0] = (!n) ? cv[0]/st->csw[0] :
921                     /* checkerboard+1 */
922       (st->mix[n]<5) ? (cv[0]*st->csw[0]+cv[n])/st->csw[n] :
923                /* checkerboard+ncol/2 */
924       (st->mix[n]<12) ? cv[0]+(cv[n]/st->csw[n]*st->ncolors/2) :
925                            /* add mix */
926       (st->mix[n]<16) ? cv[0]+(cv[n]/st->csw[n]) :
927                       /* subtract mix */
928       (st->mix[n]<18) ? cv[0]-(cv[n]/st->csw[n]) :
929                   /* r to l morph mix */
930       (st->mix[n]==18) ? ((cv[0]*x)+(cv[n]*(st->gridx-x)/st->csw[n]))/st->gridx :
931                   /* u to d morph mix */
932       ((cv[0]*y)+(cv[n]*(st->gridy-y)/st->csw[n]))/st->gridy;
933   }
934   return(cv[0]);
935 }
936
937 /* return value=line direction
938    st->olen=open space to edge or next blocking line
939    st->bln=blocking line number or -1 if edge blocks */
940 static int
941 _findopen(struct state *st, int x, int y, int z)
942 {
943   int dir, od[4], no=0;
944
945   if (((st->grid[z].hl) || (st->grid[z].hr)) &&
946     ((st->grid[z].vu) || (st->grid[z].vd)))
947     return(DIR_NONE);
948   if ((z>st->gridx) && (!st->grid[z].hl) && (!st->grid[z].hr) &&
949     (!st->grid[z-st->gridx].line)) {
950     od[no]=DIR_UP;
951     no++;
952   }
953   if ((z<st->gridn-st->gridx) && (!st->grid[z].hl) &&
954     (!st->grid[z].hr) && (!st->grid[z+st->gridx].line)) {
955     od[no]=DIR_DOWN;
956     no++;
957   }
958   if ((x) && (!st->grid[z].hl) && (!st->grid[z].hr) &&
959     (!st->grid[z-1].line)) {
960     od[no]=DIR_LEFT;
961     no++;
962   }
963   if (((z+1)%st->gridx) && (!st->grid[z].hl) && (!st->grid[z].hr) &&
964     (!st->grid[z+1].line)) {
965     od[no]=DIR_RIGHT;
966     no++;
967   }
968   if (!no)
969     return(DIR_NONE);
970   dir=od[random()%no];
971   st->olen=st->bln=0;
972   while ((st->olen<=st->maxlen) && (!st->bln)) {
973     st->olen++;
974     if (dir==DIR_UP)
975       st->bln = (y-st->olen<0) ? -1 :
976         st->grid[z-(st->olen*st->gridx)].line;
977     if (dir==DIR_DOWN)
978       st->bln = (y+st->olen>=st->gridy) ? -1 :
979         st->grid[z+(st->olen*st->gridx)].line;
980     if (dir==DIR_LEFT)
981       st->bln = (x-st->olen<0) ? -1 :
982         st->grid[z-st->olen].line;
983     if (dir==DIR_RIGHT)
984       st->bln = (x+st->olen>=st->gridx) ? -1 :
985         st->grid[z+st->olen].line;
986   }
987   st->olen--;
988   return(dir);
989 }
990
991 static void
992 _fillgrid(struct state *st)
993 {
994   unsigned int gridc, n, add;
995
996   gridc=st->gridx*st->dline[st->li].y+st->dline[st->li].x;
997   add = (st->dline[st->li].hv) ? 1 : st->gridx;
998   for (n=0;  n<=st->dline[st->li].len; n++) {
999     if (n)
1000       gridc+=add;
1001     if (!st->grid[gridc].line) {
1002       st->fi++;
1003       st->grid[gridc].line=st->li;
1004     }
1005     if (st->dline[st->li].hv) {
1006       if (n)
1007         st->grid[gridc].hr=st->li;
1008       if (n<st->dline[st->li].len)
1009         st->grid[gridc].hl=st->li;
1010     } else {
1011       if (n)
1012         st->grid[gridc].vd=st->li;
1013       if (n<st->dline[st->li].len)
1014         st->grid[gridc].vu=st->li;
1015     }
1016     if (st->fi>=st->gridn) {
1017       st->grid_full=True;
1018       return;
1019     }
1020   }
1021 }
1022
1023 static void
1024 _newline(struct state *st)
1025 {
1026   int bl, bz, dir, lt, x, y, z;
1027
1028   bl=0;
1029   z=st->zlist[st->zi];
1030   x=z%st->gridx;
1031   y=z/st->gridx;
1032   st->zi++;
1033   dir=_findopen(st,x,y,z);
1034
1035   if (!st->grid[z].line) {
1036   /* this is an empty space, make a new line unless nothing is open around it */
1037     if (dir==DIR_NONE) {
1038       /* nothing is open, force a len 1 branch in any direction */
1039       lt=LINE_FORCE;
1040       while ((dir==DIR_NONE) ||
1041         ((dir==DIR_UP) && (!y)) ||
1042         ((dir==DIR_DOWN) && (y+1==st->gridy)) ||
1043         ((dir==DIR_LEFT) && (!x)) ||
1044         ((dir==DIR_RIGHT) && (x+1==st->gridx))) {
1045           dir=random()%4;
1046       }
1047       bz = (dir==DIR_UP) ? z-st->gridx : (dir==DIR_DOWN) ? z+st->gridx : (dir==DIR_LEFT) ? z-1 : z+1;
1048       bl = st->grid[bz].line;
1049     } else if ((st->bnratio>1) && (st->bln>0) &&
1050       (st->olen<st->maxlen) && (random()%st->bnratio)) {
1051       /* branch into blocking line */
1052       lt=LINE_BRIN;
1053       bl = st->bln;
1054     } else {
1055       /* make a new line and new object */
1056       lt=LINE_NEW;
1057       st->oi++;
1058     }
1059   } else {
1060     /* this is a filled space, make a branch unless nothing is open around it */
1061     if (dir==DIR_NONE)
1062       return;
1063     /* make a branch out of this line */
1064     lt=LINE_BROUT;
1065     bl=st->grid[z].line;
1066   }
1067   st->li++;
1068   st->dline[st->li].len = (lt==LINE_FORCE) ? 1 :  (lt==LINE_BRIN) ?
1069     st->olen+1 : (!st->forcemax) ? st->olen : 1+random()%st->olen;
1070   st->dline[st->li].x=x;
1071   if (dir==DIR_LEFT)
1072     st->dline[st->li].x-=st->dline[st->li].len;
1073   st->dline[st->li].y=y;
1074   if (dir==DIR_UP)
1075     st->dline[st->li].y-=st->dline[st->li].len;
1076   st->dline[st->li].hv = ((dir==DIR_LEFT) || (dir==DIR_RIGHT)) ?
1077     True : False;
1078   st->dline[st->li].obj = (lt==LINE_NEW) ? st->oi :
1079     st->dline[bl].obj;
1080   if (lt==LINE_NEW) {
1081     int color =  (_getcolor(st,x,y))%st->ncolors;
1082     if (color < 0) color += st->ncolors;
1083     st->dline[st->li].color = color;
1084   } else {
1085     st->dline[st->li].color = st->dline[bl].color;
1086   }
1087   st->dline[st->li].deo=(_getdeo(st,x,y,st->dmap,1) +
1088     (random()%st->dvar) + (random()%st->dvar))*st->ddir;
1089   st->dline[st->li].ndol=0;
1090   _fillgrid(st);
1091 }
1092
1093 static void
1094 _create_screen(struct state *st)
1095 {
1096   while(!st->grid_full)
1097     _newline(st);
1098   qsort(st->dline, st->li+1, sizeof(struct lineStruct), _comparedeo);
1099 /*st->lpu=st->li/20/((6-st->speed)*3);
1100   Used to use a computed lpu, lines per update to control draw speed
1101   draw 1/lpu of the lines before each XSync which takes a split second
1102   the higher the lpu, the quicker the screen draws.  This worked somewhat
1103   after the 4->5 update, however with the Mac updating so much more slowly,
1104   values tuned for it draw the screen in a blink on Linux.  Therefore we
1105   draw 1/200th of the screen with each update and sleep, if necessary */
1106   st->lpu = (st->dialog) ? st->li/50 : st->li/200;
1107   if (!st->lpu) st->lpu = 1;
1108   st->bi=1;
1109   st->mode=MODE_ERASE;
1110 }
1111
1112 static void
1113 _fill_outline(struct state *st, int di)
1114 {
1115   int x, y, h, w;
1116
1117   if (!di)
1118     return;
1119   x=st->dline[di].x*st->lwid+1;
1120   y=st->dline[di].y*st->lwid+1;
1121   if (st->dline[di].hv) {
1122     w=(st->dline[di].len+1)*st->lwid-3;
1123     h=st->lwid-3;
1124   } else {
1125     w=st->lwid-3;
1126     h=(st->dline[di].len+1)*st->lwid-3;
1127   }
1128   XFillRectangle (st->display, st->window, st->bgc, x, y, w, h);
1129 }
1130
1131 static void
1132 _XFillRectangle(struct state *st, int di, int adj)
1133 {
1134   int a, b, x, y, w, h;
1135
1136   x=st->dline[di].x*st->lwid;
1137   y=st->dline[di].y*st->lwid;
1138   if (st->dline[di].hv) {
1139     w=(st->dline[di].len+1)*st->lwid-1;
1140     h=st->lwid-1;
1141   } else {
1142     w=st->lwid-1;
1143     h=(st->dline[di].len+1)*st->lwid-1;
1144   }
1145   switch (st->d3d) {
1146     case D3D_NEON:
1147       x+=adj;
1148       y+=adj;
1149       w-=adj*2;
1150       h-=adj*2;
1151     break;
1152     case D3D_BLOCK:
1153       x+=adj;
1154       y+=adj;
1155       w-=st->lwid/2-1;
1156       h-=st->lwid/2-1;
1157     break;
1158   }
1159   if (!st->round) {
1160     XFillRectangle(st->display, st->window, st->fgc, x, y, w, h);
1161   } else {
1162     if (h<st->lwid) {                                   /* horizontal */
1163       a=(h-1)/2;
1164       for (b=0; b<=a; b++)
1165         XFillRectangle(st->display, st->window, st->fgc,
1166           x+b, y+a-b, w-b*2, h-((a-b)*2));
1167     } else {                                               /* vertical */
1168       a=(w-1)/2;
1169       for (b=0; b<=a; b++)
1170         XFillRectangle(st->display, st->window, st->fgc,
1171           x+a-b, y+b, w-((a-b)*2), h-b*2);
1172     }
1173   }
1174 }
1175
1176 static void
1177 _XFillTriangle(struct state *st, int color, int x1, int y1, int x2, int y2,
1178   int x3, int y3)
1179 {
1180   XPoint points[3];
1181
1182   points[0].x=x1;
1183   points[0].y=y1;
1184   points[1].x=x2;
1185   points[1].y=y2;
1186   points[2].x=x3;
1187   points[2].y=y3;
1188   XSetForeground(st->display, st->fgc, st->colors[color].pixel);
1189   XFillPolygon (st->display, st->window, st->fgc, points, 3, Convex,
1190       CoordModeOrigin);
1191 }
1192
1193 static void
1194 _XFillPolygon4(struct state *st, int color, int x1, int y1, int x2, int y2,
1195   int x3, int y3, int x4, int y4)
1196 {
1197   XPoint points[4];
1198
1199   points[0].x=x1;
1200   points[0].y=y1;
1201   points[1].x=x2;
1202   points[1].y=y2;
1203   points[2].x=x3;
1204   points[2].y=y3;
1205   points[3].x=x4;
1206   points[3].y=y4;
1207   XSetForeground(st->display, st->fgc, st->colors[color].pixel);
1208   XFillPolygon (st->display, st->window, st->fgc, points, 4, Convex,
1209       CoordModeOrigin);
1210 }
1211
1212 static void
1213 _draw_tiled(struct state *st, int color)
1214 {
1215   int a, c, d, x, y, z, m1, m2, lr, nl, w, h;
1216   a = (st->dline[st->di].hv) ? 1 : st->gridx;
1217   z = st->dline[st->di].y*st->gridx+st->dline[st->di].x;
1218   m1 = (st->lwid-1)/2;
1219   m2 = st->lwid/2;
1220   lr = st->lwid-1;
1221   nl = st->lwid;
1222
1223   /* draw tiles one grid cell at a time */
1224   for (c=0; c<=st->dline[st->di].len; c++) {
1225     if (st->dline[st->di].hv) {
1226       x = (st->dline[st->di].x+c)*st->lwid;
1227       y = st->dline[st->di].y*st->lwid;
1228       if (c)
1229         st->grid[z].dhr=st->di;
1230       if (c<st->dline[st->di].len)
1231         st->grid[z].dhl=st->di;
1232     } else {
1233       x = st->dline[st->di].x*st->lwid;
1234       y = (st->dline[st->di].y+c)*st->lwid;
1235       if (c)
1236         st->grid[z].dvd=st->di;
1237       if (c<st->dline[st->di].len)
1238         st->grid[z].dvu=st->di;
1239     }
1240     d=0;
1241     if (st->grid[z].dhl)
1242       d+=8;
1243     if (st->grid[z].dhr)
1244       d+=4;
1245     if (st->grid[z].dvu)
1246       d+=2;
1247     if (st->grid[z].dvd)
1248       d++;
1249     /* draw line base */
1250     switch (d) {
1251       case 1:
1252       case 2:                                    /* vertical */
1253       case 3:
1254       case 5:
1255       case 6:
1256       case 7:
1257       case 11:
1258       case 15:
1259         h = ((d==1) || (d==5)) ? lr : nl;
1260         XSetForeground(st->display, st->fgc,
1261           st->colors[color].pixel);
1262         XFillRectangle (st->display, st->window, st->fgc,
1263           x, y, m2, h);
1264         XSetForeground(st->display, st->fgc,
1265            st->colors[color+3].pixel);
1266         XFillRectangle (st->display, st->window, st->fgc,
1267           x+m2, y, m1, h);
1268         break;
1269       case 4:
1270       case 8:                                     /* horizontal */
1271       case 9:
1272       case 10:
1273       case 12:
1274       case 13:
1275       case 14:
1276         w = (d==4) ? lr : nl;
1277         XSetForeground(st->display, st->fgc,
1278           st->colors[color+1].pixel);
1279         XFillRectangle (st->display, st->window, st->fgc,
1280           x, y, w, m2);
1281         XSetForeground(st->display, st->fgc,
1282            st->colors[color+2].pixel);
1283         XFillRectangle (st->display, st->window, st->fgc,
1284           x, y+m2, w, m1);
1285         break;
1286     }
1287     /* draw angles */
1288     switch(d) {
1289       case 1:                                      /* bottom end ^ */
1290         _XFillTriangle(st,color+2, x, y+lr, x+lr, y+lr, x+m2, y+m2);
1291         break;
1292       case 2:                                       /* top end \/ */
1293         _XFillTriangle(st,color+1, x, y, x+lr, y, x+m2, y+m2);
1294         break;
1295       case 4:                                       /* right end < */
1296         _XFillTriangle(st,color+3, x+lr, y, x+lr, y+lr, x+m2, y+m2);
1297         break;
1298       case 5:                                        /* LR corner */
1299         _XFillTriangle(st,color+1, x, y+m2, x+m2, y+m2, x, y);
1300         _XFillPolygon4(st,color+2, x, y+m2, x+m2, y+m2, x+lr, y+lr, x, y+lr);
1301         break;
1302       case 6:                                        /* UR corner */
1303         _XFillPolygon4(st,color+1, x, y+m2, x+m2, y+m2, x+lr, y, x, y);
1304         _XFillTriangle(st,color+2, x, y+m2, x+m2, y+m2, x, y+lr);
1305         break;
1306       case 7:                                        /* T > into line */
1307         _XFillTriangle(st,color+1, x, y+m2, x+m2, y+m2, x, y);
1308         _XFillTriangle(st,color+2, x, y+m2, x+m2, y+m2, x, y+lr);
1309         break;
1310       case 8:                                       /* left end > */
1311         _XFillTriangle(st,color, x, y, x, y+lr, x+m2, y+m2);
1312         break;
1313       case 9:                                       /* LL corner */
1314         _XFillPolygon4(st,color, x+m2, y, x+m2, y+m2, x, y+lr, x, y);
1315         _XFillTriangle(st,color+3, x+m2, y, x+m2, y+m2, x+lr, y);
1316         break;
1317       case 10:                                      /* UL corner */
1318         _XFillPolygon4(st,color, x+m2, y+nl, x+m2, y+m2, x, y, x, y+nl);
1319         _XFillPolygon4(st,color+3, x+m2, y+nl, x+m2, y+m2, x+lr, y+lr, x+lr, y+nl);
1320         break;
1321       case 11:                                       /* T < into line */
1322         _XFillPolygon4(st,color+1, x+nl, y+m2, x+m2, y+m2, x+lr, y, x+nl, y);
1323         _XFillPolygon4(st,color+2, x+nl, y+m2, x+m2, y+m2, x+lr, y+lr, x+nl, y+lr);
1324         break;
1325       case 13:                                        /* T \/ into line */
1326         _XFillTriangle(st,color, x+m2, y, x+m2, y+m2, x, y);
1327         _XFillTriangle(st,color+3, x+m2, y, x+m2, y+m2, x+lr, y);
1328         break;
1329       case 14:                                        /* T ^ into line */
1330         _XFillPolygon4(st,color, x+m2, y+nl, x+m2, y+m2, x, y+lr, x, y+nl);
1331         _XFillPolygon4(st,color+3, x+m2, y+nl, x+m2, y+m2, x+lr, y+lr, x+lr, y+nl);
1332         break;
1333       case 15:                                        /* X intersection */
1334         _XFillTriangle(st,color+1, x, y+m2, x+m2, y+m2, x, y);
1335         _XFillTriangle(st,color+2, x, y+m2, x+m2, y+m2, x, y+lr);
1336         _XFillPolygon4(st,color+1, x+nl, y+m2, x+m2, y+m2, x+lr, y, x+nl, y);
1337         _XFillPolygon4(st,color+2, x+nl, y+m2, x+m2, y+m2, x+lr, y+lr, x+nl, y+lr);
1338         break;
1339     }
1340     z+=a;
1341   }
1342 }
1343
1344 static long
1345 _mselapsed(struct state *st)
1346 {
1347   struct timeval t;
1348   gettimeofday(&t, NULL);
1349   t.tv_sec -= st->time.tv_sec;
1350   t.tv_usec -= st->time.tv_usec;
1351   return ((long)t.tv_sec*1000000+t.tv_usec);
1352 }
1353
1354 static void
1355 _draw_lines(struct state *st)
1356 {
1357   int n, z, a, color, sh, di;
1358   if (st->bi==1)
1359     for (a=0; a<=st->oi; a++)
1360       st->fdol[a]=0;
1361
1362   for (st->di=st->bi; st->di<_min(st->li+1,st->bi+st->lpu); st->di++) {
1363     color=(st->dline[st->di].color%st->ncolors)*st->shades;
1364     XSetForeground(st->display, st->fgc, st->colors[color].pixel);
1365
1366     switch (st->d3d) {
1367       case D3D_NEON:
1368         st->dline[st->di].ndol=st->fdol[st->dline[st->di].obj];
1369         st->fdol[st->dline[st->di].obj]=st->di;
1370         for (sh=0; sh<st->lwid/2; sh++) {
1371           XSetForeground(st->display, st->fgc,
1372             st->colors[color+sh].pixel);
1373           di=st->di;
1374           while(di>0) {
1375             _XFillRectangle(st,di,sh);
1376             di=st->dline[di].ndol;
1377           }
1378         }
1379         break;
1380       case D3D_BLOCK:
1381         st->dline[st->di].ndol=st->fdol[st->dline[st->di].obj];
1382         st->fdol[st->dline[st->di].obj]=st->di;
1383         for (sh=0; sh<st->lwid/2; sh++) {
1384           XSetForeground(st->display, st->fgc,
1385             st->colors[color+(st->lwid/2)-sh-1].pixel);
1386           di=st->di;
1387           while(di>0) {
1388             _XFillRectangle(st,di,sh);
1389             di=st->dline[di].ndol;
1390           }
1391         }
1392         break;
1393       case D3D_TILED:
1394         _draw_tiled(st,color);
1395         break;
1396       default:               /* D3D_NONE */
1397         _XFillRectangle(st,st->di,0);
1398         if (st->outline) {
1399           _fill_outline(st, st->di);
1400           z=st->dline[st->di].y*st->gridx+st->dline[st->di].x;
1401           a = (st->dline[st->di].hv) ? 1 : st->gridx;
1402           for (n=0; n<=st->dline[st->di].len; n++) {
1403             _fill_outline(st, st->grid[z].dhl);
1404             _fill_outline(st, st->grid[z].dhr);
1405             _fill_outline(st, st->grid[z].dvu);
1406             _fill_outline(st, st->grid[z].dvd);
1407             if (st->dline[st->di].hv) {
1408               if (n)
1409                 st->grid[z].dhr=st->di;
1410               if (n<st->dline[st->di].len)
1411                 st->grid[z].dhl=st->di;
1412             } else {
1413               if (n)
1414                 st->grid[z].dvd=st->di;
1415               if (n<st->dline[st->di].len)
1416                 st->grid[z].dvu=st->di;
1417             }
1418             z+=a;
1419           }
1420         }
1421         break;
1422     }
1423   }
1424   if (st->di>st->li) {
1425       st->bi=1;
1426       st->mode=MODE_CREATE;
1427   } else {
1428       st->bi+=st->lpu;
1429   }
1430 }
1431
1432 static void
1433 _erase_lines(struct state *st)
1434 {
1435   if (!st->ii)
1436     return;
1437   for (st->di=st->bi; st->di<_min(st->eli+1,st->bi+st->elpu); st->di++) {
1438     if (st->eline[st->di].hv) {
1439       XFillRectangle (st->display, st->window, st->bgc,
1440       st->eline[st->di].x*st->elwid,
1441       st->eline[st->di].y*st->elwid,
1442       (st->eline[st->di].len+1)*st->elwid, st->elwid);
1443     } else {
1444       XFillRectangle (st->display, st->window, st->bgc,
1445       st->eline[st->di].x*st->elwid,
1446       st->eline[st->di].y*st->elwid,
1447       st->elwid, (st->eline[st->di].len+1)*st->elwid);
1448     }
1449     if (st->di==st->eli) /* clear just in case */
1450       XFillRectangle(st->display, st->window, st->bgc, 0, 0,
1451         st->xgwa.width, st->xgwa.height);
1452   }
1453   if (st->di>st->eli) {
1454       st->bi=1;
1455       if (st->resized) {
1456          st->mode=MODE_CREATE;
1457       } else {
1458          st->mode=MODE_DRAW;
1459       }
1460   } else {
1461       st->bi+=st->elpu;
1462   }
1463 }
1464
1465 static void *
1466 abstractile_init(Display *display, Window window)
1467 {
1468   struct state *st = (struct state *) calloc (1, sizeof(*st));
1469   XGCValues gcv;
1470 /*  struct utsname os;*/
1471
1472   char *tile = get_string_resource(display, "tile", "Tile");
1473   if      (tile && !strcmp(tile, "random")) st->tile = TILE_RANDOM;
1474   else if (tile && !strcmp(tile, "flat")) st->tile = TILE_FLAT;
1475   else if (tile && !strcmp(tile, "thin")) st->tile = TILE_THIN;
1476   else if (tile && !strcmp(tile, "outline")) st->tile = TILE_OUTLINE;
1477   else if (tile && !strcmp(tile, "block")) st->tile = TILE_BLOCK;
1478   else if (tile && !strcmp(tile, "neon")) st->tile = TILE_NEON;
1479   else if (tile && !strcmp(tile, "tiled")) st->tile = TILE_TILED;
1480   else {
1481     if (tile && *tile && !!strcmp(tile, "random"))
1482       fprintf(stderr, "%s: unknown tile option %s\n", progname, tile);
1483     st->tile = TILE_RANDOM;
1484   }
1485
1486   st->speed = get_integer_resource(display, "speed", "Integer");
1487   if (st->speed < 0) st->speed = 0;
1488   if (st->speed > 5) st->speed = 5;
1489   st->sleep = get_integer_resource(display, "sleep", "Integer");
1490   if (st->sleep < 0) st->sleep = 0;
1491   if (st->sleep > 60) st->sleep = 60;
1492
1493   st->display=display;
1494   st->window=window;
1495
1496   /* get screen size and create Graphics Contexts */
1497   XGetWindowAttributes (display, window, &st->xgwa);
1498   gcv.foreground = get_pixel_resource(display, st->xgwa.colormap,
1499       "foreground", "Foreground");
1500   st->fgc = XCreateGC (display, window, GCForeground, &gcv);
1501   gcv.foreground = get_pixel_resource(display, st->xgwa.colormap,
1502       "background", "Background");
1503   st->bgc = XCreateGC (display, window, GCForeground, &gcv);
1504
1505 /* Um, no. This is obscene. -jwz.
1506   uname(&os);
1507   st->newcols=((!strcmp(os.sysname,"Linux")) || (!strcmp(os.sysname,"Darwin")))
1508       ? True : False;
1509 */
1510   st->newcols=False;
1511
1512   st->mode=MODE_CREATE;
1513   st->ii=0;
1514   st->resized=True;
1515   return st;
1516 }
1517
1518 static unsigned long
1519 abstractile_draw (Display *dpy, Window window, void *closure)
1520 {
1521   struct state *st = (struct state *) closure;
1522   int mse, usleep;
1523
1524   gettimeofday(&st->time, NULL);
1525
1526   /* If the window is too small, do nothing, sorry! */
1527   if (st->xgwa.width > 20 && st->xgwa.height > 20) {
1528     switch (st->mode) {
1529     case MODE_CREATE:
1530       _init_screen(st);
1531       _create_screen(st);
1532       break;
1533     case MODE_ERASE:
1534       _erase_lines(st);
1535       break;
1536     case MODE_DRAW:
1537       _draw_lines(st);
1538       break;
1539     }
1540   }
1541   mse=_mselapsed(st);
1542   usleep = ((!st->ii) && (st->mode==MODE_CREATE)) ?  0 :
1543       (st->mode==MODE_CREATE) ?  st->sleep*1000000-mse :
1544       /* speed=0-5, goal is 10,8,6,4,2,0 sec normal and 5,4,3,2,1,0 dialog */
1545       (5-st->speed)*(2-st->dialog)*100000/st->lpu-mse;
1546   if (usleep>=0)
1547       return usleep;
1548   return 0;
1549 }
1550
1551 static void
1552 abstractile_reshape (Display *dpy, Window window, void *closure,
1553                  unsigned int w, unsigned int h)
1554 {
1555   struct state *st = (struct state *) closure;
1556   st->xgwa.width = w;
1557   st->xgwa.height = h;
1558   if (w*h>st->max_wxh)
1559     st->resized=True;
1560 }
1561
1562 static Bool
1563 abstractile_event (Display *dpy, Window window, void *closure, XEvent *event)
1564 {
1565   struct state *st = (struct state *) closure;
1566   if (screenhack_event_helper (dpy, window, event))
1567     {
1568       st->mode=MODE_CREATE;
1569       return True;
1570     }
1571
1572   return False;
1573 }
1574
1575 static void
1576 abstractile_free (Display *dpy, Window window, void *closure)
1577 {
1578   struct state *st = (struct state *) closure;
1579   free (st);
1580 }
1581
1582 static const char *abstractile_defaults [] = {
1583   ".background:    black",
1584   ".foreground:    white",
1585   "*fpsSolid:      true",
1586   "*sleep:             3",
1587   "*speed:             3",
1588   "*tile:         random",
1589 #ifdef USE_IPHONE
1590   "*ignoreRotation: True",
1591 #endif
1592   0
1593 };
1594
1595 static XrmOptionDescRec abstractile_options [] = {
1596   { "-sleep",  ".sleep",  XrmoptionSepArg, 0 },
1597   { "-speed",  ".speed",  XrmoptionSepArg, 0 },
1598   { "-tile",   ".tile",   XrmoptionSepArg, 0 },
1599   { 0, 0, 0, 0 }
1600 };
1601
1602 XSCREENSAVER_MODULE ("Abstractile", abstractile)