1 /* -*- Mode: C; tab-width: 4 -*-
2 * swirl --- swirly color-cycling patterns.
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 swirl->fg = MI_FG_COLOR(mi);
229 swirl->bg = MI_BG_COLOR(mi);
230 swirl->fgcol.pixel = swirl->fg;
231 swirl->bgcol.pixel = swirl->bg;
232 XQueryColor(display, MI_COLORMAP(mi), &(swirl->fgcol));
233 XQueryColor(display, MI_COLORMAP(mi), &(swirl->bgcol));
234 #endif /* !STANDALONE */
237 /****************************************************************/
242 * Initialise the image for drawing to
244 * - swirl is the swirl data
247 initialise_image(ModeInfo * mi, SWIRL_P swirl)
249 Display *dpy = MI_DISPLAY(mi);
251 if (swirl->ximage != NULL)
252 XDestroyImage(swirl->ximage);
255 #ifdef HAVE_XSHM_EXTENSION
258 swirl->ximage = create_xshm_image(dpy, swirl->visual, swirl->rdepth,
259 ZPixmap, 0, &mi->shm_info,
260 swirl->width, swirl->height);
264 #endif /* HAVE_XSHM_EXTENSION */
268 swirl->ximage = XCreateImage(dpy, swirl->visual, swirl->rdepth, ZPixmap,
269 0, 0, swirl->width, swirl->height,
271 swirl->image = (unsigned char *)
272 calloc(swirl->height, swirl->ximage->bytes_per_line);
273 swirl->ximage->data = (char *) swirl->image;
277 /****************************************************************/
283 * Initialise the list of colours from which the colourmaps are derived
285 * - colours is the array to initialise
286 * - saturation is the saturation value to use 0->grey,
287 * 1.0->full saturation
290 initialise_colours(COLOUR * colours, float saturate)
294 /* start off fully saturated, medium and bright colours */
295 colours[0].r = 0xA000;
296 colours[0].g = 0x0000;
297 colours[0].b = 0x0000;
298 colours[1].r = 0xD000;
299 colours[1].g = 0x0000;
300 colours[1].b = 0x0000;
301 colours[2].r = 0x0000;
302 colours[2].g = 0x6000;
303 colours[2].b = 0x0000;
304 colours[3].r = 0x0000;
305 colours[3].g = 0x9000;
306 colours[3].b = 0x0000;
307 colours[4].r = 0x0000;
308 colours[4].g = 0x0000;
309 colours[4].b = 0xC000;
310 colours[5].r = 0x0000;
311 colours[5].g = 0x0000;
312 colours[5].b = 0xF000;
313 colours[6].r = 0xA000;
314 colours[6].g = 0x6000;
315 colours[6].b = 0x0000;
316 colours[7].r = 0xD000;
317 colours[7].g = 0x9000;
318 colours[7].b = 0x0000;
319 colours[8].r = 0xA000;
320 colours[8].g = 0x0000;
321 colours[8].b = 0xC000;
322 colours[9].r = 0xD000;
323 colours[9].g = 0x0000;
324 colours[9].b = 0xF000;
325 colours[10].r = 0x0000;
326 colours[10].g = 0x6000;
327 colours[10].b = 0xC000;
328 colours[11].r = 0x0000;
329 colours[11].g = 0x9000;
330 colours[11].b = 0xF000;
331 colours[12].r = 0xA000;
332 colours[12].g = 0xA000;
333 colours[12].b = 0xA000;
335 /* add white for low saturation */
336 for (i = 0; i < SWIRLCOLOURS - 1; i++) {
337 unsigned short max_rg, max;
339 /* what is the max intensity for this colour? */
340 max_rg = (colours[i].r > colours[i].g) ? colours[i].r : colours[i].g;
341 max = (max_rg > colours[i].b) ? max_rg : colours[i].b;
343 /* bring elements up to max as saturation approaches 0.0 */
344 colours[i].r += (unsigned short) ((float) (1.0 - saturate) *
345 ((float) max - colours[i].r));
346 colours[i].g += (unsigned short) ((float) (1.0 - saturate) *
347 ((float) max - colours[i].g));
348 colours[i].b += (unsigned short) ((float) (1.0 - saturate) *
349 ((float) max - colours[i].b));
352 #endif /* !STANDALONE */
354 /****************************************************************/
358 * set_black_and_white
360 * Set the entries for foreground & background pixels and
361 * WhitePixel & BlackPixel in an array of colour specifications.
363 * - swirl is the swirl data
364 * - values is the array of specifications
367 set_black_and_white(SWIRL_P swirl, XColor * values)
369 unsigned long white, black;
371 /* where is black and white? */
372 white = swirl->white;
373 black = swirl->black;
375 /* set black and white up */
376 values[white].flags = DoRed | DoGreen | DoBlue;
377 values[white].pixel = white;
378 values[white].red = 0xFFFF;
379 values[white].green = 0xFFFF;
380 values[white].blue = 0xFFFF;
381 values[black].flags = DoRed | DoGreen | DoBlue;
382 values[black].pixel = black;
383 values[black].red = 0;
384 values[black].green = 0;
385 values[black].blue = 0;
387 /* copy the colour specs from the original entries */
388 values[swirl->fg] = swirl->fgcol;
389 values[swirl->bg] = swirl->bgcol;
392 /****************************************************************/
397 * Set an entry in an array of XColor specifications. The given entry will be
398 * set to the given colour. If the entry corresponds to the foreground,
399 * background, WhitePixel, or BlackPixel it is ignored and the given colour
400 * is is put in the next entry.
402 * Therefore, the given colour may be placed up to four places after the
403 * specified entry in the array, if foreground, background, white, or black
406 * - swirl is the swirl data
407 * - value points to a pointer to the array entry. It gets updated to
408 * point to the next free entry.
409 * - pixel points to the current pixel number. It gets updated.
410 * - c points to the colour to add
413 set_colour(SWIRL_P swirl, XColor ** value, unsigned long *pixel, COLOUR_P c)
416 unsigned long fg, bg, white, black;
418 /* where are foreground, background, white, and black? */
421 white = swirl->white;
422 black = swirl->black;
424 /* haven't set it yet */
427 /* try and set the colour */
429 (**value).flags = DoRed | DoGreen | DoBlue;
430 (**value).pixel = *pixel;
432 /* white, black, fg, bg, or a colour? */
433 if ((*pixel != fg) && (*pixel != bg) &&
434 (*pixel != white) && (*pixel != black)) {
435 (**value).red = c->r;
436 (**value).green = c->g;
437 (**value).blue = c->b;
439 /* now we've done it */
448 /****************************************************************/
453 * Get an entry from an array of XColor specifications. The next colour from
454 * the array will be returned. Foreground, background, WhitePixel, or
455 * BlackPixel will be ignored.
457 * - swirl is the swirl data
458 * - value points the array entry. It is updated to point to the entry
459 * following the one returned.
460 * - c is set to the colour found
463 get_colour(SWIRL_P swirl, XColor ** value, COLOUR_P c)
466 unsigned long fg, bg, white, black;
468 /* where is white and black? */
471 white = swirl->white;
472 black = swirl->black;
474 /* haven't set it yet */
477 /* try and set the colour */
479 /* black, white or a colour? */
480 if (((*value)->pixel != fg) && ((*value)->pixel != bg) &&
481 ((*value)->pixel != white) && ((*value)->pixel != black)) {
482 c->r = (*value)->red;
483 c->g = (*value)->green;
484 c->b = (*value)->blue;
486 /* now we've done it */
493 #endif /* !STANDALONE */
495 /****************************************************************/
501 * Generate n colours between c1 and c2. n XColors at *value are set up with
502 * ascending pixel values.
504 * If the pixel range includes BlackPixel or WhitePixel they are set to black
505 * and white respectively but otherwise ignored. Therefore, up to n+2 colours
506 * may actually be set by this function.
508 * - swirl is the swirl data
509 * - values points a pointer to an array of XColors to update
510 * - pixel points to the pixel number to start at
511 * - k n is the number of colours to generate
512 * - c1, c2 are the colours to interpolate between
515 interpolate(SWIRL_P swirl, XColor ** values, unsigned long *pixel, int n, COLOUR_P c1, COLOUR_P c2)
524 for (i = 0; i < n / 2 && (int) *pixel < swirl->colours; i++) {
525 /* work out the colour */
526 r = c1->r + 2 * i * ((int) c2->r) / n;
527 c.r = (r > (int) maxv) ? maxv : r;
528 g = c1->g + 2 * i * ((int) c2->g) / n;
529 c.g = (g > (int) maxv) ? maxv : g;
530 b = c1->b + 2 * i * ((int) c2->b) / n;
531 c.b = (b > (int) maxv) ? maxv : b;
534 set_colour(swirl, values, pixel, &c);
536 for (i = n / 2; i >= 0 && (int) *pixel < swirl->colours; i--) {
537 r = c2->r + 2 * i * ((int) c1->r) / n;
538 c.r = (r > (int) maxv) ? maxv : r;
539 g = c2->g + 2 * i * ((int) c1->g) / n;
540 c.g = (g > (int) maxv) ? maxv : g;
541 b = c2->b + 2 * i * ((int) c1->b) / n;
542 c.b = (b > (int) maxv) ? maxv : b;
545 set_colour(swirl, values, pixel, &c);
549 /****************************************************************/
554 * Generate a `random' closed loop colourmap that occupies the whole colour
557 * - swirl is the swirl data
558 * - values is the array of colour definitions to set up
561 basic_map(SWIRL_P swirl, XColor * values)
565 unsigned short r1, g1, b1, r2, g2, b2, r3, g3, b3;
570 /* start at the beginning of the colour map */
574 /* choose 3 different basic colours at random */
575 for (i = 0; i < 3;) {
579 /* choose colour i */
580 c[i] = basic_colours[random_no(SWIRLCOLOURS - 1)];
582 /* assume different */
585 /* different from the rest? */
586 for (j = 0; j < i; j++)
587 if ((c[i].r == c[j].r) &&
588 (c[i].g == c[j].g) &&
592 /* ready for the next colour? */
597 /* extract components into variables */
608 /* work out the lengths of each side of the triangle */
609 L1 = (int) sqrt((((double) r1 - r2) * ((double) r1 - r2) +
610 ((double) g1 - g2) * ((double) g1 - g2) +
611 ((double) b1 - b2) * ((double) b1 - b2)));
613 L2 = (int) sqrt((((double) r3 - r2) * ((double) r3 - r2) +
614 ((double) g3 - g2) * ((double) g3 - g2) +
615 ((double) b3 - b2) * ((double) b3 - b2)));
617 L3 = (int) sqrt((((double) r1 - r3) * ((double) r1 - r3) +
618 ((double) g1 - g3) * ((double) g1 - g3) +
619 ((double) b1 - b3) * ((double) b1 - b3)));
623 /* allocate colours in proportion to the lengths of the sides */
624 interpolate(swirl, &value, &pixel,
625 (int) ((double) swirl->dcolours * ((double) L1 / (double) L)) + 1, c, c + 1);
626 interpolate(swirl, &value, &pixel,
627 (int) ((double) swirl->dcolours * ((double) L2 / (double) L)) + 1, c + 1, c + 2);
628 interpolate(swirl, &value, &pixel,
629 (int) ((double) swirl->dcolours * ((double) L3 / (double) L)) + 1, c + 2, c);
631 /* fill up any remaining slots (due to rounding) */
632 while ((int) pixel < swirl->colours) {
633 /* repeat the last colour */
634 set_colour(swirl, &value, &pixel, c);
637 /* ensure black and white are correct */
638 if (!swirl->fixed_colourmap)
639 set_black_and_white(swirl, values);
642 /****************************************************************/
647 * Generate pre-rotated versions of the colour specifications
649 * - swirl is the swirl data
650 * - values is an array of colour specifications
653 pre_rotate(SWIRL_P swirl, XColor * values)
660 /* how many colours to display? */
661 dcolours = swirl->dcolours;
663 /* start at the first map */
665 dest = values + swirl->colours;
667 /* generate dcolours-1 rotated maps */
668 for (i = 0; i < dcolours - 1; i++) {
671 /* start at the first pixel */
674 /* remember the first one and skip it */
675 get_colour(swirl, &src, &first);
677 /* put a rotated version of src at dest */
678 for (j = 0; j < dcolours - 1; j++) {
681 /* get the source colour */
682 get_colour(swirl, &src, &c);
685 set_colour(swirl, &dest, &pixel, &c);
688 /* put the first one at the end */
689 set_colour(swirl, &dest, &pixel, &first);
691 /* NB: src and dest should now be ready for the next table */
693 /* ensure black and white are properly set */
694 set_black_and_white(swirl, src);
698 /****************************************************************/
703 * Create a read/write colourmap to use
705 * - swirl is the swirl data
709 create_colourmap(ModeInfo * mi, SWIRL_P swirl)
711 Display *display = MI_DISPLAY(mi);
716 unsigned long redmask, greenmask, bluemask;
718 swirl->fixed_colourmap = !setupColormap(mi, &(swirl->colours),
719 &truecolor, &redmask, &greenmask, &bluemask);
720 preserve = preserveColors(swirl->fg, swirl->bg, swirl->white, swirl->black);
722 /* how many colours should we animate? */
723 swirl->dcolours = (swirl->colours > preserve + 1) ?
724 swirl->colours - preserve : swirl->colours;
726 if (MI_NPIXELS(mi) < 2)
729 /* how fast to shift the colourmap? */
730 swirl->shift = (swirl->colours > 64) ? swirl->colours / 64 : 1;
731 swirl->dshift = (swirl->shift > 1) ? swirl->shift * 2 : 1;
733 /* how may colour map rotations are there? */
734 n_rotations = (swirl->fixed_colourmap) ? 1 : swirl->dcolours;
736 /* allocate space for colour definitions (if not already there) */
737 if (swirl->rgb_values == NULL) {
738 swirl->rgb_values = (XColor *) calloc((swirl->colours + 3) * n_rotations,
741 /* create a colour map */
742 if (!swirl->fixed_colourmap)
744 XCreateColormap(display, swirl->win, swirl->visual, AllocAll);
746 /* select a set of colours for the colour map */
747 basic_map(swirl, swirl->rgb_values);
749 /* are we rotating them? */
750 if (!swirl->fixed_colourmap) {
751 /* generate rotations of the colour maps */
752 pre_rotate(swirl, swirl->rgb_values);
754 /* store the colours in the colour map */
755 XStoreColors(display, swirl->cmap, swirl->rgb_values, swirl->colours);
762 for (i = 0; (int) t > 0; i++, t >>= 1);
765 for (i = 0; (int) t > 0; i++, t >>= 1);
768 for (i = 0; (int) t > 0; i++, t >>= 1);
770 for (i = 0; i < swirl->colours; i++)
771 swirl->rgb_values[i].pixel =
772 ((rsh > 0 ? (swirl->rgb_values[i].red) >> rsh :
773 (swirl->rgb_values[i].red) << (-rsh)) & redmask) |
774 ((gsh > 0 ? (swirl->rgb_values[i].green) >> gsh :
775 (swirl->rgb_values[i].green) << (-gsh)) & greenmask) |
776 ((bsh > 0 ? (swirl->rgb_values[i].blue) >> bsh :
777 (swirl->rgb_values[i].blue) << (-bsh)) & bluemask);
779 /* lookup the colours in the fixed colour map */
780 for (i = 0; i < swirl->colours; i++)
781 (void) XAllocColor(display, MI_COLORMAP(mi),
782 &(swirl->rgb_values[i]));
787 /****************************************************************/
792 * Install a new set of colours into the colour map
794 * - dpy is the display
795 * - swirl is the swirl data
796 * - shift is the amount to rotate the colour map by
799 install_map(Display * dpy, SWIRL_P swirl, int shift)
801 if (!swirl->fixed_colourmap) {
802 /* shift the colour map */
803 swirl->current_map = (swirl->current_map + shift) %
807 XStoreColors(dpy, swirl->cmap,
809 swirl->current_map * swirl->colours,
813 #endif /* !STANDALONE */
815 /****************************************************************/
820 * Initialise the array of knot
822 * swirl is the swirl data
825 create_knots(SWIRL_P swirl)
828 Bool orbit, wheel, picasso, ray, hook;
831 /* create array for knots */
833 (void) free((void *) swirl->knots);
834 swirl->knots = (KNOT_P) calloc(swirl->n_knots, sizeof (KNOT));
837 orbit = wheel = picasso = ray = hook = False;
839 /* what types do we have? */
840 if ((int) swirl->knot_type & (int) ALL) {
841 orbit = wheel = ray = hook = True;
843 if ((int) swirl->knot_type & (int) ORBIT)
845 if ((int) swirl->knot_type & (int) WHEEL)
847 if ((int) swirl->knot_type & (int) PICASSO)
849 if ((int) swirl->knot_type & (int) RAY)
851 if ((int) swirl->knot_type & (int) HOOK)
855 /* initialise each knot */
857 for (k = 0; k < swirl->n_knots; k++) {
859 knot->x = random_no((unsigned int) swirl->width);
860 knot->y = random_no((unsigned int) swirl->height);
863 knot->m = random_no(MASS) + 1;
865 /* can be negative */
866 if (random_no(100) > 50)
871 while (knot->t == NONE) {
872 /* choose a random one from the types available */
873 switch (random_no(4)) {
897 /* if two planes, do same for second plane */
898 if (swirl->two_plane) {
900 while (knot->T == NONE || knot->T == knot->t) {
901 /* choose a different type */
902 switch (random_no(4)) {
931 /****************************************************************/
936 * Work out the pixel value at i, j. Ensure it does not clash with BlackPixel
939 * - swirl is the swirl data
940 * - i, j is the point to calculate
942 * Returns the value of the point
945 do_point(SWIRL_P swirl, int i, int j)
947 int tT, k, value, add;
948 double dx, dy, theta, dist;
949 int dcolours, qcolours;
953 /* how many colours? */
954 dcolours = swirl->dcolours;
955 qcolours = dcolours / 4;
957 /* colour step round a circle */
958 rads = (double) dcolours / (2.0 * M_PI);
963 /* go through all the knots */
965 for (k = 0; k < swirl->n_knots; k++) {
969 /* in two_plane mode get the appropriate knot type */
970 if (swirl->two_plane)
971 tT = (int) ((swirl->first_plane) ? knot->t : knot->T);
975 /* distance from knot */
976 dist = sqrt(dx * dx + dy * dy);
978 /* nothing to add at first */
981 /* work out the contribution (if close enough) */
985 add = (int) (dcolours / (1.0 + 0.01 * abs(knot->m) * dist));
988 /* Avoid atan2: DOMAIN error message */
989 if (dy == 0.0 && dx == 0.0)
992 theta = (atan2(dy, dx) + M_PI) / M_PI;
994 add = (int) (dcolours * theta +
995 sin(0.1 * knot->m * dist) *
996 qcolours * exp(-0.01 * dist));
998 add = (int) (dcolours * (theta - 1.0) +
999 sin(0.1 * knot->m * dist) *
1000 qcolours * exp(-0.01 * dist));
1003 add = (int) (dcolours *
1004 fabs(cos(0.002 * knot->m * dist)));
1007 /* Avoid atan2: DOMAIN error message */
1008 if (dy == 0.0 && dx == 0.0)
1011 add = (int) (dcolours * fabs(sin(2.0 * atan2(dy, dx))));
1015 /* Avoid atan2: DOMAIN error message */
1016 if (dy == 0.0 && dx == 0.0)
1017 add = (int) (0.05 * (abs(knot->m) - 1) * dist);
1019 add = (int) (rads * atan2(dy, dx) +
1020 0.05 * (abs(knot->m) - 1) * dist);
1023 /* for a +ve mass add on the contribution else take it off */
1034 swirl->first_plane = (!swirl->first_plane);
1036 /* make sure we handle -ve values properly */
1038 value = (value % dcolours) + 2;
1040 value = dcolours - (abs(value) % (dcolours - 1));
1043 /* if fg and bg are 1 and 0 we should be OK, but just in case */
1044 while ((dcolours > 2) &&
1045 (((value % swirl->colours) == (int) swirl->fg) ||
1046 ((value % swirl->colours) == (int) swirl->bg) ||
1047 ((value % swirl->colours) == (int) swirl->white) ||
1048 ((value % swirl->colours) == (int) swirl->black))) {
1051 #endif /* !STANDALONE */
1053 /* definitely make sure it is in range */
1054 value = value % swirl->colours;
1056 /* lookup the pixel value if necessary */
1058 if (swirl->fixed_colourmap && swirl->dcolours > 2)
1060 value = swirl->rgb_values[value].pixel;
1063 return ((unsigned long) value);
1066 /****************************************************************/
1071 * Draw a square block of points with the same value.
1073 * - ximage is the XImage to draw on.
1074 * - x, y is the top left corner
1075 * - s is the length of each side
1079 draw_block(XImage * ximage, int x, int y, int s, unsigned long v)
1083 for (a = 0; a < s; a++)
1084 for (b = 0; b < s; b++) {
1085 XPutPixel(ximage, x + b, y + a, v);
1089 /****************************************************************/
1092 * draw_point Draw the current point in a swirl pattern onto the XImage
1094 * - swirl is the swirl
1095 * - win is the window to update
1098 draw_point(ModeInfo * mi, SWIRL_P swirl)
1103 /* get current point coordinates and resolution */
1108 /* check we are within the window */
1109 if ((x < 0) || (x > swirl->width - r) || (y < 0) || (y > swirl->height - r))
1112 /* what style are we drawing? */
1113 if (swirl->two_plane) {
1116 /* halve the block size */
1119 /* interleave blocks at half r */
1120 draw_block(swirl->ximage, x, y, r2, do_point(swirl, x, y));
1121 draw_block(swirl->ximage, x + r2, y, r2, do_point(swirl, x + r2, y));
1122 draw_block(swirl->ximage, x + r2, y + r2, r2, do_point(swirl,
1124 draw_block(swirl->ximage, x, y + r2, r2, do_point(swirl, x, y + r2));
1126 draw_block(swirl->ximage, x, y, r, do_point(swirl, x, y));
1128 /* update the screen */
1130 #ifdef HAVE_XSHM_EXTENSION
1132 XShmPutImage(MI_DISPLAY(mi), MI_WINDOW(mi), MI_GC(mi), swirl->ximage,
1133 x, y, x, y, r, r, False);
1135 #endif /* !HAVE_XSHM_EXTENSION */
1136 /* PURIFY 4.0.1 on SunOS4 and on Solaris 2 reports a 256 byte memory
1137 leak on the next line. */
1138 XPutImage(MI_DISPLAY(mi), MI_WINDOW(mi), MI_GC(mi), swirl->ximage,
1142 /****************************************************************/
1145 * next_point Move to the next point in the spiral pattern
1146 * - swirl is the swirl
1147 * - win is the window to update
1150 next_point(SWIRL_P swirl)
1152 /* more to do in this direction? */
1153 if (swirl->dir_done < swirl->dir_todo) {
1154 /* move in the current direction */
1155 switch (swirl->direction) {
1157 swirl->x += swirl->r;
1160 swirl->y += swirl->r;
1163 swirl->x -= swirl->r;
1166 swirl->y -= swirl->r;
1170 /* done another point */
1173 /* none drawn yet */
1174 swirl->dir_done = 0;
1176 /* change direction - check and record if off screen */
1177 switch (swirl->direction) {
1179 swirl->direction = DRAW_DOWN;
1180 if (swirl->x > swirl->width - swirl->r) {
1181 /* skip these points */
1182 swirl->dir_done = swirl->dir_todo;
1183 swirl->y += (swirl->dir_todo * swirl->r);
1185 /* check for finish */
1186 if (swirl->off_screen)
1187 swirl->drawing = False;
1188 swirl->off_screen = True;
1190 swirl->off_screen = False;
1193 swirl->direction = DRAW_LEFT;
1195 if (swirl->y > swirl->height - swirl->r) {
1196 /* skip these points */
1197 swirl->dir_done = swirl->dir_todo;
1198 swirl->x -= (swirl->dir_todo * swirl->r);
1200 /* check for finish */
1201 if (swirl->off_screen)
1202 swirl->drawing = False;
1203 swirl->off_screen = True;
1205 swirl->off_screen = False;
1208 swirl->direction = DRAW_UP;
1210 /* skip these points */
1211 swirl->dir_done = swirl->dir_todo;
1212 swirl->y -= (swirl->dir_todo * swirl->r);
1214 /* check for finish */
1215 if (swirl->off_screen)
1216 swirl->drawing = False;
1217 swirl->off_screen = True;
1219 swirl->off_screen = False;
1222 swirl->direction = DRAW_RIGHT;
1225 /* skip these points */
1226 swirl->dir_done = swirl->dir_todo;
1227 swirl->x += (swirl->dir_todo * swirl->r);
1229 /* check for finish */
1230 if (swirl->off_screen)
1231 swirl->drawing = False;
1232 swirl->off_screen = True;
1234 swirl->off_screen = False;
1240 /****************************************************************/
1245 * Initialise things for swirling
1247 * - win is the window to draw in
1250 init_swirl(ModeInfo * mi)
1252 Display *display = MI_DISPLAY(mi);
1253 Window window = MI_WINDOW(mi);
1256 /* does the swirls array exist? */
1257 if (swirls == NULL) {
1260 /* allocate an array, one entry for each screen */
1261 swirls = (SWIRL_P) calloc(ScreenCount(display), sizeof (SWIRL));
1263 /* initialise them all */
1264 for (i = 0; i < ScreenCount(display); i++)
1265 initialise_swirl(mi, &swirls[i]);
1267 /* get a pointer to this swirl */
1268 swirl = &(swirls[MI_SCREEN(mi)]);
1270 /* get window parameters */
1271 swirl->win = window;
1272 swirl->width = MI_WIN_WIDTH(mi);
1273 swirl->height = MI_WIN_HEIGHT(mi);
1274 swirl->depth = MI_WIN_DEPTH(mi);
1275 swirl->rdepth = swirl->depth;
1276 swirl->visual = MI_VISUAL(mi);
1278 if (swirl->depth > 16)
1281 /* initialise image for speeding up drawing */
1282 initialise_image(mi, swirl);
1284 /* clear the window (before setting the colourmap) */
1285 XClearWindow(display, MI_WINDOW(mi));
1289 swirl->rgb_values = mi->colors;
1290 swirl->colours = mi->npixels;
1291 swirl->dcolours = swirl->colours;
1292 /* swirl->fixed_colourmap = !mi->writable_p;*/
1294 #else /* !STANDALONE */
1296 /* initialise the colours from which the colourmap is derived */
1297 initialise_colours(basic_colours, MI_SATURATION(mi));
1299 /* set up the colour map */
1300 create_colourmap(mi, swirl);
1302 /* attach the colour map to the window (if we have one) */
1303 if (!swirl->fixed_colourmap) {
1305 setColormap(display, window, swirl->cmap, MI_WIN_IS_INWINDOW(mi));
1307 XSetWindowColormap(display, window, swirl->cmap);
1308 (void) XSetWMColormapWindows(display, window, &window, 1);
1309 XInstallColormap(display, swirl->cmap);
1312 #endif /* STANDALONE */
1314 /* resolution starts off chunky */
1315 swirl->resolution = MIN_RES + 1;
1317 /* calculate the pixel step for this resulution */
1318 swirl->r = (1 << (swirl->resolution - 1));
1320 /* how many knots? */
1321 swirl->n_knots = random_no((unsigned int) MI_BATCHCOUNT(mi) / 2) +
1322 MI_BATCHCOUNT(mi) + 1;
1324 /* what type of knots? */
1325 swirl->knot_type = ALL; /* for now */
1327 /* use two_plane mode occaisionally */
1328 if (random_no(100) <= TWO_PLANE_PCNT) {
1329 swirl->two_plane = swirl->first_plane = True;
1330 swirl->max_resolution = 2;
1332 swirl->two_plane = False;
1334 /* fix the knot values */
1335 create_knots(swirl);
1338 swirl->started = True;
1339 swirl->drawing = False;
1342 /****************************************************************/
1347 * Draw one iteration of swirling
1349 * - win is the window to draw in
1352 draw_swirl(ModeInfo * mi)
1354 SWIRL_P swirl = &(swirls[MI_SCREEN(mi)]);
1357 if (swirl->started) {
1358 /* in the middle of drawing? */
1359 if (swirl->drawing) {
1362 rotate_colors(MI_DISPLAY(mi), MI_COLORMAP(mi),
1363 swirl->rgb_values, swirl->colours, 1);
1364 #else /* !STANDALONE */
1365 /* rotate the colours */
1366 install_map(MI_DISPLAY(mi), swirl, swirl->dshift);
1367 #endif /* !STANDALONE */
1369 /* draw a batch of points */
1370 swirl->batch_todo = BATCH_DRAW;
1371 while ((swirl->batch_todo > 0) && swirl->drawing) {
1373 draw_point(mi, swirl);
1375 /* move to the next point */
1379 swirl->batch_todo--;
1384 rotate_colors(MI_DISPLAY(mi), MI_COLORMAP(mi),
1385 swirl->rgb_values, swirl->colours, 1);
1386 #else /* !STANDALONE */
1387 /* rotate the colours */
1388 install_map(MI_DISPLAY(mi), swirl, swirl->shift);
1389 #endif /* !STANDALONE */
1391 /* time for a higher resolution? */
1392 if (swirl->resolution > swirl->max_resolution) {
1393 /* move to higher resolution */
1394 swirl->resolution--;
1396 /* calculate the pixel step for this resulution */
1397 swirl->r = (1 << (swirl->resolution - 1));
1399 /* start drawing again */
1400 swirl->drawing = True;
1402 /* start in the middle of the screen */
1403 swirl->x = (swirl->width - swirl->r) / 2;
1404 swirl->y = (swirl->height - swirl->r) / 2;
1406 /* initialise spiral drawing parameters */
1407 swirl->direction = DRAW_RIGHT;
1408 swirl->dir_todo = 1;
1409 swirl->dir_done = 0;
1411 /* all done, decide when to restart */
1412 if (swirl->start_again == -1) {
1413 /* start the counter */
1414 swirl->start_again = RESTART;
1415 } else if (swirl->start_again == 0) {
1416 /* reset the counter */
1417 swirl->start_again = -1;
1420 /* Pick a new colormap! */
1421 XClearWindow (MI_DISPLAY(mi), MI_WINDOW(mi));
1422 free_colors (MI_DISPLAY(mi), MI_COLORMAP(mi),
1423 mi->colors, mi->npixels);
1424 make_smooth_colormap (MI_DISPLAY(mi),
1427 mi->colors, &mi->npixels, True,
1428 &mi->writable_p, True);
1429 swirl->colours = mi->npixels;
1430 #endif /* STANDALONE */
1435 /* decrement the counter */
1436 swirl->start_again--;
1442 /****************************************************************/
1445 release_swirl(ModeInfo * mi)
1447 /* does the swirls array exist? */
1448 if (swirls != NULL) {
1452 for (i = 0; i < MI_NUM_SCREENS(mi); i++) {
1453 SWIRL_P swirl = &(swirls[i]);
1455 if (swirl->cmap != (Colormap) NULL)
1456 XFreeColormap(MI_DISPLAY(mi), swirl->cmap);
1457 if (swirl->rgb_values != NULL)
1458 XFree((void *) swirl->rgb_values);
1459 if (swirl->ximage != NULL)
1460 XDestroyImage(swirl->ximage);
1462 (void) free((void *) swirl->knots);
1464 /* deallocate an array, one entry for each screen */
1465 (void) free((void *) swirls);
1470 /****************************************************************/
1473 refresh_swirl(ModeInfo * mi)
1475 SWIRL_P swirl = &(swirls[MI_SCREEN(mi)]);
1477 if (swirl->started) {
1479 swirl->resolution = swirl->resolution + 1;
1480 swirl->drawing = False;