ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/hd/rhd_glx.c
(Generate patch)

Comparing ray/src/hd/rhd_glx.c (file contents):
Revision 3.3 by gregl, Wed Dec 24 15:55:43 1997 UTC vs.
Revision 3.30 by schorsch, Mon Jul 21 22:30:18 2003 UTC

# Line 1 | Line 1
1 /* Copyright (c) 1997 Silicon Graphics, Inc. */
2
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ SGI";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6
4   /*
5   * OpenGL GLX driver for holodeck display.
6 < * Based on x11 driver.
6 > * Based on old GLX driver using cones.
7 > *
8 > * Define symbol STEREO for stereo viewing.
9 > * Define symbol DOBJ for display object viewing.
10   */
11  
12 + #ifdef NOSTEREO
13 + #ifdef STEREO
14 + #undef STEREO
15 + #else
16 + #undef NOSTEREO
17 + #endif
18 + #endif
19 +
20   #include "standard.h"
13 #include "tonemap.h"
14 #include "rhdriver.h"
21  
22 < #include  <GL/glx.h>
23 <
24 < #include  "x11icon.h"
25 <
20 < #ifndef FEQ
21 < #define FEQ(a,b)        ((a)-(b) <= FTINY && (a)-(b) >= -FTINY)
22 > #include <GL/glx.h>
23 > #include <GL/glu.h>
24 > #ifdef STEREO
25 > #include <X11/extensions/SGIStereo.h>
26   #endif
27  
28 < #ifndef int4
29 < #define int4    int
28 > #include "rhd_sample.h"
29 > #ifdef DOBJ
30 > #include "rhdobj.h"
31   #endif
32  
33 < #ifndef FREEPCT
34 < #define FREEPCT         10              /* percentage of values to free */
33 > #include "x11icon.h"
34 >
35 > #ifndef RAYQLEN
36 > #define RAYQLEN         1024            /* max. rays to queue before flush */
37   #endif
38  
39 < #ifndef NCONEV
40 < #define NCONEV          7               /* number of cone base vertices */
39 > #ifndef FEQ
40 > #define FEQ(a,b)        ((a)-(b) <= FTINY && (a)-(b) >= -FTINY)
41   #endif
35 #ifndef CONEH
36 #define CONEH           3.              /* cone height (fraction of depth) */
37 #endif
38 #ifndef CONEW
39 #define CONEW           0.05            /* cone width (fraction of screen) */
40 #endif
41 #ifndef DIRPEN
42 #define DIRPEN          0.001           /* direction penalty factor */
43 #endif
44 #ifndef VALUA
45 #define VALUA           16              /* target value area (pixels) */
46 #endif
42  
43   #define GAMMA           1.4             /* default gamma correction */
44  
45 + #define FRAMESTATE(s)   (((s)&(ShiftMask|ControlMask))==(ShiftMask|ControlMask))
46 +
47   #define MOVPCT          7               /* percent distance to move /frame */
48   #define MOVDIR(b)       ((b)==Button1 ? 1 : (b)==Button2 ? 0 : -1)
49   #define MOVDEG          (-5)            /* degrees to orbit CW/down /frame */
50   #define MOVORB(s)       ((s)&ShiftMask ? 1 : (s)&ControlMask ? -1 : 0)
51  
52 + #ifndef TARGETFPS
53 + #define TARGETFPS       2.0             /* target frames/sec during motion */
54 + #endif
55 +
56   #define MINWIDTH        480             /* minimum graphics window width */
57   #define MINHEIGHT       400             /* minimum graphics window height */
58  
# Line 59 | Line 60 | static char SCCSid[] = "$SunId$ SGI";
60  
61   #define BORWIDTH        5               /* border width */
62  
63 + #define setstereobuf(bid)       (glXWaitGL(), \
64 +                                XSGISetStereoBuffer(ourdisplay, gwind, bid), \
65 +                                glXWaitX())
66 +
67   #define  ourscreen      DefaultScreen(ourdisplay)
68   #define  ourroot        RootWindow(ourdisplay,ourscreen)
69   #define  ourmask        (StructureNotifyMask|ExposureMask|KeyPressMask|\
# Line 68 | Line 73 | static char SCCSid[] = "$SunId$ SGI";
73  
74   struct driver   odev;                   /* global device driver structure */
75  
76 + char odev_args[64];                     /* command arguments */
77 +
78 + #ifdef STEREO
79 + static VIEW     vwright;                /* right eye view */
80 + #endif
81 +
82 + static int      rayqleft = 0;           /* rays left to queue before flush */
83 +
84   static XEvent  currentevent;            /* current event */
85  
86   static int  mapped = 0;                 /* window is mapped? */
# Line 78 | Line 91 | static XVisualInfo  *ourvinf;          /* our visual informati
91   static Window  gwind = 0;               /* our graphics window */
92   static GLXContext       gctx;           /* our GLX context */
93  
81 static double   mindepth = FHUGE;       /* minimum depth value so far */
82 static double   maxdepth = 0.;          /* maximum depth value so far */
83
94   static double   pwidth, pheight;        /* pixel dimensions (mm) */
95  
96 < static FVECT    conev[NCONEV];          /* drawing cone */
87 < static double   coneh;                  /* cone height */
96 > static double   mindpth, maxdpth;       /* min. and max. depth */
97  
98 + double  dev_zmin, dev_zmax;             /* fore and aft clipping plane dist. */
99 +
100   static int      inpresflags;            /* input result flags */
101  
102   static int      headlocked = 0;         /* lock vertical motion */
103  
104 < static int      quicken = 0;            /* quicker, sloppier update rate? */
104 > static int  resizewindow(), getevent(), getkey(), moveview(), wipeclean(),
105 >                setglpersp(), getframe(), getmove(), fixwindow(), mytmflags();
106  
107 < static struct {
108 <        float           (*wp)[3];       /* world intersection point array */
109 <        int4            *wd;            /* world direction array */
110 <        TMbright        *brt;           /* encoded brightness array */
111 <        BYTE            (*chr)[3];      /* encoded chrominance array */
112 <        BYTE            (*rgb)[3];      /* tone-mapped color array */
101 <        BYTE            *alpha;         /* alpha values */
102 <        int             nl;             /* count of values */
103 <        int             bl, tl;         /* bottom and top (next) value index */
104 <        int             tml;            /* next value needing tone-mapping */
105 <        int             drl;            /* next value in need of drawing */
106 <        char            *base;          /* base of allocated memory */
107 < }       rV;                     /* our collection of values */
107 > #ifdef STEREO
108 > static int  pushright(), popright();
109 > #endif
110 > static double   getdistance();
111 > #define mapdepth(d) ((d)> 0.9995 ? FHUGE: dev_zmin/ \
112 >                     (1.-(d)*(1.-dev_zmin/dev_zmax)))
113  
114 < static int      *valmap = NULL;         /* sorted map of screen values */
110 < static int      vmaplen = 0;            /* value map length */
114 > extern time_t   time();
115  
112 #define redraw()        (rV.drl = rV.bl)
116  
117 < static int  resizewindow(), getevent(), getkey(), moveview(),
115 <                setGLview(), getmove(), fixwindow(), mytmflags(),
116 <                drawvalue(), valcmp(), clralphas(), setalphas(), mergalphas(),
117 <                IndexValue(), Compost(), FindValue(), TMapValues(),
118 <                AllocValues(), FreeValues();
119 <
120 < extern int4     encodedir();
121 < extern double   fdir2diff(), dir2diff();
122 <
123 <
124 < dev_open(id)                    /* initialize X11 driver */
117 > dev_open(id)                    /* initialize GLX driver */
118   char  *id;
119   {
120          extern char     *getenv();
121 +        static RGBPRIMS myprims = STDPRIMS;
122          static int      atlBest[] = {GLX_RGBA, GLX_RED_SIZE,8,
123                                  GLX_GREEN_SIZE,8, GLX_BLUE_SIZE,8,
124 <                                GLX_ALPHA_SIZE,8, GLX_DEPTH_SIZE,15,
125 <                                None};
132 <        char    *gv;
124 >                                GLX_DEPTH_SIZE,15, None};
125 >        char    *ev;
126          double  gamval = GAMMA;
127 +        RGBPRIMP        dpri = stdprims;
128          XSetWindowAttributes    ourwinattr;
129          XWMHints        ourxwmhints;
130          XSizeHints      oursizhints;
131 +                                        /* check for unsupported stereo */
132 + #ifdef NOSTEREO
133 +        error(INTERNAL, "stereo display driver unavailable");
134 + #endif
135                                          /* open display server */
136          ourdisplay = XOpenDisplay(NULL);
137          if (ourdisplay == NULL)
138                  error(USER, "cannot open X-windows; DISPLAY variable set?\n");
139 + #ifdef STEREO
140 +        switch (XSGIQueryStereoMode(ourdisplay, ourroot)) {
141 +        case STEREO_TOP:
142 +        case STEREO_BOTTOM:
143 +                break;
144 +        case STEREO_OFF:
145 +                error(USER,
146 +        "wrong video mode: run \"/usr/gfx/setmon -n STR_TOP\" first");
147 +        case X_STEREO_UNSUPPORTED:
148 +                error(USER, "stereo mode not supported on this screen");
149 +        default:
150 +                error(INTERNAL, "unknown stereo mode");
151 +        }
152 + #endif
153                                          /* find a usable visual */
154          ourvinf = glXChooseVisual(ourdisplay, ourscreen, atlBest);
155          if (ourvinf == NULL)
# Line 145 | Line 157 | char  *id;
157                                          /* get a context */
158          gctx = glXCreateContext(ourdisplay, ourvinf, NULL, GL_TRUE);
159                                          /* set gamma and tone mapping */
160 <        if ((gv = XGetDefault(ourdisplay, "radiance", "gamma")) != NULL
161 <                        || (gv = getenv("DISPLAY_GAMMA")) != NULL)
162 <                gamval = atof(gv);
163 <        if (tmInit(mytmflags(), stdprims, gamval) == NULL)
160 >        if ((ev = XGetDefault(ourdisplay, "radiance", "gamma")) != NULL
161 >                        || (ev = getenv("DISPLAY_GAMMA")) != NULL)
162 >                gamval = atof(ev);
163 >        if ((ev = getenv("DISPLAY_PRIMARIES")) != NULL &&
164 >                        sscanf(ev, "%f %f %f %f %f %f %f %f",
165 >                                &myprims[RED][CIEX],&myprims[RED][CIEY],
166 >                                &myprims[GRN][CIEX],&myprims[GRN][CIEY],
167 >                                &myprims[BLU][CIEX],&myprims[BLU][CIEY],
168 >                                &myprims[WHT][CIEX],&myprims[WHT][CIEY]) >= 6)
169 >                dpri = myprims;
170 >        if (tmInit(mytmflags(), dpri, gamval) == NULL)
171                  error(SYSTEM, "not enough memory in dev_open");
172                                          /* open window */
173          ourwinattr.background_pixel = ourblack;
# Line 159 | Line 178 | char  *id;
178                                  ourvinf->visual, AllocNone);
179          gwind = XCreateWindow(ourdisplay, ourroot, 0, 0,
180                  DisplayWidth(ourdisplay,ourscreen)-2*BORWIDTH,
181 + #ifdef STEREO
182 +                (DisplayHeight(ourdisplay,ourscreen)-2*BORWIDTH)/2,
183 + #else
184                  DisplayHeight(ourdisplay,ourscreen)-2*BORWIDTH,
185 + #endif
186                  BORWIDTH, ourvinf->depth, InputOutput, ourvinf->visual,
187                  CWBackPixel|CWBorderPixel|CWColormap|CWEventMask, &ourwinattr);
188          if (gwind == 0)
# Line 172 | Line 195 | char  *id;
195                          gwind, x11icon_bits, x11icon_width, x11icon_height);
196          XSetWMHints(ourdisplay, gwind, &ourxwmhints);
197          oursizhints.min_width = MINWIDTH;
198 + #ifdef STEREO
199 +        oursizhints.min_height = MINHEIGHT/2;
200 +        oursizhints.max_width = DisplayWidth(ourdisplay,ourscreen)-2*BORWIDTH;
201 +        oursizhints.max_height = (DisplayHeight(ourdisplay,ourscreen) -
202 +                                        2*BORWIDTH)/2;
203 +        oursizhints.flags = PMinSize|PMaxSize;
204 + #else
205          oursizhints.min_height = MINHEIGHT;
206          oursizhints.flags = PMinSize;
207 + #endif
208          XSetNormalHints(ourdisplay, gwind, &oursizhints);
209                                          /* set GLX context */
210          glXMakeCurrent(ourdisplay, gwind, gctx);
211          glEnable(GL_DEPTH_TEST);
212          glDepthFunc(GL_LEQUAL);
213 <        glShadeModel(GL_FLAT);
213 >        glShadeModel(GL_SMOOTH);
214          glDisable(GL_DITHER);
215 +        glDisable(GL_CULL_FACE);
216                                          /* figure out sensible view */
217          pwidth = (double)DisplayWidthMM(ourdisplay, ourscreen) /
218                          DisplayWidth(ourdisplay, ourscreen);
219          pheight = (double)DisplayHeightMM(ourdisplay, ourscreen) /
220                          DisplayHeight(ourdisplay, ourscreen);
221 <        copystruct(&odev.v, &stdview);
221 > #ifdef STEREO
222 >        pheight *= 2.;
223 >        setstereobuf(STEREO_BUFFER_LEFT);
224 > #endif
225 >        checkglerr("setting rendering parameters");
226 >        odev.v = stdview;
227          odev.v.type = VT_PER;
228                                          /* map the window */
229          XMapWindow(ourdisplay, gwind);
230          dev_input();                    /* sets size and view angles */
231 <        glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
232 <                                        /* allocate our value list */
233 <        if (!AllocValues(DisplayWidth(ourdisplay,ourscreen) *
197 <                        DisplayHeight(ourdisplay,ourscreen) / VALUA))
231 >                                        /* allocate our samples */
232 >        if (!smInit(DisplayWidth(ourdisplay,ourscreen) *
233 >                        DisplayHeight(ourdisplay,ourscreen) / 10))
234                  error(SYSTEM, "insufficient memory for value storage");
235 +        mindpth = FHUGE; maxdpth = FTINY;
236          odev.name = id;
237          odev.ifd = ConnectionNumber(ourdisplay);
238   }
# Line 203 | Line 240 | char  *id;
240  
241   dev_close()                     /* close our display and free resources */
242   {
243 +        smInit(0);
244 + #ifdef DOBJ
245 +        dobj_cleanup();
246 + #endif
247          glXMakeCurrent(ourdisplay, None, NULL);
248          glXDestroyContext(ourdisplay, gctx);
249          XDestroyWindow(ourdisplay, gwind);
# Line 210 | Line 251 | dev_close()                    /* close our display and free resources
251          XCloseDisplay(ourdisplay);
252          ourdisplay = NULL;
253          tmDone(NULL);
213        FreeValues();
254          odev.v.type = 0;
255          odev.hres = odev.vres = 0;
256          odev.ifd = -1;
257   }
258  
259  
260 + dev_clear()                     /* clear our representation */
261 + {
262 +        smInit(rsL.max_samp);
263 +        wipeclean(1);
264 +        rayqleft = 0;                   /* hold off update */
265 + }
266 +
267 +
268   int
269   dev_view(nv)                    /* assign new driver view */
270   register VIEW   *nv;
271   {
272 +        double  d;
273 +
274          if (nv->type != VT_PER ||               /* check view legality */
275 <                        nv->horiz > 120. || nv->vert > 120.) {
275 >                        nv->horiz > 160. || nv->vert > 160.) {
276                  error(COMMAND, "illegal view type/angle");
277 <                nv->type = VT_PER;
277 >                nv->type = odev.v.type;
278                  nv->horiz = odev.v.horiz;
279                  nv->vert = odev.v.vert;
280                  return(0);
281          }
282          if (nv != &odev.v) {
283 <                if (!FEQ(nv->horiz,odev.v.horiz) ||     /* resize window? */
283 >                                                /* resize window? */
284 >                if (!FEQ(nv->horiz,odev.v.horiz) ||
285                                  !FEQ(nv->vert,odev.v.vert)) {
286                          int     dw = DisplayWidth(ourdisplay,ourscreen);
287                          int     dh = DisplayHeight(ourdisplay,ourscreen);
288  
289                          dw -= 25;       /* for window frame */
290                          dh -= 50;
291 + #ifdef STEREO
292 +                        dh /= 2;
293 + #endif
294                          odev.hres = 2.*VIEWDIST/pwidth *
295                                          tan(PI/180./2.*nv->horiz);
296                          odev.vres = 2.*VIEWDIST/pheight *
# Line 252 | Line 306 | register VIEW  *nv;
306                          XResizeWindow(ourdisplay, gwind, odev.hres, odev.vres);
307                          dev_input();    /* get resize event */
308                  }
309 <                copystruct(&odev.v, nv);
256 <                setGLview();
257 <                glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
258 <                redraw();
309 >                odev.v = *nv;   /* setview() already called */
310          }
311 + #ifdef STEREO
312 +        vwright = *nv;
313 +        d = eyesepdist / sqrt(nv->hn2);
314 +        VSUM(vwright.vp, nv->vp, nv->hvec, d);
315 +        /* setview(&vwright);   -- Unnecessary */
316 + #endif
317 +        wipeclean(0);
318          return(1);
319   }
320  
321  
322 < int
323 < dev_input()                     /* get X11 input */
322 > dev_section(gfn, pfn)           /* add octree for geometry rendering */
323 > char    *gfn, *pfn;
324   {
325 <        inpresflags = 0;
268 <
269 <        do
270 <                getevent();
271 <
272 <        while (XQLength(ourdisplay) > 0);
273 <
274 <        return(inpresflags);
325 >        /* unimplemented */
326   }
327  
328  
329 < dev_value(c, p, v)              /* add a pixel value to our list */
330 < COLR    c;
280 < FVECT   p, v;
329 > dev_auxcom(cmd, args)           /* process an auxiliary command */
330 > char    *cmd, *args;
331   {
332 <        register int    li;
333 <
284 <        li = rV.tl++;
285 <        if (rV.tl >= rV.nl)     /* get next leaf in ring */
286 <                rV.tl = 0;
287 <        if (rV.tl == rV.bl)     /* need to shake some free */
288 <                Compost(FREEPCT);
289 <        VCOPY(rV.wp[li], p);
290 <        rV.wd[li] = encodedir(v);
291 <        tmCvColrs(&rV.brt[li], rV.chr[li], c, 1);
292 < }
293 <
294 <
295 < int
296 < dev_flush()                     /* flush output */
297 < {
298 <        if (mapped) {
299 <                TMapValues(0);
300 <                while (rV.drl != rV.tl) {
301 <                        drawvalue(rV.drl);
302 <                        if (++rV.drl >= rV.nl)
303 <                                rV.drl = 0;
304 <                }
305 <                glFlush();
306 <        }
307 <        return(XPending(ourdisplay));
308 < }
309 <
310 <
311 < static int
312 < mytmflags()                     /* figure out tone mapping flags */
313 < {
314 <        extern char     *progname;
315 <        register char   *cp, *tail;
316 <                                        /* find basic name */
317 <        for (cp = tail = progname; *cp; cp++)
318 <                if (*cp == '/')
319 <                        tail = cp+1;
320 <        for (cp = tail; *cp && *cp != '.'; cp++)
321 <                ;
322 <        if (cp-tail == 3 && !strncmp(tail, "glx", 3))
323 <                return(TM_F_CAMERA);
324 <        if (cp-tail == 4 && !strncmp(tail, "glxh", 4))
325 <                return(TM_F_HUMAN);
326 <        error(USER, "illegal driver name");
327 < }
328 <
329 <
330 < static
331 < setGLview()                     /* set our GL view */
332 < {
333 <        double  xmin, xmax, ymin, ymax, zmin, zmax;
334 <        double  d, cx, sx, crad;
335 <        FVECT   vx, vy;
336 <        register int    i, j;
337 <                                        /* compute view frustum */
338 <        if (normalize(odev.v.vdir) == 0.0)
332 > #ifdef DOBJ
333 >        if (dobj_command(cmd, args) >= 0)
334                  return;
335 <        if (mindepth < maxdepth) {
336 <                zmin = 0.25*mindepth;
337 <                zmax = 4.0*(1.+CONEH)*maxdepth;
343 <        } else {
344 <                zmin = 0.01;
345 <                zmax = 1000.;
346 <        }
347 <        if (odev.v.vfore > FTINY)
348 <                zmin = odev.v.vfore;
349 <        if (odev.v.vaft > FTINY)
350 <                zmax = odev.v.vaft;
351 <        xmax = zmin * tan(PI/180./2. * odev.v.horiz);
352 <        xmin = -xmax;
353 <        d = odev.v.hoff * (xmax - xmin);
354 <        xmin += d; xmax += d;
355 <        ymax = zmin * tan(PI/180./2. * odev.v.vert);
356 <        ymin = -ymax;
357 <        d = odev.v.voff * (ymax - ymin);
358 <        ymin += d; ymax += d;
359 <                                        /* set view matrix */
360 <        glMatrixMode(GL_PROJECTION);
361 <        glLoadIdentity();
362 <        glFrustum(xmin, xmax, ymin, ymax, zmin, zmax);
363 <        gluLookAt(odev.v.vp[0], odev.v.vp[1], odev.v.vp[2],
364 <                odev.v.vp[0] + odev.v.vdir[0],
365 <                odev.v.vp[1] + odev.v.vdir[1],
366 <                odev.v.vp[2] + odev.v.vdir[2],
367 <                odev.v.vup[0], odev.v.vup[1], odev.v.vup[2]);
368 <                                        /* set viewport */
369 <        glViewport(0, 0, odev.hres, odev.vres);
370 <                                        /* initialize cone for Vornoi polys */
371 <        coneh = CONEH*(zmax - zmin);
372 <        crad = 0.5 * CONEW * 0.5*(xmax-xmin + ymax-ymin) * (zmin+coneh)/zmin;
373 <        vy[0] = vy[1] = vy[2] = 0.;
374 <        for (i = 0; i < 3; i++)
375 <                if (odev.v.vdir[i] < 0.6 && odev.v.vdir[i] > -0.6)
376 <                        break;
377 <        vy[i] = 1.;
378 <        fcross(vx, vy, odev.v.vdir);
379 <        normalize(vx);
380 <        fcross(vy, odev.v.vdir, vx);
381 <        for (j = 0, d = 0.; j < NCONEV; j++, d += 2.*PI/NCONEV) {
382 <                cx = crad*cos(d); sx = crad*sin(d);
383 <                for (i = 0; i < 3; i++)
384 <                        conev[j][i] = coneh*odev.v.vdir[i] +
385 <                                        cx*vx[i] + sx*vy[i];
386 <        }
335 > #endif
336 >        sprintf(errmsg, "%s: unknown command", cmd);
337 >        error(COMMAND, errmsg);
338   }
339  
340  
341 < #define SUCCSTEP        8       /* skip step when successful */
342 < #define MAXSTEP         64
343 <
344 < static
394 < drawvalue(li)                   /* draw a pixel value as a cone */
395 < register int    li;
341 > VIEW *
342 > dev_auxview(n, hvres)           /* return nth auxiliary view */
343 > int     n;
344 > int     hvres[2];
345   {
346 <        static int      skipstep = 1;
347 <        static FVECT    disp;
348 <        FVECT   apex;
349 <        double  d, dorg, dnew, h, v;
350 <        register int    i;
351 <                                /* check for quicker update */
352 <        if (quicken) {
353 <                if (li % skipstep)
405 <                        return;
406 <                if (skipstep < MAXSTEP)
407 <                        skipstep++;
408 <        }
409 <                                /* compute cone coordinates */
410 <        disp[0] = rV.wp[li][0] - odev.v.vp[0];
411 <        disp[1] = rV.wp[li][1] - odev.v.vp[1];
412 <        disp[2] = rV.wp[li][2] - odev.v.vp[2];
413 <        dorg = DOT(disp,odev.v.vdir);
414 <        if (dorg <= odev.v.vfore)
415 <                return;         /* clipped too near */
416 <        if (odev.v.vaft > FTINY && dorg > odev.v.vaft)
417 <                return;         /* clipped too far */
418 <        if (dorg > 1e5) {       /* background pixel */
419 <                dnew = maxdepth;
420 <                d = dnew/dorg;
421 <                dorg = maxdepth;
422 <        } else {                /* foreground pixel, compute penalty */
423 <                normalize(disp);
424 <                d = dnew = dorg + coneh*fdir2diff(rV.wd[li],disp)*DIRPEN;
425 <        }
426 <                                /* compute adjusted apex position */
427 <        disp[0] *= d; disp[1] *= d; disp[2] *= d;
428 <        apex[0] = odev.v.vp[0] + disp[0];
429 <        apex[1] = odev.v.vp[1] + disp[1];
430 <        apex[2] = odev.v.vp[2] + disp[2];
431 <                                /* compute view position and base offset */
432 <        h = DOT(disp,odev.v.hvec)/(dnew*odev.v.hn2);
433 <        v = DOT(disp,odev.v.vvec)/(dnew*odev.v.vn2);
434 <        if (fabs(h - odev.v.hoff) > 0.5 || fabs(v - odev.v.voff) > 0.5)
435 <                return;         /* clipped off screen */
436 <        if (dorg < mindepth)
437 <                mindepth = dorg;
438 <        if (dorg > maxdepth)
439 <                maxdepth = dorg;
440 <        for (i = 0; i < 3; i++)
441 <                disp[i] = apex[i] + coneh*(h*odev.v.hvec[i] + v*odev.v.vvec[i]);
442 <                                /* draw cone (pyramid approx.) */
443 <        glColor4ub(rV.rgb[li][0], rV.rgb[li][1], rV.rgb[li][2], rV.alpha[li]);
444 <        glBegin(GL_TRIANGLE_FAN);
445 <        glVertex3f(apex[0], apex[1], apex[2]);
446 <        for (i = 0; i < NCONEV; i++)
447 <                glVertex3d(conev[i][0] + disp[0], conev[i][1] + disp[1],
448 <                                conev[i][2] + disp[2]);
449 <                                /* connect last face to first */
450 <        glVertex3d(conev[0][0] + disp[0], conev[0][1] + disp[1],
451 <                        conev[0][2] + disp[2]);
452 <        glEnd();                /* done */
453 <        skipstep = SUCCSTEP;
346 >        hvres[0] = odev.hres; hvres[1] = odev.vres;
347 >        if (n == 0)
348 >                return(&odev.v);
349 > #ifdef STEREO
350 >        if (n == 1)
351 >                return(&vwright);
352 > #endif
353 >        return(NULL);
354   }
355  
456 #undef SUCCSTEP
457 #undef MAXSTEP
356  
357 <
358 < #define LEAFSIZ         (3*sizeof(float)+sizeof(int4)+\
461 <                        sizeof(TMbright)+7*sizeof(BYTE))
462 <
463 < static
464 < AllocValues(n)                  /* allocate space for n values */
465 < register int    n;
357 > int
358 > dev_input()                     /* get X11 input */
359   {
360 <        unsigned        nbytes;
468 <        register unsigned       i;
360 >        inpresflags = 0;
361  
362 <        if (n <= 0)
363 <                return(0);
472 <        if (rV.nl >= n)
473 <                return(rV.nl);
474 <        else if (rV.nl > 0)
475 <                free(rV.base);
476 <                                /* round space up to nearest power of 2 */
477 <        nbytes = n*LEAFSIZ + 8;
478 <        for (i = 1024; nbytes > i; i <<= 1)
479 <                ;
480 <        n = (i - 8) / LEAFSIZ;  /* should we make sure n is even? */
481 <        rV.base = (char *)malloc(n*LEAFSIZ);
482 <        if (rV.base == NULL)
483 <                return(0);
484 <                                /* assign larger alignment types earlier */
485 <        rV.wp = (float (*)[3])rV.base;
486 <        rV.wd = (int4 *)(rV.wp + n);
487 <        rV.brt = (TMbright *)(rV.wd + n);
488 <        rV.chr = (BYTE (*)[3])(rV.brt + n);
489 <        rV.rgb = (BYTE (*)[3])(rV.chr + n);
490 <        rV.alpha = (BYTE *)(rV.rgb + n);
491 <        rV.nl = n;
492 <        rV.drl = rV.tml = rV.bl = rV.tl = 0;
493 <        return(n);
494 < }
362 >        do
363 >                getevent();
364  
365 < #undef  LEAFSIZ
365 >        while (XPending(ourdisplay) > 0);
366  
367 +        odev.inpready = 0;
368  
369 < static
500 < FreeValues()                    /* free our allocated values */
501 < {
502 <        if (rV.nl <= 0)
503 <                return;
504 <        free(rV.base);
505 <        rV.base = NULL;
506 <        rV.nl = 0;
369 >        return(inpresflags);
370   }
371  
372  
373 < static
374 < clralphas()                     /* prepare for new alpha values */
373 > dev_value(c, d, p)              /* add a pixel value to our mesh */
374 > COLR    c;
375 > FVECT   d, p;
376   {
377 <        if (!vmaplen)
377 >        double  depth;
378 > #ifdef DOBJ
379 >        if (dobj_lightsamp != NULL) {   /* in light source sampling */
380 >                (*dobj_lightsamp)(c, d, p);
381                  return;
382 <        free((char *)valmap);
383 <        valmap = NULL;
384 <        vmaplen = 0;
382 >        }
383 > #endif
384 >        if (p != NULL) {                /* add depth to our range */
385 >                depth = (p[0] - odev.v.vp[0])*d[0] +
386 >                        (p[1] - odev.v.vp[1])*d[1] +
387 >                        (p[2] - odev.v.vp[2])*d[2];
388 >                if (depth > FTINY) {
389 >                        if (depth < mindpth)
390 >                                mindpth = depth;
391 >                        if (depth > maxdpth)
392 >                                maxdpth = depth;
393 >                }
394 >        }
395 >        smNewSamp(c, d, p);             /* add to display representation */
396 >        if (!--rayqleft)
397 >                dev_flush();            /* flush output */
398   }
399  
400  
401 < static int
402 < valcmp(v1p, v2p)                /* compare two pixel values */
523 < int     *v1p, *v2p;
401 > int
402 > dev_flush()                     /* flush output */
403   {
404 <        register int    v1 = *v1p, v2 = *v2p;
405 <        register int    c;
406 <
407 <        if ((c = rV.rgb[v1][0] - rV.rgb[v2][0])) return(c);
408 <        if ((c = rV.rgb[v1][1] - rV.rgb[v2][1])) return(c);
409 <        if ((c = rV.rgb[v1][2] - rV.rgb[v2][2])) return(c);
410 <        return(rV.alpha[v1] - rV.alpha[v2]);
404 >        if (mapped) {
405 > #ifdef STEREO
406 >                pushright();                    /* update right eye */
407 >                smUpdate(&vwright, MAXQUALITY);
408 > #ifdef DOBJ
409 >                dobj_render();                  /* usually in foreground */
410 > #endif
411 >                popright();                     /* update left eye */
412 > #endif
413 >                smUpdate(&odev.v, MAXQUALITY);
414 >                checkglerr("rendering mesh");
415 > #ifdef DOBJ
416 >                dobj_render();
417 > #endif
418 >                glFlush();                      /* flush OGL */
419 >        }
420 >        rayqleft = RAYQLEN;
421 >                                        /* flush X11 and return # pending */
422 >        return(odev.inpready = XPending(ourdisplay));
423   }
424  
425  
426 < static
427 < mergalphas(adest, al1, n1, al2, n2)     /* merge two sorted alpha lists */
537 < register int    *adest, *al1, *al2;
538 < int     n1, n2;
426 > checkglerr(where)               /* check for GL or GLU error */
427 > char    *where;
428   {
429 <        register int    cmp;
429 >        register GLenum errcode;
430  
431 <        while (n1 | n2) {
432 <                if (!n1) cmp = 1;
433 <                else if (!n2) cmp = -1;
434 <                else cmp = valcmp(al1, al2);
546 <                if (cmp > 0) {
547 <                        *adest++ = *al2++;
548 <                        n2--;
549 <                } else {
550 <                        *adest++ = *al1++;
551 <                        n1--;
552 <                }
431 >        while ((errcode = glGetError()) != GL_NO_ERROR) {
432 >                sprintf(errmsg, "OpenGL error %s: %s",
433 >                                where, gluErrorString(errcode));
434 >                error(WARNING, errmsg);
435          }
436   }
437  
438  
439 + #ifdef STEREO
440   static
441 < setalphas(vbeg, nvals)          /* add values to our map and set alphas */
559 < int     vbeg, nvals;
441 > pushright()                     /* push on right view */
442   {
443 <        register int    *newmap;
562 <        short   ccmp[3], lastalpha;
563 <        int     newmaplen;
443 >        double  d;
444  
445 <        if (nvals <= 0)
446 <                return;
447 <        newmaplen = vmaplen + nvals;    /* allocate new map */
448 <        newmap = (int *)malloc(newmaplen*sizeof(int));
449 <        if (newmap == NULL)
450 <                error(SYSTEM, "out of memory in setalphas");
571 <        while (nvals--) {               /* add new values to end */
572 <                rV.alpha[vbeg] = 255;
573 <                newmap[vmaplen+nvals] = vbeg++;
574 <        }
575 <        if (nvals >= 3*vmaplen) {       /* resort the combined array */
576 <                while (vmaplen--)
577 <                        newmap[vmaplen] = valmap[vmaplen];
578 <                qsort((char *)newmap, newmaplen, sizeof(int), valcmp);
579 <        } else {                        /* perform merge sort */
580 <                qsort((char *)(newmap+vmaplen), newmaplen-vmaplen,
581 <                                sizeof(int), valcmp);
582 <                mergalphas(newmap, valmap, vmaplen,
583 <                                newmap+vmaplen, newmaplen-vmaplen);
584 <        }
585 <        if (valmap != NULL)             /* free old map and assign new one */
586 <                free((char *)valmap);
587 <        valmap = newmap;
588 <        vmaplen = newmaplen;
589 <        lastalpha = 0;                  /* set new alpha values */
590 <        ccmp[0] = ccmp[1] = ccmp[2] = 256;
591 <        while (newmaplen--)
592 <                if (rV.rgb[*newmap][0] == ccmp[0] &&
593 <                                rV.rgb[*newmap][1] == ccmp[1] &&
594 <                                rV.rgb[*newmap][2] == ccmp[2]) {
595 <                        if (lastalpha >= 255)
596 <                                newmap++;
597 <                        else if (rV.alpha[*newmap] < 255)
598 <                                lastalpha = rV.alpha[*newmap++];
599 <                        else
600 <                                rV.alpha[*newmap++] = ++lastalpha;
601 <                } else {
602 <                        ccmp[0] = rV.rgb[*newmap][0];
603 <                        ccmp[1] = rV.rgb[*newmap][1];
604 <                        ccmp[2] = rV.rgb[*newmap][2];
605 <                        if (rV.alpha[*newmap] < 255)
606 <                                lastalpha = rV.alpha[*newmap++];
607 <                        else
608 <                                rV.alpha[*newmap++] = lastalpha = 1;
609 <                }
445 >        setstereobuf(STEREO_BUFFER_RIGHT);
446 >        glMatrixMode(GL_MODELVIEW);
447 >        glPushMatrix();
448 >        d = -eyesepdist / sqrt(odev.v.hn2);
449 >        glTranslated(d*odev.v.hvec[0], d*odev.v.hvec[1], d*odev.v.hvec[2]);
450 >        checkglerr("setting right view");
451   }
452  
453  
454   static
455 < TMapValues(redo)                /* map our values to RGB */
615 < int     redo;
455 > popright()                      /* pop off right view */
456   {
457 <        int     aorg, alen, borg, blen;
458 <                                        /* recompute mapping? */
459 <        if (redo)
620 <                rV.tml = rV.bl;
621 <                                        /* already done? */
622 <        if (rV.tml == rV.tl)
623 <                return(1);
624 <                                        /* compute segments */
625 <        aorg = rV.tml;
626 <        if (rV.tl >= aorg) {
627 <                alen = rV.tl - aorg;
628 <                blen = 0;
629 <        } else {
630 <                alen = rV.nl - aorg;
631 <                borg = 0;
632 <                blen = rV.tl;
633 <        }
634 <                                        /* (re)compute tone mapping? */
635 <        if (rV.tml == rV.bl) {
636 <                tmClearHisto();
637 <                tmAddHisto(rV.brt+aorg, alen, 1);
638 <                if (blen > 0)
639 <                        tmAddHisto(rV.brt+borg, blen, 1);
640 <                if (tmComputeMapping(0., 0., 0.) != TM_E_OK)
641 <                        return(0);
642 <                clralphas();            /* restart value list */
643 <                rV.drl = rV.bl;         /* need to redraw */
644 <        }
645 <        if (tmMapPixels(rV.rgb+aorg, rV.brt+aorg,
646 <                        rV.chr+aorg, alen) != TM_E_OK)
647 <                return(0);
648 <        if (blen > 0)
649 <                tmMapPixels(rV.rgb+borg, rV.brt+borg,
650 <                                rV.chr+borg, blen);
651 <        setalphas(aorg, alen);          /* compute add'l alpha values */
652 <        if (blen > 0)
653 <                setalphas(borg, blen);
654 <        rV.tml = rV.tl;                 /* we're all up to date */
655 <        return(1);
457 >        glMatrixMode(GL_MODELVIEW);
458 >        glPopMatrix();
459 >        setstereobuf(STEREO_BUFFER_LEFT);
460   }
461 + #endif
462  
463  
464   static int
465 < Compost(pct)                    /* free up some values */
661 < int     pct;
465 > mytmflags()                     /* figure out tone mapping flags */
466   {
467 <        int     nused, nclear, nmapped, ndrawn;
468 <                                /* figure out how many values to clear */
469 <        nclear = rV.nl * pct / 100;
470 <        nused = rV.tl - rV.bl;
471 <        if (nused <= 0) nused += rV.nl;
472 <        nclear -= rV.nl - nused;
473 <        if (nclear <= 0)
474 <                return(0);
475 <        if (nclear >= nused) {  /* clear them all? */
476 <                rV.drl = rV.tml = rV.bl = rV.tl = 0;
477 <                return(nused);
478 <        }
675 <                                /* else clear values from bottom */
676 <        ndrawn = rV.drl - rV.bl;
677 <        if (ndrawn < 0) ndrawn += rV.nl;
678 <        nmapped = rV.tml - rV.bl;
679 <        if (nmapped < 0) nmapped += rV.nl;
680 <        rV.bl += nclear;
681 <        if (rV.bl >= rV.nl) rV.bl -= rV.nl;
682 <        if (ndrawn < nclear) rV.drl = rV.bl;
683 <        if (nmapped < nclear) rV.tml = rV.bl;
684 <        return(nclear);
467 >        extern char     *progname;
468 >        register char   *cp, *tail;
469 >                                        /* find basic name */
470 >        for (cp = tail = progname; *cp; cp++)
471 >                if (*cp == '/')
472 >                        tail = cp+1;
473 >        for (cp = tail; *cp && *cp != '.'; cp++)
474 >                ;
475 >        if (cp > tail && cp[-1] == 'h')
476 >                return(TM_F_HUMAN|TM_F_NOSTDERR);
477 >        else
478 >                return(TM_F_CAMERA|TM_F_NOSTDERR);
479   }
480  
481  
# Line 706 | Line 500 | getevent()                     /* get next event */
500                  getkey(levptr(XKeyPressedEvent));
501                  break;
502          case ButtonPress:
503 <                getmove(levptr(XButtonPressedEvent));
503 >                if (FRAMESTATE(levptr(XButtonPressedEvent)->state))
504 >                        getframe(levptr(XButtonPressedEvent));
505 >                else
506 >                        getmove(levptr(XButtonPressedEvent));
507                  break;
508          }
509   }
# Line 722 | Line 519 | register FVECT wp[2];
519  
520  
521   static
522 < draw_grids()                    /* draw holodeck section grids */
522 > draw_grids(fore)                /* draw holodeck section grids */
523 > int     fore;
524   {
525 <        static BYTE     gridrgba[4] = {0x0, 0xff, 0xff, 0x00};
526 <
527 <        if (!mapped)
528 <                return;
529 <        glColor4ub(gridrgba[0], gridrgba[1], gridrgba[2], gridrgba[3]);
732 <                                        /* draw each grid line */
733 <        glBegin(GL_LINES);
525 >        if (fore)
526 >                glColor4ub(0, 255, 255, 0);
527 >        else
528 >                glColor4ub(0, 0, 0, 0);
529 >        glBegin(GL_LINES);              /* draw each grid line */
530          gridlines(draw3dline);
531          glEnd();
532 +        checkglerr("drawing grid lines");
533   }
534  
535 <
536 < static int
740 < IndexValue(rgba)                /* locate a pixel by it's framebuffer value */
741 < register BYTE   rgba[4];
742 < {
743 <        register int    *vp;
744 <                                        /* check legality */
745 <        if (rgba[3] == 0 || rgba[3] == 255)
746 <                return(-1);
747 <                                        /* borrow a value slot */
748 <        rV.rgb[rV.tl][0] = rgba[0];
749 <        rV.rgb[rV.tl][1] = rgba[1];
750 <        rV.rgb[rV.tl][2] = rgba[2];
751 <        rV.alpha[rV.tl] = rgba[3];
752 <                                        /* find it */
753 <        vp = (int *)bsearch((char *)&rV.tl, (char *)valmap, vmaplen,
754 <                        sizeof(int), valcmp);
755 <        if (vp == NULL)
756 <                return(-1);
757 <        return(*vp);
758 < }
759 <
760 <
761 < static int
762 < FindValue(dx, dy)               /* find a value on the display */
535 > static double
536 > getdistance(dx, dy, direc)      /* distance from fore plane along view ray */
537   int     dx, dy;
538 + FVECT   direc;
539   {
540 <        BYTE    rgba[4];
540 >        GLfloat gldepth;
541 >        double  dist;
542  
543 <        if (dx < 0 || dy < 0 || dx >= odev.hres || dy >= odev.vres)
544 <                return(-1);
545 <        glReadPixels(dx, dy, 1, 1, GL_RGBA, GL_UNSIGNED_BYTE, rgba);
546 <        return(IndexValue(rgba));
543 >        if (dx<0 | dx>=odev.hres | dy<0 | dy>=odev.vres)
544 >                return(FHUGE);
545 >        glReadPixels(dx,dy, 1,1, GL_DEPTH_COMPONENT,GL_FLOAT, &gldepth);
546 >        dist = mapdepth(gldepth);
547 >        if (dist >= .99*FHUGE)
548 >                return(FHUGE);
549 >        return((dist-odev.v.vfore)/DOT(direc,odev.v.vdir));
550   }
551  
552  
# Line 776 | Line 555 | moveview(dx, dy, mov, orb)     /* move our view */
555   int     dx, dy, mov, orb;
556   {
557          VIEW    nv;
558 <        FVECT   odir, v1;
559 <        double  d;
558 >        FVECT   odir, v1, wip;
559 >        double  d,d1;
560          register int    li;
561                                  /* start with old view */
562 <        copystruct(&nv, &odev.v);
563 <                                /* change view direction */
564 <        if (mov | orb) {
565 <                if ((li = FindValue(dx, dy)) < 0)
566 <                        return(0);      /* not on window */
567 <                VSUM(odir, rV.wp[li], nv.vp, -1.);
568 <        } else {
569 <                if (viewray(nv.vp, nv.vdir, &odev.v,
570 <                                (dx+.5)/odev.hres, (dy+.5)/odev.vres) < -FTINY)
571 <                        return(0);      /* outside view */
572 <        }
562 >        nv = odev.v;
563 >                                /* orient our motion */
564 >        if (viewray(v1, odir, &odev.v,
565 >                        (dx+.5)/odev.hres, (dy+.5)/odev.vres) < -FTINY)
566 >                return(0);              /* outside view */
567 >        if (mov | orb) {        /* moving relative to geometry */
568 >                d = getdistance(dx, dy, odir);/*distance from front plane */
569 > #ifdef DOBJ
570 >                d1 = dobj_trace(NULL, v1, odir);        /* check objects */
571 >                                                /* check holodeck */
572 >                if (d1 < d)
573 >                        d = d1;
574 > #endif
575 >                if (d >= .99*FHUGE)
576 >                        d = 0.5*(dev_zmax+dev_zmin);    /* just guess */
577 >                VSUM(wip, v1, odir, d);
578 >                VSUB(odir, wip, odev.v.vp);
579 > #if 0
580 >                fprintf(stderr, "moveview: hit %s at (%f,%f,%f) (t=%f)\n",
581 >                                li < 0 ? "object" : "mesh",
582 >                                wip[0], wip[1], wip[2],
583 >                                (wip[0]-odev.v.vp[0])*odir[0] +
584 >                                (wip[1]-odev.v.vp[1])*odir[1] +
585 >                                (wip[2]-odev.v.vp[2])*odir[2]);
586 > #endif
587 >        } else                  /* panning with constant viewpoint */
588 >                VCOPY(nv.vdir, odir);
589          if (orb && mov) {               /* orbit left/right */
590                  spinvector(odir, odir, nv.vup, d=MOVDEG*PI/180.*mov);
591 <                VSUM(nv.vp, rV.wp[li], odir, -1.);
591 >                VSUM(nv.vp, wip, odir, -1.);
592                  spinvector(nv.vdir, nv.vdir, nv.vup, d);
593          } else if (orb) {               /* orbit up/down */
594                  fcross(v1, odir, nv.vup);
595                  if (normalize(v1) == 0.)
596                          return(0);
597                  spinvector(odir, odir, v1, d=MOVDEG*PI/180.*orb);
598 <                VSUM(nv.vp, rV.wp[li], odir, -1.);
598 >                VSUM(nv.vp, wip, odir, -1.);
599                  spinvector(nv.vdir, nv.vdir, v1, d);
600          } else if (mov) {               /* move forward/backward */
601                  d = MOVPCT/100. * mov;
# Line 808 | Line 603 | int    dx, dy, mov, orb;
603          }
604          if (!mov ^ !orb && headlocked) {        /* restore head height */
605                  VSUM(v1, odev.v.vp, nv.vp, -1.);
606 <                d = DOT(v1, odev.v.vup);
606 >                d = DOT(v1, nv.vup);
607                  VSUM(nv.vp, nv.vp, odev.v.vup, d);
608          }
609          if (setview(&nv) != NULL)
# Line 820 | Line 615 | int    dx, dy, mov, orb;
615  
616  
617   static
618 + getframe(ebut)                          /* get focus frame */
619 + XButtonPressedEvent     *ebut;
620 + {
621 +        int     startx = ebut->x, starty = ebut->y;
622 +        int     endx, endy;
623 +
624 +        XMaskEvent(ourdisplay, ButtonReleaseMask, levptr(XEvent));
625 +        endx = levptr(XButtonReleasedEvent)->x;
626 +        endy = levptr(XButtonReleasedEvent)->y;
627 +        if (endx == startx | endy == starty) {
628 +                XBell(ourdisplay, 0);
629 +                return;
630 +        }
631 +        if (endx < startx) {register int c = endx; endx = startx; startx = c;}
632 +        if (endy < starty) {register int c = endy; endy = starty; starty = c;}
633 +        sprintf(odev_args, "%.3f %.3f %.3f %.3f",
634 +                        (startx+.5)/odev.hres, 1.-(endy+.5)/odev.vres,
635 +                        (endx+.5)/odev.hres, 1.-(starty+.5)/odev.vres);
636 +        inpresflags |= DFL(DC_FOCUS);
637 + }
638 +
639 +
640 + static
641   getmove(ebut)                           /* get view change */
642   XButtonPressedEvent     *ebut;
643   {
644          int     movdir = MOVDIR(ebut->button);
645          int     movorb = MOVORB(ebut->state);
646 +        int     qlevel = MAXQUALITY-1;
647 +        time_t  lasttime, thistime;
648 +        int     nframes;
649          Window  rootw, childw;
650          int     rootx, rooty, wx, wy;
651          unsigned int    statemask;
652  
653 <        quicken = 1;                    /* accelerate update rate */
833 <        XNoOp(ourdisplay);
653 >        XNoOp(ourdisplay);              /* makes sure we're not idle */
654  
655 +        nframes = 0;
656          while (!XCheckMaskEvent(ourdisplay,
657                          ButtonReleaseMask, levptr(XEvent))) {
658 <
658 >                                        /* get cursor position */
659                  if (!XQueryPointer(ourdisplay, gwind, &rootw, &childw,
660                                  &rootx, &rooty, &wx, &wy, &statemask))
661                          break;          /* on another screen */
662  
663 +                draw_grids(0);          /* clear old grid lines */
664 + #ifdef STEREO
665 +                pushright(); draw_grids(0); popright();
666 + #endif
667 +                                        /* compute view motion */
668                  if (!moveview(wx, odev.vres-1-wy, movdir, movorb)) {
669                          sleep(1);
670 +                        lasttime++;
671 +                        continue;       /* cursor in bad place */
672 +                }
673 +                draw_grids(1);          /* redraw grid */
674 + #ifdef STEREO
675 +                pushright();
676 +                draw_grids(1);
677 +                smUpdate(&vwright, qlevel);
678 + #ifdef DOBJ
679 +                dobj_render();
680 + #endif
681 +                popright();
682 + #endif
683 +                                        /* redraw mesh */
684 +                smUpdate(&odev.v, qlevel);
685 + #ifdef DOBJ
686 +                dobj_render();          /* redraw object */
687 + #endif
688 +                glFlush();
689 +                                        /* figure out good quality level */
690 +                if (!nframes++) {               /* ignore first frame */
691 +                        lasttime = time(0);
692                          continue;
693                  }
694 <                draw_grids();
695 <                dev_flush();
694 >                thistime = time(0);
695 >                if (thistime - lasttime >= 6 ||
696 >                                nframes > (int)(3*3*TARGETFPS)) {
697 >                        qlevel = thistime<=lasttime ? 3*MAXQUALITY :
698 >                                (int)((double)nframes/(thistime-lasttime)
699 >                                        / TARGETFPS * qlevel + 0.5);
700 >                        nframes = 0;
701 >                        if (qlevel >= MAXQUALITY) {
702 >                                if (qlevel >= 3*MAXQUALITY) {   /* brakes!! */
703 >                                        sleep(1);
704 >                                        lasttime++;
705 >                                }
706 >                                qlevel = MAXQUALITY-1;
707 >                        } else if (qlevel < 1)
708 >                                qlevel = 1;
709 >                }
710          }
711          if (!(inpresflags & DFL(DC_SETVIEW))) { /* do final motion */
712                  movdir = MOVDIR(levptr(XButtonReleasedEvent)->button);
713                  wx = levptr(XButtonReleasedEvent)->x;
714                  wy = levptr(XButtonReleasedEvent)->y;
715                  moveview(wx, odev.vres-1-wy, movdir, movorb);
716 +        }
717 + }
718 +
719 +
720 + static
721 + setglpersp(vp)                  /* set perspective view in GL */
722 + register VIEW   *vp;
723 + {
724 +        double  d, xmin, xmax, ymin, ymax;
725 +
726 +        if (mindpth >= maxdpth) {
727 +                dev_zmin = 1.;
728 +                dev_zmax = 100.;
729          } else {
730 <                glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
731 <                redraw();
730 >                dev_zmin = 0.1*mindpth;
731 >                dev_zmax = 5.0*maxdpth;
732          }
733 <        quicken = 0;
734 <        dev_flush();
733 >        if (odev.v.vfore > FTINY)
734 >                dev_zmin = odev.v.vfore;
735 >        if (odev.v.vaft > FTINY)
736 >                dev_zmax = odev.v.vaft;
737 >        if (dev_zmin*500. < dev_zmax)
738 >                dev_zmax = dev_zmin*500.;
739 >        xmax = dev_zmin * tan(PI/180./2. * odev.v.horiz);
740 >        xmin = -xmax;
741 >        d = odev.v.hoff * (xmax - xmin);
742 >        xmin += d; xmax += d;
743 >        ymax = dev_zmin * tan(PI/180./2. * odev.v.vert);
744 >        ymin = -ymax;
745 >        d = odev.v.voff * (ymax - ymin);
746 >        ymin += d; ymax += d;
747 >                                        /* set view matrix */
748 >        glMatrixMode(GL_PROJECTION);
749 >        glLoadIdentity();
750 >        glFrustum(xmin, xmax, ymin, ymax, dev_zmin, dev_zmax);
751 >        gluLookAt(odev.v.vp[0], odev.v.vp[1], odev.v.vp[2],
752 >                odev.v.vp[0] + odev.v.vdir[0],
753 >                odev.v.vp[1] + odev.v.vdir[1],
754 >                odev.v.vp[2] + odev.v.vdir[2],
755 >                odev.v.vup[0], odev.v.vup[1], odev.v.vup[2]);
756 >        checkglerr("setting perspective view");
757   }
758  
759  
760   static
761 + wipeclean(tmflag)                       /* prepare for redraw */
762 + int     tmflag;
763 + {
764 +                                        /* clear depth buffer */
765 + #ifdef STEREO
766 +        setstereobuf(STEREO_BUFFER_RIGHT);
767 +        glClear(GL_DEPTH_BUFFER_BIT);
768 +        setstereobuf(STEREO_BUFFER_LEFT);
769 + #endif
770 +        glClear(GL_DEPTH_BUFFER_BIT);
771 +        smClean(tmflag);                /* reset drawing routines */
772 +        setglpersp(&odev.v);            /* reset view & clipping planes */
773 + }
774 +
775 +
776 + static
777   getkey(ekey)                            /* get input key */
778   register XKeyPressedEvent  *ekey;
779   {
780 +        Window  rootw, childw;
781 +        int     rootx, rooty, wx, wy;
782 +        unsigned int    statemask;
783          int  n;
784          char    buf[8];
785  
# Line 886 | Line 802 | register XKeyPressedEvent  *ekey;
802          case 'v':                       /* spit out view */
803                  inpresflags |= DFL(DC_GETVIEW);
804                  return;
805 +        case 'f':                       /* frame view position */
806 +                if (!XQueryPointer(ourdisplay, gwind, &rootw, &childw,
807 +                                &rootx, &rooty, &wx, &wy, &statemask))
808 +                        return;         /* on another screen */
809 +                sprintf(odev_args, "%.4f %.4f", (wx+.5)/odev.hres,
810 +                                1.-(wy+.5)/odev.vres);
811 +                inpresflags |= DFL(DC_FOCUS);
812 +                return;
813 +        case 'F':                       /* unfocus */
814 +                odev_args[0] = '\0';
815 +                inpresflags |= DFL(DC_FOCUS);
816 +                return;
817          case '\n':
818          case '\r':                      /* resume computation */
819                  inpresflags |= DFL(DC_RESUME);
820                  return;
821          case CTRL('R'):                 /* redraw screen */
822 <                TMapValues(1);
895 <                glClear(GL_DEPTH_BUFFER_BIT);
896 <                redraw();
822 >                wipeclean(1);
823                  return;
824          case CTRL('L'):                 /* refresh from server */
825                  if (inpresflags & DFL(DC_REDRAW))
826                          return;
827 +                setglpersp(&odev.v);            /* reset clipping planes */
828                  glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
829 <                glDisable(GL_DEPTH_TEST);
830 <                draw_grids();
829 >                glDisable(GL_DEPTH_TEST);       /* so grids will clear */
830 >                draw_grids(1);
831 > #ifdef STEREO
832 >                pushright();
833 >                glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
834 >                draw_grids(1);
835 >                popright();
836 > #endif
837                  glEnable(GL_DEPTH_TEST);
838                  glFlush();
839 <                Compost(100);                   /* get rid of old values */
839 >                smInit(rsL.max_samp);           /* get rid of old values */
840                  inpresflags |= DFL(DC_REDRAW);  /* resend values from server */
841 +                rayqleft = 0;                   /* hold off update */
842                  return;
843          case 'K':                       /* kill rtrace process(es) */
844                  inpresflags |= DFL(DC_KILL);
# Line 929 | Line 863 | static
863   fixwindow(eexp)                         /* repair damage to window */
864   register XExposeEvent  *eexp;
865   {
866 <        if (odev.hres == 0 || odev.vres == 0)   /* first exposure */
866 >        if (odev.hres == 0 | odev.vres == 0) {  /* first exposure */
867                  resizewindow((XConfigureEvent *)eexp);
868 <        if (eexp->width == odev.hres && eexp->height == odev.vres)
935 <                TMapValues(1);
936 <        if (!eexp->count) {
937 <                glClear(GL_DEPTH_BUFFER_BIT);
938 <                redraw();
868 >                return;
869          }
870 +        if (eexp->count)                /* wait for final exposure */
871 +                return;
872 +        wipeclean(0);                   /* clear depth */
873   }
874  
875  
# Line 944 | Line 877 | static
877   resizewindow(ersz)                      /* resize window */
878   register XConfigureEvent  *ersz;
879   {
880 +        glViewport(0, 0, ersz->width, ersz->height);
881 +
882          if (ersz->width == odev.hres && ersz->height == odev.vres)
883                  return;
884  
# Line 952 | Line 887 | register XConfigureEvent  *ersz;
887  
888          odev.v.horiz = 2.*180./PI * atan(0.5/VIEWDIST*pwidth*odev.hres);
889          odev.v.vert = 2.*180./PI * atan(0.5/VIEWDIST*pheight*odev.vres);
955
956        setGLview();
890  
891          inpresflags |= DFL(DC_SETVIEW);
892   }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines