ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/hd/rhdobj.c
Revision: 3.15
Committed: Mon Jul 21 22:30:18 2003 UTC (20 years, 9 months ago) by schorsch
Content type: text/plain
Branch: MAIN
Changes since 3.14: +3 -3 lines
Log Message:
Eliminated copystruct() macro, which is unnecessary in ANSI.
Reduced ambiguity warnings for nested if/if/else clauses.

File Contents

# Content
1 #ifndef lint
2 static const char RCSid[] = "$Id: rhdobj.c,v 3.14 2003/06/30 14:59:11 schorsch Exp $";
3 #endif
4 /*
5 * Routines for loading and displaying Radiance objects in rholo with GLX.
6 */
7
8 #include <string.h>
9
10 #include "radogl.h"
11 #include "tonemap.h"
12 #include "rhdisp.h"
13 #include "rhdriver.h"
14 #include "rhdobj.h"
15 #include "rtprocess.h"
16
17 extern FILE *sstdout; /* user standard output */
18
19 char rhdcmd[DO_NCMDS][8] = DO_INIT; /* user command list */
20
21 int (*dobj_lightsamp)() = NULL; /* pointer to function to get lights */
22
23 #define AVGREFL 0.5 /* assumed average reflectance */
24
25 #define MAXAC 512 /* maximum number of args */
26
27 #ifndef MINTHRESH
28 #define MINTHRESH 5.0 /* source threshold w.r.t. mean */
29 #endif
30
31 #ifndef NALT
32 #define NALT 11 /* # sampling altitude angles */
33 #endif
34 #ifndef NAZI
35 #define NAZI ((int)(PI/2.*NALT+.5))
36 #endif
37
38 typedef struct dlights {
39 struct dlights *next; /* next in lighting set list */
40 FVECT lcent; /* computed lighting center */
41 double ravg; /* harmonic mean free range radius */
42 TMbright larb; /* average reflected brightness */
43 COLOR lamb; /* local ambient value */
44 short nl; /* number of lights in this set */
45 struct lsource {
46 FVECT direc; /* source direction */
47 double omega; /* source solid angle */
48 COLOR val; /* source color */
49 } li[MAXLIGHTS]; /* light sources */
50 } DLIGHTS; /* a light source set */
51
52 typedef struct dobject {
53 struct dobject *next; /* next object in list */
54 char name[64]; /* object name */
55 FVECT center; /* orig. object center */
56 RREAL radius; /* orig. object radius */
57 int listid; /* GL display list identifier */
58 int nlists; /* number of lists allocated */
59 SUBPROC rtp; /* associated rtrace process */
60 DLIGHTS *ol; /* object lights */
61 FULLXF xfb; /* coordinate transform */
62 short drawcode; /* drawing code */
63 short xfac; /* transform argument count */
64 char *xfav[MAXAC+1]; /* transform args */
65 } DOBJECT; /* a displayable object */
66
67 static DLIGHTS *dlightsets; /* lighting sets */
68 static DOBJECT *dobjects; /* display object list */
69 static DOBJECT *curobj; /* current (last referred) object */
70 static int lastxfac; /* last number of transform args */
71 static char *lastxfav[MAXAC+1]; /* saved transform arguments */
72
73 #define getdcent(c,op) multp3(c,(op)->center,(op)->xfb.f.xfm)
74 #define getdrad(op) ((op)->radius*(op)->xfb.f.sca)
75
76 #define RTARGC 8
77 static char *rtargv[RTARGC+1] = {"rtrace", "-h-", "-w-", "-fdd",
78 "-x", "1", "-oL"};
79
80 static struct {
81 int nsamp; /* number of ray samples */
82 COLOR val; /* value (sum) */
83 } ssamp[NALT][NAZI]; /* current sphere samples */
84
85 #define curname (curobj==NULL ? (char *)NULL : curobj->name)
86
87
88 static DOBJECT *
89 getdobj(nm) /* get object from list by name */
90 char *nm;
91 {
92 register DOBJECT *op;
93
94 if (nm == NULL)
95 return(NULL);
96 if (nm == curname)
97 return(curobj);
98 for (op = dobjects; op != NULL; op = op->next)
99 if (!strcmp(op->name, nm))
100 break;
101 return(op);
102 }
103
104
105 static
106 freedobj(op) /* free resources and memory assoc. with op */
107 register DOBJECT *op;
108 {
109 int foundlink = 0;
110 DOBJECT ohead;
111 register DOBJECT *opl;
112
113 if (op == NULL)
114 return(0);
115 ohead.next = dobjects;
116 for (opl = &ohead; opl->next != NULL; opl = opl->next) {
117 if (opl->next == op && (opl->next = op->next) == NULL)
118 break;
119 foundlink += opl->next->listid == op->listid;
120 }
121 dobjects = ohead.next;
122 if (!foundlink) {
123 glDeleteLists(op->listid, op->nlists);
124 close_process(&(op->rtp));
125 }
126 while (op->xfac)
127 freestr(op->xfav[--op->xfac]);
128 free((void *)op);
129 return(1);
130 }
131
132
133 static
134 savedxf(op) /* save transform for display object */
135 register DOBJECT *op;
136 {
137 /* free old */
138 while (lastxfac)
139 freestr(lastxfav[--lastxfac]);
140 /* nothing to save? */
141 if (op == NULL) {
142 lastxfav[0] = NULL;
143 return(0);
144 }
145 /* else save new */
146 for (lastxfac = 0; lastxfac < op->xfac; lastxfac++)
147 lastxfav[lastxfac] = savestr(op->xfav[lastxfac]);
148 lastxfav[lastxfac] = NULL;
149 return(1);
150 }
151
152
153 static
154 ssph_sample(clr, direc, pos) /* add sample to current source sphere */
155 COLR clr;
156 FVECT direc, pos;
157 {
158 COLOR col;
159 double d;
160 register int alt, azi;
161
162 if (dlightsets == NULL)
163 return;
164 if (pos == NULL)
165 d = FHUGE; /* sample is at infinity */
166 else if ((d = (pos[0] - dlightsets->lcent[0])*direc[0] +
167 (pos[1] - dlightsets->lcent[1])*direc[1] +
168 (pos[2] - dlightsets->lcent[2])*direc[2]) > FTINY)
169 dlightsets->ravg += 1./d;
170 else
171 return; /* sample is behind us */
172 alt = NALT*(1.-FTINY)*(.5*(1.+FTINY) + direc[2]*.5);
173 azi = NAZI*(1.-FTINY)*(.5*(1.+FTINY) + .5/PI*atan2(direc[1],direc[0]));
174 colr_color(col, clr);
175 addcolor(ssamp[alt][azi].val, col);
176 ssamp[alt][azi].nsamp++;
177 }
178
179
180 static
181 ssph_direc(direc, alt, azi) /* compute sphere sampling direction */
182 FVECT direc;
183 int alt, azi;
184 {
185 double phi, d;
186
187 direc[2] = 2./NALT*(alt+.5) - 1.;
188 d = sqrt(1. - direc[2]*direc[2]);
189 phi = 2.*PI/NAZI*(azi+.5) - PI;
190 direc[0] = d*cos(phi);
191 direc[1] = d*sin(phi);
192 }
193
194
195 static int
196 ssph_neigh(sp, next) /* neighbor counter on sphere */
197 register int sp[2];
198 int next;
199 {
200 static short nneigh = 0; /* neighbor count */
201 static short neighlist[NAZI+6][2]; /* neighbor list (0 is home) */
202 register int i;
203
204 if (next) {
205 if (nneigh <= 0)
206 return(0);
207 sp[0] = neighlist[--nneigh][0];
208 sp[1] = neighlist[nneigh][1];
209 return(1);
210 }
211 if (sp[0] < 0 | sp[0] >= NALT | sp[1] < 0 | sp[1] >= NAZI)
212 return(nneigh=0);
213 neighlist[0][0] = sp[0]; neighlist[0][1] = sp[1];
214 nneigh = 1;
215 if (sp[0] == 0) {
216 neighlist[nneigh][0] = 1;
217 neighlist[nneigh++][1] = (sp[1]+1)%NAZI;
218 neighlist[nneigh][0] = 1;
219 neighlist[nneigh++][1] = sp[1];
220 neighlist[nneigh][0] = 1;
221 neighlist[nneigh++][1] = (sp[1]+(NAZI-1))%NAZI;
222 for (i = 0; i < NAZI; i++)
223 if (i != sp[1]) {
224 neighlist[nneigh][0] = 0;
225 neighlist[nneigh++][1] = i;
226 }
227 } else if (sp[0] == NALT-1) {
228 neighlist[nneigh][0] = NALT-2;
229 neighlist[nneigh++][1] = (sp[1]+1)%NAZI;
230 neighlist[nneigh][0] = NALT-2;
231 neighlist[nneigh++][1] = sp[1];
232 neighlist[nneigh][0] = NALT-2;
233 neighlist[nneigh++][1] = (sp[1]+(NAZI-1))%NAZI;
234 for (i = 0; i < NAZI; i++)
235 if (i != sp[1]) {
236 neighlist[nneigh][0] = NALT-1;
237 neighlist[nneigh++][1] = i;
238 }
239 } else {
240 neighlist[nneigh][0] = sp[0]-1;
241 neighlist[nneigh++][1] = (sp[1]+1)%NAZI;
242 neighlist[nneigh][0] = sp[0]-1;
243 neighlist[nneigh++][1] = sp[1];
244 neighlist[nneigh][0] = sp[0]-1;
245 neighlist[nneigh++][1] = (sp[1]+(NAZI-1))%NAZI;
246 neighlist[nneigh][0] = sp[0];
247 neighlist[nneigh++][1] = (sp[1]+1)%NAZI;
248 neighlist[nneigh][0] = sp[0];
249 neighlist[nneigh++][1] = (sp[1]+(NAZI-1))%NAZI;
250 neighlist[nneigh][0] = sp[0]+1;
251 neighlist[nneigh++][1] = (sp[1]+1)%NAZI;
252 neighlist[nneigh][0] = sp[0]+1;
253 neighlist[nneigh++][1] = sp[1];
254 neighlist[nneigh][0] = sp[0]+1;
255 neighlist[nneigh++][1] = (sp[1]+(NAZI-1))%NAZI;
256 }
257 return(nneigh);
258 }
259
260
261 static
262 ssph_compute() /* compute source set from sphere samples */
263 {
264 int ncells, nsamps;
265 COLOR csum;
266 FVECT v;
267 double d, thresh, maxbr;
268 int maxalt, maxazi, spos[2];
269 register int alt, azi;
270 register struct lsource *ls;
271 /* count & average sampled cells */
272 setcolor(csum, 0., 0., 0.);
273 ncells = nsamps = 0;
274 for (alt = 0; alt < NALT; alt++)
275 for (azi = 0; azi < NAZI; azi++)
276 if (ssamp[alt][azi].nsamp) {
277 if (ssamp[alt][azi].nsamp > 1) {
278 d = 1.0/ssamp[alt][azi].nsamp;
279 scalecolor(ssamp[alt][azi].val, d);
280 }
281 addcolor(csum, ssamp[alt][azi].val);
282 nsamps += ssamp[alt][azi].nsamp;
283 ncells++;
284 }
285 if (dlightsets == NULL | ncells < NALT*NAZI/4) {
286 ncells = 0;
287 goto done;
288 }
289 /* harmonic mean distance */
290 if (dlightsets->ravg > FTINY)
291 dlightsets->ravg = nsamps / dlightsets->ravg;
292 else
293 dlightsets->ravg = FHUGE;
294 /* light source threshold */
295 thresh = MINTHRESH*bright(csum)/ncells;
296 if (thresh <= FTINY) {
297 ncells = 0;
298 goto done;
299 }
300 /* avg. reflected brightness */
301 d = AVGREFL / (double)ncells;
302 scalecolor(csum, d);
303 if (tmCvColors(&dlightsets->larb, TM_NOCHROM, &csum, 1) != TM_E_OK)
304 error(CONSISTENCY, "tone mapping problem in ssph_compute");
305 /* greedy light source clustering */
306 while (dlightsets->nl < MAXLIGHTS) {
307 maxbr = 0.; /* find brightest cell */
308 for (alt = 0; alt < NALT; alt++)
309 for (azi = 0; azi < NAZI; azi++)
310 if ((d = bright(ssamp[alt][azi].val)) > maxbr) {
311 maxalt = alt; maxazi = azi;
312 maxbr = d;
313 }
314 if (maxbr < thresh) /* below threshold? */
315 break;
316 ls = dlightsets->li + dlightsets->nl++;
317 spos[0] = maxalt; spos[1] = maxazi; /* cluster */
318 for (ssph_neigh(spos, 0); ssph_neigh(spos, 1); ) {
319 alt = spos[0]; azi = spos[1];
320 if ((d = bright(ssamp[alt][azi].val)) < .75*thresh)
321 continue; /* too dim */
322 ssph_direc(v, alt, azi); /* else add it in */
323 VSUM(ls->direc, ls->direc, v, d);
324 ls->omega += 1.;
325 addcolor(ls->val, ssamp[alt][azi].val);
326 /* remove from list */
327 setcolor(ssamp[alt][azi].val, 0., 0., 0.);
328 ssamp[alt][azi].nsamp = 0;
329 }
330 d = 1./ls->omega; /* avg. brightness */
331 scalecolor(ls->val, d);
332 ls->omega *= 4.*PI/(NALT*NAZI); /* solid angle */
333 normalize(ls->direc); /* direction */
334 }
335 /* compute ambient remainder */
336 for (alt = 0; alt < NALT; alt++)
337 for (azi = 0; azi < NAZI; azi++)
338 if (ssamp[alt][azi].nsamp)
339 addcolor(dlightsets->lamb, ssamp[alt][azi].val);
340 d = 1.0/ncells;
341 scalecolor(dlightsets->lamb, d);
342 done: /* clear sphere sample array */
343 memset((void *)ssamp, '\0', sizeof(ssamp));
344 return(ncells);
345 }
346
347
348 static
349 getdlights(op, force) /* get lights for display object */
350 register DOBJECT *op;
351 int force;
352 {
353 double d2, mind2 = FHUGE*FHUGE;
354 FVECT ocent;
355 VIEW cvw;
356 register DLIGHTS *dl;
357
358 op->ol = NULL;
359 if (op->drawcode != DO_LIGHT)
360 return(0);
361 /* check for usable light set */
362 getdcent(ocent, op);
363 for (dl = dlightsets; dl != NULL; dl = dl->next)
364 if ((d2 = dist2(dl->lcent, ocent)) < mind2) {
365 op->ol = dl;
366 mind2 = d2;
367 }
368 /* the following is heuristic */
369 d2 = 2.*getdrad(op); d2 *= d2;
370 if ((dl = op->ol) != NULL && (mind2 < 0.0625*dl->ravg*dl->ravg ||
371 mind2 < 4.*getdrad(op)*getdrad(op)))
372 return(1);
373 if (!force)
374 return(0);
375 /* need to compute new light set */
376 cvw = stdview;
377 cvw.type = VT_PER;
378 VCOPY(cvw.vp, ocent);
379 cvw.vup[0] = 1.; cvw.vup[1] = cvw.vup[2] = 0.;
380 cvw.horiz = 90; cvw.vert = 90.;
381 beam_init(1); /* query beams through center */
382 cvw.vdir[0] = cvw.vdir[1] = 0.; cvw.vdir[2] = 1.;
383 setview(&cvw); beam_view(&cvw, 0, 0);
384 cvw.vdir[0] = cvw.vdir[1] = 0.; cvw.vdir[2] = -1.;
385 setview(&cvw); beam_view(&cvw, 0, 0);
386 cvw.vup[0] = cvw.vup[1] = 0.; cvw.vup[2] = 1.;
387 cvw.vdir[0] = cvw.vdir[2] = 0.; cvw.vdir[1] = 1.;
388 setview(&cvw); beam_view(&cvw, 0, 0);
389 cvw.vdir[0] = cvw.vdir[2] = 0.; cvw.vdir[1] = -1.;
390 setview(&cvw); beam_view(&cvw, 0, 0);
391 cvw.vdir[1] = cvw.vdir[2] = 0.; cvw.vdir[0] = 1.;
392 setview(&cvw); beam_view(&cvw, 0, 0);
393 cvw.vdir[1] = cvw.vdir[2] = 0.; cvw.vdir[0] = -1.;
394 setview(&cvw); beam_view(&cvw, 0, 0);
395 /* allocate new light set */
396 dl = (DLIGHTS *)calloc(1, sizeof(DLIGHTS));
397 if (dl == NULL)
398 goto memerr;
399 VCOPY(dl->lcent, ocent);
400 /* push onto our light set list */
401 dl->next = dlightsets;
402 dlightsets = dl;
403 dobj_lightsamp = ssph_sample; /* get beams from server */
404 imm_mode = beam_sync(-1) > 0;
405 while (imm_mode)
406 if (serv_result() == DS_SHUTDOWN)
407 quit(0);
408 if (!ssph_compute()) { /* compute light sources from sphere */
409 dlightsets = dl->next;
410 free((void *)dl);
411 return(0);
412 }
413 op->ol = dl;
414 return(1);
415 memerr:
416 error(SYSTEM, "out of memory in getdlights");
417 }
418
419
420 static
421 cmderror(cn, err) /* report command error */
422 int cn;
423 char *err;
424 {
425 sprintf(errmsg, "%s: %s", rhdcmd[cn], err);
426 error(COMMAND, errmsg);
427 }
428
429
430 int
431 dobj_command(cmd, args) /* run object display command */
432 char *cmd;
433 register char *args;
434 {
435 int somechange = 0;
436 int cn, na, doxfm;
437 register int nn;
438 char *alist[MAXAC+1], *nm;
439 /* find command */
440 for (cn = 0; cn < DO_NCMDS; cn++)
441 if (!strcmp(cmd, rhdcmd[cn]))
442 break;
443 if (cn >= DO_NCMDS)
444 return(-1); /* not in our list */
445 /* make argument list */
446 for (na = 0; *args; na++) {
447 if (na > MAXAC)
448 goto toomany;
449 alist[na] = args;
450 while (*args && !isspace(*args))
451 args++;
452 while (isspace(*args))
453 *args++ = '\0';
454 }
455 alist[na] = NULL;
456 /* execute command */
457 switch (cn) {
458 case DO_LOAD: /* load an octree */
459 if (na == 1)
460 dobj_load(alist[0], alist[0]);
461 else if (na == 2)
462 dobj_load(alist[0], alist[1]);
463 else {
464 cmderror(cn, "need octree [name]");
465 return(0);
466 }
467 break;
468 case DO_UNLOAD: /* clear an object */
469 if (na > 1) goto toomany;
470 if (na && alist[0][0] == '*')
471 somechange += dobj_cleanup();
472 else
473 somechange += dobj_unload(na ? alist[0] : curname);
474 break;
475 case DO_XFORM: /* transform object */
476 case DO_MOVE:
477 if (na && alist[0][0] != '-') {
478 nm = alist[0]; nn = 1;
479 } else {
480 nm = curname; nn = 0;
481 }
482 if (cn == DO_MOVE && nn >= na) {
483 cmderror(cn, "missing transform");
484 return(0);
485 }
486 somechange += dobj_xform(nm, cn==DO_MOVE, na-nn, alist+nn);
487 break;
488 case DO_UNMOVE: /* undo last transform */
489 somechange += dobj_unmove();
490 break;
491 case DO_OBJECT: /* print object statistics */
492 if (dobj_putstats(na ? alist[0] : curname, sstdout))
493 if (na && alist[0][0] != '*' && (curobj == NULL ||
494 strcmp(alist[0], curobj->name)))
495 savedxf(curobj = getdobj(alist[0]));
496 break;
497 case DO_DUP: /* duplicate object */
498 for (nn = 0; nn < na; nn++)
499 if (alist[nn][0] == '-')
500 break;
501 switch (nn) {
502 case 0:
503 cmderror(cn, "need new object name");
504 return(0);
505 case 1:
506 nm = curname;
507 break;
508 case 2:
509 nm = alist[0];
510 break;
511 default:
512 goto toomany;
513 }
514 if (!dobj_dup(nm, alist[nn-1]))
515 break;
516 if (na > nn)
517 somechange += dobj_xform(curname, 1, na-nn, alist+nn);
518 else
519 curobj->drawcode = DO_HIDE;
520 savedxf(curobj);
521 break;
522 case DO_SHOW: /* change rendering option */
523 case DO_LIGHT:
524 case DO_HIDE:
525 if (na > 1) goto toomany;
526 dobj_lighting(na ? alist[0] : curname, cn);
527 somechange++;
528 break;
529 default:
530 error(CONSISTENCY, "bad command id in dobj_command");
531 }
532 return(somechange);
533 toomany:
534 return(cmderror(cn, "too many arguments"));
535 }
536
537
538 dobj_load(oct, nam) /* create/load an octree object */
539 char *oct, *nam;
540 {
541 char *fpp, fpath[128];
542 register DOBJECT *op;
543 /* check arguments */
544 if (oct == NULL) {
545 error(COMMAND, "missing octree");
546 return(0);
547 }
548 if (nam == NULL) {
549 error(COMMAND, "missing name");
550 return(0);
551 }
552 if (*nam == '*' | *nam == '-') {
553 error(COMMAND, "illegal name");
554 return(0);
555 }
556 if (getdobj(nam) != NULL) {
557 error(COMMAND, "name already taken (clear first)");
558 return(0);
559 }
560 /* get octree path */
561 if ((fpp = getpath(oct, getrlibpath(), R_OK)) == NULL) {
562 sprintf(errmsg, "cannot find octree \"%s\"", oct);
563 error(COMMAND, errmsg);
564 return(0);
565 }
566 strcpy(fpath, fpp);
567 op = (DOBJECT *)malloc(sizeof(DOBJECT));
568 if (op == NULL)
569 error(SYSTEM, "out of memory in dobj_load");
570 /* set struct fields */
571 strcpy(op->name, nam);
572 op->ol = NULL;
573 op->drawcode = DO_HIDE;
574 setident4(op->xfb.f.xfm); op->xfb.f.sca = 1.;
575 setident4(op->xfb.b.xfm); op->xfb.b.sca = 1.;
576 op->xfav[op->xfac=0] = NULL;
577 /* load octree into display list */
578 dolights = 0;
579 domats = 1;
580 op->listid = rgl_octlist(fpath, op->center, &op->radius, &op->nlists);
581 /* start rtrace */
582 rtargv[RTARGC-1] = fpath;
583 rtargv[RTARGC] = NULL;
584 open_process(&(op->rtp), rtargv);
585 /* insert into main list */
586 op->next = dobjects;
587 curobj = dobjects = op;
588 savedxf(NULL);
589 return(1);
590 }
591
592
593 dobj_unload(nam) /* free the named object */
594 char *nam;
595 {
596 register DOBJECT *op;
597
598 if ((op = getdobj(nam)) == NULL) {
599 error(COMMAND, "no object");
600 return(0);
601 }
602 freedobj(op);
603 savedxf(curobj = NULL);
604 return(1);
605 }
606
607
608 dobj_cleanup() /* free all resources */
609 {
610 register DLIGHTS *lp;
611
612 while (dobjects != NULL)
613 freedobj(dobjects);
614 savedxf(curobj = NULL);
615 while ((lp = dlightsets) != NULL) {
616 dlightsets = lp->next;
617 free((void *)lp);
618 }
619 return(1);
620 }
621
622
623 dobj_xform(nam, rel, ac, av) /* set/add transform for nam */
624 char *nam;
625 int rel, ac;
626 char **av;
627 {
628 register DOBJECT *op;
629 FVECT cent;
630 double rad;
631 char scoord[16];
632 int i;
633
634 if ((op = getdobj(nam)) == NULL) {
635 error(COMMAND, "no object");
636 return(0);
637 }
638 if (rel)
639 rel = op->xfac + 8;
640 if (ac + rel > MAXAC) {
641 error(COMMAND, "too many transform arguments");
642 return(0);
643 }
644 savedxf(curobj = op); /* remember current transform */
645 if (rel && ac == 4 && !strcmp(av[0], "-t"))
646 rel = -1; /* don't move for translate */
647 else {
648 getdcent(cent, op); /* don't move if near orig. */
649 rad = getdrad(op);
650 if (DOT(cent,cent) < rad*rad)
651 rel = -1;
652 }
653 if (!rel) { /* remove old transform */
654 while (op->xfac)
655 freestr(op->xfav[--op->xfac]);
656 } else if (rel > 0) { /* relative move */
657 op->xfav[op->xfac++] = savestr("-t");
658 for (i = 0; i < 3; i++) {
659 sprintf(scoord, "%.4e", -cent[i]);
660 op->xfav[op->xfac++] = savestr(scoord);
661 }
662 }
663 while (ac--)
664 op->xfav[op->xfac++] = savestr(*av++);
665 if (rel > 0) { /* move back */
666 op->xfav[op->xfac++] = savestr("-t");
667 for (i = 0; i < 3; i++) {
668 sprintf(scoord, "%.4e", cent[i]);
669 op->xfav[op->xfac++] = savestr(scoord);
670 }
671 }
672 op->xfav[op->xfac] = NULL;
673 if (fullxf(&op->xfb, op->xfac, op->xfav) != op->xfac) {
674 error(COMMAND, "bad transform arguments");
675 dobj_unmove();
676 savedxf(op); /* save current transform instead */
677 return(0);
678 }
679 /* don't know local lights anymore */
680 getdlights(op, 0);
681 return(1);
682 }
683
684
685 dobj_putstats(nam, fp) /* put out statistics for nam */
686 char *nam;
687 FILE *fp;
688 {
689 FVECT ocent;
690 register DOBJECT *op;
691 register int i;
692
693 if (nam == NULL) {
694 error(COMMAND, "no current object");
695 return(0);
696 }
697 if (nam[0] == '*') {
698 i = 0;
699 for (op = dobjects; op != NULL; op = op->next)
700 i += dobj_putstats(op->name, fp);
701 return(i);
702 }
703 if ((op = getdobj(nam)) == NULL) {
704 error(COMMAND, "unknown object");
705 return(0);
706 }
707 getdcent(ocent, op);
708 fprintf(fp, "%s: %s, center [%g %g %g], radius %g", op->name,
709 op->drawcode==DO_HIDE ? "hidden" :
710 op->drawcode==DO_LIGHT && op->ol!=NULL ? "lighted" :
711 "shown",
712 ocent[0],ocent[1],ocent[2], getdrad(op));
713 if (op->xfac)
714 fputs(", (xform", fp);
715 for (i = 0; i < op->xfac; i++) {
716 putc(' ', fp);
717 fputs(op->xfav[i], fp);
718 }
719 if (op->xfac)
720 fputc(')', fp);
721 fputc('\n', fp);
722 return(1);
723 }
724
725
726 dobj_unmove() /* undo last transform change */
727 {
728 int txfac;
729 char *txfav[MAXAC+1];
730
731 if (curobj == NULL) {
732 error(COMMAND, "no current object");
733 return(0);
734 }
735 /* hold last transform */
736 memcpy((void *)txfav, (void *)lastxfav,
737 (txfac=lastxfac)*sizeof(char *));
738 /* save this transform */
739 memcpy((void *)lastxfav, (void *)curobj->xfav,
740 (lastxfac=curobj->xfac)*sizeof(char *));
741 /* copy back last transform */
742 memcpy((void *)curobj->xfav, (void *)txfav,
743 (curobj->xfac=txfac)*sizeof(char *));
744 /* set matrices */
745 fullxf(&curobj->xfb, curobj->xfac, curobj->xfav);
746 /* don't know local lights anymore */
747 getdlights(curobj, 0);
748 return(1);
749 }
750
751
752 dobj_dup(oldnm, nam) /* duplicate object oldnm as nam */
753 char *oldnm, *nam;
754 {
755 register DOBJECT *op, *opdup;
756 /* check arguments */
757 if ((op = getdobj(oldnm)) == NULL) {
758 error(COMMAND, "no object");
759 return(0);
760 }
761 if (nam == NULL) {
762 error(COMMAND, "missing name");
763 return(0);
764 }
765 if (*nam == '*' | *nam == '-') {
766 error(COMMAND, "illegal name");
767 return(0);
768 }
769 if (getdobj(nam) != NULL) {
770 error(COMMAND, "name already taken (clear first)");
771 return(0);
772 }
773 /* allocate and copy struct */
774 opdup = (DOBJECT *)malloc(sizeof(DOBJECT));
775 if (opdup == NULL)
776 error(SYSTEM, "out of memory in dobj_dup");
777 *opdup = *op;
778 /* rename */
779 strcpy(opdup->name, nam);
780 /* get our own copy of transform */
781 for (opdup->xfac = 0; opdup->xfac < op->xfac; opdup->xfac++)
782 opdup->xfav[opdup->xfac] = savestr(op->xfav[opdup->xfac]);
783 opdup->xfav[opdup->xfac] = NULL;
784 /* insert it into our list */
785 opdup->next = dobjects;
786 curobj = dobjects = opdup;
787 return(1);
788 }
789
790
791 dobj_lighting(nam, cn) /* set up lighting for display object */
792 char *nam;
793 int cn;
794 {
795 int i, res[2];
796 VIEW *dv;
797 register DOBJECT *op;
798
799 if (nam == NULL) {
800 error(COMMAND, "no current object");
801 return(0);
802 }
803 if (nam[0] == '*') {
804 for (op = dobjects; op != NULL; op = op->next)
805 if ((op->drawcode = cn) == DO_LIGHT)
806 getdlights(op, 1);
807 else
808 op->ol = NULL;
809 } else if ((op = getdobj(nam)) == NULL) {
810 error(COMMAND, "unknown object");
811 return(0);
812 } else if ((op->drawcode = cn) == DO_LIGHT) {
813 if (!getdlights(op, 1))
814 error(COMMAND, "insufficient samples to light object");
815 } else
816 op->ol = NULL;
817
818 if (dobj_lightsamp != NULL) { /* restore beam set */
819 dobj_lightsamp = NULL;
820 beam_init(1);
821 for (i = 0; (dv = dev_auxview(i, res)) != NULL; i++)
822 beam_view(dv, res[0], res[1]);
823 beam_sync(1); /* update server */
824 }
825 }
826
827
828 double
829 dobj_trace(nm, rorg, rdir) /* check for ray intersection with object(s) */
830 char nm[];
831 FVECT rorg, rdir;
832 {
833 register DOBJECT *op;
834 FVECT xorg, xdir;
835 double darr[6];
836 /* check each visible object? */
837 if (nm == NULL || *nm == '*') {
838 double dist, mindist = 1.01*FHUGE;
839
840 if (nm != NULL) nm[0] = '\0';
841 for (op = dobjects; op != NULL; op = op->next) {
842 if (op->drawcode == DO_HIDE)
843 continue;
844 dist = dobj_trace(op->name, rorg, rdir);
845 if (dist < mindist) {
846 if (nm != NULL) strcpy(nm, op->name);
847 mindist = dist;
848 }
849 }
850 return(mindist);
851 }
852 /* else check particular object */
853 if ((op = getdobj(nm)) == NULL) {
854 error(COMMAND, "unknown object");
855 return(FHUGE);
856 }
857 if (op->xfac) { /* put ray in local coordinates */
858 multp3(xorg, rorg, op->xfb.b.xfm);
859 multv3(xdir, rdir, op->xfb.b.xfm);
860 VCOPY(darr, xorg); VCOPY(darr+3, xdir);
861 } else {
862 VCOPY(darr, rorg); VCOPY(darr+3, rdir);
863 }
864 /* trace it */
865 if (process(&(op->rtp), (char *)darr, (char *)darr, sizeof(double),
866 6*sizeof(double)) != sizeof(double))
867 error(SYSTEM, "rtrace communication error");
868 /* return distance */
869 if (darr[0] >= .99*FHUGE)
870 return(FHUGE);
871 return(darr[0]*op->xfb.f.sca);
872 }
873
874
875 int
876 dobj_render() /* render our objects in OpenGL */
877 {
878 int nrendered = 0;
879 GLboolean normalizing;
880 GLfloat vec[4];
881 FVECT v1;
882 register DOBJECT *op;
883 register int i;
884 /* anything to render? */
885 for (op = dobjects; op != NULL; op = op->next)
886 if (op->drawcode != DO_HIDE)
887 break;
888 if (op == NULL)
889 return(0);
890 /* set up general rendering params */
891 glGetBooleanv(GL_NORMALIZE, &normalizing);
892 glPushAttrib(GL_LIGHTING_BIT|GL_TRANSFORM_BIT|GL_ENABLE_BIT|
893 GL_COLOR_BUFFER_BIT|GL_DEPTH_BUFFER_BIT|GL_POLYGON_BIT);
894 glDepthFunc(GL_LESS);
895 glEnable(GL_DEPTH_TEST);
896 glLightModeli(GL_LIGHT_MODEL_TWO_SIDE, GL_TRUE);
897 glPolygonMode(GL_FRONT_AND_BACK, GL_FILL);
898 glShadeModel(GL_SMOOTH);
899 glFrontFace(GL_CCW);
900 glDisable(GL_CULL_FACE);
901 for (i = MAXLIGHTS; i--; )
902 glDisable(glightid[i]);
903 glEnable(GL_LIGHTING);
904 rgl_checkerr("setting rendering modes in dobj_render");
905 /* render each object */
906 for (op = dobjects; op != NULL; op = op->next) {
907 if (op->drawcode == DO_HIDE)
908 continue;
909 /* set up lighting */
910 if (op->drawcode == DO_LIGHT && op->ol != NULL) {
911 BYTE pval;
912 double expval, d;
913 /* use computed sources */
914 if (tmMapPixels(&pval, &op->ol->larb, TM_NOCHROM, 1)
915 != TM_E_OK)
916 error(CONSISTENCY, "dobj_render w/o tone map");
917 expval = pval * (WHTEFFICACY/256.) /
918 tmLuminance(op->ol->larb);
919 vec[0] = expval * colval(op->ol->lamb,RED);
920 vec[1] = expval * colval(op->ol->lamb,GRN);
921 vec[2] = expval * colval(op->ol->lamb,BLU);
922 vec[3] = 1.;
923 glLightModelfv(GL_LIGHT_MODEL_AMBIENT, vec);
924 for (i = op->ol->nl; i--; ) {
925 VCOPY(vec, op->ol->li[i].direc);
926 vec[3] = 0.;
927 glLightfv(glightid[i], GL_POSITION, vec);
928 d = expval * op->ol->li[i].omega;
929 vec[0] = d * colval(op->ol->li[i].val,RED);
930 vec[1] = d * colval(op->ol->li[i].val,GRN);
931 vec[2] = d * colval(op->ol->li[i].val,BLU);
932 vec[3] = 1.;
933 glLightfv(glightid[i], GL_SPECULAR, vec);
934 glLightfv(glightid[i], GL_DIFFUSE, vec);
935 vec[0] = vec[1] = vec[2] = 0.; vec[3] = 1.;
936 glLightfv(glightid[i], GL_AMBIENT, vec);
937 glEnable(glightid[i]);
938 }
939 } else { /* fake lighting */
940 vec[0] = vec[1] = vec[2] = 0.; vec[3] = 1.;
941 glLightModelfv(GL_LIGHT_MODEL_AMBIENT, vec);
942 getdcent(v1, op);
943 VSUB(v1, odev.v.vp, v1);
944 if (normalize(v1) <= getdrad(op)) {
945 vec[0] = -odev.v.vdir[0];
946 vec[1] = -odev.v.vdir[1];
947 vec[2] = -odev.v.vdir[2];
948 } else
949 VCOPY(vec, v1);
950 vec[3] = 0.;
951 glLightfv(GL_LIGHT0, GL_POSITION, vec);
952 vec[0] = vec[1] = vec[2] = .7; vec[3] = 1.;
953 glLightfv(GL_LIGHT0, GL_SPECULAR, vec);
954 glLightfv(GL_LIGHT0, GL_DIFFUSE, vec);
955 vec[0] = vec[1] = vec[2] = .3; vec[3] = 1.;
956 glLightfv(GL_LIGHT0, GL_AMBIENT, vec);
957 glEnable(GL_LIGHT0);
958 }
959 /* set up object transform */
960 if (op->xfac) {
961 if (!normalizing && op->xfb.f.sca < 1.-FTINY |
962 op->xfb.f.sca > 1.+FTINY)
963 glEnable(GL_NORMALIZE);
964 glMatrixMode(GL_MODELVIEW);
965 glPushMatrix();
966 /* matrix order works out to same */
967 #ifdef SMLFLT
968 glMultMatrixf((GLfloat *)op->xfb.f.xfm);
969 #else
970 glMultMatrixd((GLdouble *)op->xfb.f.xfm);
971 #endif
972 }
973 /* render the display list */
974 glCallList(op->listid);
975 nrendered++;
976 /* restore matrix */
977 if (op->xfac) {
978 glMatrixMode(GL_MODELVIEW);
979 glPopMatrix();
980 if (!normalizing)
981 glDisable(GL_NORMALIZE);
982 }
983 /* restore lighting */
984 if (op->drawcode == DO_LIGHT && op->ol != NULL)
985 for (i = op->ol->nl; i--; )
986 glDisable(glightid[i]);
987 else
988 glDisable(GL_LIGHT0);
989 /* check errors */
990 }
991 glPopAttrib(); /* restore rendering params */
992 rgl_checkerr("rendering objects in dobj_render");
993 return(nrendered);
994 }