1 /* -*- Mode: C; tab-width: 4 -*-
2 * swirl --- swirly color-cycling patterns.
4 #if !defined( lint ) && !defined( SABER )
5 static const char sccsid[] = "@(#)swirl.c 4.00 97/01/01 xlockmore";
8 /* Copyright (c) 1994 M.Dobie <mrd@ecs.soton.ac.uk>
10 * Permission to use, copy, modify, and distribute this software and its
11 * documentation for any purpose and without fee is hereby granted,
12 * provided that the above copyright notice appear in all copies and that
13 * both that copyright notice and this permission notice appear in
14 * supporting documentation.
16 * This file is provided AS IS with no warranties of any kind. The author
17 * shall have no liability with respect to the infringement of copyrights,
18 * trade secrets or any patents by this file or any part thereof. In no
19 * event will the author be liable for any lost revenue or profits or
20 * other special, indirect and consequential damages.
22 * 13-May-97: jwz@jwz.org: turned into a standalone program.
23 * 21-Apr-95: improved startup time for TrueColour displays
24 * (limited to 16bpp to save memory) S.Early <sde1000@cam.ac.uk>
25 * 09-Jan-95: fixed colour maps (more colourful) and the image now spirals
26 * outwards from the centre with a fixed number of points drawn
27 * every iteration. Thanks to M.Dobie <mrd@ecs.soton.ac.uk>.
28 * 1994: written. Copyright (c) 1994 M.Dobie <mrd@ecs.soton.ac.uk>
29 * based on original code by R.Taylor
33 # define PROGCLASS "Swirl"
34 # define HACK_INIT init_swirl
35 # define HACK_DRAW draw_swirl
36 # define swirl_opts xlockmore_opts
37 # define DEFAULTS "*count: 5 \n" \
41 # define SMOOTH_COLORS
42 # define WRITABLE_COLORS
43 # include "xlockmore.h" /* from the xscreensaver distribution */
44 # include <X11/Xutil.h>
45 # ifdef HAVE_XSHM_EXTENSION
47 # endif /* HAVE_XSHM_EXTENSION */
48 #else /* !STANDALONE */
49 # include "xlock.h" /* from the xlockmore distribution */
50 # undef HAVE_XSHM_EXTENSION
51 #endif /* !STANDALONE */
53 ModeSpecOpt swirl_opts = {
54 0, NULL, 0, NULL, NULL };
58 /****************************************************************/
60 #define MASS 4 /* maximum mass of a knot */
61 #define MIN_RES 5 /* minimim resolution (>= MIN_RES) */
62 #define MAX_RES 1 /* maximum resolution (>0) */
63 #define TWO_PLANE_PCNT 30 /* probability for two plane mode (0-100) */
64 #define RESTART 2500 /* number of cycles before restart */
65 #define BATCH_DRAW 100 /* points to draw per iteration */
80 int x, y; /* position */
82 KNOT_T t; /* type in the first (or only) plane */
83 KNOT_T T; /* type in second plane if there is one */
84 int M; /* mass in second plane if there is one */
87 /* a colour specification */
88 typedef struct Colour {
89 unsigned short r, g, b;
92 /* drawing direction */
94 DRAW_RIGHT, DRAW_DOWN, DRAW_LEFT, DRAW_UP
97 /****************************************************************/
99 /* data associated with a swirl window */
100 typedef struct swirl_data {
101 /* window paramaters */
102 Window win; /* the window */
103 int width, height; /* window size */
104 int depth; /* depth */
105 int rdepth; /* real depth (for XImage) */
106 Visual *visual; /* visual */
108 /* swirl drawing parameters */
109 int n_knots; /* number of knots */
110 KNOT_P knots; /* knot details */
111 KNOT_T knot_type; /* general type of knots */
112 int resolution; /* drawing resolution, 1..5 */
113 int max_resolution; /* maximum resolution, MAX_RES */
114 int r; /* pixel step */
115 Bool two_plane; /* two plane mode? */
116 Bool first_plane; /* doing first plane? */
117 int start_again; /* when to restart */
119 /* spiral drawing parameters */
120 int x, y; /* current point */
121 DIR_T direction; /* current direction */
122 int dir_todo, dir_done; /* how many points in current direction? */
123 int batch_todo, batch_done; /* how many points in this batch */
124 Bool started, drawing; /* are we drawing? */
127 unsigned char *image; /* image data */
131 int colours; /* how many colours possible */
132 int dcolours; /* how many colours for shading */
134 Bool fixed_colourmap; /* fixed colourmap? */
135 #endif /* !STANDALONE */
136 Bool monochrome; /* monochrome? */
137 Colormap cmap; /* colour map for the window */
138 XColor *rgb_values; /* colour definitions array */
140 int current_map; /* current colour map, 0..dcolours-1 */
141 unsigned long fg, bg, white, black; /* black and white pixel values */
142 int shift; /* colourmap shift */
143 int dshift; /* colourmap shift while drawing */
144 XColor fgcol, bgcol; /* foreground and background colour specs */
145 #endif /* !STANDALONE */
149 #define SWIRLCOLOURS 13
153 static COLOUR basic_colours[SWIRLCOLOURS];
154 #endif /* !STANDALONE */
156 /* an array of swirls for each screen */
157 static SWIRL_P swirls = NULL;
162 Return a random integer between 0 and n inclusive
164 - n is the maximum number
166 Returns a random integer */
169 random_no(unsigned int n)
171 return ((int) ((n + 1) * (double) LRAND() / MAXRAND));
174 /****************************************************************/
179 Initialise all the swirl data
181 - swirl is the swirl data */
184 initialise_swirl(ModeInfo * mi, SWIRL_P swirl)
187 Display *display = MI_DISPLAY(mi);
188 #endif /* !STANDALONE */
190 swirl->width = 0; /* width and height of window */
194 swirl->visual = NULL;
195 swirl->resolution = MIN_RES + 1; /* current resolution */
196 swirl->max_resolution = MAX_RES; /* maximum resolution */
197 swirl->n_knots = 0; /* number of knots */
198 swirl->knot_type = ALL; /* general type of knots */
199 swirl->two_plane = False; /* two plane mode? */
200 swirl->first_plane = False; /* doing first plane? */
201 swirl->start_again = -1; /* restart counter */
203 /* drawing parameters */
206 swirl->started = False;
207 swirl->drawing = False;
210 swirl->image = NULL; /* image data */
211 swirl->ximage = NULL;
214 swirl->colours = 0; /* how many colours possible */
215 swirl->dcolours = 0; /* how many colours for shading */
216 swirl->cmap = (Colormap) NULL;
217 swirl->rgb_values = NULL; /* colour definitions array */
219 swirl->current_map = 0; /* current colour map, 0..dcolours-1 */
221 /* set up fg fb colour specs */
222 swirl->white = MI_WIN_WHITE_PIXEL(mi);
223 swirl->black = MI_WIN_BLACK_PIXEL(mi);
224 #endif /* !STANDALONE */
228 # define MI_COLORMAP MI_WIN_COLORMAP
229 #else /* !STANDALONE */
230 swirl->fg = MI_FG_COLOR(mi);
231 swirl->bg = MI_BG_COLOR(mi);
232 swirl->fgcol.pixel = swirl->fg;
233 swirl->bgcol.pixel = swirl->bg;
234 XQueryColor(display, MI_COLORMAP(mi), &(swirl->fgcol));
235 XQueryColor(display, MI_COLORMAP(mi), &(swirl->bgcol));
236 #endif /* !STANDALONE */
239 /****************************************************************/
244 * Initialise the image for drawing to
246 * - swirl is the swirl data
249 initialise_image(ModeInfo * mi, SWIRL_P swirl)
251 Display *dpy = MI_DISPLAY(mi);
253 if (swirl->ximage != NULL)
254 XDestroyImage(swirl->ximage);
257 #ifdef HAVE_XSHM_EXTENSION
260 swirl->ximage = create_xshm_image(dpy, swirl->visual, swirl->rdepth,
261 ZPixmap, 0, &mi->shm_info,
262 swirl->width, swirl->height);
266 #endif /* HAVE_XSHM_EXTENSION */
270 swirl->ximage = XCreateImage(dpy, swirl->visual, swirl->rdepth, ZPixmap,
271 0, 0, swirl->width, swirl->height,
273 swirl->image = (unsigned char *)
274 calloc(swirl->height, swirl->ximage->bytes_per_line);
275 swirl->ximage->data = (char *) swirl->image;
279 /****************************************************************/
285 * Initialise the list of colours from which the colourmaps are derived
287 * - colours is the array to initialise
288 * - saturation is the saturation value to use 0->grey,
289 * 1.0->full saturation
292 initialise_colours(COLOUR * colours, float saturate)
296 /* start off fully saturated, medium and bright colours */
297 colours[0].r = 0xA000;
298 colours[0].g = 0x0000;
299 colours[0].b = 0x0000;
300 colours[1].r = 0xD000;
301 colours[1].g = 0x0000;
302 colours[1].b = 0x0000;
303 colours[2].r = 0x0000;
304 colours[2].g = 0x6000;
305 colours[2].b = 0x0000;
306 colours[3].r = 0x0000;
307 colours[3].g = 0x9000;
308 colours[3].b = 0x0000;
309 colours[4].r = 0x0000;
310 colours[4].g = 0x0000;
311 colours[4].b = 0xC000;
312 colours[5].r = 0x0000;
313 colours[5].g = 0x0000;
314 colours[5].b = 0xF000;
315 colours[6].r = 0xA000;
316 colours[6].g = 0x6000;
317 colours[6].b = 0x0000;
318 colours[7].r = 0xD000;
319 colours[7].g = 0x9000;
320 colours[7].b = 0x0000;
321 colours[8].r = 0xA000;
322 colours[8].g = 0x0000;
323 colours[8].b = 0xC000;
324 colours[9].r = 0xD000;
325 colours[9].g = 0x0000;
326 colours[9].b = 0xF000;
327 colours[10].r = 0x0000;
328 colours[10].g = 0x6000;
329 colours[10].b = 0xC000;
330 colours[11].r = 0x0000;
331 colours[11].g = 0x9000;
332 colours[11].b = 0xF000;
333 colours[12].r = 0xA000;
334 colours[12].g = 0xA000;
335 colours[12].b = 0xA000;
337 /* add white for low saturation */
338 for (i = 0; i < SWIRLCOLOURS - 1; i++) {
339 unsigned short max_rg, max;
341 /* what is the max intensity for this colour? */
342 max_rg = (colours[i].r > colours[i].g) ? colours[i].r : colours[i].g;
343 max = (max_rg > colours[i].b) ? max_rg : colours[i].b;
345 /* bring elements up to max as saturation approaches 0.0 */
346 colours[i].r += (unsigned short) ((float) (1.0 - saturate) *
347 ((float) max - colours[i].r));
348 colours[i].g += (unsigned short) ((float) (1.0 - saturate) *
349 ((float) max - colours[i].g));
350 colours[i].b += (unsigned short) ((float) (1.0 - saturate) *
351 ((float) max - colours[i].b));
354 #endif /* !STANDALONE */
356 /****************************************************************/
360 * set_black_and_white
362 * Set the entries for foreground & background pixels and
363 * WhitePixel & BlackPixel in an array of colour specifications.
365 * - swirl is the swirl data
366 * - values is the array of specifications
369 set_black_and_white(SWIRL_P swirl, XColor * values)
371 unsigned long white, black;
373 /* where is black and white? */
374 white = swirl->white;
375 black = swirl->black;
377 /* set black and white up */
378 values[white].flags = DoRed | DoGreen | DoBlue;
379 values[white].pixel = white;
380 values[white].red = 0xFFFF;
381 values[white].green = 0xFFFF;
382 values[white].blue = 0xFFFF;
383 values[black].flags = DoRed | DoGreen | DoBlue;
384 values[black].pixel = black;
385 values[black].red = 0;
386 values[black].green = 0;
387 values[black].blue = 0;
389 /* copy the colour specs from the original entries */
390 values[swirl->fg] = swirl->fgcol;
391 values[swirl->bg] = swirl->bgcol;
394 /****************************************************************/
399 * Set an entry in an array of XColor specifications. The given entry will be
400 * set to the given colour. If the entry corresponds to the foreground,
401 * background, WhitePixel, or BlackPixel it is ignored and the given colour
402 * is is put in the next entry.
404 * Therefore, the given colour may be placed up to four places after the
405 * specified entry in the array, if foreground, background, white, or black
408 * - swirl is the swirl data
409 * - value points to a pointer to the array entry. It gets updated to
410 * point to the next free entry.
411 * - pixel points to the current pixel number. It gets updated.
412 * - c points to the colour to add
415 set_colour(SWIRL_P swirl, XColor ** value, unsigned long *pixel, COLOUR_P c)
418 unsigned long fg, bg, white, black;
420 /* where are foreground, background, white, and black? */
423 white = swirl->white;
424 black = swirl->black;
426 /* haven't set it yet */
429 /* try and set the colour */
431 (**value).flags = DoRed | DoGreen | DoBlue;
432 (**value).pixel = *pixel;
434 /* white, black, fg, bg, or a colour? */
435 if ((*pixel != fg) && (*pixel != bg) &&
436 (*pixel != white) && (*pixel != black)) {
437 (**value).red = c->r;
438 (**value).green = c->g;
439 (**value).blue = c->b;
441 /* now we've done it */
450 /****************************************************************/
455 * Get an entry from an array of XColor specifications. The next colour from
456 * the array will be returned. Foreground, background, WhitePixel, or
457 * BlackPixel will be ignored.
459 * - swirl is the swirl data
460 * - value points the array entry. It is updated to point to the entry
461 * following the one returned.
462 * - c is set to the colour found
465 get_colour(SWIRL_P swirl, XColor ** value, COLOUR_P c)
468 unsigned long fg, bg, white, black;
470 /* where is white and black? */
473 white = swirl->white;
474 black = swirl->black;
476 /* haven't set it yet */
479 /* try and set the colour */
481 /* black, white or a colour? */
482 if (((*value)->pixel != fg) && ((*value)->pixel != bg) &&
483 ((*value)->pixel != white) && ((*value)->pixel != black)) {
484 c->r = (*value)->red;
485 c->g = (*value)->green;
486 c->b = (*value)->blue;
488 /* now we've done it */
495 #endif /* !STANDALONE */
497 /****************************************************************/
503 * Generate n colours between c1 and c2. n XColors at *value are set up with
504 * ascending pixel values.
506 * If the pixel range includes BlackPixel or WhitePixel they are set to black
507 * and white respectively but otherwise ignored. Therefore, up to n+2 colours
508 * may actually be set by this function.
510 * - swirl is the swirl data
511 * - values points a pointer to an array of XColors to update
512 * - pixel points to the pixel number to start at
513 * - k n is the number of colours to generate
514 * - c1, c2 are the colours to interpolate between
517 interpolate(SWIRL_P swirl, XColor ** values, unsigned long *pixel, int n, COLOUR_P c1, COLOUR_P c2)
526 for (i = 0; i < n / 2 && (int) *pixel < swirl->colours; i++) {
527 /* work out the colour */
528 r = c1->r + 2 * i * ((int) c2->r) / n;
529 c.r = (r > (int) maxv) ? maxv : r;
530 g = c1->g + 2 * i * ((int) c2->g) / n;
531 c.g = (g > (int) maxv) ? maxv : g;
532 b = c1->b + 2 * i * ((int) c2->b) / n;
533 c.b = (b > (int) maxv) ? maxv : b;
536 set_colour(swirl, values, pixel, &c);
538 for (i = n / 2; i >= 0 && (int) *pixel < swirl->colours; i--) {
539 r = c2->r + 2 * i * ((int) c1->r) / n;
540 c.r = (r > (int) maxv) ? maxv : r;
541 g = c2->g + 2 * i * ((int) c1->g) / n;
542 c.g = (g > (int) maxv) ? maxv : g;
543 b = c2->b + 2 * i * ((int) c1->b) / n;
544 c.b = (b > (int) maxv) ? maxv : b;
547 set_colour(swirl, values, pixel, &c);
551 /****************************************************************/
556 * Generate a `random' closed loop colourmap that occupies the whole colour
559 * - swirl is the swirl data
560 * - values is the array of colour definitions to set up
563 basic_map(SWIRL_P swirl, XColor * values)
567 unsigned short r1, g1, b1, r2, g2, b2, r3, g3, b3;
572 /* start at the beginning of the colour map */
576 /* choose 3 different basic colours at random */
577 for (i = 0; i < 3;) {
581 /* choose colour i */
582 c[i] = basic_colours[random_no(SWIRLCOLOURS - 1)];
584 /* assume different */
587 /* different from the rest? */
588 for (j = 0; j < i; j++)
589 if ((c[i].r == c[j].r) &&
590 (c[i].g == c[j].g) &&
594 /* ready for the next colour? */
599 /* extract components into variables */
610 /* work out the lengths of each side of the triangle */
611 L1 = (int) sqrt((((double) r1 - r2) * ((double) r1 - r2) +
612 ((double) g1 - g2) * ((double) g1 - g2) +
613 ((double) b1 - b2) * ((double) b1 - b2)));
615 L2 = (int) sqrt((((double) r3 - r2) * ((double) r3 - r2) +
616 ((double) g3 - g2) * ((double) g3 - g2) +
617 ((double) b3 - b2) * ((double) b3 - b2)));
619 L3 = (int) sqrt((((double) r1 - r3) * ((double) r1 - r3) +
620 ((double) g1 - g3) * ((double) g1 - g3) +
621 ((double) b1 - b3) * ((double) b1 - b3)));
625 /* allocate colours in proportion to the lengths of the sides */
626 interpolate(swirl, &value, &pixel,
627 (int) ((double) swirl->dcolours * ((double) L1 / (double) L)) + 1, c, c + 1);
628 interpolate(swirl, &value, &pixel,
629 (int) ((double) swirl->dcolours * ((double) L2 / (double) L)) + 1, c + 1, c + 2);
630 interpolate(swirl, &value, &pixel,
631 (int) ((double) swirl->dcolours * ((double) L3 / (double) L)) + 1, c + 2, c);
633 /* fill up any remaining slots (due to rounding) */
634 while ((int) pixel < swirl->colours) {
635 /* repeat the last colour */
636 set_colour(swirl, &value, &pixel, c);
639 /* ensure black and white are correct */
640 if (!swirl->fixed_colourmap)
641 set_black_and_white(swirl, values);
644 /****************************************************************/
649 * Generate pre-rotated versions of the colour specifications
651 * - swirl is the swirl data
652 * - values is an array of colour specifications
655 pre_rotate(SWIRL_P swirl, XColor * values)
662 /* how many colours to display? */
663 dcolours = swirl->dcolours;
665 /* start at the first map */
667 dest = values + swirl->colours;
669 /* generate dcolours-1 rotated maps */
670 for (i = 0; i < dcolours - 1; i++) {
673 /* start at the first pixel */
676 /* remember the first one and skip it */
677 get_colour(swirl, &src, &first);
679 /* put a rotated version of src at dest */
680 for (j = 0; j < dcolours - 1; j++) {
683 /* get the source colour */
684 get_colour(swirl, &src, &c);
687 set_colour(swirl, &dest, &pixel, &c);
690 /* put the first one at the end */
691 set_colour(swirl, &dest, &pixel, &first);
693 /* NB: src and dest should now be ready for the next table */
695 /* ensure black and white are properly set */
696 set_black_and_white(swirl, src);
700 /****************************************************************/
705 * Create a read/write colourmap to use
707 * - swirl is the swirl data
711 create_colourmap(ModeInfo * mi, SWIRL_P swirl)
713 Display *display = MI_DISPLAY(mi);
718 unsigned long redmask, greenmask, bluemask;
720 swirl->fixed_colourmap = !setupColormap(mi, &(swirl->colours),
721 &truecolor, &redmask, &greenmask, &bluemask);
722 preserve = preserveColors(swirl->fg, swirl->bg, swirl->white, swirl->black);
724 /* how many colours should we animate? */
725 swirl->dcolours = (swirl->colours > preserve + 1) ?
726 swirl->colours - preserve : swirl->colours;
728 if (MI_NPIXELS(mi) < 2)
731 /* how fast to shift the colourmap? */
732 swirl->shift = (swirl->colours > 64) ? swirl->colours / 64 : 1;
733 swirl->dshift = (swirl->shift > 1) ? swirl->shift * 2 : 1;
735 /* how may colour map rotations are there? */
736 n_rotations = (swirl->fixed_colourmap) ? 1 : swirl->dcolours;
738 /* allocate space for colour definitions (if not already there) */
739 if (swirl->rgb_values == NULL) {
740 swirl->rgb_values = (XColor *) calloc((swirl->colours + 3) * n_rotations,
743 /* create a colour map */
744 if (!swirl->fixed_colourmap)
746 XCreateColormap(display, swirl->win, swirl->visual, AllocAll);
748 /* select a set of colours for the colour map */
749 basic_map(swirl, swirl->rgb_values);
751 /* are we rotating them? */
752 if (!swirl->fixed_colourmap) {
753 /* generate rotations of the colour maps */
754 pre_rotate(swirl, swirl->rgb_values);
756 /* store the colours in the colour map */
757 XStoreColors(display, swirl->cmap, swirl->rgb_values, swirl->colours);
764 for (i = 0; (int) t > 0; i++, t >>= 1);
767 for (i = 0; (int) t > 0; i++, t >>= 1);
770 for (i = 0; (int) t > 0; i++, t >>= 1);
772 for (i = 0; i < swirl->colours; i++)
773 swirl->rgb_values[i].pixel =
774 ((rsh > 0 ? (swirl->rgb_values[i].red) >> rsh :
775 (swirl->rgb_values[i].red) << (-rsh)) & redmask) |
776 ((gsh > 0 ? (swirl->rgb_values[i].green) >> gsh :
777 (swirl->rgb_values[i].green) << (-gsh)) & greenmask) |
778 ((bsh > 0 ? (swirl->rgb_values[i].blue) >> bsh :
779 (swirl->rgb_values[i].blue) << (-bsh)) & bluemask);
781 /* lookup the colours in the fixed colour map */
782 for (i = 0; i < swirl->colours; i++)
783 (void) XAllocColor(display, MI_COLORMAP(mi),
784 &(swirl->rgb_values[i]));
789 /****************************************************************/
794 * Install a new set of colours into the colour map
796 * - dpy is the display
797 * - swirl is the swirl data
798 * - shift is the amount to rotate the colour map by
801 install_map(Display * dpy, SWIRL_P swirl, int shift)
803 if (!swirl->fixed_colourmap) {
804 /* shift the colour map */
805 swirl->current_map = (swirl->current_map + shift) %
809 XStoreColors(dpy, swirl->cmap,
811 swirl->current_map * swirl->colours,
815 #endif /* !STANDALONE */
817 /****************************************************************/
822 * Initialise the array of knot
824 * swirl is the swirl data
827 create_knots(SWIRL_P swirl)
830 Bool orbit, wheel, picasso, ray, hook;
833 /* create array for knots */
835 (void) free((void *) swirl->knots);
836 swirl->knots = (KNOT_P) calloc(swirl->n_knots, sizeof (KNOT));
839 orbit = wheel = picasso = ray = hook = False;
841 /* what types do we have? */
842 if ((int) swirl->knot_type & (int) ALL) {
843 orbit = wheel = ray = hook = True;
845 if ((int) swirl->knot_type & (int) ORBIT)
847 if ((int) swirl->knot_type & (int) WHEEL)
849 if ((int) swirl->knot_type & (int) PICASSO)
851 if ((int) swirl->knot_type & (int) RAY)
853 if ((int) swirl->knot_type & (int) HOOK)
857 /* initialise each knot */
859 for (k = 0; k < swirl->n_knots; k++) {
861 knot->x = random_no((unsigned int) swirl->width);
862 knot->y = random_no((unsigned int) swirl->height);
865 knot->m = random_no(MASS) + 1;
867 /* can be negative */
868 if (random_no(100) > 50)
873 while (knot->t == NONE) {
874 /* choose a random one from the types available */
875 switch (random_no(4)) {
899 /* if two planes, do same for second plane */
900 if (swirl->two_plane) {
902 while (knot->T == NONE || knot->T == knot->t) {
903 /* choose a different type */
904 switch (random_no(4)) {
933 /****************************************************************/
938 * Work out the pixel value at i, j. Ensure it does not clash with BlackPixel
941 * - swirl is the swirl data
942 * - i, j is the point to calculate
944 * Returns the value of the point
947 do_point(SWIRL_P swirl, int i, int j)
949 int tT, k, value, add;
950 double dx, dy, theta, dist;
951 int dcolours, qcolours;
955 /* how many colours? */
956 dcolours = swirl->dcolours;
957 qcolours = dcolours / 4;
959 /* colour step round a circle */
960 rads = (double) dcolours / (2.0 * M_PI);
965 /* go through all the knots */
967 for (k = 0; k < swirl->n_knots; k++) {
971 /* in two_plane mode get the appropriate knot type */
972 if (swirl->two_plane)
973 tT = (int) ((swirl->first_plane) ? knot->t : knot->T);
977 /* distance from knot */
978 dist = sqrt(dx * dx + dy * dy);
980 /* nothing to add at first */
983 /* work out the contribution (if close enough) */
987 add = (int) (dcolours / (1.0 + 0.01 * abs(knot->m) * dist));
990 /* Avoid atan2: DOMAIN error message */
991 if (dy == 0.0 && dx == 0.0)
994 theta = (atan2(dy, dx) + M_PI) / M_PI;
996 add = (int) (dcolours * theta +
997 sin(0.1 * knot->m * dist) *
998 qcolours * exp(-0.01 * dist));
1000 add = (int) (dcolours * (theta - 1.0) +
1001 sin(0.1 * knot->m * dist) *
1002 qcolours * exp(-0.01 * dist));
1005 add = (int) (dcolours *
1006 fabs(cos(0.002 * knot->m * dist)));
1009 /* Avoid atan2: DOMAIN error message */
1010 if (dy == 0.0 && dx == 0.0)
1013 add = (int) (dcolours * fabs(sin(2.0 * atan2(dy, dx))));
1017 /* Avoid atan2: DOMAIN error message */
1018 if (dy == 0.0 && dx == 0.0)
1019 add = (int) (0.05 * (abs(knot->m) - 1) * dist);
1021 add = (int) (rads * atan2(dy, dx) +
1022 0.05 * (abs(knot->m) - 1) * dist);
1025 /* for a +ve mass add on the contribution else take it off */
1036 swirl->first_plane = (!swirl->first_plane);
1038 /* make sure we handle -ve values properly */
1040 value = (value % dcolours) + 2;
1042 value = dcolours - (abs(value) % (dcolours - 1));
1045 /* if fg and bg are 1 and 0 we should be OK, but just in case */
1046 while ((dcolours > 2) &&
1047 (((value % swirl->colours) == (int) swirl->fg) ||
1048 ((value % swirl->colours) == (int) swirl->bg) ||
1049 ((value % swirl->colours) == (int) swirl->white) ||
1050 ((value % swirl->colours) == (int) swirl->black))) {
1053 #endif /* !STANDALONE */
1055 /* definitely make sure it is in range */
1056 value = value % swirl->colours;
1058 /* lookup the pixel value if necessary */
1060 if (swirl->fixed_colourmap && swirl->dcolours > 2)
1062 value = swirl->rgb_values[value].pixel;
1065 return ((unsigned long) value);
1068 /****************************************************************/
1073 * Draw a square block of points with the same value.
1075 * - ximage is the XImage to draw on.
1076 * - x, y is the top left corner
1077 * - s is the length of each side
1081 draw_block(XImage * ximage, int x, int y, int s, unsigned long v)
1085 for (a = 0; a < s; a++)
1086 for (b = 0; b < s; b++) {
1087 XPutPixel(ximage, x + b, y + a, v);
1091 /****************************************************************/
1094 * draw_point Draw the current point in a swirl pattern onto the XImage
1096 * - swirl is the swirl
1097 * - win is the window to update
1100 draw_point(ModeInfo * mi, SWIRL_P swirl)
1105 /* get current point coordinates and resolution */
1110 /* check we are within the window */
1111 if ((x < 0) || (x > swirl->width - r) || (y < 0) || (y > swirl->height - r))
1114 /* what style are we drawing? */
1115 if (swirl->two_plane) {
1118 /* halve the block size */
1121 /* interleave blocks at half r */
1122 draw_block(swirl->ximage, x, y, r2, do_point(swirl, x, y));
1123 draw_block(swirl->ximage, x + r2, y, r2, do_point(swirl, x + r2, y));
1124 draw_block(swirl->ximage, x + r2, y + r2, r2, do_point(swirl,
1126 draw_block(swirl->ximage, x, y + r2, r2, do_point(swirl, x, y + r2));
1128 draw_block(swirl->ximage, x, y, r, do_point(swirl, x, y));
1130 /* update the screen */
1132 #ifdef HAVE_XSHM_EXTENSION
1134 XShmPutImage(MI_DISPLAY(mi), MI_WINDOW(mi), MI_GC(mi), swirl->ximage,
1135 x, y, x, y, r, r, False);
1137 #endif /* !HAVE_XSHM_EXTENSION */
1138 /* PURIFY 4.0.1 on SunOS4 and on Solaris 2 reports a 256 byte memory
1139 leak on the next line. */
1140 XPutImage(MI_DISPLAY(mi), MI_WINDOW(mi), MI_GC(mi), swirl->ximage,
1144 /****************************************************************/
1147 * next_point Move to the next point in the spiral pattern
1148 * - swirl is the swirl
1149 * - win is the window to update
1152 next_point(SWIRL_P swirl)
1154 /* more to do in this direction? */
1155 if (swirl->dir_done < swirl->dir_todo) {
1156 /* move in the current direction */
1157 switch (swirl->direction) {
1159 swirl->x += swirl->r;
1162 swirl->y += swirl->r;
1165 swirl->x -= swirl->r;
1168 swirl->y -= swirl->r;
1172 /* done another point */
1175 /* none drawn yet */
1176 swirl->dir_done = 0;
1178 /* change direction - check and record if off screen */
1179 switch (swirl->direction) {
1181 swirl->direction = DRAW_DOWN;
1182 if (swirl->x > swirl->width - swirl->r) {
1183 /* skip these points */
1184 swirl->dir_done = swirl->dir_todo;
1185 swirl->y += (swirl->dir_todo * swirl->r);
1187 /* check for finish */
1188 if (swirl->off_screen)
1189 swirl->drawing = False;
1190 swirl->off_screen = True;
1192 swirl->off_screen = False;
1195 swirl->direction = DRAW_LEFT;
1197 if (swirl->y > swirl->height - swirl->r) {
1198 /* skip these points */
1199 swirl->dir_done = swirl->dir_todo;
1200 swirl->x -= (swirl->dir_todo * swirl->r);
1202 /* check for finish */
1203 if (swirl->off_screen)
1204 swirl->drawing = False;
1205 swirl->off_screen = True;
1207 swirl->off_screen = False;
1210 swirl->direction = DRAW_UP;
1212 /* skip these points */
1213 swirl->dir_done = swirl->dir_todo;
1214 swirl->y -= (swirl->dir_todo * swirl->r);
1216 /* check for finish */
1217 if (swirl->off_screen)
1218 swirl->drawing = False;
1219 swirl->off_screen = True;
1221 swirl->off_screen = False;
1224 swirl->direction = DRAW_RIGHT;
1227 /* skip these points */
1228 swirl->dir_done = swirl->dir_todo;
1229 swirl->x += (swirl->dir_todo * swirl->r);
1231 /* check for finish */
1232 if (swirl->off_screen)
1233 swirl->drawing = False;
1234 swirl->off_screen = True;
1236 swirl->off_screen = False;
1242 /****************************************************************/
1247 * Initialise things for swirling
1249 * - win is the window to draw in
1252 init_swirl(ModeInfo * mi)
1254 Display *display = MI_DISPLAY(mi);
1255 Window window = MI_WINDOW(mi);
1258 /* does the swirls array exist? */
1259 if (swirls == NULL) {
1262 /* allocate an array, one entry for each screen */
1263 swirls = (SWIRL_P) calloc(ScreenCount(display), sizeof (SWIRL));
1265 /* initialise them all */
1266 for (i = 0; i < ScreenCount(display); i++)
1267 initialise_swirl(mi, &swirls[i]);
1269 /* get a pointer to this swirl */
1270 swirl = &(swirls[MI_SCREEN(mi)]);
1272 /* get window parameters */
1273 swirl->win = window;
1274 swirl->width = MI_WIN_WIDTH(mi);
1275 swirl->height = MI_WIN_HEIGHT(mi);
1276 swirl->depth = MI_WIN_DEPTH(mi);
1277 swirl->rdepth = swirl->depth;
1278 swirl->visual = MI_VISUAL(mi);
1280 if (swirl->depth > 16)
1283 /* initialise image for speeding up drawing */
1284 initialise_image(mi, swirl);
1286 /* clear the window (before setting the colourmap) */
1287 XClearWindow(display, MI_WINDOW(mi));
1291 swirl->rgb_values = mi->colors;
1292 swirl->colours = mi->npixels;
1293 swirl->dcolours = swirl->colours;
1294 /* swirl->fixed_colourmap = !mi->writable_p;*/
1296 #else /* !STANDALONE */
1298 /* initialise the colours from which the colourmap is derived */
1299 initialise_colours(basic_colours, MI_SATURATION(mi));
1301 /* set up the colour map */
1302 create_colourmap(mi, swirl);
1304 /* attach the colour map to the window (if we have one) */
1305 if (!swirl->fixed_colourmap) {
1307 setColormap(display, window, swirl->cmap, MI_WIN_IS_INWINDOW(mi));
1309 XSetWindowColormap(display, window, swirl->cmap);
1310 (void) XSetWMColormapWindows(display, window, &window, 1);
1311 XInstallColormap(display, swirl->cmap);
1314 #endif /* STANDALONE */
1316 /* resolution starts off chunky */
1317 swirl->resolution = MIN_RES + 1;
1319 /* calculate the pixel step for this resulution */
1320 swirl->r = (1 << (swirl->resolution - 1));
1322 /* how many knots? */
1323 swirl->n_knots = random_no((unsigned int) MI_BATCHCOUNT(mi) / 2) +
1324 MI_BATCHCOUNT(mi) + 1;
1326 /* what type of knots? */
1327 swirl->knot_type = ALL; /* for now */
1329 /* use two_plane mode occaisionally */
1330 if (random_no(100) <= TWO_PLANE_PCNT) {
1331 swirl->two_plane = swirl->first_plane = True;
1332 swirl->max_resolution = 2;
1334 swirl->two_plane = False;
1336 /* fix the knot values */
1337 create_knots(swirl);
1340 swirl->started = True;
1341 swirl->drawing = False;
1344 /****************************************************************/
1349 * Draw one iteration of swirling
1351 * - win is the window to draw in
1354 draw_swirl(ModeInfo * mi)
1356 SWIRL_P swirl = &(swirls[MI_SCREEN(mi)]);
1359 if (swirl->started) {
1360 /* in the middle of drawing? */
1361 if (swirl->drawing) {
1364 rotate_colors(MI_DISPLAY(mi), MI_COLORMAP(mi),
1365 swirl->rgb_values, swirl->colours, 1);
1366 #else /* !STANDALONE */
1367 /* rotate the colours */
1368 install_map(MI_DISPLAY(mi), swirl, swirl->dshift);
1369 #endif /* !STANDALONE */
1371 /* draw a batch of points */
1372 swirl->batch_todo = BATCH_DRAW;
1373 while ((swirl->batch_todo > 0) && swirl->drawing) {
1375 draw_point(mi, swirl);
1377 /* move to the next point */
1381 swirl->batch_todo--;
1386 rotate_colors(MI_DISPLAY(mi), MI_COLORMAP(mi),
1387 swirl->rgb_values, swirl->colours, 1);
1388 #else /* !STANDALONE */
1389 /* rotate the colours */
1390 install_map(MI_DISPLAY(mi), swirl, swirl->shift);
1391 #endif /* !STANDALONE */
1393 /* time for a higher resolution? */
1394 if (swirl->resolution > swirl->max_resolution) {
1395 /* move to higher resolution */
1396 swirl->resolution--;
1398 /* calculate the pixel step for this resulution */
1399 swirl->r = (1 << (swirl->resolution - 1));
1401 /* start drawing again */
1402 swirl->drawing = True;
1404 /* start in the middle of the screen */
1405 swirl->x = (swirl->width - swirl->r) / 2;
1406 swirl->y = (swirl->height - swirl->r) / 2;
1408 /* initialise spiral drawing parameters */
1409 swirl->direction = DRAW_RIGHT;
1410 swirl->dir_todo = 1;
1411 swirl->dir_done = 0;
1413 /* all done, decide when to restart */
1414 if (swirl->start_again == -1) {
1415 /* start the counter */
1416 swirl->start_again = RESTART;
1417 } else if (swirl->start_again == 0) {
1418 /* reset the counter */
1419 swirl->start_again = -1;
1422 /* Pick a new colormap! */
1423 XClearWindow (MI_DISPLAY(mi), MI_WINDOW(mi));
1424 free_colors (MI_DISPLAY(mi), MI_COLORMAP(mi),
1425 mi->colors, mi->npixels);
1426 make_smooth_colormap (MI_DISPLAY(mi),
1429 mi->colors, &mi->npixels, True,
1430 &mi->writable_p, True);
1431 swirl->colours = mi->npixels;
1432 #endif /* STANDALONE */
1437 /* decrement the counter */
1438 swirl->start_again--;
1444 /****************************************************************/
1447 release_swirl(ModeInfo * mi)
1449 /* does the swirls array exist? */
1450 if (swirls != NULL) {
1454 for (i = 0; i < MI_NUM_SCREENS(mi); i++) {
1455 SWIRL_P swirl = &(swirls[i]);
1457 if (swirl->cmap != (Colormap) NULL)
1458 XFreeColormap(MI_DISPLAY(mi), swirl->cmap);
1459 if (swirl->rgb_values != NULL)
1460 XFree((void *) swirl->rgb_values);
1461 if (swirl->ximage != NULL)
1462 XDestroyImage(swirl->ximage);
1464 (void) free((void *) swirl->knots);
1466 /* deallocate an array, one entry for each screen */
1467 (void) free((void *) swirls);
1472 /****************************************************************/
1475 refresh_swirl(ModeInfo * mi)
1477 SWIRL_P swirl = &(swirls[MI_SCREEN(mi)]);
1479 if (swirl->started) {
1481 swirl->resolution = swirl->resolution + 1;
1482 swirl->drawing = False;