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

# User Rev Content
1 gwlarson 3.1 #ifndef lint
2 schorsch 3.15 static const char RCSid[] = "$Id: rhdobj.c,v 3.14 2003/06/30 14:59:11 schorsch Exp $";
3 gwlarson 3.1 #endif
4     /*
5 gwlarson 3.5 * Routines for loading and displaying Radiance objects in rholo with GLX.
6 gwlarson 3.1 */
7    
8 schorsch 3.14 #include <string.h>
9    
10 gwlarson 3.1 #include "radogl.h"
11     #include "tonemap.h"
12     #include "rhdisp.h"
13     #include "rhdriver.h"
14     #include "rhdobj.h"
15 schorsch 3.13 #include "rtprocess.h"
16 gwlarson 3.1
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 gwlarson 3.6 #define MAXAC 512 /* maximum number of args */
26 gwlarson 3.1
27     #ifndef MINTHRESH
28 gwlarson 3.2 #define MINTHRESH 5.0 /* source threshold w.r.t. mean */
29 gwlarson 3.1 #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 schorsch 3.13 RREAL radius; /* orig. object radius */
57 gwlarson 3.1 int listid; /* GL display list identifier */
58 gwlarson 3.10 int nlists; /* number of lists allocated */
59 schorsch 3.13 SUBPROC rtp; /* associated rtrace process */
60 gwlarson 3.1 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 gwlarson 3.2 #define getdcent(c,op) multp3(c,(op)->center,(op)->xfb.f.xfm)
74     #define getdrad(op) ((op)->radius*(op)->xfb.f.sca)
75    
76 gwlarson 3.1 #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 gwlarson 3.10 glDeleteLists(op->listid, op->nlists);
124 schorsch 3.13 close_process(&(op->rtp));
125 gwlarson 3.1 }
126     while (op->xfac)
127     freestr(op->xfav[--op->xfac]);
128 greg 3.11 free((void *)op);
129 gwlarson 3.1 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 gwlarson 3.2 ncells = 0;
287     goto done;
288 gwlarson 3.1 }
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 gwlarson 3.2 if (thresh <= FTINY) {
297     ncells = 0;
298     goto done;
299     }
300 gwlarson 3.1 /* avg. reflected brightness */
301     d = AVGREFL / (double)ncells;
302     scalecolor(csum, d);
303 greg 3.11 if (tmCvColors(&dlightsets->larb, TM_NOCHROM, &csum, 1) != TM_E_OK)
304 gwlarson 3.2 error(CONSISTENCY, "tone mapping problem in ssph_compute");
305 gwlarson 3.1 /* 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 gwlarson 3.6 ls->omega += 1.;
325 gwlarson 3.1 addcolor(ls->val, ssamp[alt][azi].val);
326 gwlarson 3.6 /* remove from list */
327 gwlarson 3.1 setcolor(ssamp[alt][azi].val, 0., 0., 0.);
328 gwlarson 3.6 ssamp[alt][azi].nsamp = 0;
329 gwlarson 3.1 }
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 gwlarson 3.2 done: /* clear sphere sample array */
343 schorsch 3.14 memset((void *)ssamp, '\0', sizeof(ssamp));
344 gwlarson 3.1 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 gwlarson 3.2 getdcent(ocent, op);
363 gwlarson 3.1 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 gwlarson 3.2 d2 = 2.*getdrad(op); d2 *= d2;
370 gwlarson 3.1 if ((dl = op->ol) != NULL && (mind2 < 0.0625*dl->ravg*dl->ravg ||
371 gwlarson 3.2 mind2 < 4.*getdrad(op)*getdrad(op)))
372 gwlarson 3.1 return(1);
373     if (!force)
374     return(0);
375     /* need to compute new light set */
376 schorsch 3.15 cvw = stdview;
377 gwlarson 3.1 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 greg 3.11 free((void *)dl);
411 gwlarson 3.1 return(0);
412     }
413     op->ol = dl;
414     return(1);
415     memerr:
416     error(SYSTEM, "out of memory in getdlights");
417     }
418    
419    
420 gwlarson 3.9 static
421 gwlarson 3.1 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 gwlarson 3.8 int somechange = 0;
436 gwlarson 3.1 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 gwlarson 3.9 else {
464     cmderror(cn, "need octree [name]");
465     return(0);
466     }
467 gwlarson 3.1 break;
468 gwlarson 3.5 case DO_UNLOAD: /* clear an object */
469 gwlarson 3.1 if (na > 1) goto toomany;
470     if (na && alist[0][0] == '*')
471 gwlarson 3.8 somechange += dobj_cleanup();
472 gwlarson 3.1 else
473 gwlarson 3.8 somechange += dobj_unload(na ? alist[0] : curname);
474 gwlarson 3.1 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 gwlarson 3.9 if (cn == DO_MOVE && nn >= na) {
483     cmderror(cn, "missing transform");
484     return(0);
485     }
486 gwlarson 3.8 somechange += dobj_xform(nm, cn==DO_MOVE, na-nn, alist+nn);
487 gwlarson 3.1 break;
488     case DO_UNMOVE: /* undo last transform */
489 gwlarson 3.8 somechange += dobj_unmove();
490 gwlarson 3.1 break;
491     case DO_OBJECT: /* print object statistics */
492 gwlarson 3.3 if (dobj_putstats(na ? alist[0] : curname, sstdout))
493 gwlarson 3.5 if (na && alist[0][0] != '*' && (curobj == NULL ||
494     strcmp(alist[0], curobj->name)))
495 gwlarson 3.3 savedxf(curobj = getdobj(alist[0]));
496 gwlarson 3.1 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 gwlarson 3.9 cmderror(cn, "need new object name");
504     return(0);
505 gwlarson 3.1 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 gwlarson 3.8 somechange += dobj_xform(curname, 1, na-nn, alist+nn);
518 gwlarson 3.1 else
519     curobj->drawcode = DO_HIDE;
520 gwlarson 3.3 savedxf(curobj);
521 gwlarson 3.1 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 gwlarson 3.8 somechange++;
528 gwlarson 3.1 break;
529     default:
530     error(CONSISTENCY, "bad command id in dobj_command");
531     }
532 gwlarson 3.9 return(somechange);
533 gwlarson 3.1 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 gwlarson 3.3 if (getdobj(nam) != NULL) {
557 gwlarson 3.5 error(COMMAND, "name already taken (clear first)");
558 gwlarson 3.3 return(0);
559     }
560 gwlarson 3.1 /* get octree path */
561 greg 3.12 if ((fpp = getpath(oct, getrlibpath(), R_OK)) == NULL) {
562 gwlarson 3.1 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 gwlarson 3.8 domats = 1;
580 gwlarson 3.10 op->listid = rgl_octlist(fpath, op->center, &op->radius, &op->nlists);
581 gwlarson 3.1 /* start rtrace */
582     rtargv[RTARGC-1] = fpath;
583     rtargv[RTARGC] = NULL;
584 schorsch 3.13 open_process(&(op->rtp), rtargv);
585 gwlarson 3.1 /* 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 greg 3.11 free((void *)lp);
618 gwlarson 3.1 }
619     return(1);
620     }
621    
622    
623 gwlarson 3.6 dobj_xform(nam, rel, ac, av) /* set/add transform for nam */
624 gwlarson 3.1 char *nam;
625 gwlarson 3.6 int rel, ac;
626 gwlarson 3.1 char **av;
627     {
628     register DOBJECT *op;
629 gwlarson 3.6 FVECT cent;
630     double rad;
631     char scoord[16];
632     int i;
633 gwlarson 3.1
634     if ((op = getdobj(nam)) == NULL) {
635     error(COMMAND, "no object");
636     return(0);
637     }
638 gwlarson 3.6 if (rel)
639     rel = op->xfac + 8;
640     if (ac + rel > MAXAC) {
641 gwlarson 3.1 error(COMMAND, "too many transform arguments");
642     return(0);
643     }
644 gwlarson 3.6 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 gwlarson 3.1 while (op->xfac)
655     freestr(op->xfav[--op->xfac]);
656 gwlarson 3.6 } 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 gwlarson 3.1 while (ac--)
664     op->xfav[op->xfac++] = savestr(*av++);
665 gwlarson 3.6 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 gwlarson 3.1 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 gwlarson 3.2 getdcent(ocent, op);
708 gwlarson 3.6 fprintf(fp, "%s: %s, center [%g %g %g], radius %g", op->name,
709 gwlarson 3.3 op->drawcode==DO_HIDE ? "hidden" :
710     op->drawcode==DO_LIGHT && op->ol!=NULL ? "lighted" :
711 gwlarson 3.1 "shown",
712 gwlarson 3.2 ocent[0],ocent[1],ocent[2], getdrad(op));
713 gwlarson 3.1 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 schorsch 3.14 memcpy((void *)txfav, (void *)lastxfav,
737 gwlarson 3.1 (txfac=lastxfac)*sizeof(char *));
738     /* save this transform */
739 schorsch 3.14 memcpy((void *)lastxfav, (void *)curobj->xfav,
740 gwlarson 3.1 (lastxfac=curobj->xfac)*sizeof(char *));
741     /* copy back last transform */
742 schorsch 3.14 memcpy((void *)curobj->xfav, (void *)txfav,
743 gwlarson 3.1 (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 gwlarson 3.3 if (getdobj(nam) != NULL) {
770 gwlarson 3.5 error(COMMAND, "name already taken (clear first)");
771 gwlarson 3.3 return(0);
772     }
773 gwlarson 3.1 /* allocate and copy struct */
774     opdup = (DOBJECT *)malloc(sizeof(DOBJECT));
775     if (opdup == NULL)
776     error(SYSTEM, "out of memory in dobj_dup");
777 schorsch 3.15 *opdup = *op;
778 gwlarson 3.1 /* 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 gwlarson 3.4 } else if ((op->drawcode = cn) == DO_LIGHT) {
813     if (!getdlights(op, 1))
814     error(COMMAND, "insufficient samples to light object");
815     } else
816 gwlarson 3.1 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 gwlarson 3.2 dobj_trace(nm, rorg, rdir) /* check for ray intersection with object(s) */
830     char nm[];
831 gwlarson 3.1 FVECT rorg, rdir;
832     {
833     register DOBJECT *op;
834     FVECT xorg, xdir;
835 gwlarson 3.2 double darr[6];
836 gwlarson 3.3 /* check each visible object? */
837 gwlarson 3.2 if (nm == NULL || *nm == '*') {
838     double dist, mindist = 1.01*FHUGE;
839 gwlarson 3.3
840     if (nm != NULL) nm[0] = '\0';
841 gwlarson 3.2 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 gwlarson 3.3 if (nm != NULL) strcpy(nm, op->name);
847     mindist = dist;
848 gwlarson 3.2 }
849 gwlarson 3.1 }
850 gwlarson 3.2 return(mindist);
851 gwlarson 3.1 }
852 gwlarson 3.2 /* else check particular object */
853     if ((op = getdobj(nm)) == NULL) {
854     error(COMMAND, "unknown object");
855     return(FHUGE);
856     }
857 gwlarson 3.3 if (op->xfac) { /* put ray in local coordinates */
858 gwlarson 3.2 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 schorsch 3.13 if (process(&(op->rtp), (char *)darr, (char *)darr, sizeof(double),
866 gwlarson 3.2 6*sizeof(double)) != sizeof(double))
867     error(SYSTEM, "rtrace communication error");
868     /* return distance */
869     if (darr[0] >= .99*FHUGE)
870     return(FHUGE);
871 gwlarson 3.3 return(darr[0]*op->xfb.f.sca);
872 gwlarson 3.1 }
873    
874    
875 gwlarson 3.7 int
876 gwlarson 3.1 dobj_render() /* render our objects in OpenGL */
877     {
878 gwlarson 3.7 int nrendered = 0;
879 gwlarson 3.1 GLboolean normalizing;
880     GLfloat vec[4];
881 gwlarson 3.2 FVECT v1;
882 gwlarson 3.1 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 gwlarson 3.7 return(0);
890 gwlarson 3.1 /* set up general rendering params */
891     glGetBooleanv(GL_NORMALIZE, &normalizing);
892 gwlarson 3.4 glPushAttrib(GL_LIGHTING_BIT|GL_TRANSFORM_BIT|GL_ENABLE_BIT|
893     GL_COLOR_BUFFER_BIT|GL_DEPTH_BUFFER_BIT|GL_POLYGON_BIT);
894 gwlarson 3.1 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 gwlarson 3.2 /* use computed sources */
914 gwlarson 3.1 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 gwlarson 3.3 } else { /* fake lighting */
940 gwlarson 3.2 vec[0] = vec[1] = vec[2] = 0.; vec[3] = 1.;
941 gwlarson 3.1 glLightModelfv(GL_LIGHT_MODEL_AMBIENT, vec);
942 gwlarson 3.2 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 gwlarson 3.1 }
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 gwlarson 3.7 nrendered++;
976 gwlarson 3.1 /* 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 gwlarson 3.2 else
988     glDisable(GL_LIGHT0);
989 gwlarson 3.1 /* check errors */
990     }
991     glPopAttrib(); /* restore rendering params */
992 gwlarson 3.3 rgl_checkerr("rendering objects in dobj_render");
993 gwlarson 3.7 return(nrendered);
994 gwlarson 3.1 }