ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/hd/rhd_ogl.c
Revision: 3.1
Committed: Fri Dec 18 11:55:19 1998 UTC (25 years, 4 months ago) by gwlarson
Content type: text/plain
Branch: MAIN
Log Message:
Initial revision

File Contents

# User Rev Content
1 gwlarson 3.1 /* Copyright (c) 1998 Silicon Graphics, Inc. */
2    
3     #ifndef lint
4     static char SCCSid[] = "$SunId$ SGI";
5     #endif
6    
7     /*
8     * OpenGL driver for holodeck display.
9     * Based on GLX driver using T-mesh.
10     *
11     * Define symbol STEREO for stereo viewing.
12     * Define symbol DOBJ for display object viewing.
13     */
14    
15     #ifdef NOSTEREO
16     #ifdef STEREO
17     #undef STEREO
18     #else
19     #undef NOSTEREO
20     #endif
21     #endif
22    
23     #include "standard.h"
24    
25     #include <sys/types.h>
26     #include <GL/glx.h>
27     #include <GL/glu.h>
28     #ifdef STEREO
29     #include <X11/extensions/SGIStereo.h>
30     #endif
31    
32     #include "rhd_odraw.h"
33     #ifdef DOBJ
34     #include "rhdobj.h"
35     #endif
36    
37     #include "x11icon.h"
38    
39     #ifndef RAYQLEN
40     #define RAYQLEN 10240 /* max. rays to queue before flush */
41     #endif
42    
43     #ifndef FEQ
44     #define FEQ(a,b) ((a)-(b) <= FTINY && (a)-(b) >= -FTINY)
45     #endif
46    
47     #define GAMMA 1.4 /* default gamma correction */
48    
49     #define FRAMESTATE(s) (((s)&(ShiftMask|ControlMask))==(ShiftMask|ControlMask))
50    
51     #define MOVPCT 7 /* percent distance to move /frame */
52     #define MOVDIR(b) ((b)==Button1 ? 1 : (b)==Button2 ? 0 : -1)
53     #define MOVDEG (-5) /* degrees to orbit CW/down /frame */
54     #define MOVORB(s) ((s)&ShiftMask ? 1 : (s)&ControlMask ? -1 : 0)
55    
56     #define MINWIDTH 480 /* minimum graphics window width */
57     #define MINHEIGHT 400 /* minimum graphics window height */
58    
59     #define VIEWDIST 356 /* assumed viewing distance (mm) */
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|\
70     ButtonPressMask|ButtonReleaseMask)
71    
72     #define levptr(etype) ((etype *)&currentevent)
73    
74     struct driver odev; /* global device driver structure */
75    
76     char odev_args[64]; /* command arguments */
77    
78     static GLfloat *depthbuffer = NULL; /* depth buffer */
79    
80     #ifdef STEREO
81     static VIEW vwright; /* right eye view */
82     static GLfloat *depthright = NULL; /* right depth buffer */
83     #endif
84    
85     static int rayqleft = 0; /* rays left to queue before flush */
86    
87     static XEvent currentevent; /* current event */
88    
89     static int mapped = 0; /* window is mapped? */
90     static unsigned long ourblack=0, ourwhite=~0;
91    
92     static Display *ourdisplay = NULL; /* our display */
93     static XVisualInfo *ourvinf; /* our visual information */
94     static Window gwind = 0; /* our graphics window */
95     static GLXContext gctx; /* our GLX context */
96    
97     static double pwidth, pheight; /* pixel dimensions (mm) */
98    
99     static double dev_zmin, dev_zmax; /* fore and aft clipping plane dist. */
100     static double dev_zrat; /* (1. - dev_zmin/dev_zmax) */
101    
102     #define setzrat() (dev_zrat = 1. - dev_zmin/dev_zmax)
103     #define mapdepth(d) ((d)>0.9995 ? FHUGE : dev_zmin/(1.-(d)*dev_zrat))
104    
105     static int inpresflags; /* input result flags */
106    
107     static int headlocked; /* lock vertical motion */
108    
109     static int isperspective; /* perspective/ortho view */
110    
111     static int viewsteady; /* is view steady? */
112    
113     static int resizewindow(), getevent(), getkey(), moveview(), wipeclean(),
114     xferdepth(), setglortho(),
115     setglpersp(), getframe(), getmove(), fixwindow(), mytmflags();
116    
117     static double getdistance();
118    
119     #ifdef STEREO
120     static int pushright(), popright();
121     #endif
122    
123     extern time_t time();
124    
125    
126     dev_open(id) /* initialize GLX driver */
127     char *id;
128     {
129     extern char *getenv();
130     static RGBPRIMS myprims = STDPRIMS;
131     static int atlBest[] = {GLX_RGBA, GLX_RED_SIZE,8,
132     GLX_GREEN_SIZE,8, GLX_BLUE_SIZE,8,
133     GLX_DEPTH_SIZE,15, None};
134     char *ev;
135     double gamval = GAMMA;
136     RGBPRIMP dpri = stdprims;
137     XSetWindowAttributes ourwinattr;
138     XWMHints ourxwmhints;
139     XSizeHints oursizhints;
140     /* check for unsupported stereo */
141     #ifdef NOSTEREO
142     error(INTERNAL, "stereo display driver unavailable");
143     #endif
144     /* open display server */
145     ourdisplay = XOpenDisplay(NULL);
146     if (ourdisplay == NULL)
147     error(USER, "cannot open X-windows; DISPLAY variable set?\n");
148     #ifdef STEREO
149     switch (XSGIQueryStereoMode(ourdisplay, ourroot)) {
150     case STEREO_TOP:
151     case STEREO_BOTTOM:
152     break;
153     case STEREO_OFF:
154     error(USER,
155     "wrong video mode: run \"/usr/gfx/setmon -n STR_TOP\" first");
156     case X_STEREO_UNSUPPORTED:
157     error(USER, "stereo mode not supported on this screen");
158     default:
159     error(INTERNAL, "unknown stereo mode");
160     }
161     #endif
162     /* find a usable visual */
163     ourvinf = glXChooseVisual(ourdisplay, ourscreen, atlBest);
164     if (ourvinf == NULL)
165     error(USER, "no suitable visuals available");
166     /* get a context */
167     gctx = glXCreateContext(ourdisplay, ourvinf, NULL, GL_TRUE);
168     /* set gamma and tone mapping */
169     if ((ev = XGetDefault(ourdisplay, "radiance", "gamma")) != NULL
170     || (ev = getenv("DISPLAY_GAMMA")) != NULL)
171     gamval = atof(ev);
172     if ((ev = getenv("DISPLAY_PRIMARIES")) != NULL &&
173     sscanf(ev, "%f %f %f %f %f %f %f %f",
174     &myprims[RED][CIEX],&myprims[RED][CIEY],
175     &myprims[GRN][CIEX],&myprims[GRN][CIEY],
176     &myprims[BLU][CIEX],&myprims[BLU][CIEY],
177     &myprims[WHT][CIEX],&myprims[WHT][CIEY]) >= 6)
178     dpri = myprims;
179     if (tmInit(mytmflags(), dpri, gamval) == NULL)
180     error(SYSTEM, "not enough memory in dev_open");
181     /* open window */
182     ourwinattr.background_pixel = ourblack;
183     ourwinattr.border_pixel = ourblack;
184     ourwinattr.event_mask = ourmask;
185     /* this is stupid */
186     ourwinattr.colormap = XCreateColormap(ourdisplay, ourroot,
187     ourvinf->visual, AllocNone);
188     gwind = XCreateWindow(ourdisplay, ourroot, 0, 0,
189     DisplayWidth(ourdisplay,ourscreen)-2*BORWIDTH,
190     #ifdef STEREO
191     (DisplayHeight(ourdisplay,ourscreen)-2*BORWIDTH)/2,
192     #else
193     DisplayHeight(ourdisplay,ourscreen)-2*BORWIDTH,
194     #endif
195     BORWIDTH, ourvinf->depth, InputOutput, ourvinf->visual,
196     CWBackPixel|CWBorderPixel|CWColormap|CWEventMask, &ourwinattr);
197     if (gwind == 0)
198     error(SYSTEM, "cannot create window\n");
199     XStoreName(ourdisplay, gwind, id);
200     /* set window manager hints */
201     ourxwmhints.flags = InputHint|IconPixmapHint;
202     ourxwmhints.input = True;
203     ourxwmhints.icon_pixmap = XCreateBitmapFromData(ourdisplay,
204     gwind, x11icon_bits, x11icon_width, x11icon_height);
205     XSetWMHints(ourdisplay, gwind, &ourxwmhints);
206     oursizhints.min_width = MINWIDTH;
207     #ifdef STEREO
208     oursizhints.min_height = MINHEIGHT/2;
209     oursizhints.max_width = DisplayWidth(ourdisplay,ourscreen)-2*BORWIDTH;
210     oursizhints.max_height = (DisplayHeight(ourdisplay,ourscreen) -
211     2*BORWIDTH)/2;
212     oursizhints.flags = PMinSize|PMaxSize;
213     #else
214     oursizhints.min_height = MINHEIGHT;
215     oursizhints.flags = PMinSize;
216     #endif
217     XSetNormalHints(ourdisplay, gwind, &oursizhints);
218     /* set GLX context */
219     glXMakeCurrent(ourdisplay, gwind, gctx);
220     glEnable(GL_DEPTH_TEST);
221     glDepthFunc(GL_LEQUAL);
222     glDisable(GL_DITHER);
223     glFrontFace(GL_CCW);
224     glDisable(GL_CULL_FACE);
225     glLightModeli(GL_LIGHT_MODEL_TWO_SIDE, GL_TRUE);
226     glPolygonMode(GL_FRONT_AND_BACK, GL_FILL);
227     /* figure out sensible view */
228     pwidth = (double)DisplayWidthMM(ourdisplay, ourscreen) /
229     DisplayWidth(ourdisplay, ourscreen);
230     pheight = (double)DisplayHeightMM(ourdisplay, ourscreen) /
231     DisplayHeight(ourdisplay, ourscreen);
232     #ifdef STEREO
233     pheight *= 2.;
234     setstereobuf(STEREO_BUFFER_LEFT);
235     #endif
236     checkglerr("setting rendering parameters");
237     copystruct(&odev.v, &stdview);
238     odev.v.type = VT_PER;
239     headlocked = 0; /* free up head movement */
240     viewsteady = 1; /* view starts static */
241     isperspective = -1; /* but no view set, yet */
242     /* map the window */
243     XMapWindow(ourdisplay, gwind);
244     dev_input(); /* sets size and view angles */
245     if (!odInit(DisplayWidth(ourdisplay,ourscreen) *
246     DisplayHeight(ourdisplay,ourscreen) / 10))
247     error(SYSTEM, "insufficient memory for value storage");
248     odev.name = id;
249     odev.firstuse = 1; /* can't recycle samples */
250     odev.ifd = ConnectionNumber(ourdisplay);
251     }
252    
253    
254     dev_close() /* close our display and free resources */
255     {
256     #ifdef DOBJ
257     dobj_cleanup();
258     #endif
259     odDone();
260     glXMakeCurrent(ourdisplay, None, NULL);
261     glXDestroyContext(ourdisplay, gctx);
262     XDestroyWindow(ourdisplay, gwind);
263     gwind = 0;
264     XCloseDisplay(ourdisplay);
265     ourdisplay = NULL;
266     tmDone(NULL);
267     odev.v.type = 0;
268     odev.hres = odev.vres = 0;
269     odev.ifd = -1;
270     }
271    
272    
273     dev_clear() /* clear our representation */
274     {
275     wipeclean();
276     rayqleft = 0; /* hold off update */
277     }
278    
279    
280     int
281     dev_view(nv) /* assign new driver view */
282     register VIEW *nv;
283     {
284     double d;
285    
286     if (nv->type != VT_PER || /* check view legality */
287     nv->horiz > 160. || nv->vert > 160.) {
288     error(COMMAND, "illegal view type/angle");
289     nv->type = odev.v.type;
290     nv->horiz = odev.v.horiz;
291     nv->vert = odev.v.vert;
292     return(0);
293     }
294     if (nv != &odev.v) {
295     /* resize window? */
296     if (!FEQ(nv->horiz,odev.v.horiz) ||
297     !FEQ(nv->vert,odev.v.vert)) {
298     int dw = DisplayWidth(ourdisplay,ourscreen);
299     int dh = DisplayHeight(ourdisplay,ourscreen);
300    
301     dw -= 25; /* for window frame */
302     dh -= 50;
303     #ifdef STEREO
304     dh /= 2;
305     #endif
306     odev.hres = 2.*VIEWDIST/pwidth *
307     tan(PI/180./2.*nv->horiz);
308     odev.vres = 2.*VIEWDIST/pheight *
309     tan(PI/180./2.*nv->vert);
310     if (odev.hres > dw) {
311     odev.vres = dw * odev.vres / odev.hres;
312     odev.hres = dw;
313     }
314     if (odev.vres > dh) {
315     odev.hres = dh * odev.hres / odev.vres;
316     odev.vres = dh;
317     }
318     XResizeWindow(ourdisplay, gwind, odev.hres, odev.vres);
319     dev_input(); /* get resize event */
320     }
321     copystruct(&odev.v, nv); /* setview() already called */
322     #ifdef STEREO
323     copystruct(&vwright, nv);
324     d = eyesepdist / sqrt(nv->hn2);
325     VSUM(vwright.vp, nv->vp, nv->hvec, d);
326     /* setview(&vwright); -- Unnecessary */
327     #endif
328     } else
329     viewsteady = 1;
330     wipeclean();
331     return(1);
332     }
333    
334    
335     dev_section(ofn) /* add octree for geometry rendering */
336     char *ofn;
337     {
338     extern char *index();
339     char *cp;
340    
341     if (ofn == NULL)
342     otEndOctree();
343     else if (access(ofn, R_OK) == 0)
344     otNewOctree(ofn);
345     #ifdef DEBUG
346     else {
347     sprintf(errmsg, "cannot load octree \"%s\"", ofn);
348     error(WARNING, errmsg);
349     }
350     #endif
351     }
352    
353    
354     dev_auxcom(cmd, args) /* process an auxiliary command */
355     char *cmd, *args;
356     {
357     #ifdef DOBJ
358     if (dobj_command(cmd, args) >= 0)
359     return;
360     #endif
361     sprintf(errmsg, "%s: unknown command", cmd);
362     error(COMMAND, errmsg);
363     }
364    
365    
366     VIEW *
367     dev_auxview(n, hvres) /* return nth auxiliary view */
368     int n;
369     int hvres[2];
370     {
371     hvres[0] = odev.hres; hvres[1] = odev.vres;
372     if (n == 0)
373     return(&odev.v);
374     #ifdef STEREO
375     if (n == 1)
376     return(&vwright);
377     #endif
378     return(NULL);
379     }
380    
381    
382     int
383     dev_input() /* get X11 input */
384     {
385     inpresflags = 0;
386    
387     do
388     getevent();
389    
390     while (XPending(ourdisplay) > 0);
391    
392     odev.inpready = 0;
393    
394     return(inpresflags);
395     }
396    
397    
398     dev_value(c, d, p) /* add a pixel value to our texture */
399     COLR c;
400     FVECT d, p;
401     {
402     #ifdef DOBJ
403     if (dobj_lightsamp != NULL) { /* in light source sampling */
404     (*dobj_lightsamp)(c, d, p);
405     return;
406     }
407     #endif
408     odSample(c, d, p); /* add to display representation */
409     if (!--rayqleft)
410     dev_flush(); /* flush output */
411     }
412    
413    
414     int
415     dev_flush() /* flush output as appropriate */
416     {
417     int ndrawn;
418    
419     if (mapped && isperspective > 0) {
420     #ifdef STEREO
421     pushright(); /* draw right eye */
422     otDrawOctrees();
423     #ifdef DOBJ
424     dobj_render();
425     #endif
426     checkglerr("rendering right eye");
427     popright(); /* draw left eye */
428     #endif
429     ndrawn = otDrawOctrees();
430     #ifdef DOBJ
431     ndrawn += dobj_render();
432     #endif
433     checkglerr("rendering base view");
434     }
435     if (mapped && viewsteady) {
436     if (isperspective) { /* first time after steady */
437     #ifdef STEREO
438     pushright();
439     popright();
440     #endif
441     if (ndrawn)
442     xferdepth(); /* transfer and clear depth */
443     setglortho(); /* set orthographic view */
444     } else {
445     #ifdef STEREO
446     pushright();
447     odUpdate(1); /* draw right eye */
448     popright();
449     #endif
450     odUpdate(0); /* draw left eye */
451     }
452     }
453     glFlush(); /* flush OpenGL */
454     rayqleft = RAYQLEN;
455     /* flush X11 and return # pending */
456     return(odev.inpready = XPending(ourdisplay));
457     }
458    
459    
460     checkglerr(where) /* check for GL or GLU error */
461     char *where;
462     {
463     register GLenum errcode;
464    
465     while ((errcode = glGetError()) != GL_NO_ERROR) {
466     sprintf(errmsg, "OpenGL error %s: %s",
467     where, gluErrorString(errcode));
468     error(WARNING, errmsg);
469     }
470     }
471    
472    
473     static
474     xferdepth() /* load and clear depth buffer */
475     {
476     register GLfloat *dbp;
477    
478     if (depthbuffer == NULL) {
479     #ifdef STEREO
480     depthright = (GLfloat *)malloc(
481     odev.hres*odev.vres*sizeof(GLfloat));
482     #endif
483     depthbuffer = (GLfloat *)malloc(
484     odev.hres*odev.vres*sizeof(GLfloat));
485     if (depthbuffer == NULL)
486     error(SYSTEM, "out of memory in xferdepth");
487     }
488     #ifdef STEREO
489     setstereobuf(STEREO_BUFFER_RIGHT);
490     glReadPixels(0, 0, odev.hres, odev.vres,
491     GL_DEPTH_COMPONENT, GL_FLOAT, depthright);
492     setstereobuf(STEREO_BUFFER_LEFT);
493     for (dbp = depthright + odev.hres*odev.vres; dbp-- > depthright; )
494     *dbp = mapdepth(*dbp);
495     odDepthMap(1, depthright);
496     glClear(GL_DEPTH_BUFFER_BIT);
497     #endif
498     glReadPixels(0, 0, odev.hres, odev.vres,
499     GL_DEPTH_COMPONENT, GL_FLOAT, depthbuffer);
500     for (dbp = depthbuffer + odev.hres*odev.vres; dbp-- > depthbuffer; )
501     *dbp = mapdepth(*dbp);
502     odDepthMap(0, depthbuffer); /* transfer depth data */
503     glClear(GL_DEPTH_BUFFER_BIT); /* clear system buffer */
504     }
505    
506    
507     static double
508     getdistance(dx, dy, direc) /* distance from fore plane along view ray */
509     int dx, dy;
510     FVECT direc;
511     {
512     GLfloat gldepth;
513     double dist;
514    
515     if (dx<0 | dx>=odev.hres | dy<0 | dy>=odev.vres)
516     return(FHUGE);
517     if (depthbuffer != NULL)
518     dist = depthbuffer[dy*odev.hres + dx];
519     else {
520     glReadPixels(dx,dy, 1,1, GL_DEPTH_COMPONENT,
521     GL_FLOAT, &gldepth);
522     dist = mapdepth(gldepth);
523     }
524     if (dist >= .99*FHUGE)
525     return(FHUGE);
526     return((dist-odev.v.vfore)/DOT(direc,odev.v.vdir));
527     }
528    
529    
530     #ifdef STEREO
531     static
532     pushright() /* push on right view & buffer */
533     {
534     double d;
535    
536     setstereobuf(STEREO_BUFFER_RIGHT);
537     if (isperspective > 0) {
538     glMatrixMode(GL_MODELVIEW);
539     glPushMatrix();
540     d = -eyesepdist / sqrt(odev.v.hn2);
541     glTranslated(d*odev.v.hvec[0], d*odev.v.hvec[1],
542     d*odev.v.hvec[2]);
543     checkglerr("setting right view");
544     }
545     }
546    
547    
548     static
549     popright() /* pop off right view & buffer */
550     {
551     if (isperspective > 0) {
552     glMatrixMode(GL_MODELVIEW);
553     glPopMatrix();
554     }
555     setstereobuf(STEREO_BUFFER_LEFT);
556     }
557     #endif
558    
559    
560     static int
561     mytmflags() /* figure out tone mapping flags */
562     {
563     extern char *progname;
564     register char *cp, *tail;
565     /* find basic name */
566     for (cp = tail = progname; *cp; cp++)
567     if (*cp == '/')
568     tail = cp+1;
569     for (cp = tail; *cp && *cp != '.'; cp++)
570     ;
571     #ifdef DEBUG
572     if (cp > tail && cp[-1] == 'h')
573     return(TM_F_HUMAN);
574     else
575     return(TM_F_CAMERA);
576     #else
577     if (cp > tail && cp[-1] == 'h')
578     return(TM_F_HUMAN|TM_F_NOSTDERR);
579     else
580     return(TM_F_CAMERA|TM_F_NOSTDERR);
581     #endif
582     }
583    
584    
585     static
586     getevent() /* get next event */
587     {
588     XNextEvent(ourdisplay, levptr(XEvent));
589     switch (levptr(XEvent)->type) {
590     case ConfigureNotify:
591     resizewindow(levptr(XConfigureEvent));
592     break;
593     case UnmapNotify:
594     mapped = 0;
595     break;
596     case MapNotify:
597     mapped = 1;
598     break;
599     case Expose:
600     fixwindow(levptr(XExposeEvent));
601     break;
602     case KeyPress:
603     getkey(levptr(XKeyPressedEvent));
604     break;
605     case ButtonPress:
606     if (FRAMESTATE(levptr(XButtonPressedEvent)->state))
607     getframe(levptr(XButtonPressedEvent));
608     else
609     getmove(levptr(XButtonPressedEvent));
610     break;
611     }
612     }
613    
614    
615     static
616     draw3dline(wp) /* draw 3d line in world coordinates */
617     register FVECT wp[2];
618     {
619     glVertex3d(wp[0][0], wp[0][1], wp[0][2]);
620     glVertex3d(wp[1][0], wp[1][1], wp[1][2]);
621     }
622    
623    
624     static
625     draw_grids(fore) /* draw holodeck section grids */
626     int fore;
627     {
628     glPushAttrib(GL_LIGHTING_BIT|GL_ENABLE_BIT);
629     glDisable(GL_LIGHTING);
630     if (fore)
631     glColor3ub(0, 255, 255);
632     else
633     glColor3ub(0, 0, 0);
634     glBegin(GL_LINES); /* draw each grid line */
635     gridlines(draw3dline);
636     glEnd();
637     checkglerr("drawing grid lines");
638     glPopAttrib();
639     }
640    
641    
642     static
643     moveview(dx, dy, mov, orb) /* move our view */
644     int dx, dy, mov, orb;
645     {
646     VIEW nv;
647     FVECT odir, v1, wip;
648     double d, d1;
649     register int li;
650     /* start with old view */
651     copystruct(&nv, &odev.v);
652     /* orient our motion */
653     if (viewray(v1, odir, &odev.v,
654     (dx+.5)/odev.hres, (dy+.5)/odev.vres) < -FTINY)
655     return(0); /* outside view */
656     if (mov | orb) { /* moving relative to geometry */
657     d = getdistance(dx, dy, odir); /* distance from front plane */
658     #ifdef DOBJ
659     d1 = dobj_trace(NULL, v1, odir);
660     if (d1 < d)
661     d = d1;
662     #endif
663     if (d >= .99*FHUGE)
664     d = 0.5*(dev_zmax+dev_zmin); /* just guess */
665     VSUM(wip, v1, odir, d);
666     VSUB(odir, wip, odev.v.vp);
667     } else /* panning with constant viewpoint */
668     VCOPY(nv.vdir, odir);
669     if (orb && mov) { /* orbit left/right */
670     spinvector(odir, odir, nv.vup, d=MOVDEG*PI/180.*mov);
671     VSUM(nv.vp, wip, odir, -1.);
672     spinvector(nv.vdir, nv.vdir, nv.vup, d);
673     } else if (orb) { /* orbit up/down */
674     fcross(v1, odir, nv.vup);
675     if (normalize(v1) == 0.)
676     return(0);
677     spinvector(odir, odir, v1, d=MOVDEG*PI/180.*orb);
678     VSUM(nv.vp, wip, odir, -1.);
679     spinvector(nv.vdir, nv.vdir, v1, d);
680     } else if (mov) { /* move forward/backward */
681     d = MOVPCT/100. * mov;
682     VSUM(nv.vp, nv.vp, odir, d);
683     }
684     if (!mov ^ !orb && headlocked) { /* restore head height */
685     VSUM(v1, odev.v.vp, nv.vp, -1.);
686     d = DOT(v1, nv.vup);
687     VSUM(nv.vp, nv.vp, odev.v.vup, d);
688     }
689     if (setview(&nv) != NULL)
690     return(0); /* illegal view */
691     dev_view(&nv);
692     inpresflags |= DFL(DC_SETVIEW)|DFL(DC_REDRAW);
693     return(1);
694     }
695    
696    
697     static
698     getframe(ebut) /* get focus frame */
699     XButtonPressedEvent *ebut;
700     {
701     int startx = ebut->x, starty = ebut->y;
702     int endx, endy;
703    
704     XMaskEvent(ourdisplay, ButtonReleaseMask, levptr(XEvent));
705     endx = levptr(XButtonReleasedEvent)->x;
706     endy = levptr(XButtonReleasedEvent)->y;
707     if (endx == startx | endy == starty) {
708     XBell(ourdisplay, 0);
709     return;
710     }
711     if (endx < startx) {register int c = endx; endx = startx; startx = c;}
712     if (endy < starty) {register int c = endy; endy = starty; starty = c;}
713     sprintf(odev_args, "%.3f %.3f %.3f %.3f",
714     (startx+.5)/odev.hres, 1.-(endy+.5)/odev.vres,
715     (endx+.5)/odev.hres, 1.-(starty+.5)/odev.vres);
716     inpresflags |= DFL(DC_FOCUS);
717     }
718    
719    
720     static
721     getmove(ebut) /* get view change */
722     XButtonPressedEvent *ebut;
723     {
724     int movdir = MOVDIR(ebut->button);
725     int movorb = MOVORB(ebut->state);
726     Window rootw, childw;
727     int rootx, rooty, wx, wy;
728     unsigned int statemask;
729    
730     XNoOp(ourdisplay); /* makes sure we're not idle */
731    
732     viewsteady = 0; /* flag moving view */
733     setglpersp(&odev.v); /* start us off in perspective */
734     while (!XCheckMaskEvent(ourdisplay,
735     ButtonReleaseMask, levptr(XEvent))) {
736     /* get cursor position */
737     if (!XQueryPointer(ourdisplay, gwind, &rootw, &childw,
738     &rootx, &rooty, &wx, &wy, &statemask))
739     break; /* on another screen */
740     /* compute view motion */
741     if (!moveview(wx, odev.vres-1-wy, movdir, movorb)) {
742     sleep(1);
743     continue; /* cursor in bad place */
744     }
745     draw_grids(1); /* redraw grid */
746     #ifdef STEREO
747     pushright();
748     draw_grids(1);
749     otDrawOctrees();
750     #ifdef DOBJ
751     dobj_render();
752     #endif
753     popright();
754     #endif
755     /* redraw octrees */
756     otDrawOctrees();
757     #ifdef DOBJ
758     dobj_render(); /* redraw objects */
759     #endif
760     glFlush();
761     }
762     if (!(inpresflags & DFL(DC_SETVIEW))) { /* do final motion */
763     movdir = MOVDIR(levptr(XButtonReleasedEvent)->button);
764     wx = levptr(XButtonReleasedEvent)->x;
765     wy = levptr(XButtonReleasedEvent)->y;
766     moveview(wx, odev.vres-1-wy, movdir, movorb);
767     }
768     }
769    
770    
771     static
772     setglpersp(vp) /* set perspective view in GL */
773     register VIEW *vp;
774     {
775     double d, xmin, xmax, ymin, ymax;
776     GLfloat vec[4];
777     double depthlim[2];
778     /* set depth limits */
779     otDepthLimit(depthlim, odev.v.vp, odev.v.vdir);
780     if (depthlim[0] >= depthlim[1]) {
781     dev_zmin = 1.;
782     dev_zmax = 100.;
783     } else {
784     dev_zmin = 0.5*depthlim[0];
785     dev_zmax = 1.5*depthlim[1];
786     if (dev_zmin > dev_zmax/5.)
787     dev_zmin = dev_zmax/5.;
788     }
789     if (odev.v.vfore > FTINY)
790     dev_zmin = odev.v.vfore;
791     if (odev.v.vaft > FTINY)
792     dev_zmax = odev.v.vaft;
793     if (dev_zmin < dev_zmax/100.)
794     dev_zmin = dev_zmax/100.;
795     setzrat();
796     xmax = dev_zmin * tan(PI/180./2. * odev.v.horiz);
797     xmin = -xmax;
798     d = odev.v.hoff * (xmax - xmin);
799     xmin += d; xmax += d;
800     ymax = dev_zmin * tan(PI/180./2. * odev.v.vert);
801     ymin = -ymax;
802     d = odev.v.voff * (ymax - ymin);
803     ymin += d; ymax += d;
804     /* set view matrix */
805     glMatrixMode(GL_PROJECTION);
806     glLoadIdentity();
807     glFrustum(xmin, xmax, ymin, ymax, dev_zmin, dev_zmax);
808     gluLookAt(odev.v.vp[0], odev.v.vp[1], odev.v.vp[2],
809     odev.v.vp[0] + odev.v.vdir[0],
810     odev.v.vp[1] + odev.v.vdir[1],
811     odev.v.vp[2] + odev.v.vdir[2],
812     odev.v.vup[0], odev.v.vup[1], odev.v.vup[2]);
813     checkglerr("setting perspective view");
814     isperspective = 1;
815     vec[0] = vec[1] = vec[2] = 0.; vec[3] = 1.;
816     glLightModelfv(GL_LIGHT_MODEL_AMBIENT, vec);
817     vec[0] = -odev.v.vdir[0];
818     vec[1] = -odev.v.vdir[1];
819     vec[2] = -odev.v.vdir[2];
820     vec[3] = 0.;
821     glLightfv(GL_LIGHT0, GL_POSITION, vec);
822     vec[0] = vec[1] = vec[2] = .7; vec[3] = 1.;
823     glLightfv(GL_LIGHT0, GL_SPECULAR, vec);
824     glLightfv(GL_LIGHT0, GL_DIFFUSE, vec);
825     vec[0] = vec[1] = vec[2] = .3; vec[3] = 1.;
826     glLightfv(GL_LIGHT0, GL_AMBIENT, vec);
827     glEnable(GL_LIGHT0);
828     glEnable(GL_LIGHTING); /* light our GL objects */
829     glShadeModel(GL_SMOOTH);
830     }
831    
832    
833     static
834     setglortho() /* set up orthographic view for cone drawing */
835     {
836     /* set view matrix */
837     glMatrixMode(GL_PROJECTION);
838     glLoadIdentity();
839     glOrtho(0., (double)odev.hres, 0., (double)odev.vres,
840     0.001*OMAXDEPTH, 1.001*(-OMAXDEPTH));
841     checkglerr("setting orthographic view");
842     isperspective = 0;
843     glDisable(GL_LIGHTING); /* cones are constant color */
844     glShadeModel(GL_FLAT);
845     }
846    
847    
848     static
849     wipeclean() /* prepare for redraw */
850     {
851     /* clear depth buffer */
852     #ifdef STEREO
853     setstereobuf(STEREO_BUFFER_RIGHT);
854     glClear(GL_DEPTH_BUFFER_BIT);
855     setstereobuf(STEREO_BUFFER_LEFT);
856     #endif
857     glClear(GL_DEPTH_BUFFER_BIT);
858     if (viewsteady) /* clear samples if steady */
859     odClean();
860     if (depthbuffer != NULL) { /* free recorded depth buffer */
861     #ifdef STEREO
862     odDepthMap(1, NULL);
863     free((char *)depthright);
864     depthright = NULL;
865     #endif
866     odDepthMap(0, NULL);
867     free((char *)depthbuffer);
868     depthbuffer = NULL;
869     }
870     setglpersp(&odev.v); /* reset view & clipping planes */
871     }
872    
873    
874     static
875     getkey(ekey) /* get input key */
876     register XKeyPressedEvent *ekey;
877     {
878     Window rootw, childw;
879     int rootx, rooty, wx, wy;
880     unsigned int statemask;
881     int n;
882     char buf[8];
883    
884     n = XLookupString(ekey, buf, sizeof(buf), NULL, NULL);
885     if (n != 1)
886     return;
887     switch (buf[0]) {
888     case 'h': /* turn on height motion lock */
889     headlocked = 1;
890     return;
891     case 'H': /* turn off height motion lock */
892     headlocked = 0;
893     return;
894     case 'l': /* retrieve last view */
895     inpresflags |= DFL(DC_LASTVIEW);
896     return;
897     case 'p': /* pause computation */
898     inpresflags |= DFL(DC_PAUSE);
899     return;
900     case 'v': /* spit out view */
901     inpresflags |= DFL(DC_GETVIEW);
902     return;
903     case 'f': /* frame view position */
904     if (!XQueryPointer(ourdisplay, gwind, &rootw, &childw,
905     &rootx, &rooty, &wx, &wy, &statemask))
906     return; /* on another screen */
907     sprintf(odev_args, "%.4f %.4f", (wx+.5)/odev.hres,
908     1.-(wy+.5)/odev.vres);
909     inpresflags |= DFL(DC_FOCUS);
910     return;
911     case 'F': /* unfocus */
912     odev_args[0] = '\0';
913     inpresflags |= DFL(DC_FOCUS);
914     return;
915     case '\n':
916     case '\r': /* resume computation */
917     inpresflags |= DFL(DC_RESUME);
918     return;
919     case CTRL('R'): /* redraw screen */
920     odRemap();
921     glClear(GL_COLOR_BUFFER_BIT|GL_DEPTH_BUFFER_BIT);
922     #ifdef STEREO
923     setstereobuf(STEREO_BUFFER_RIGHT);
924     glClear(GL_COLOR_BUFFER_BIT|GL_DEPTH_BUFFER_BIT);
925     setstereobuf(STEREO_BUFFER_LEFT);
926     #endif
927     return;
928     case CTRL('L'): /* refresh from server */
929     if (inpresflags & DFL(DC_REDRAW))
930     return;
931     XRaiseWindow(ourdisplay, gwind);
932     XFlush(ourdisplay);
933     sleep(1);
934     wipeclean(); /* fresh display */
935     odRemap(); /* new tone mapping */
936     dev_flush(); /* draw octrees */
937     inpresflags |= DFL(DC_REDRAW); /* resend values from server */
938     rayqleft = 0; /* hold off update */
939     return;
940     case 'K': /* kill rtrace process(es) */
941     inpresflags |= DFL(DC_KILL);
942     break;
943     case 'R': /* restart rtrace */
944     inpresflags |= DFL(DC_RESTART);
945     break;
946     case 'C': /* clobber holodeck */
947     inpresflags |= DFL(DC_CLOBBER);
948     break;
949     case 'q': /* quit the program */
950     inpresflags |= DFL(DC_QUIT);
951     return;
952     default:
953     XBell(ourdisplay, 0);
954     return;
955     }
956     }
957    
958    
959     static
960     fixwindow(eexp) /* repair damage to window */
961     register XExposeEvent *eexp;
962     {
963     int xmin, ymin, xmax, ymax;
964    
965     if (odev.hres == 0 | odev.vres == 0) { /* first exposure */
966     resizewindow((XConfigureEvent *)eexp);
967     return;
968     }
969     xmin = eexp->x; xmax = eexp->x + eexp->width;
970     ymin = odev.vres - eexp->y - eexp->height; ymax = odev.vres - eexp->y;
971     /* clear portion of depth */
972     glColorMask(GL_FALSE, GL_FALSE, GL_FALSE, GL_FALSE);
973     glDepthFunc(GL_ALWAYS);
974     glBegin(GL_POLYGON);
975     glVertex3i(xmin, ymin, OMAXDEPTH);
976     glVertex3i(xmax, ymin, OMAXDEPTH);
977     glVertex3i(xmax, ymax, OMAXDEPTH);
978     glVertex3i(xmin, ymax, OMAXDEPTH);
979     glEnd();
980     #ifdef STEREO
981     setstereobuf(STEREO_BUFFER_RIGHT);
982     glBegin(GL_POLYGON);
983     glVertex3i(xmin, ymin, OMAXDEPTH);
984     glVertex3i(xmax, ymin, OMAXDEPTH);
985     glVertex3i(xmax, ymax, OMAXDEPTH);
986     glVertex3i(xmin, ymax, OMAXDEPTH);
987     glEnd();
988     odRedraw(1, xmin, ymin, xmax, ymax);
989     setstereobuf(STEREO_BUFFER_LEFT);
990     #endif
991     glColorMask(GL_TRUE, GL_TRUE, GL_TRUE, GL_TRUE);
992     glDepthFunc(GL_LEQUAL);
993     odRedraw(0, xmin, ymin, xmax, ymax);
994     }
995    
996    
997     static
998     resizewindow(ersz) /* resize window */
999     register XConfigureEvent *ersz;
1000     {
1001     glViewport(0, 0, ersz->width, ersz->height);
1002    
1003     if (ersz->width == odev.hres && ersz->height == odev.vres)
1004     return;
1005    
1006     odev.hres = ersz->width;
1007     odev.vres = ersz->height;
1008    
1009     odev.v.horiz = 2.*180./PI * atan(0.5/VIEWDIST*pwidth*odev.hres);
1010     odev.v.vert = 2.*180./PI * atan(0.5/VIEWDIST*pheight*odev.vres);
1011    
1012     inpresflags |= DFL(DC_SETVIEW)|DFL(DC_REDRAW);
1013     }