http://ftp.x.org/contrib/applications/xscreensaver-3.07.tar.gz
[xscreensaver] / hacks / hopalong.c
index c762a162f71b044f7a6c4d43538c7998cb0b05ec..7b252e87d77c0e2cb8e39830e3e952b68c1bf61c 100644 (file)
@@ -21,6 +21,10 @@ static const char sccsid[] = "@(#)hop.c      4.02 97/04/01 xlockmore";
  *
  * Revision History:
  * Changes of David Bagley <bagleyd@bigfoot.com>
+ * 24-Jun-97: EJK and RR functions stolen from xmartin2.2
+ *            Ed Kubaitis <ejk@ux2.cso.uiuc.edu> ejk functions and xmartin 
+ *            Renaldo Recuerdo rr function, generalized exponent version
+ *              of the Barry Martin's square root function
  * 10-May-97: jwz@jwz.org: ported from xlockmore 4.03a10 to be a 
  *                       standalone program and thus usable with xscreensaver
  *                       (I threw away my 1992 port and started over.)
@@ -28,7 +32,7 @@ static const char sccsid[] = "@(#)hop.c       4.02 97/04/01 xlockmore";
  *            July 87 p. 111.  Sometimes they are amazing but there are a
  *            few duds (I did not see a pattern in the parameters).
  * 29-Mar-95: changed name from hopalong to hop
- * 09-Dec-94: added sine hop
+ * 09-Dec-94: added Barry Martin's sine hop
  *
  * (12-Aug-92: jwz@lucid.com: made xlock version run standalone.)
  *
@@ -44,6 +48,8 @@ static const char sccsid[] = "@(#)hop.c       4.02 97/04/01 xlockmore";
  * 20-Sep-89: Lint.
  * 31-Aug-88: Forked from xlock.c for modularity.
  * 23-Mar-88: Coded HOPALONG routines from Scientific American Sept. 86 p. 14.
+ *            Hopalong was attributed to Barry Martin of Aston University
+ *            (Birmingham, England)
  */
 
 #ifdef STANDALONE
@@ -52,9 +58,9 @@ static const char sccsid[] = "@(#)hop.c       4.02 97/04/01 xlockmore";
 # define HACK_DRAW                                     draw_hop
 # define HACK_FREE                                     release_hop
 # define hop_opts                                      xlockmore_opts
-# define DEFAULTS      "*count:                1000    \n"                     \
+# define DEFAULTS      "*delay:                10000   \n"                     \
+                                       "*count:                1000    \n"                     \
                                        "*cycles:               2500    \n"                     \
-                                       "*delay:                10000   \n"                     \
                                        "*ncolors:              200     \n"
 # define SMOOTH_COLORS
 # include "xlockmore.h"                                /* from the xscreensaver distribution */
@@ -63,14 +69,50 @@ static const char sccsid[] = "@(#)hop.c     4.02 97/04/01 xlockmore";
 # include "xlock.h"                                    /* from the xlockmore distribution */
 #endif /* !STANDALONE */
 
-static Bool jong;
-static Bool sine;
-
+#define DEF_MARTIN "False"
+#define DEF_POPCORN "False"
+#define DEF_EJK1 "False"
+#define DEF_EJK2 "False"
+#define DEF_EJK3 "False"
+#define DEF_EJK4 "False"
+#define DEF_EJK5 "False"
+#define DEF_EJK6 "False"
+#define DEF_RR "False"
 #define DEF_JONG "False"
 #define DEF_SINE "False"
 
+static Bool martin;
+static Bool popcorn;
+static Bool ejk1;
+static Bool ejk2;
+static Bool ejk3;
+static Bool ejk4;
+static Bool ejk5;
+static Bool ejk6;
+static Bool rr;
+static Bool jong;
+static Bool sine;
+
 static XrmOptionDescRec opts[] =
 {
+       {"-martin", ".hop.martin", XrmoptionNoArg, (caddr_t) "on"},
+       {"+martin", ".hop.martin", XrmoptionNoArg, (caddr_t) "off"},
+       {"-popcorn", ".hop.popcorn", XrmoptionNoArg, (caddr_t) "on"},
+       {"+popcorn", ".hop.popcorn", XrmoptionNoArg, (caddr_t) "off"},
+       {"-ejk1", ".hop.ejk1", XrmoptionNoArg, (caddr_t) "on"},
+       {"+ejk1", ".hop.ejk1", XrmoptionNoArg, (caddr_t) "off"},
+       {"-ejk2", ".hop.ejk2", XrmoptionNoArg, (caddr_t) "on"},
+       {"+ejk2", ".hop.ejk2", XrmoptionNoArg, (caddr_t) "off"},
+       {"-ejk3", ".hop.ejk3", XrmoptionNoArg, (caddr_t) "on"},
+       {"+ejk3", ".hop.ejk3", XrmoptionNoArg, (caddr_t) "off"},
+       {"-ejk4", ".hop.ejk4", XrmoptionNoArg, (caddr_t) "on"},
+       {"+ejk4", ".hop.ejk4", XrmoptionNoArg, (caddr_t) "off"},
+       {"-ejk5", ".hop.ejk5", XrmoptionNoArg, (caddr_t) "on"},
+       {"+ejk5", ".hop.ejk5", XrmoptionNoArg, (caddr_t) "off"},
+       {"-ejk6", ".hop.ejk6", XrmoptionNoArg, (caddr_t) "on"},
+       {"+ejk6", ".hop.ejk6", XrmoptionNoArg, (caddr_t) "off"},
+       {"-rr", ".hop.rr", XrmoptionNoArg, (caddr_t) "on"},
+       {"+rr", ".hop.rr", XrmoptionNoArg, (caddr_t) "off"},
        {"-jong", ".hop.jong", XrmoptionNoArg, (caddr_t) "on"},
        {"+jong", ".hop.jong", XrmoptionNoArg, (caddr_t) "off"},
        {"-sine", ".hop.sine", XrmoptionNoArg, (caddr_t) "on"},
@@ -78,45 +120,75 @@ static XrmOptionDescRec opts[] =
 };
 static argtype vars[] =
 {
+       {(caddr_t *) & martin, "martin", "Martin", DEF_MARTIN, t_Bool},
+       {(caddr_t *) & popcorn, "popcorn", "Popcorn", DEF_POPCORN, t_Bool},
+       {(caddr_t *) & ejk1, "ejk1", "EJK1", DEF_EJK1, t_Bool},
+       {(caddr_t *) & ejk2, "ejk2", "EJK2", DEF_EJK2, t_Bool},
+       {(caddr_t *) & ejk3, "ejk3", "EJK3", DEF_EJK3, t_Bool},
+       {(caddr_t *) & ejk4, "ejk4", "EJK4", DEF_EJK4, t_Bool},
+       {(caddr_t *) & ejk5, "ejk5", "EJK5", DEF_EJK5, t_Bool},
+       {(caddr_t *) & ejk6, "ejk6", "EJK6", DEF_EJK6, t_Bool},
+       {(caddr_t *) & rr, "rr", "RR", DEF_RR, t_Bool},
        {(caddr_t *) & jong, "jong", "Jong", DEF_JONG, t_Bool},
        {(caddr_t *) & sine, "sine", "Sine", DEF_SINE, t_Bool}
 };
 static OptionStruct desc[] =
 {
+       {"-/+martin", "turn on/off sqrt format"},
+       {"-/+popcorn", "turn on/off Clifford A. Pickover's popcorn format"},
+       {"-/+ejk1", "turn on/off ejk1 format"},
+       {"-/+ejk2", "turn on/off ejk2 format"},
+       {"-/+ejk3", "turn on/off ejk3 format"},
+       {"-/+ejk4", "turn on/off ejk4 format"},
+       {"-/+ejk5", "turn on/off ejk5 format"},
+       {"-/+ejk6", "turn on/off ejk6 format"},
+       {"-/+rr", "turn on/off rr format"},
        {"-/+jong", "turn on/off jong format"},
        {"-/+sine", "turn on/off sine format"}
 };
 
-ModeSpecOpt hop_opts = { 4, opts, 2, vars, desc };
+ModeSpecOpt hop_opts =
+{sizeof opts / sizeof opts[0], opts, sizeof vars / sizeof vars[0], vars, desc};
 
+#ifdef USE_MODULES
+ModStruct   hop_description =
+{"hop", "init_hop", "draw_hop", "release_hop",
+ "refresh_hop", "init_hop", NULL, &hop_opts,
+ 10000, 1000, 2500, 1, 64, 1.0, "",
+ "Shows real plane iterated fractals", 0, NULL};
 
-#define SQRT 0
-#define JONG 1
-#define SINE 2
+#endif
+
+#define MARTIN 0
+#define POPCORN 7
+#define SINE 8
+#define EJK1 1
+#define EJK2 2
+#define EJK3 9
+#define EJK4 3
+#define EJK5 4
+#define EJK6 10
+#define RR 5
+#define JONG 6
 #ifdef OFFENDING
-#define OPS 2                  /* Sine might be too close to a Swastika for some... */
+#define OPS 8                  /* 8, 9, 10 might be too close to a swastika for some... */
 #else
-#define OPS 3
+#define OPS 11
 #endif
 
 typedef struct {
-       int         centerx;
-       int         centery;    /* center of the screen */
-       double      a;
-       double      b;
-       double      c;
-       double      d;
-       double      i;
-       double      j;          /* hopalong parameters */
+       int         centerx, centery;   /* center of the screen */
+       double      a, b, c, d;
+       double      i, j;       /* hopalong parameters */
        int         inc;
        int         pix;
        int         op;
        int         count;
        int         bufsize;
+       XPoint     *pointBuffer;        /* pointer for XDrawPoints */
 } hopstruct;
 
 static hopstruct *hops = NULL;
-static XPoint *pointBuffer = 0;        /* pointer for XDrawPoints */
 
 void
 init_hop(ModeInfo * mi)
@@ -133,53 +205,186 @@ init_hop(ModeInfo * mi)
        }
        hp = &hops[MI_SCREEN(mi)];
 
-       hp->centerx = MI_WIN_WIDTH(mi) / 2;
-       hp->centery = MI_WIN_HEIGHT(mi) / 2;
+       hp->centerx = MI_WIDTH(mi) / 2;
+       hp->centery = MI_HEIGHT(mi) / 2;
        /* Make the other operations less common since they are less interesting */
-       if (MI_WIN_IS_FULLRANDOM(mi)) {
-         hp->op = NRAND(OPS+2); /* jwz: make the others a bit more likely. */
-         if (hp->op >= OPS)
-               hp->op = SQRT;
+       if (MI_IS_FULLRANDOM(mi)) {
+               hp->op = NRAND(OPS);
        } else {
-               hp->op = SQRT;
-               if (jong)
+               if (martin)
+                       hp->op = MARTIN;
+               else if (popcorn)
+                       hp->op = POPCORN;
+               else if (ejk1)
+                       hp->op = EJK1;
+               else if (ejk2)
+                       hp->op = EJK2;
+               else if (ejk3)
+                       hp->op = EJK3;
+               else if (ejk4)
+                       hp->op = EJK4;
+               else if (ejk5)
+                       hp->op = EJK5;
+               else if (ejk6)
+                       hp->op = EJK6;
+               else if (rr)
+                       hp->op = RR;
+               else if (jong)
                        hp->op = JONG;
                else if (sine)
                        hp->op = SINE;
+               else
+                       hp->op = NRAND(OPS);
        }
-       switch (hp->op) {
-               case SQRT:
-                       range = sqrt((double) hp->centerx * hp->centerx +
-                                    (double) hp->centery * hp->centery) / (10.0 + NRAND(10));
 
-                       hp->a = (LRAND() / MAXRAND) * range - range / 2.0;
-                       hp->b = (LRAND() / MAXRAND) * range - range / 2.0;
-                       hp->c = (LRAND() / MAXRAND) * range - range / 2.0;
+       range = sqrt((double) hp->centerx * hp->centerx +
+            (double) hp->centery * hp->centery) / (1.0 + LRAND() / MAXRAND);
+       hp->i = hp->j = 0.0;
+       hp->inc = (int) ((LRAND() / MAXRAND) * 200) - 100;
+#undef XMARTIN
+       switch (hp->op) {
+               case MARTIN:
+#ifdef XMARTIN
+                       hp->a = (LRAND() / MAXRAND) * 1500.0 + 40.0;
+                       hp->b = (LRAND() / MAXRAND) * 17.0 + 3.0;
+                       hp->c = (LRAND() / MAXRAND) * 3000.0 + 100.0;
+#else
+                       hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 20.0;
+                       hp->b = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 20.0;
                        if (LRAND() & 1)
+                               hp->c = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 20.0;
+                       else
                                hp->c = 0.0;
+#endif
+                       if (MI_IS_VERBOSE(mi))
+                               (void) fprintf(stdout, "sqrt a=%g, b=%g, c=%g\n", hp->a, hp->b, hp->c);
+                       break;
+               case EJK1:
+#ifdef XMARTIN
+                       hp->a = (LRAND() / MAXRAND) * 500.0;
+                       hp->c = (LRAND() / MAXRAND) * 100.0 + 10.0;
+#else
+                       hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 30.0;
+                       hp->c = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 40.0;
+#endif
+                       hp->b = (LRAND() / MAXRAND) * 0.4;
+                       if (MI_IS_VERBOSE(mi))
+                               (void) fprintf(stdout, "ejk1 a=%g, b=%g, c=%g\n", hp->a, hp->b, hp->c);
+                       break;
+               case EJK2:
+#ifdef XMARTIN
+                       hp->a = (LRAND() / MAXRAND) * 500.0;
+#else
+                       hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 30.0;
+#endif
+                       hp->b = pow(10.0, 6.0 + (LRAND() / MAXRAND) * 24.0);
+                       if (LRAND() & 1)
+                               hp->b = -hp->b;
+                       hp->c = pow(10.0, (LRAND() / MAXRAND) * 9.0);
+                       if (LRAND() & 1)
+                               hp->c = -hp->c;
+                       if (MI_IS_VERBOSE(mi))
+                               (void) fprintf(stdout, "ejk2 a=%g, b=%g, c=%g\n", hp->a, hp->b, hp->c);
+                       break;
+               case EJK3:
+#ifdef XMARTIN
+                       hp->a = (LRAND() / MAXRAND) * 500.0;
+                       hp->c = (LRAND() / MAXRAND) * 80.0 + 30.0;
+#else
+                       hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 30.0;
+                       hp->c = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 70.0;
+#endif
+                       hp->b = (LRAND() / MAXRAND) * 0.35 + 0.5;
+                       if (MI_IS_VERBOSE(mi))
+                               (void) fprintf(stdout, "ejk3 a=%g, b=%g, c=%g\n", hp->a, hp->b, hp->c);
+                       break;
+               case EJK4:
+#ifdef XMARTIN
+                       hp->a = (LRAND() / MAXRAND) * 1000.0;
+                       hp->c = (LRAND() / MAXRAND) * 40.0 + 30.0;
+#else
+                       hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 2.0;
+                       hp->c = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 200.0;
+#endif
+                       hp->b = (LRAND() / MAXRAND) * 9.0 + 1.0;
+                       if (MI_IS_VERBOSE(mi))
+                               (void) fprintf(stdout, "ejk4 a=%g, b=%g, c=%g\n", hp->a, hp->b, hp->c);
+                       break;
+               case EJK5:
+#ifdef XMARTIN
+                       hp->a = (LRAND() / MAXRAND) * 600.0;
+                       hp->c = (LRAND() / MAXRAND) * 90.0 + 20.0;
+#else
+                       hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 2.0;
+                       hp->c = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 200.0;
+#endif
+                       hp->b = (LRAND() / MAXRAND) * 0.3 + 0.1;
+                       if (MI_IS_VERBOSE(mi))
+                               (void) fprintf(stdout, "ejk5 a=%g, b=%g, c=%g\n", hp->a, hp->b, hp->c);
+                       break;
+               case EJK6:
+#ifdef XMARTIN
+                       hp->a = (LRAND() / MAXRAND) * 100.0 + 550.0;
+#else
+                       hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 30.0;
+#endif
+                       hp->b = (LRAND() / MAXRAND) + 0.5;
+                       if (MI_IS_VERBOSE(mi))
+                               (void) fprintf(stdout, "ejk6 a=%g, b=%g\n", hp->a, hp->b);
+                       break;
+               case RR:
+#ifdef XMARTIN
+                       hp->a = (LRAND() / MAXRAND) * 100.0;
+                       hp->b = (LRAND() / MAXRAND) * 20.0;
+                       hp->c = (LRAND() / MAXRAND) * 200.0;
+#else
+                       hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 40.0;
+                       hp->b = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 200.0;
+                       hp->c = ((LRAND() / MAXRAND) * 2.0 - 1.0) * range / 20.0;
+#endif
+                       hp->d = (LRAND() / MAXRAND) * 0.9;
+                       if (MI_IS_VERBOSE(mi))
+                               (void) fprintf(stdout, "rr a=%g, b=%g, c=%g, d=%g\n",
+                                              hp->a, hp->b, hp->c, hp->d);
+                       break;
+               case POPCORN:
+                       hp->a = 0.0;
+                       hp->b = 0.0;
+                       hp->c = ((LRAND() / MAXRAND) * 2.0 - 1.0) * 0.24 + 0.25;
+                       hp->inc = 100;
+                       if (MI_IS_VERBOSE(mi))
+                               (void) fprintf(stdout, "popcorn a=%g, b=%g, c=%g, d=%g\n",
+                                              hp->a, hp->b, hp->c, hp->d);
                        break;
                case JONG:
-                       hp->a = (LRAND() / MAXRAND) * 2.0 * M_PI - M_PI;
-                       hp->b = (LRAND() / MAXRAND) * 2.0 * M_PI - M_PI;
-                       hp->c = (LRAND() / MAXRAND) * 2.0 * M_PI - M_PI;
-                       hp->d = (LRAND() / MAXRAND) * 2.0 * M_PI - M_PI;
+                       hp->a = ((LRAND() / MAXRAND) * 2.0 - 1.0) * M_PI;
+                       hp->b = ((LRAND() / MAXRAND) * 2.0 - 1.0) * M_PI;
+                       hp->c = ((LRAND() / MAXRAND) * 2.0 - 1.0) * M_PI;
+                       hp->d = ((LRAND() / MAXRAND) * 2.0 - 1.0) * M_PI;
+                       if (MI_IS_VERBOSE(mi))
+                               (void) fprintf(stdout, "jong a=%g, b=%g, c=%g, d=%g\n",
+                                              hp->a, hp->b, hp->c, hp->d);
                        break;
-               case SINE:
+               case SINE:      /* MARTIN2 */
+#ifdef XMARTIN
+                       hp->a = M_PI + ((LRAND() / MAXRAND) * 2.0 - 1.0) * 0.07;
+#else
                        hp->a = M_PI + ((LRAND() / MAXRAND) * 2.0 - 1.0) * 0.7;
+#endif
+                       if (MI_IS_VERBOSE(mi))
+                               (void) fprintf(stdout, "sine a=%g\n", hp->a);
                        break;
        }
        if (MI_NPIXELS(mi) > 2)
                hp->pix = NRAND(MI_NPIXELS(mi));
-       hp->i = hp->j = 0.0;
-       hp->inc = (int) ((LRAND() / MAXRAND) * 200) - 100;
-       hp->bufsize = MI_BATCHCOUNT(mi);
+       hp->bufsize = MI_COUNT(mi);
 
-       if (!pointBuffer)
-               pointBuffer = (XPoint *) malloc(hp->bufsize * sizeof (XPoint));
+       if (!hp->pointBuffer)
+               hp->pointBuffer = (XPoint *) malloc(hp->bufsize * sizeof (XPoint));
 
-       XClearWindow(display, MI_WINDOW(mi));
+       MI_CLEARWINDOW(mi);
 
-       XSetForeground(display, gc, MI_WIN_WHITE_PIXEL(mi));
+       XSetForeground(display, gc, MI_WHITE_PIXEL(mi));
        hp->count = 0;
 }
 
@@ -189,9 +394,11 @@ draw_hop(ModeInfo * mi)
 {
        hopstruct  *hp = &hops[MI_SCREEN(mi)];
        double      oldj, oldi;
-       XPoint     *xp = pointBuffer;
+       XPoint     *xp = hp->pointBuffer;
        int         k = hp->bufsize;
 
+       MI_IS_DRAWN(mi) = True;
+
        hp->inc++;
        if (MI_NPIXELS(mi) > 2) {
                XSetForeground(MI_DISPLAY(mi), MI_GC(mi), MI_PIXEL(mi, hp->pix));
@@ -201,7 +408,7 @@ draw_hop(ModeInfo * mi)
        while (k--) {
                oldj = hp->j;
                switch (hp->op) {
-                       case SQRT:
+                       case MARTIN:    /* SQRT, MARTIN1 */
                                oldi = hp->i + hp->inc;
                                hp->j = hp->a - hp->i;
                                hp->i = oldj + ((hp->i < 0)
@@ -210,6 +417,86 @@ draw_hop(ModeInfo * mi)
                                xp->x = hp->centerx + (int) (hp->i + hp->j);
                                xp->y = hp->centery - (int) (hp->i - hp->j);
                                break;
+                       case EJK1:
+                               oldi = hp->i + hp->inc;
+                               hp->j = hp->a - hp->i;
+                               hp->i = oldj - ((hp->i > 0) ? (hp->b * oldi - hp->c) :
+                                               -(hp->b * oldi - hp->c));
+                               xp->x = hp->centerx + (int) (hp->i + hp->j);
+                               xp->y = hp->centery - (int) (hp->i - hp->j);
+                               break;
+                       case EJK2:
+                               oldi = hp->i + hp->inc;
+                               hp->j = hp->a - hp->i;
+                               hp->i = oldj - ((hp->i < 0) ? log(fabs(hp->b * oldi - hp->c)) :
+                                          -log(fabs(hp->b * oldi - hp->c)));
+                               xp->x = hp->centerx + (int) (hp->i + hp->j);
+                               xp->y = hp->centery - (int) (hp->i - hp->j);
+                               break;
+                       case EJK3:
+                               oldi = hp->i + hp->inc;
+                               hp->j = hp->a - hp->i;
+                               hp->i = oldj - ((hp->i > 0) ? sin(hp->b * oldi) - hp->c :
+                                               -sin(hp->b * oldi) - hp->c);
+                               xp->x = hp->centerx + (int) (hp->i + hp->j);
+                               xp->y = hp->centery - (int) (hp->i - hp->j);
+                               break;
+                       case EJK4:
+                               oldi = hp->i + hp->inc;
+                               hp->j = hp->a - hp->i;
+                               hp->i = oldj - ((hp->i > 0) ? sin(hp->b * oldi) - hp->c :
+                                         -sqrt(fabs(hp->b * oldi - hp->c)));
+                               xp->x = hp->centerx + (int) (hp->i + hp->j);
+                               xp->y = hp->centery - (int) (hp->i - hp->j);
+                               break;
+                       case EJK5:
+                               oldi = hp->i + hp->inc;
+                               hp->j = hp->a - hp->i;
+                               hp->i = oldj - ((hp->i > 0) ? sin(hp->b * oldi) - hp->c :
+                                               -(hp->b * oldi - hp->c));
+                               xp->x = hp->centerx + (int) (hp->i + hp->j);
+                               xp->y = hp->centery - (int) (hp->i - hp->j);
+                               break;
+                       case EJK6:
+                               oldi = hp->i + hp->inc;
+                               hp->j = hp->a - hp->i;
+                               hp->i = oldj - asin((hp->b * oldi) - (long) (hp->b * oldi));
+                               xp->x = hp->centerx + (int) (hp->i + hp->j);
+                               xp->y = hp->centery - (int) (hp->i - hp->j);
+                               break;
+                       case RR:        /* RR1 */
+                               oldi = hp->i + hp->inc;
+                               hp->j = hp->a - hp->i;
+                               hp->i = oldj - ((hp->i < 0) ? -pow(fabs(hp->b * oldi - hp->c), hp->d) :
+                                    pow(fabs(hp->b * oldi - hp->c), hp->d));
+                               xp->x = hp->centerx + (int) (hp->i + hp->j);
+                               xp->y = hp->centery - (int) (hp->i - hp->j);
+                               break;
+                       case POPCORN:
+#define HVAL 0.05
+#define INCVAL 50
+                               {
+                                       double      tempi, tempj;
+
+                                       if (hp->inc >= 100)
+                                               hp->inc = 0;
+                                       if (hp->inc == 0) {
+                                               if (hp->a++ >= INCVAL) {
+                                                       hp->a = 0;
+                                                       if (hp->b++ >= INCVAL)
+                                                               hp->b = 0;
+                                               }
+                                               hp->i = (-hp->c * INCVAL / 2 + hp->c * hp->a) * M_PI / 180.0;
+                                               hp->j = (-hp->c * INCVAL / 2 + hp->c * hp->b) * M_PI / 180.0;
+                                       }
+                                       tempi = hp->i - HVAL * sin(hp->j + tan(3.0 * hp->j));
+                                       tempj = hp->j - HVAL * sin(hp->i + tan(3.0 * hp->i));
+                                       xp->x = hp->centerx + (int) (MI_WIDTH(mi) / 40 * tempi);
+                                       xp->y = hp->centery + (int) (MI_HEIGHT(mi) / 40 * tempj);
+                                       hp->i = tempi;
+                                       hp->j = tempj;
+                               }
+                               break;
                        case JONG:
                                if (hp->centerx > 0)
                                        oldi = hp->i + 4 * hp->inc / hp->centerx;
@@ -220,7 +507,7 @@ draw_hop(ModeInfo * mi)
                                xp->x = hp->centerx + (int) (hp->centerx * (hp->i + hp->j) / 4.0);
                                xp->y = hp->centery - (int) (hp->centery * (hp->i - hp->j) / 4.0);
                                break;
-                       case SINE:
+                       case SINE:      /* MARTIN2 */
                                oldi = hp->i + hp->inc;
                                hp->j = hp->a - hp->i;
                                hp->i = oldj - sin(oldi);
@@ -231,12 +518,12 @@ draw_hop(ModeInfo * mi)
                xp++;
        }
        XDrawPoints(MI_DISPLAY(mi), MI_WINDOW(mi), MI_GC(mi),
-                   pointBuffer, hp->bufsize, CoordModeOrigin);
+                   hp->pointBuffer, hp->bufsize, CoordModeOrigin);
        if (++hp->count > MI_CYCLES(mi)) {
 #ifdef STANDALONE
-         erase_full_window(MI_DISPLAY(mi), MI_WINDOW(mi));
+           erase_full_window(MI_DISPLAY(mi), MI_WINDOW(mi));
 #endif /* STANDALONE */
-         init_hop(mi);
+               init_hop(mi);
        }
 }
 
@@ -244,17 +531,22 @@ void
 release_hop(ModeInfo * mi)
 {
        if (hops != NULL) {
+               int         screen;
+
+               for (screen = 0; screen < MI_NUM_SCREENS(mi); screen++) {
+                       hopstruct  *hp = &hops[screen];
+
+                       if (hp->pointBuffer)
+                               (void) free((void *) hp->pointBuffer);
+               }
+
                (void) free((void *) hops);
                hops = NULL;
        }
-       if (pointBuffer) {
-               (void) free((void *) pointBuffer);
-               pointBuffer = NULL;
-       }
 }
 
 void
 refresh_hop(ModeInfo * mi)
 {
-       XClearWindow(MI_DISPLAY(mi), MI_WINDOW(mi));
+       MI_CLEARWINDOW(mi);
 }