From http://www.jwz.org/xscreensaver/xscreensaver-5.22.tar.gz
[xscreensaver] / hacks / hopalong.c
1 /* -*- Mode: C; tab-width: 4 -*- */
2 /* hop --- real plane fractals */
3
4 #if 0
5 static const char sccsid[] = "@(#)hop.c 5.00 2000/11/01 xlockmore";
6 #endif
7
8 /*-
9  * Copyright (c) 1991 by Patrick J. Naughton.
10  *
11  * Permission to use, copy, modify, and distribute this software and its
12  * documentation for any purpose and without fee is hereby granted,
13  * provided that the above copyright notice appear in all copies and that
14  * both that copyright notice and this permission notice appear in
15  * supporting documentation.
16  *
17  * This file is provided AS IS with no warranties of any kind.  The author
18  * shall have no liability with respect to the infringement of copyrights,
19  * trade secrets or any patents by this file or any part thereof.  In no
20  * event will the author be liable for any lost revenue or profits or
21  * other special, indirect and consequential damages.
22  *
23  * Revision History:
24  * Changes in xlockmore distribution
25  * 01-Nov-2000: Allocation checks
26  * 24-Jun-1997: EJK and RR functions stolen from xmartin2.2
27  *              Ed Kubaitis <ejk@ux2.cso.uiuc.edu> ejk functions and xmartin
28  *              Renaldo Recuerdo rr function, generalized exponent version
29  *              of the Barry Martin's square root function
30  * 10-May-1997: Compatible with xscreensaver
31  * 27-Jul-1995: added Peter de Jong's hop from Scientific American
32  *              July 87 p. 111.  Sometimes they are amazing but there are a
33  *              few duds (I did not see a pattern in the parameters).
34  * 29-Mar-1995: changed name from hopalong to hop
35  * 09-Dec-1994: added Barry Martin's sine hop
36  * Changes in original xlock
37  * 29-Oct-1990: fix bad (int) cast.
38  * 29-Jul-1990: support for multiple screens.
39  * 08-Jul-1990: new timing and colors and new algorithm for fractals.
40  * 15-Dec-1989: Fix for proper skipping of {White,Black}Pixel() in colors.
41  * 08-Oct-1989: Fixed long standing typo bug in RandomInitHop();
42  *                  Fixed bug in memory allocation in init_hop();
43  *                  Moved seconds() to an extern.
44  *                  Got rid of the % mod since .mod is slow on a sparc.
45  * 20-Sep-1989: Lint.
46  * 31-Aug-1988: Forked from xlock.c for modularity.
47  * 23-Mar-1988: Coded HOPALONG routines from Scientific American Sept. 86 p. 14.
48  *              Hopalong was attributed to Barry Martin of Aston University
49  *              (Birmingham, England)
50  */
51
52
53 #ifdef STANDALONE
54 #define MODE_hop
55 #define DEFAULTS        "*delay: 10000 \n" \
56                                         "*count: 1000 \n" \
57                                         "*cycles: 2500 \n" \
58                                         "*ncolors: 200 \n" \
59                                         "*fpsSolid: true \n" \
60                                         "*ignoreRotation: True \n" \
61
62 # define SMOOTH_COLORS
63 # define reshape_hop 0
64 # define hop_handle_event 0
65 # include "xlockmore.h"         /* in xscreensaver distribution */
66 # include "erase.h"
67 #else /* STANDALONE */
68 # include "xlock.h"             /* in xlockmore distribution */
69 #endif /* STANDALONE */
70
71 #ifdef MODE_hop
72
73 #define DEF_MARTIN "False"
74 #define DEF_POPCORN "False"
75 #define DEF_EJK1 "False"
76 #define DEF_EJK2 "False"
77 #define DEF_EJK3 "False"
78 #define DEF_EJK4 "False"
79 #define DEF_EJK5 "False"
80 #define DEF_EJK6 "False"
81 #define DEF_RR "False"
82 #define DEF_JONG "False"
83 #define DEF_SINE "False"
84
85 static Bool martin;
86 static Bool popcorn;
87 static Bool ejk1;
88 static Bool ejk2;
89 static Bool ejk3;
90 static Bool ejk4;
91 static Bool ejk5;
92 static Bool ejk6;
93 static Bool rr;
94 static Bool jong;
95 static Bool sine;
96
97 static XrmOptionDescRec opts[] =
98 {
99         {"-martin", ".hop.martin", XrmoptionNoArg, "on"},
100         {"+martin", ".hop.martin", XrmoptionNoArg, "off"},
101         {"-popcorn", ".hop.popcorn", XrmoptionNoArg, "on"},
102         {"+popcorn", ".hop.popcorn", XrmoptionNoArg, "off"},
103         {"-ejk1", ".hop.ejk1", XrmoptionNoArg, "on"},
104         {"+ejk1", ".hop.ejk1", XrmoptionNoArg, "off"},
105         {"-ejk2", ".hop.ejk2", XrmoptionNoArg, "on"},
106         {"+ejk2", ".hop.ejk2", XrmoptionNoArg, "off"},
107         {"-ejk3", ".hop.ejk3", XrmoptionNoArg, "on"},
108         {"+ejk3", ".hop.ejk3", XrmoptionNoArg, "off"},
109         {"-ejk4", ".hop.ejk4", XrmoptionNoArg, "on"},
110         {"+ejk4", ".hop.ejk4", XrmoptionNoArg, "off"},
111         {"-ejk5", ".hop.ejk5", XrmoptionNoArg, "on"},
112         {"+ejk5", ".hop.ejk5", XrmoptionNoArg, "off"},
113         {"-ejk6", ".hop.ejk6", XrmoptionNoArg, "on"},
114         {"+ejk6", ".hop.ejk6", XrmoptionNoArg, "off"},
115         {"-rr", ".hop.rr", XrmoptionNoArg, "on"},
116         {"+rr", ".hop.rr", XrmoptionNoArg, "off"},
117         {"-jong", ".hop.jong", XrmoptionNoArg, "on"},
118         {"+jong", ".hop.jong", XrmoptionNoArg, "off"},
119         {"-sine", ".hop.sine", XrmoptionNoArg, "on"},
120         {"+sine", ".hop.sine", XrmoptionNoArg, "off"}
121 };
122 static argtype vars[] =
123 {
124         {&martin,  "martin",  "Martin",  DEF_MARTIN,  t_Bool},
125         {&popcorn, "popcorn", "Popcorn", DEF_POPCORN, t_Bool},
126         {&ejk1, "ejk1", "EJK1", DEF_EJK1, t_Bool},
127         {&ejk2, "ejk2", "EJK2", DEF_EJK2, t_Bool},
128         {&ejk3, "ejk3", "EJK3", DEF_EJK3, t_Bool},
129         {&ejk4, "ejk4", "EJK4", DEF_EJK4, t_Bool},
130         {&ejk5, "ejk5", "EJK5", DEF_EJK5, t_Bool},
131         {&ejk6, "ejk6", "EJK6", DEF_EJK6, t_Bool},
132         {&rr,   "rr",   "RR",   DEF_RR,   t_Bool},
133         {&jong, "jong", "Jong", DEF_JONG, t_Bool},
134         {&sine, "sine", "Sine", DEF_SINE, t_Bool}
135 };
136 static OptionStruct desc[] =
137 {
138         {"-/+martin", "turn on/off sqrt format"},
139         {"-/+popcorn", "turn on/off Clifford A. Pickover's popcorn format"},
140         {"-/+ejk1", "turn on/off ejk1 format"},
141         {"-/+ejk2", "turn on/off ejk2 format"},
142         {"-/+ejk3", "turn on/off ejk3 format"},
143         {"-/+ejk4", "turn on/off ejk4 format"},
144         {"-/+ejk5", "turn on/off ejk5 format"},
145         {"-/+ejk6", "turn on/off ejk6 format"},
146         {"-/+rr",   "turn on/off rr format"},
147         {"-/+jong", "turn on/off jong format"},
148         {"-/+sine", "turn on/off sine format"}
149 };
150
151 ENTRYPOINT ModeSpecOpt hop_opts =
152 {sizeof opts / sizeof opts[0], opts, sizeof vars / sizeof vars[0], vars, desc};
153
154 #ifdef USE_MODULES
155 ModStruct   hop_description =
156 {"hop", "init_hop", "draw_hop", "release_hop",
157  "refresh_hop", "init_hop", (char *) NULL, &hop_opts,
158  10000, 1000, 2500, 1, 64, 1.0, "",
159  "Shows real plane iterated fractals", 0, NULL};
160
161 #endif
162
163 #define MARTIN 0
164 #define POPCORN 7
165 #define SINE 8
166 #define EJK1 1
167 #define EJK2 2
168 #define EJK3 9
169 #define EJK4 3
170 #define EJK5 4
171 #define EJK6 10
172 #define RR 5
173 #define JONG 6
174 #ifdef OFFENDING
175 #define OPS 8                   /* 8, 9, 10 might be too close to a swastika for some... */
176 #else
177 #define OPS 11
178 #endif
179
180 typedef struct {
181         int         centerx, centery;   /* center of the screen */
182         double      a, b, c, d;
183         double      i, j;       /* hopalong parameters */
184         int         inc;
185         int         pix;
186         int         op;
187         int         count;
188         int         bufsize;
189         XPoint     *pointBuffer;        /* pointer for XDrawPoints */
190 #ifdef STANDALONE
191   eraser_state *eraser;
192 #endif
193 } hopstruct;
194
195 static hopstruct *hops = (hopstruct *) NULL;
196
197 ENTRYPOINT void
198 init_hop(ModeInfo * mi)
199 {
200         Display    *display = MI_DISPLAY(mi);
201         GC          gc = MI_GC(mi);
202         double      range;
203         hopstruct  *hp;
204
205         if (hops == NULL) {
206                 if ((hops = (hopstruct *) calloc(MI_NUM_SCREENS(mi),
207                                                  sizeof (hopstruct))) == NULL)
208                         return;
209         }
210         hp = &hops[MI_SCREEN(mi)];
211
212         hp->centerx = MI_WIDTH(mi) / 2;
213         hp->centery = MI_HEIGHT(mi) / 2;
214         /* Make the other operations less common since they are less interesting */
215         if (MI_IS_FULLRANDOM(mi)) {
216                 hp->op = NRAND(OPS);
217         } else {
218                 if (martin)
219                         hp->op = MARTIN;
220                 else if (popcorn)
221                         hp->op = POPCORN;
222                 else if (ejk1)
223                         hp->op = EJK1;
224                 else if (ejk2)
225                         hp->op = EJK2;
226                 else if (ejk3)
227                         hp->op = EJK3;
228                 else if (ejk4)
229                         hp->op = EJK4;
230                 else if (ejk5)
231                         hp->op = EJK5;
232                 else if (ejk6)
233                         hp->op = EJK6;
234                 else if (rr)
235                         hp->op = RR;
236                 else if (jong)
237                         hp->op = JONG;
238                 else if (sine)
239                         hp->op = SINE;
240                 else
241                         hp->op = NRAND(OPS);
242         }
243
244         range = sqrt((double) hp->centerx * hp->centerx +
245              (double) hp->centery * hp->centery) / (1.0 + LRAND() / MAXRAND);
246         hp->i = hp->j = 0.0;
247         hp->inc = (int) ((LRAND() / MAXRAND) * 200) - 100;
248 #undef XMARTIN
249         switch (hp->op) {
250                 case MARTIN:
251 #ifdef XMARTIN
252                         hp->a = (LRAND() / MAXRAND) * 1500.0 + 40.0;
253                         hp->b = (LRAND() / MAXRAND) * 17.0 + 3.0;
254                         hp->c = (LRAND() / MAXRAND) * 3000.0 + 100.0;
255 #else
256                         hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 20.0;
257                         hp->b = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 20.0;
258                         if (LRAND() & 1)
259                                 hp->c = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 20.0;
260                         else
261                                 hp->c = 0.0;
262 #endif
263                         if (MI_IS_VERBOSE(mi))
264                                 (void) fprintf(stdout, "sqrt a=%g, b=%g, c=%g\n", hp->a, hp->b, hp->c);
265                         break;
266                 case EJK1:
267 #ifdef XMARTIN
268                         hp->a = (LRAND() / MAXRAND) * 500.0;
269                         hp->c = (LRAND() / MAXRAND) * 100.0 + 10.0;
270 #else
271                         hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 30.0;
272                         hp->c = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 40.0;
273 #endif
274                         hp->b = (LRAND() / MAXRAND) * 0.4;
275                         if (MI_IS_VERBOSE(mi))
276                                 (void) fprintf(stdout, "ejk1 a=%g, b=%g, c=%g\n", hp->a, hp->b, hp->c);
277                         break;
278                 case EJK2:
279 #ifdef XMARTIN
280                         hp->a = (LRAND() / MAXRAND) * 500.0;
281 #else
282                         hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 30.0;
283 #endif
284                         hp->b = pow(10.0, 6.0 + (LRAND() / MAXRAND) * 24.0);
285                         if (LRAND() & 1)
286                                 hp->b = -hp->b;
287                         hp->c = pow(10.0, (LRAND() / MAXRAND) * 9.0);
288                         if (LRAND() & 1)
289                                 hp->c = -hp->c;
290                         if (MI_IS_VERBOSE(mi))
291                                 (void) fprintf(stdout, "ejk2 a=%g, b=%g, c=%g\n", hp->a, hp->b, hp->c);
292                         break;
293                 case EJK3:
294 #ifdef XMARTIN
295                         hp->a = (LRAND() / MAXRAND) * 500.0;
296                         hp->c = (LRAND() / MAXRAND) * 80.0 + 30.0;
297 #else
298                         hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 30.0;
299                         hp->c = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 70.0;
300 #endif
301                         hp->b = (LRAND() / MAXRAND) * 0.35 + 0.5;
302                         if (MI_IS_VERBOSE(mi))
303                                 (void) fprintf(stdout, "ejk3 a=%g, b=%g, c=%g\n", hp->a, hp->b, hp->c);
304                         break;
305                 case EJK4:
306 #ifdef XMARTIN
307                         hp->a = (LRAND() / MAXRAND) * 1000.0;
308                         hp->c = (LRAND() / MAXRAND) * 40.0 + 30.0;
309 #else
310                         hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 2.0;
311                         hp->c = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 200.0;
312 #endif
313                         hp->b = (LRAND() / MAXRAND) * 9.0 + 1.0;
314                         if (MI_IS_VERBOSE(mi))
315                                 (void) fprintf(stdout, "ejk4 a=%g, b=%g, c=%g\n", hp->a, hp->b, hp->c);
316                         break;
317                 case EJK5:
318 #ifdef XMARTIN
319                         hp->a = (LRAND() / MAXRAND) * 600.0;
320                         hp->c = (LRAND() / MAXRAND) * 90.0 + 20.0;
321 #else
322                         hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 2.0;
323                         hp->c = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 200.0;
324 #endif
325                         hp->b = (LRAND() / MAXRAND) * 0.3 + 0.1;
326                         if (MI_IS_VERBOSE(mi))
327                                 (void) fprintf(stdout, "ejk5 a=%g, b=%g, c=%g\n", hp->a, hp->b, hp->c);
328                         break;
329                 case EJK6:
330 #ifdef XMARTIN
331                         hp->a = (LRAND() / MAXRAND) * 100.0 + 550.0;
332 #else
333                         hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 30.0;
334 #endif
335                         hp->b = (LRAND() / MAXRAND) + 0.5;
336                         if (MI_IS_VERBOSE(mi))
337                                 (void) fprintf(stdout, "ejk6 a=%g, b=%g\n", hp->a, hp->b);
338                         break;
339                 case RR:
340 #ifdef XMARTIN
341                         hp->a = (LRAND() / MAXRAND) * 100.0;
342                         hp->b = (LRAND() / MAXRAND) * 20.0;
343                         hp->c = (LRAND() / MAXRAND) * 200.0;
344 #else
345                         hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 40.0;
346                         hp->b = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 200.0;
347                         hp->c = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 20.0;
348 #endif
349                         hp->d = (LRAND() / MAXRAND) * 0.9;
350                         if (MI_IS_VERBOSE(mi))
351                                 (void) fprintf(stdout, "rr a=%g, b=%g, c=%g, d=%g\n",
352                                                hp->a, hp->b, hp->c, hp->d);
353                         break;
354                 case POPCORN:
355                         hp->a = 0.0;
356                         hp->b = 0.0;
357                         hp->c = ((LRAND() / MAXRAND) * 2.0 - 1.0) * 0.24 + 0.25;
358                         hp->inc = 100;
359                         if (MI_IS_VERBOSE(mi))
360                                 (void) fprintf(stdout, "popcorn a=%g, b=%g, c=%g, d=%g\n",
361                                                hp->a, hp->b, hp->c, hp->d);
362                         break;
363                 case JONG:
364                         hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * M_PI;
365                         hp->b = ((LRAND() / MAXRAND) * 2.0 - 1.0) * M_PI;
366                         hp->c = ((LRAND() / MAXRAND) * 2.0 - 1.0) * M_PI;
367                         hp->d = ((LRAND() / MAXRAND) * 2.0 - 1.0) * M_PI;
368                         if (MI_IS_VERBOSE(mi))
369                                 (void) fprintf(stdout, "jong a=%g, b=%g, c=%g, d=%g\n",
370                                                hp->a, hp->b, hp->c, hp->d);
371                         break;
372                 case SINE:      /* MARTIN2 */
373 #ifdef XMARTIN
374                         hp->a = M_PI + ((LRAND() / MAXRAND) * 2.0 - 1.0) * 0.07;
375 #else
376                         hp->a = M_PI + ((LRAND() / MAXRAND) * 2.0 - 1.0) * 0.7;
377 #endif
378                         if (MI_IS_VERBOSE(mi))
379                                 (void) fprintf(stdout, "sine a=%g\n", hp->a);
380                         break;
381         }
382         if (MI_NPIXELS(mi) > 2)
383                 hp->pix = NRAND(MI_NPIXELS(mi));
384         hp->bufsize = MI_COUNT(mi);
385
386         if (hp->pointBuffer == NULL) {
387                 if ((hp->pointBuffer = (XPoint *) malloc(hp->bufsize *
388                                 sizeof (XPoint))) == NULL)
389                         return;
390         }
391
392 #ifndef STANDALONE
393         MI_CLEARWINDOW(mi);
394 #endif
395
396         XSetForeground(display, gc, MI_WHITE_PIXEL(mi));
397         hp->count = 0;
398 }
399
400
401 ENTRYPOINT void
402 draw_hop(ModeInfo * mi)
403 {
404         double      oldj, oldi;
405         XPoint     *xp;
406         int         k;
407         hopstruct  *hp;
408
409         if (hops == NULL)
410                 return;
411         hp = &hops[MI_SCREEN(mi)];
412
413 #ifdef STANDALONE
414     if (hp->eraser) {
415       hp->eraser = erase_window (MI_DISPLAY(mi), MI_WINDOW(mi), hp->eraser);
416       return;
417     }
418 #endif
419
420
421         if (hp->pointBuffer == NULL)
422                 return;
423         xp = hp->pointBuffer;
424         k = hp->bufsize;
425
426         MI_IS_DRAWN(mi) = True;
427         hp->inc++;
428         if (MI_NPIXELS(mi) > 2) {
429                 XSetForeground(MI_DISPLAY(mi), MI_GC(mi), MI_PIXEL(mi, hp->pix));
430                 if (++hp->pix >= MI_NPIXELS(mi))
431                         hp->pix = 0;
432         }
433         while (k--) {
434                 oldj = hp->j;
435                 switch (hp->op) {
436                         case MARTIN:    /* SQRT, MARTIN1 */
437                                 oldi = hp->i + hp->inc;
438                                 hp->j = hp->a - hp->i;
439                                 hp->i = oldj + ((hp->i < 0)
440                                            ? sqrt(fabs(hp->b * oldi - hp->c))
441                                         : -sqrt(fabs(hp->b * oldi - hp->c)));
442                                 xp->x = hp->centerx + (int) (hp->i + hp->j);
443                                 xp->y = hp->centery - (int) (hp->i - hp->j);
444                                 break;
445                         case EJK1:
446                                 oldi = hp->i + hp->inc;
447                                 hp->j = hp->a - hp->i;
448                                 hp->i = oldj - ((hp->i > 0) ? (hp->b * oldi - hp->c) :
449                                                 -(hp->b * oldi - hp->c));
450                                 xp->x = hp->centerx + (int) (hp->i + hp->j);
451                                 xp->y = hp->centery - (int) (hp->i - hp->j);
452                                 break;
453                         case EJK2:
454                                 oldi = hp->i + hp->inc;
455                                 hp->j = hp->a - hp->i;
456                                 hp->i = oldj - ((hp->i < 0) ? log(fabs(hp->b * oldi - hp->c)) :
457                                            -log(fabs(hp->b * oldi - hp->c)));
458                                 xp->x = hp->centerx + (int) (hp->i + hp->j);
459                                 xp->y = hp->centery - (int) (hp->i - hp->j);
460                                 break;
461                         case EJK3:
462                                 oldi = hp->i + hp->inc;
463                                 hp->j = hp->a - hp->i;
464                                 hp->i = oldj - ((hp->i > 0) ? sin(hp->b * oldi) - hp->c :
465                                                 -sin(hp->b * oldi) - hp->c);
466                                 xp->x = hp->centerx + (int) (hp->i + hp->j);
467                                 xp->y = hp->centery - (int) (hp->i - hp->j);
468                                 break;
469                         case EJK4:
470                                 oldi = hp->i + hp->inc;
471                                 hp->j = hp->a - hp->i;
472                                 hp->i = oldj - ((hp->i > 0) ? sin(hp->b * oldi) - hp->c :
473                                           -sqrt(fabs(hp->b * oldi - hp->c)));
474                                 xp->x = hp->centerx + (int) (hp->i + hp->j);
475                                 xp->y = hp->centery - (int) (hp->i - hp->j);
476                                 break;
477                         case EJK5:
478                                 oldi = hp->i + hp->inc;
479                                 hp->j = hp->a - hp->i;
480                                 hp->i = oldj - ((hp->i > 0) ? sin(hp->b * oldi) - hp->c :
481                                                 -(hp->b * oldi - hp->c));
482                                 xp->x = hp->centerx + (int) (hp->i + hp->j);
483                                 xp->y = hp->centery - (int) (hp->i - hp->j);
484                                 break;
485                         case EJK6:
486                                 oldi = hp->i + hp->inc;
487                                 hp->j = hp->a - hp->i;
488                                 hp->i = oldj - asin((hp->b * oldi) - (long) (hp->b * oldi));
489                                 xp->x = hp->centerx + (int) (hp->i + hp->j);
490                                 xp->y = hp->centery - (int) (hp->i - hp->j);
491                                 break;
492                         case RR:        /* RR1 */
493                                 oldi = hp->i + hp->inc;
494                                 hp->j = hp->a - hp->i;
495                                 hp->i = oldj - ((hp->i < 0) ? -pow(fabs(hp->b * oldi - hp->c), hp->d) :
496                                      pow(fabs(hp->b * oldi - hp->c), hp->d));
497                                 xp->x = hp->centerx + (int) (hp->i + hp->j);
498                                 xp->y = hp->centery - (int) (hp->i - hp->j);
499                                 break;
500                         case POPCORN:
501 #define HVAL 0.05
502 #define INCVAL 50
503                                 {
504                                         double      tempi, tempj;
505
506                                         if (hp->inc >= 100)
507                                                 hp->inc = 0;
508                                         if (hp->inc == 0) {
509                                                 if (hp->a++ >= INCVAL) {
510                                                         hp->a = 0;
511                                                         if (hp->b++ >= INCVAL)
512                                                                 hp->b = 0;
513                                                 }
514                                                 hp->i = (-hp->c * INCVAL / 2 + hp->c * hp->a) * M_PI / 180.0;
515                                                 hp->j = (-hp->c * INCVAL / 2 + hp->c * hp->b) * M_PI / 180.0;
516                                         }
517                                         tempi = hp->i - HVAL * sin(hp->j + tan(3.0 * hp->j));
518                                         tempj = hp->j - HVAL * sin(hp->i + tan(3.0 * hp->i));
519                                         xp->x = hp->centerx + (int) (MI_WIDTH(mi) / 40 * tempi);
520                                         xp->y = hp->centery + (int) (MI_HEIGHT(mi) / 40 * tempj);
521                                         hp->i = tempi;
522                                         hp->j = tempj;
523                                 }
524                                 break;
525                         case JONG:
526                                 if (hp->centerx > 0)
527                                         oldi = hp->i + 4 * hp->inc / hp->centerx;
528                                 else
529                                         oldi = hp->i;
530                                 hp->j = sin(hp->c * hp->i) - cos(hp->d * hp->j);
531                                 hp->i = sin(hp->a * oldj) - cos(hp->b * oldi);
532                                 xp->x = hp->centerx + (int) (hp->centerx * (hp->i + hp->j) / 4.0);
533                                 xp->y = hp->centery - (int) (hp->centery * (hp->i - hp->j) / 4.0);
534                                 break;
535                         case SINE:      /* MARTIN2 */
536                                 oldi = hp->i + hp->inc;
537                                 hp->j = hp->a - hp->i;
538                                 hp->i = oldj - sin(oldi);
539                                 xp->x = hp->centerx + (int) (hp->i + hp->j);
540                                 xp->y = hp->centery - (int) (hp->i - hp->j);
541                                 break;
542                 }
543                 xp++;
544         }
545         XDrawPoints(MI_DISPLAY(mi), MI_WINDOW(mi), MI_GC(mi),
546                     hp->pointBuffer, hp->bufsize, CoordModeOrigin);
547         if (++hp->count > MI_CYCLES(mi)) {
548 #ifdef STANDALONE
549       hp->eraser = erase_window (MI_DISPLAY(mi), MI_WINDOW(mi), hp->eraser);
550 #endif /* STANDALONE */
551                 init_hop(mi);
552         }
553 }
554
555 ENTRYPOINT void
556 release_hop(ModeInfo * mi)
557 {
558         if (hops != NULL) {
559                 int         screen;
560
561                 for (screen = 0; screen < MI_NUM_SCREENS(mi); screen++) {
562                         hopstruct  *hp = &hops[screen];
563
564                         if (hp->pointBuffer != NULL)
565                                 (void) free((void *) hp->pointBuffer);
566                 }
567                 (void) free((void *) hops);
568                 hops = (hopstruct *) NULL;
569         }
570 }
571
572 ENTRYPOINT void
573 refresh_hop(ModeInfo * mi)
574 {
575         MI_CLEARWINDOW(mi);
576 }
577
578 XSCREENSAVER_MODULE_2 ("Hopalong", hopalong, hop)
579
580 #endif /* MODE_hop */