diff options
Diffstat (limited to 'hacks/glx/hilbert.c')
| -rw-r--r-- | hacks/glx/hilbert.c | 1164 |
1 files changed, 0 insertions, 1164 deletions
diff --git a/hacks/glx/hilbert.c b/hacks/glx/hilbert.c deleted file mode 100644 index 5585f63..0000000 --- a/hacks/glx/hilbert.c +++ /dev/null @@ -1,1164 +0,0 @@ -/* hilbert, Copyright (c) 2011-2014 Jamie Zawinski <jwz@jwz.org> - * - * Permission to use, copy, modify, distribute, and sell this software and its - * documentation for any purpose is hereby granted without fee, provided that - * the above copyright notice appear in all copies and that both that - * copyright notice and this permission notice appear in supporting - * documentation. No representations are made about the suitability of this - * software for any purpose. It is provided "as is" without express or - * implied warranty. - * - * 2D and 3D Hilbert space-filling curves. - * - * Inspired by "Visualizing Hilbert Curves" by Nelson Max, 1998: - * https://e-reports-ext.llnl.gov/pdf/234149.pdf - */ - -#define DEFAULTS "*delay: 30000 \n" \ - "*count: 30 \n" \ - "*showFPS: False \n" \ - "*wireframe: False \n" \ - "*geometry: 800x800\n" \ - "*suppressRotationAnimation: True\n" \ - -# define release_hilbert 0 - -#include "xlockmore.h" -#include "colors.h" -#include "sphere.h" -#include "tube.h" -#include "rotator.h" -#include "gltrackball.h" -#include <ctype.h> - -#ifdef USE_GL /* whole file */ - - -#define DEF_SPIN "True" -#define DEF_WANDER "False" -#define DEF_SPEED "1.0" -#define DEF_MODE "random" -#define DEF_ENDS "random" -#define DEF_MAX_DEPTH "5" -#define DEF_THICKNESS "0.25" - -#define PAUSE_TICKS 180 - -typedef struct { - double x,y,z; -} XYZ; - -typedef struct { - unsigned short x,y,z; -} XYZb; - -typedef struct { - int size; - XYZb *values; /* assume max depth of 20 (2^16 on each side) */ -} hilbert_cache; - - -static int dlist_faces[] = { 20, 10, 8, 4, 3 }; - - -typedef struct { - GLXContext *glx_context; - rotator *rot, *rot2; - trackball_state *trackball; - Bool button_down_p; - Bool twodee_p; - Bool closed_p; - int ncolors; - XColor *colors; - - int depth; - int depth_tick; - - GLfloat path_start, path_end; - int path_tick; - int pause; - GLfloat diam; - - hilbert_cache **caches; /* filled lazily */ - - GLuint dlists [40][2]; /* we don't actually alloc all of these */ - int dlist_polys [40][2]; - -} hilbert_configuration; - -static hilbert_configuration *bps = NULL; - -static Bool do_spin; -static GLfloat speed; -static Bool do_wander; -static char *mode_str; -static char *ends_str; -static int max_depth; -static GLfloat thickness; - -static XrmOptionDescRec opts[] = { - { "-spin", ".spin", XrmoptionNoArg, "True" }, - { "+spin", ".spin", XrmoptionNoArg, "False" }, - { "-speed", ".speed", XrmoptionSepArg, 0 }, - { "-wander", ".wander", XrmoptionNoArg, "True" }, - { "+wander", ".wander", XrmoptionNoArg, "False" }, - { "-mode", ".mode", XrmoptionSepArg, 0 }, - { "-2d", ".mode", XrmoptionNoArg, "2D" }, - { "-3d", ".mode", XrmoptionNoArg, "3D" }, - { "-ends", ".ends", XrmoptionSepArg, 0 }, - { "-closed", ".ends", XrmoptionNoArg, "closed" }, - { "-open", ".ends", XrmoptionNoArg, "open" }, - { "-max-depth", ".maxDepth", XrmoptionSepArg, 0 }, - { "-thickness", ".thickness",XrmoptionSepArg, 0 }, -}; - -static argtype vars[] = { - {&do_spin, "spin", "Spin", DEF_SPIN, t_Bool}, - {&do_wander, "wander", "Wander", DEF_WANDER, t_Bool}, - {&speed, "speed", "Speed", DEF_SPEED, t_Float}, - {&mode_str, "mode", "Mode", DEF_MODE, t_String}, - {&ends_str, "ends", "Ends", DEF_ENDS, t_String}, - {&max_depth, "maxDepth", "MaxDepth", DEF_MAX_DEPTH, t_Int}, - {&thickness, "thickness","Thickness",DEF_THICKNESS, t_Float}, -}; - -ENTRYPOINT ModeSpecOpt hilbert_opts = {countof(opts), opts, countof(vars), vars, NULL}; - - -/* 2D Hilbert, and closed-loop variant. - */ - -/* These arrays specify which bits of the T index contribute to the - X, Y and Z values. Higher order bits are defined recursively. - */ -static const int xbit2d[4] = { 0, 0, 1, 1 }; -static const int ybit2d[4] = { 0, 1, 1, 0 }; - -static const int xbit3d[8] = { 0,0,0,0,1,1,1,1 }; -static const int ybit3d[8] = { 0,1,1,0,0,1,1,0 }; -static const int zbit3d[8] = { 0,0,1,1,1,1,0,0 }; - -/* These matrixes encapsulate the necessary reflection and translation - of each 4 sub-squares or 8 sub-cubes. The r2d and r3d arrays are - the normal Hilbert descent, and the s2d and s3d arrays are the - modified curve used for only level 0 of a closed-form path. - */ - -static int r2d[4][2][2] = { - /* _ _ - | |..| | - :_ _: - _| |_ - - */ - {{ 0, 1}, - { 1, 0}}, - {{ 1, 0}, - { 0, 1}}, - {{ 1, 0}, - { 0, 1}}, - {{ 0,-1}, - {-1, 0}}, -}; - -static int s2d[4][2][2] = { - /* __ __ - | |..| | Used for outermost iteration only, in "closed" mode - : .. : - |__| |__| - - */ - {{-1, 0}, - { 0,-1}}, - {{ 1, 0}, - { 0, 1}}, - {{ 1, 0}, - { 0, 1}}, - {{-1, 0}, - { 0,-1}}, -}; - - -static int r3d[8][3][3] = { - /* - /| /| - / | / | - /__|____/ | - | | - / / - / / - */ - {{ 0, 1, 0}, - { 1, 0, 0}, - { 0, 0, 1}}, - {{ 0, 0, 1}, - { 0, 1, 0}, - { 1, 0, 0}}, - {{ 1, 0, 0}, - { 0, 1, 0}, - { 0, 0, 1}}, - {{ 0, 0,-1}, - {-1, 0, 0}, - { 0, 1, 0}}, - {{ 0, 0, 1}, - { 1, 0, 0}, - { 0, 1, 0}}, - {{ 1, 0, 0}, - { 0, 1, 0}, - { 0, 0, 1}}, - {{ 0, 0,-1}, - { 0, 1, 0}, - {-1, 0, 0}}, - {{ 0,-1, 0}, - {-1, 0, 0}, - { 0, 0, 1}}, -}; - - -static int s3d[8][3][3] = { - /* - /| /| Used for outermost iteration only, in "closed" mode - / | / | - /__|____/ | - | | - / / - /_______/ - */ - {{-1, 0, 0}, - { 0, 0,-1}, - { 0, 1, 0}}, - {{ 0, 0, 1}, - { 0, 1, 0}, - { 1, 0, 0}}, - {{ 1, 0, 0}, - { 0, 1, 0}, - { 0, 0, 1}}, - {{ 0, 0,-1}, - {-1, 0, 0}, - { 0, 1, 0}}, - {{ 0, 0, 1}, - { 1, 0, 0}, - { 0, 1, 0}}, - {{ 1, 0, 0}, - { 0, 1, 0}, - { 0, 0, 1}}, - {{ 0, 0,-1}, - { 0, 1, 0}, - {-1, 0, 0}}, - - {{-1, 0, 0}, - { 0, 0,-1}, - { 0, 1, 0}}, -}; - - - -/* Matrix utilities - */ - -static void -matrix_times_vector2d (int m[2][2], int v[2], int dest[2]) -{ - dest[0] = m[0][0]*v[0] + m[0][1]*v[1]; - dest[1] = m[1][0]*v[0] + m[1][1]*v[1]; -} - -static void -matrix_times_vector3d (int m[3][3], int v[3], int dest[3]) -{ - dest[0] = m[0][0]*v[0] + m[0][1]*v[1] + m[0][2]*v[2]; - dest[1] = m[1][0]*v[0] + m[1][1]*v[1] + m[1][2]*v[2]; - dest[2] = m[2][0]*v[0] + m[2][1]*v[1] + m[2][2]*v[2]; -} - - -static void -matrix_multiply2d (int m1[2][2], int m2[2][2], int dest[2][2]) -{ - dest[0][0] = m1[0][0] * m2[0][0] + m1[0][1] * m2[1][0]; - dest[0][1] = m1[0][0] * m2[0][1] + m1[0][1] * m2[1][1]; - dest[1][0] = m1[1][0] * m2[0][0] + m1[1][1] * m2[1][0]; - dest[1][1] = m1[1][0] * m2[0][1] + m1[1][1] * m2[1][1]; -} - -static void -matrix_multiply3d (int m1[3][3], int m2[3][3], int dest[3][3]) -{ - int i, j, k; - for (i = 0; i < 3; i++) - for (j = 0; j < 3; j++) - { - dest[i][j] = 0; - for (k = 0; k < 3; k++) - dest[i][j] += m1[i][k] * m2[k][j]; - } -} - -static void -identity_matrix2d (int m[2][2]) -{ - m[0][0] = m[1][1] = 1; - m[0][1] = m[1][0] = 0; -} - -static void -identity_matrix3d (int m[3][3]) -{ - m[0][0] = m[1][1] = m[2][2] = 1; - m[0][1] = m[0][2] = m[1][0] = m[1][2] = m[2][0] = m[2][1] = 0; -} - - -static void -matrix_copy2d (int src[2][2], int dest[2][2]) -{ - int i, j; - for (i = 0; i < 2; i++) - for (j = 0; j < 2; j++) - dest[i][j] = src[i][j]; -} - - -static void -matrix_copy3d (int src[3][3], int dest[3][3]) -{ - int i, j; - for (i = 0; i < 3; i++) - for (j = 0; j < 3; j++) - dest[i][j] = src[i][j]; -} - - -/* 2D and 3D Hilbert: - Given an index T along the curve, return the XY or XYZ coordinates. - N is depth of the curve. - */ - -static void -t_to_xy (int n, int t, int *x, int *y, Bool closed_p) -{ - int k, rt[2][2], rq[2][2], va[2], vb[2]; - identity_matrix2d(rt); - *x = *y = 0; - - for (k = n-1; k >= 0; k--) - { - int j = 3 & (t >> (2*k)); - va[0] = 2 * xbit2d[j] - 1; - va[1] = 2 * ybit2d[j] - 1; - matrix_times_vector2d (rt, va, vb); - *x += ((vb[0] + 1) / 2) << k; - *y += ((vb[1] + 1) / 2) << k; - if (k > 0) - { - matrix_copy2d (rt, rq); - if (k == n-1 && closed_p) - matrix_multiply2d (rq, s2d[j], rt); - else - matrix_multiply2d (rq, r2d[j], rt); - } - } -} - - -static void -t_to_xyz (int n, int t, int *x, int *y, int *z, Bool closed_p) -{ - int k, rt[3][3], rq[3][3], va[3], vb[3]; - identity_matrix3d(rt); - *x = *y = *z = 0; - - for (k = n-1; k >= 0; k--) - { - int j = 7 & (t >> (3*k)); - va[0] = 2 * xbit3d[j] - 1; - va[1] = 2 * ybit3d[j] - 1; - va[2] = 2 * zbit3d[j] - 1; - matrix_times_vector3d (rt, va, vb); - *x += ((vb[0] + 1) / 2) << k; - *y += ((vb[1] + 1) / 2) << k; - *z += ((vb[2] + 1) / 2) << k; - if (k > 0) - { - matrix_copy3d (rt, rq); - if (k == n-1 && closed_p) - matrix_multiply3d (rq, s3d[j], rt); - else - matrix_multiply3d (rq, r3d[j], rt); - } - } -} - - -/* Rendering the curve - */ - - -/* Draw a sphere at the intersection of two tubes, to round off - the connection between them. Use one of our cache display lists - of unit spheres in various sizes. - */ -static int -draw_joint (ModeInfo *mi, XYZ p, GLfloat diam, int faces, int wire) -{ - hilbert_configuration *bp = &bps[MI_SCREEN(mi)]; - int i; - diam *= 0.99; /* try to clean up the edges a bit */ - - if (faces <= 4) return 0; /* too small to see */ - - glPushMatrix(); - glTranslatef (p.x, p.y, p.z); - glScalef (diam, diam, diam); - - /* i = unit_sphere (faces, faces, wire);*/ - - /* if (!bp->dlists[faces][0]) abort(); */ - /* if (!bp->dlist_polys[faces][0]) abort(); */ - - glCallList (bp->dlists[faces][0]); - i = bp->dlist_polys[faces][0]; - glPopMatrix(); - return i; -} - - -/* Draw a tube, and associated end caps. Use one of our cache display lists - of unit tubes in various sizes. Pick the resolution of the tube based - on how large it's being drawn. It's ok to get chunkier if the thing is - only a few pixels wide on the screen. - */ -static Bool -draw_segment (ModeInfo *mi, - XYZ p0, XYZ p1, /* from, to */ - int t, int end_t, /* value and range */ - GLfloat path_start, GLfloat path_end, /* clipping */ - Bool head_cap_p, - int *last_colorP) -{ - hilbert_configuration *bp = &bps[MI_SCREEN(mi)]; - - double t0 = (double) (t-1) / end_t; /* ratio of p[01] along curve */ - double t1 = (double) t / end_t; - - int wire = MI_IS_WIREFRAME(mi); - int owire = wire; - GLfloat dd = bp->diam; - int faces; - - if (path_start >= t1) /* whole segment is before clipping region */ - return False; - if (path_end < t0) /* whole segment is after clipping region */ - return False; - - - if (bp->twodee_p) dd *= 2; /* more polys in 2D mode */ - - faces = (dd > 0.040 ? dlist_faces[0] : - dd > 0.020 ? dlist_faces[1] : - dd > 0.010 ? dlist_faces[2] : - dd > 0.005 ? dlist_faces[3] : - dd > 0.002 ? dlist_faces[4] : - 1); - - /* In 2d mode, we can drop into wireframe mode and it still looks ok... */ - if (faces == 1) - { - if (bp->twodee_p) - wire = True; - else - faces = 3; - } - - if (wire && !owire) - { - glDisable (GL_DEPTH_TEST); - glDisable (GL_CULL_FACE); - glDisable (GL_LIGHTING); - } - - if (t0 < path_start) - { - XYZ p; - GLfloat seg_range = t1 - t0; - GLfloat clip_range = path_start - t0; - GLfloat ratio = clip_range / seg_range; - p.x = p0.x + ((p1.x - p0.x) * ratio); - p.y = p0.y + ((p1.y - p0.y) * ratio); - p.z = p0.z + ((p1.z - p0.z) * ratio); - p0 = p; - } - - if (t1 > path_end) - { - XYZ p; - GLfloat seg_range = t1 - t0; - GLfloat clip_range = path_end - t0; - GLfloat ratio = clip_range / seg_range; - p.x = p0.x + ((p1.x - p0.x) * ratio); - p.y = p0.y + ((p1.y - p0.y) * ratio); - p.z = p0.z + ((p1.z - p0.z) * ratio); - p1 = p; - } - - if (p0.x == p1.x && - p0.y == p1.y && - p0.z == p1.z) - return False; - - { - int segs = bp->ncolors * (t1 - t0); - int i; - XYZ p0b, p1b; - - if (segs < 1) segs = 1; - - for (i = 0; i < segs; i++) - { - int j = i + 1; - GLfloat color[4] = {0.0, 0.0, 0.0, 1.0}; - GLfloat t0b; - int c; - - p0b.x = p0.x + i * ((p1.x - p0.x) / segs); - p0b.y = p0.y + i * ((p1.y - p0.y) / segs); - p0b.z = p0.z + i * ((p1.z - p0.z) / segs); - - p1b.x = p0.x + j * ((p1.x - p0.x) / segs); - p1b.y = p0.y + j * ((p1.y - p0.y) / segs); - p1b.z = p0.z + j * ((p1.z - p0.z) / segs); - - - - /* #### this isn't quite right */ - t0b = t0 + i * (t1 - t0) / segs; - - c = bp->ncolors * t0b; - if (c >= bp->ncolors) c = bp->ncolors-1; - - /* Above depth 6, was spending 5% of the time in glMateralfv, - so only set the color if it's different. */ - - if (c != *last_colorP) - { - color[0] = bp->colors[c].red / 65536.0; - color[1] = bp->colors[c].green / 65536.0; - color[2] = bp->colors[c].blue / 65536.0; - if (wire) - glColor3fv (color); - else - glMaterialfv (GL_FRONT, GL_AMBIENT_AND_DIFFUSE, color); - *last_colorP = c; - } - - - if (wire) - { - glBegin (GL_LINES); - glVertex3f (p0b.x, p0b.y, p0b.z); - glVertex3f (p1b.x, p1b.y, p1b.z); - glEnd (); - mi->polygon_count++; - } - else - { - /* mi->polygon_count += tube (p0b.x, p0b.y, p0b.z, - p1b.x, p1b.y, p1b.z, - bp->diam, 0, faces, True, - 0, wire); - */ - - /* Render a rotated and scaled prefab unit tube. - - There's probably a slightly more concise way to do this, - but since they're all at right angles at least we don't - have to use atan(). - */ - GLfloat s; - glPushMatrix(); - glTranslatef (p0b.x, p0b.y, p0b.z); - - if (p1b.x > p0b.x) - { - s = p1b.x - p0b.x; - } - else if (p1b.x < p0b.x) - { - glRotatef (180, 0, 0, 1); - s = p0b.x - p1b.x; - } - else if (p1b.y > p0b.y) { - glRotatef (90, 0, 0, 1); - s = p1b.y - p0b.y; - } - else if (p1b.y < p0b.y) - { - glRotatef (-90, 0, 0, 1); - s = p0b.y - p1b.y; - } - else if (p1b.z > p0b.z) - { - glRotatef (-90, 0, 1, 0); - s = p1b.z - p0b.z; - } - else /* (p1b.z < p0b.z) */ - { - glRotatef (90, 0, 1, 0); - s = p0b.z - p1b.z; - } - - glScalef (s, bp->diam, bp->diam); - glCallList (bp->dlists[faces][1]); - mi->polygon_count += bp->dlist_polys[faces][1]; - glPopMatrix(); - - - /* If this is the bleeding edge, cap it too. */ - if (head_cap_p) { - mi->polygon_count += draw_joint (mi, p0b, bp->diam, faces, wire); - head_cap_p = False; - } - } - } - - /* If the path is closed, close it. This segment doesn't animate - like the others, but, oh well. - */ - if (! wire) - mi->polygon_count += draw_joint (mi, p1b, bp->diam, faces, wire); - } - - return True; -} - - -static void -mem(void) -{ - fprintf (stderr, "%s: out of memory\n", progname); - exit (1); -} - - -/* Render the whole thing, but omit segments that lie outside of - the path_start and path_end ratios. - */ -static void -hilbert (ModeInfo *mi, int size, GLfloat path_start, GLfloat path_end) -{ - hilbert_configuration *bp = &bps[MI_SCREEN(mi)]; - int wire = MI_IS_WIREFRAME(mi); - - GLfloat w = pow(2, size); - int end_t = (bp->twodee_p ? w * w : w * w * w); - - XYZ prev = { 0, }; - XYZ first = { 0, }; - Bool first_p = False; - Bool any_p = False; - int t; - Bool fill_cache_p = False; - hilbert_cache *cc; - int last_color = -1; - - /* Do this here since at higher resolutions, we turn wireframe on - and off. */ - - if (! wire) - { - glEnable (GL_NORMALIZE); - glEnable (GL_DEPTH_TEST); - glEnable (GL_CULL_FACE); - glEnable (GL_LIGHTING); - glEnable (GL_POLYGON_OFFSET_FILL); - } - - - if (!bp->caches) - { - bp->caches = (hilbert_cache **) - calloc (max_depth + 2, sizeof (*bp->caches)); - if (!bp->caches) mem(); - fill_cache_p = True; - } - - cc = bp->caches[size]; - if (! cc) - { - cc = (hilbert_cache *) calloc (1, sizeof (*cc)); - cc->values = (XYZb *) calloc (end_t + 1, sizeof (*cc->values)); - if (!cc->values) mem(); - cc->size = end_t; - bp->caches[size] = cc; - fill_cache_p = True; - } - - for (t = 0; t < end_t; t++) - { - int x, y, z; - XYZ c; - XYZb *cb; - - if (!fill_cache_p) - { - cb = &cc->values[t]; - x = cb->x; - y = cb->y; - z = cb->z; - } - else - { - if (!bp->twodee_p) - t_to_xyz (size, t, &x, &y, &z, bp->closed_p); - else - { - t_to_xy (size, t, &x, &y, bp->closed_p); - z = w/2; - } - cb = &cc->values[t]; - cb->x = x; - cb->y = y; - cb->z = z; - } - - c.x = (GLfloat) x / w - 0.5; - c.y = (GLfloat) y / w - 0.5; - c.z = (GLfloat) z / w - 0.5; - - /* #### We could save some polygons by not drawing the spheres - between colinear segments. */ - - if (t > 0) { - if (draw_segment (mi, prev, c, t, end_t, path_start, path_end, !any_p, - &last_color)) - any_p = True; - } - prev = c; - if (! first_p) { - first = c; - first_p = True; - } - } - - if (bp->closed_p && path_end >= 1.0) { - draw_segment (mi, prev, first, t, end_t, path_start, path_end, !any_p, - &last_color); - } -} - - - -/* Window management, etc - */ -ENTRYPOINT void -reshape_hilbert (ModeInfo *mi, int width, int height) -{ - GLfloat h = (GLfloat) height / (GLfloat) width; - int y = 0; - - if (width > height * 5) { /* tiny window: show middle */ - height = width * 9/16; - y = -height/2; - h = height / (GLfloat) width; - } - - glViewport (0, y, (GLint) width, (GLint) height); - - glMatrixMode(GL_PROJECTION); - glLoadIdentity(); - gluPerspective (30.0, 1/h, 1.0, 100.0); - - glMatrixMode(GL_MODELVIEW); - glLoadIdentity(); - gluLookAt( 0.0, 0.0, 30.0, - 0.0, 0.0, 0.0, - 0.0, 1.0, 0.0); - - { - GLfloat s = (MI_WIDTH(mi) < MI_HEIGHT(mi) - ? (MI_WIDTH(mi) / (GLfloat) MI_HEIGHT(mi)) - : 1); - glScalef (s, s, s); - } - - glClear(GL_COLOR_BUFFER_BIT); -} - - -ENTRYPOINT Bool -hilbert_handle_event (ModeInfo *mi, XEvent *event) -{ - hilbert_configuration *bp = &bps[MI_SCREEN(mi)]; - - if (gltrackball_event_handler (event, bp->trackball, - MI_WIDTH (mi), MI_HEIGHT (mi), - &bp->button_down_p)) - return True; - else if (event->xany.type == KeyPress) - { - KeySym keysym; - char c = 0; - XLookupString (&event->xkey, &c, 1, &keysym, 0); - if (c == '+' || c == '=' || - keysym == XK_Up || keysym == XK_Right || keysym == XK_Next) - { - bp->depth++; - if (bp->depth > max_depth) bp->depth = max_depth; - return True; - } - else if (c == '-' || c == '_' || - keysym == XK_Down || keysym == XK_Left || keysym == XK_Prior) - { - bp->depth--; - if (bp->depth < 1) bp->depth = 1; - return True; - } - else if (screenhack_event_helper (MI_DISPLAY(mi), MI_WINDOW(mi), event)) - { - bp->depth += bp->depth_tick; - if (bp->depth > max_depth-1) - { - bp->depth = max_depth; - bp->depth_tick = -1; - } - else if (bp->depth <= 1) - { - bp->depth = 1; - bp->depth_tick = 1; - } - return True; - } - } - - return False; -} - - -ENTRYPOINT void -init_hilbert (ModeInfo *mi) -{ - hilbert_configuration *bp; - int wire = MI_IS_WIREFRAME(mi); - int i; - - MI_INIT (mi, bps); - - bp = &bps[MI_SCREEN(mi)]; - - bp->depth = 2; - bp->depth_tick = 1; - bp->path_start = 0.0; - bp->path_end = 0.0; - bp->path_tick = 1; - - if (thickness < 0.01) thickness = 0.01; - if (thickness > 1.0) thickness = 1.0; - - if (speed <= 0) speed = 0.0001; - if (max_depth < 2) max_depth = 2; - if (max_depth > 20) max_depth = 20; /* hard limit due to hilbert_cache */ - - if (bp->depth > max_depth-1) bp->depth = max_depth-1; - - bp->glx_context = init_GL(mi); - - reshape_hilbert (mi, MI_WIDTH(mi), MI_HEIGHT(mi)); - - if (!wire) - { - GLfloat pos[4] = {1.0, 1.0, 1.0, 0.0}; - GLfloat amb[4] = {0.0, 0.0, 0.0, 1.0}; - GLfloat dif[4] = {1.0, 1.0, 1.0, 1.0}; - GLfloat spc[4] = {0.0, 1.0, 1.0, 1.0}; - - glEnable (GL_LIGHTING); - glEnable (GL_LIGHT0); - - glLightfv(GL_LIGHT0, GL_POSITION, pos); - glLightfv(GL_LIGHT0, GL_AMBIENT, amb); - glLightfv(GL_LIGHT0, GL_DIFFUSE, dif); - glLightfv(GL_LIGHT0, GL_SPECULAR, spc); - } - - { - double spin_speed = 0.04; - double tilt_speed = spin_speed / 10; - double wander_speed = 0.008; - double spin_accel = 0.01; - - bp->rot = make_rotator (do_spin ? spin_speed : 0, - do_spin ? spin_speed : 0, - do_spin ? spin_speed : 0, - spin_accel, - do_wander ? wander_speed : 0, - do_spin); - bp->rot2 = make_rotator (0, 0, 0, 0, tilt_speed, True); - bp->trackball = gltrackball_init (True); - } - - if (mode_str && !strcasecmp(mode_str, "2d")) - bp->twodee_p = True; - else if (mode_str && (!strcasecmp(mode_str, "3d"))) - bp->twodee_p = False; - else - { - if (! (mode_str && !strcasecmp(mode_str, "random"))) - fprintf (stderr, "%s: 'mode' must be '2d', '3d', or 'random'\n", - progname); - bp->twodee_p = ((random() % 3) == 0); - } - - - if (ends_str && (!strcasecmp(ends_str, "closed"))) - bp->closed_p = True; - else if (ends_str && (!strcasecmp(ends_str, "open"))) - bp->closed_p = False; - else - { - if (! (ends_str && !strcasecmp(ends_str, "random"))) - fprintf (stderr, "%s: 'ends' must be 'open', 'closed', or 'random'\n", - progname); - bp->closed_p = ((random() % 3) != 0); - } - - - /* More colors results in more polygons (more tube segments) so keeping - this small is worthwhile. */ - bp->ncolors = (bp->twodee_p ? 1024 : 128); - - if (bp->closed_p) - { - /* Since the path is closed, colors must also be a closed loop */ - bp->colors = (XColor *) calloc(bp->ncolors, sizeof(XColor)); - make_smooth_colormap (0, 0, 0, - bp->colors, &bp->ncolors, - False, 0, False); - } - else - { - /* Since the path is open, first and last colors should differ. */ - bp->ncolors *= 2; - bp->colors = (XColor *) calloc(bp->ncolors, sizeof(XColor)); - make_uniform_colormap (0, 0, 0, - bp->colors, &bp->ncolors, - False, 0, False); - bp->ncolors /= 2; - } - - /* Generate display lists for several different resolutions of - a unit tube and a unit sphere. - */ - for (i = 0; i < countof(dlist_faces); i++) - { - int faces = dlist_faces[i]; - bp->dlists[faces][0] = glGenLists (1); - - glNewList (bp->dlists[faces][0], GL_COMPILE); - bp->dlist_polys[faces][0] = unit_sphere (faces, faces, wire); - glEndList (); - - bp->dlists[faces][1] = glGenLists (1); - - glNewList (bp->dlists[faces][1], GL_COMPILE); - bp->dlist_polys[faces][1] = - tube (0, 0, 0, 1, 0, 0, - 1, 0, faces, True, 0, wire); - glEndList (); - } -} - - -ENTRYPOINT void -draw_hilbert (ModeInfo *mi) -{ - hilbert_configuration *bp = &bps[MI_SCREEN(mi)]; - Display *dpy = MI_DISPLAY(mi); - Window window = MI_WINDOW(mi); - int wire = MI_IS_WIREFRAME(mi); - - static const GLfloat bspec[4] = {1.0, 1.0, 1.0, 1.0}; - static const GLfloat bshiny = 128.0; - GLfloat bcolor[4] = {1.0, 1.0, 1.0, 1.0}; - - if (!bp->glx_context) - return; - - glXMakeCurrent(MI_DISPLAY(mi), MI_WINDOW(mi), *bp->glx_context); - - glShadeModel(GL_SMOOTH); - - if (! wire) - { - glEnable (GL_NORMALIZE); - glEnable (GL_DEPTH_TEST); - glEnable (GL_CULL_FACE); - glEnable (GL_LIGHTING); - glEnable (GL_LIGHT0); - glEnable (GL_POLYGON_OFFSET_FILL); - } - - glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT); - - glPushMatrix (); - - glScalef(1.1, 1.1, 1.1); - - { - double x, y, z, z2; - get_position (bp->rot, &x, &y, &z, !bp->button_down_p); - glTranslatef((x - 0.5) * 8, - (y - 0.5) * 8, - (z - 0.5) * 15); - - gltrackball_rotate (bp->trackball); - - get_rotation (bp->rot, &x, &y, &z, !bp->button_down_p); - - if (bp->twodee_p && do_spin) /* face front */ - { - double max = 70; - get_position (bp->rot2, &x, &y, &z2, !bp->button_down_p); - glRotatef (max/2 - x*max, 1, 0, 0); - glRotatef (max/2 - y*max, 0, 1, 0); - glRotatef (z * 360, 0, 0, 1); /* but upside down is ok */ - } - else - { - glRotatef (x * 360, 1, 0, 0); - glRotatef (y * 360, 0, 1, 0); - glRotatef (z * 360, 0, 0, 1); - } - } - - mi->polygon_count = 0; - - glMaterialfv (GL_FRONT, GL_SPECULAR, bspec); - glMateriali (GL_FRONT, GL_SHININESS, bshiny); - glMaterialfv (GL_FRONT, GL_AMBIENT_AND_DIFFUSE, bcolor); - - glScalef (8,8,8); - glTranslatef (0.1, 0.1, 0); - - if (!do_spin && !bp->twodee_p) - { - /* tilt the cube a little, and make the start and end visible */ - glTranslatef (-0.1, 0, 0); - glRotatef (140, 0, 1, 0); - glRotatef (30, 1, 0, 0); - } - - if (wire) - glLineWidth (bp->depth > 4 ? 1.0 : - bp->depth > 3 ? 2.0 : - 3.0); - - if (bp->path_tick > 0) /* advancing the end point, [0 - N) */ - { /* drawing 1 partial path, 1st time only. */ - - if (!bp->button_down_p) - bp->path_end += 0.01 * speed; - - if (bp->path_end >= 1.0) - { - bp->path_start = 0.0; - bp->path_end = 1.0; - bp->path_tick = -1; - bp->pause = PAUSE_TICKS; - } - - bp->diam = thickness / pow (2, bp->depth); - hilbert (mi, bp->depth, bp->path_start, bp->path_end); - mi->recursion_depth = bp->depth + bp->path_start; - } - - else /* advancing the start point, (N - 1] */ - { /* drawing 2 paths at different rez. */ - if (bp->pause) - { - bp->pause--; - } - else - { - if (!bp->button_down_p) - bp->path_start += 0.01 * speed; - - if (bp->path_start > 1.0) - { - bp->path_start = 0.0; - bp->depth += bp->depth_tick; - bp->pause = PAUSE_TICKS; - if (bp->depth > max_depth-1) - { - bp->depth = max_depth; - bp->depth_tick = -1; - } - else if (bp->depth <= 1) - { - bp->depth = 1; - bp->depth_tick = 1; - } - } - } - - { - GLfloat d1 = thickness / pow (2, bp->depth); - GLfloat d2 = thickness / pow (2, bp->depth + bp->depth_tick); - bp->diam = (d1 * (1 - bp->path_start) + - d2 * bp->path_start); - } - - /* First time around, start is 0, and end animates forward. - Then, to display the next higher resolution, we render - depth=N while increasing start and leaving end=1.0; - and simultaneously animationg depth=N+1 with start=0 and - end increasing by the same amount. - - The two paths aren't actually connected together at the - resolution-change interface, and sometimes they overlap, - but things go fast enough that it's hard to spot those - glitches, so it looks ok. - */ - glPolygonOffset (0, 0); - hilbert (mi, bp->depth, bp->path_start, bp->path_end); - - glPolygonOffset (1.0, 1.0); - hilbert (mi, bp->depth + bp->depth_tick, 0.0, bp->path_start); - - mi->recursion_depth = bp->depth + (bp->depth_tick * bp->path_start); - } - - glPopMatrix (); - - if (mi->fps_p) do_fps (mi); - glFinish(); - - glXSwapBuffers(dpy, window); -} - - -ENTRYPOINT void -free_hilbert (ModeInfo *mi) -{ - hilbert_configuration *bp = &bps[MI_SCREEN(mi)]; - int i; - - if (!bp->glx_context) return; - glXMakeCurrent(MI_DISPLAY(mi), MI_WINDOW(mi), *bp->glx_context); - - if (bp->trackball) gltrackball_free (bp->trackball); - if (bp->rot) free_rotator (bp->rot); - if (bp->rot2) free_rotator (bp->rot2); - if (bp->colors) free (bp->colors); - if (bp->caches) - for (i = 0; i < max_depth + 2; i++) - if (bp->caches[i]) - { - free (bp->caches[i]->values); - free (bp->caches[i]); - } - /* #### free dlists */ -} - -XSCREENSAVER_MODULE ("Hilbert", hilbert) - -#endif /* USE_GL */ |