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

Comparing ray/src/hd/rhdobj.c (file contents):
Revision 3.1 by gwlarson, Wed Aug 19 17:45:24 1998 UTC vs.
Revision 3.11 by greg, Sat Feb 22 02:07:24 2003 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1998 Silicon Graphics, Inc. */
2 –
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ SGI";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6 –
4   /*
5 < * Routines for loading and displaying Radiance objects under OpenGL in rholo.
5 > * Routines for loading and displaying Radiance objects in rholo with GLX.
6   */
7  
8   #include "radogl.h"
# Line 22 | Line 19 | int    (*dobj_lightsamp)() = NULL;     /* pointer to function
19  
20   #define AVGREFL         0.5             /* assumed average reflectance */
21  
22 < #define MAXAC           64              /* maximum number of args */
22 > #define MAXAC           512             /* maximum number of args */
23  
24   #ifndef MINTHRESH
25 < #define MINTHRESH       7.0             /* source threshold w.r.t. mean */
25 > #define MINTHRESH       5.0             /* source threshold w.r.t. mean */
26   #endif
27  
28   #ifndef NALT
# Line 55 | Line 52 | typedef struct dobject {
52          FVECT   center;                 /* orig. object center */
53          FLOAT   radius;                 /* orig. object radius */
54          int     listid;                 /* GL display list identifier */
55 +        int     nlists;                 /* number of lists allocated */
56          int     rtp[3];                 /* associated rtrace process */
57          DLIGHTS *ol;                    /* object lights */
58          FULLXF  xfb;                    /* coordinate transform */
# Line 69 | Line 67 | static DOBJECT *curobj;                /* current (last referred) ob
67   static int      lastxfac;               /* last number of transform args */
68   static char     *lastxfav[MAXAC+1];     /* saved transform arguments */
69  
70 + #define getdcent(c,op)  multp3(c,(op)->center,(op)->xfb.f.xfm)
71 + #define getdrad(op)     ((op)->radius*(op)->xfb.f.sca)
72 +
73   #define RTARGC  8
74   static char     *rtargv[RTARGC+1] = {"rtrace", "-h-", "-w-", "-fdd",
75                                          "-x", "1", "-oL"};
# Line 116 | Line 117 | register DOBJECT       *op;
117          }
118          dobjects = ohead.next;
119          if (!foundlink) {
120 <                glDeleteLists(op->listid, 1);
120 >                glDeleteLists(op->listid, op->nlists);
121                  close_process(op->rtp);
122          }
123          while (op->xfac)
124                  freestr(op->xfav[--op->xfac]);
125 <        free((char *)op);
125 >        free((void *)op);
126          return(1);
127   }
128  
# Line 279 | Line 280 | ssph_compute()                 /* compute source set from sphere sam
280                                  ncells++;
281                          }
282          if (dlightsets == NULL | ncells < NALT*NAZI/4) {
283 <                bzero((char *)ssamp, sizeof(ssamp));
284 <                return(0);
283 >                ncells = 0;
284 >                goto done;
285          }
286                                                  /* harmonic mean distance */
287          if (dlightsets->ravg > FTINY)
# Line 289 | Line 290 | ssph_compute()                 /* compute source set from sphere sam
290                  dlightsets->ravg = FHUGE;
291                                                  /* light source threshold */
292          thresh = MINTHRESH*bright(csum)/ncells;
293 +        if (thresh <= FTINY) {
294 +                ncells = 0;
295 +                goto done;
296 +        }
297                                                  /* avg. reflected brightness */
298          d = AVGREFL / (double)ncells;  
299          scalecolor(csum, d);
300 <        if (tmCvColors(&dlightsets->larb, TM_NOCHROM, csum, 1) != TM_E_OK)
301 <                error(CONSISTENCY, "bad tone mapping in ssph_compute");
300 >        if (tmCvColors(&dlightsets->larb, TM_NOCHROM, &csum, 1) != TM_E_OK)
301 >                error(CONSISTENCY, "tone mapping problem in ssph_compute");
302                                          /* greedy light source clustering */
303          while (dlightsets->nl < MAXLIGHTS) {
304                  maxbr = 0.;                     /* find brightest cell */
# Line 313 | Line 318 | ssph_compute()                 /* compute source set from sphere sam
318                                  continue;               /* too dim */
319                          ssph_direc(v, alt, azi);        /* else add it in */
320                          VSUM(ls->direc, ls->direc, v, d);
321 <                        ls->omega++;
321 >                        ls->omega += 1.;
322                          addcolor(ls->val, ssamp[alt][azi].val);
323 +                                                        /* remove from list */
324                          setcolor(ssamp[alt][azi].val, 0., 0., 0.);
325 +                        ssamp[alt][azi].nsamp = 0;
326                  }
327                  d = 1./ls->omega;                       /* avg. brightness */
328                  scalecolor(ls->val, d);
# Line 329 | Line 336 | ssph_compute()                 /* compute source set from sphere sam
336                                  addcolor(dlightsets->lamb, ssamp[alt][azi].val);
337          d = 1.0/ncells;
338          scalecolor(dlightsets->lamb, d);
339 <                                        /* clear sphere sample array */
339 > done:                                   /* clear sphere sample array */
340          bzero((char *)ssamp, sizeof(ssamp));
341          return(ncells);
342   }
# Line 349 | Line 356 | int    force;
356          if (op->drawcode != DO_LIGHT)
357                  return(0);
358                                          /* check for usable light set */
359 <        multp3(ocent, op->center, op->xfb.f.xfm);
359 >        getdcent(ocent, op);
360          for (dl = dlightsets; dl != NULL; dl = dl->next)
361                  if ((d2 = dist2(dl->lcent, ocent)) < mind2) {
362                          op->ol = dl;
363                          mind2 = d2;
364                  }
365                                          /* the following is heuristic */
366 <        d2 = 2.*op->radius*op->xfb.f.sca; d2 *= d2;
366 >        d2 = 2.*getdrad(op); d2 *= d2;
367          if ((dl = op->ol) != NULL && (mind2 < 0.0625*dl->ravg*dl->ravg ||
368 <                mind2 < 4.*op->radius*op->xfb.f.sca*op->radius*op->xfb.f.sca))
368 >                        mind2 < 4.*getdrad(op)*getdrad(op)))
369                  return(1);
370          if (!force)
371                  return(0);
# Line 397 | Line 404 | int    force;
404                          quit(0);
405          if (!ssph_compute()) {          /* compute light sources from sphere */
406                  dlightsets = dl->next;
407 <                free((char *)dl);
407 >                free((void *)dl);
408                  return(0);
409          }
410          op->ol = dl;
# Line 407 | Line 414 | memerr:
414   }
415  
416  
417 < static int
417 > static
418   cmderror(cn, err)               /* report command error */
419   int     cn;
420   char    *err;
421   {
422          sprintf(errmsg, "%s: %s", rhdcmd[cn], err);
423          error(COMMAND, errmsg);
417 –        return(cn);
424   }
425  
426  
# Line 423 | Line 429 | dobj_command(cmd, args)                /* run object display command
429   char    *cmd;
430   register char   *args;
431   {
432 +        int     somechange = 0;
433          int     cn, na, doxfm;
434          register int    nn;
435          char    *alist[MAXAC+1], *nm;
# Line 450 | Line 457 | register char  *args;
457                          dobj_load(alist[0], alist[0]);
458                  else if (na == 2)
459                          dobj_load(alist[0], alist[1]);
460 <                else
461 <                        return(cmderror(cn, "need octree [name]"));
460 >                else {
461 >                        cmderror(cn, "need octree [name]");
462 >                        return(0);
463 >                }
464                  break;
465 <        case DO_UNLOAD:                         /* unload an object */
465 >        case DO_UNLOAD:                         /* clear an object */
466                  if (na > 1) goto toomany;
467                  if (na && alist[0][0] == '*')
468 <                        dobj_cleanup();
468 >                        somechange += dobj_cleanup();
469                  else
470 <                        dobj_unload(na ? alist[0] : curname);
470 >                        somechange += dobj_unload(na ? alist[0] : curname);
471                  break;
472          case DO_XFORM:                          /* transform object */
473          case DO_MOVE:
# Line 467 | Line 476 | register char  *args;
476                  } else {
477                          nm = curname; nn = 0;
478                  }
479 <                if (cn == DO_MOVE && nn >= na)
480 <                        return(cmderror(cn, "missing transform"));
481 <                dobj_xform(nm, cn==DO_MOVE, na-nn, alist+nn);
479 >                if (cn == DO_MOVE && nn >= na) {
480 >                        cmderror(cn, "missing transform");
481 >                        return(0);
482 >                }
483 >                somechange += dobj_xform(nm, cn==DO_MOVE, na-nn, alist+nn);
484                  break;
485          case DO_UNMOVE:                         /* undo last transform */
486 <                dobj_unmove();
486 >                somechange += dobj_unmove();
487                  break;
488          case DO_OBJECT:                         /* print object statistics */
489 <                dobj_putstats(na ? alist[0] : curname, sstdout);
489 >                if (dobj_putstats(na ? alist[0] : curname, sstdout))
490 >                        if (na && alist[0][0] != '*' && (curobj == NULL ||
491 >                                        strcmp(alist[0], curobj->name)))
492 >                                savedxf(curobj = getdobj(alist[0]));
493                  break;
494          case DO_DUP:                            /* duplicate object */
495                  for (nn = 0; nn < na; nn++)
# Line 483 | Line 497 | register char  *args;
497                                  break;
498                  switch (nn) {
499                  case 0:
500 <                        return(cmderror(cn, "need new object name"));
500 >                        cmderror(cn, "need new object name");
501 >                        return(0);
502                  case 1:
503                          nm = curname;
504                          break;
# Line 496 | Line 511 | register char  *args;
511                  if (!dobj_dup(nm, alist[nn-1]))
512                          break;
513                  if (na > nn)
514 <                        dobj_xform(curname, 1, na-nn, alist+nn);
514 >                        somechange += dobj_xform(curname, 1, na-nn, alist+nn);
515                  else
516                          curobj->drawcode = DO_HIDE;
517 +                savedxf(curobj);
518                  break;
519          case DO_SHOW:                           /* change rendering option */
520          case DO_LIGHT:
521          case DO_HIDE:
522                  if (na > 1) goto toomany;
523                  dobj_lighting(na ? alist[0] : curname, cn);
524 +                somechange++;
525                  break;
526          default:
527                  error(CONSISTENCY, "bad command id in dobj_command");
528          }
529 <        dev_view(&odev.v);                      /* redraw */
513 <        return(cn);
529 >        return(somechange);
530   toomany:
531          return(cmderror(cn, "too many arguments"));
532   }
# Line 535 | Line 551 | char   *oct, *nam;
551                  error(COMMAND, "illegal name");
552                  return(0);
553          }
554 +        if (getdobj(nam) != NULL) {
555 +                error(COMMAND, "name already taken (clear first)");
556 +                return(0);
557 +        }
558                                          /* get octree path */
559          if ((fpp = getpath(oct, getlibpath(), R_OK)) == NULL) {
560                  sprintf(errmsg, "cannot find octree \"%s\"", oct);
# Line 542 | Line 562 | char   *oct, *nam;
562                  return(0);
563          }
564          strcpy(fpath, fpp);
545 –        freedobj(getdobj(nam));         /* free previous use of nam */
565          op = (DOBJECT *)malloc(sizeof(DOBJECT));
566          if (op == NULL)
567                  error(SYSTEM, "out of memory in dobj_load");
# Line 555 | Line 574 | char   *oct, *nam;
574          op->xfav[op->xfac=0] = NULL;
575                                          /* load octree into display list */
576          dolights = 0;
577 <        op->listid = rgl_octlist(fpath, op->center, &op->radius);
577 >        domats = 1;
578 >        op->listid = rgl_octlist(fpath, op->center, &op->radius, &op->nlists);
579                                          /* start rtrace */
580          rtargv[RTARGC-1] = fpath;
581          rtargv[RTARGC] = NULL;
# Line 592 | Line 612 | dobj_cleanup()                         /* free all resources */
612          savedxf(curobj = NULL);
613          while ((lp = dlightsets) != NULL) {
614                  dlightsets = lp->next;
615 <                free((char *)lp);
615 >                free((void *)lp);
616          }
617          return(1);
618   }
619  
620  
621 < dobj_xform(nam, add, ac, av)            /* set/add transform for nam */
621 > dobj_xform(nam, rel, ac, av)            /* set/add transform for nam */
622   char    *nam;
623 < int     add, ac;
623 > int     rel, ac;
624   char    **av;
625   {
626          register DOBJECT        *op;
627 +        FVECT   cent;
628 +        double  rad;
629 +        char    scoord[16];
630 +        int     i;
631  
632          if ((op = getdobj(nam)) == NULL) {
633                  error(COMMAND, "no object");
634                  return(0);
635          }
636 <        if (add) add = op->xfac;
637 <        if (ac + add > MAXAC) {
636 >        if (rel)
637 >                rel = op->xfac + 8;
638 >        if (ac + rel > MAXAC) {
639                  error(COMMAND, "too many transform arguments");
640                  return(0);
641          }
642 <        savedxf(curobj = op);
643 <        if (!add)
642 >        savedxf(curobj = op);           /* remember current transform */
643 >        if (rel && ac == 4 && !strcmp(av[0], "-t"))
644 >                rel = -1;                       /* don't move for translate */
645 >        else {
646 >                getdcent(cent, op);             /* don't move if near orig. */
647 >                rad = getdrad(op);
648 >                if (DOT(cent,cent) < rad*rad)
649 >                        rel = -1;
650 >        }
651 >        if (!rel) {                             /* remove old transform */
652                  while (op->xfac)
653                          freestr(op->xfav[--op->xfac]);
654 +        } else if (rel > 0) {                   /* relative move */
655 +                op->xfav[op->xfac++] = savestr("-t");
656 +                for (i = 0; i < 3; i++) {
657 +                        sprintf(scoord, "%.4e", -cent[i]);
658 +                        op->xfav[op->xfac++] = savestr(scoord);
659 +                }
660 +        }
661          while (ac--)
662                  op->xfav[op->xfac++] = savestr(*av++);
663 +        if (rel > 0) {                          /* move back */
664 +                op->xfav[op->xfac++] = savestr("-t");
665 +                for (i = 0; i < 3; i++) {
666 +                        sprintf(scoord, "%.4e", cent[i]);
667 +                        op->xfav[op->xfac++] = savestr(scoord);
668 +                }
669 +        }
670          op->xfav[op->xfac] = NULL;
671          if (fullxf(&op->xfb, op->xfac, op->xfav) != op->xfac) {
672                  error(COMMAND, "bad transform arguments");
# Line 655 | Line 702 | FILE   *fp;
702                  error(COMMAND, "unknown object");
703                  return(0);
704          }
705 <        multp3(ocent, op->center, op->xfb.f.xfm);
706 <        fprintf(fp, "%s: %s, center [%f %f %f], radius %f", op->name,
707 <                        op->drawcode==DO_HIDE ? "hid" :
708 <                        op->drawcode==DO_LIGHT && op->ol!=NULL ? "lit" :
705 >        getdcent(ocent, op);
706 >        fprintf(fp, "%s: %s, center [%g %g %g], radius %g", op->name,
707 >                        op->drawcode==DO_HIDE ? "hidden" :
708 >                        op->drawcode==DO_LIGHT && op->ol!=NULL ? "lighted" :
709                          "shown",
710 <                        ocent[0],ocent[1],ocent[2], op->radius*op->xfb.f.sca);
710 >                        ocent[0],ocent[1],ocent[2], getdrad(op));
711          if (op->xfac)
712                  fputs(", (xform", fp);
713          for (i = 0; i < op->xfac; i++) {
# Line 717 | Line 764 | char   *oldnm, *nam;
764                  error(COMMAND, "illegal name");
765                  return(0);
766          }
767 +        if (getdobj(nam) != NULL) {
768 +                error(COMMAND, "name already taken (clear first)");
769 +                return(0);
770 +        }
771                                          /* allocate and copy struct */
772          opdup = (DOBJECT *)malloc(sizeof(DOBJECT));
773          if (opdup == NULL)
# Line 756 | Line 807 | int    cn;
807          } else if ((op = getdobj(nam)) == NULL) {
808                  error(COMMAND, "unknown object");
809                  return(0);
810 <        } else if ((op->drawcode = cn) == DO_LIGHT)
811 <                getdlights(op, 1);
812 <        else
810 >        } else if ((op->drawcode = cn) == DO_LIGHT) {
811 >                if (!getdlights(op, 1))
812 >                        error(COMMAND, "insufficient samples to light object");
813 >        } else
814                  op->ol = NULL;
815  
816          if (dobj_lightsamp != NULL) {           /* restore beam set */
# Line 772 | Line 824 | int    cn;
824  
825  
826   double
827 < dobj_trace(rorg, rdir)          /* check for ray intersection with objects */
827 > dobj_trace(nm, rorg, rdir)      /* check for ray intersection with object(s) */
828 > char    nm[];
829   FVECT   rorg, rdir;
830   {
831          register DOBJECT        *op;
832          FVECT   xorg, xdir;
833 <        double  darr[6], mindist = FHUGE;
834 <                                        /* check each visible object */
835 <        for (op = dobjects; op != NULL; op = op->next) {
836 <                if (op->drawcode == DO_HIDE)
837 <                        continue;
838 <                if (op->xfac) {         /* transform ray */
839 <                        multp3(xorg, rorg, op->xfb.b.xfm);
840 <                        multv3(xdir, rdir, op->xfb.b.xfm);
841 <                        VCOPY(darr, xorg); VCOPY(darr+3, xdir);
842 <                } else {
843 <                        VCOPY(darr, rorg); VCOPY(darr+3, rdir);
833 >        double  darr[6];
834 >                                        /* check each visible object? */
835 >        if (nm == NULL || *nm == '*') {
836 >                double  dist, mindist = 1.01*FHUGE;
837 >
838 >                if (nm != NULL) nm[0] = '\0';
839 >                for (op = dobjects; op != NULL; op = op->next) {
840 >                        if (op->drawcode == DO_HIDE)
841 >                                continue;
842 >                        dist = dobj_trace(op->name, rorg, rdir);
843 >                        if (dist < mindist) {
844 >                                if (nm != NULL) strcpy(nm, op->name);
845 >                                mindist = dist;
846 >                        }
847                  }
848 <                                        /* trace it */
793 <                if (process(op->rtp, darr, darr, sizeof(double),
794 <                                6*sizeof(double)) != sizeof(double))
795 <                        error(SYSTEM, "rtrace communication error");
796 <                                        /* get closest */
797 <                if ((darr[0] *= op->xfb.f.sca) < mindist)
798 <                        mindist = darr[0];
848 >                return(mindist);
849          }
850 <        return(mindist);
850 >                                        /* else check particular object */
851 >        if ((op = getdobj(nm)) == NULL) {
852 >                error(COMMAND, "unknown object");
853 >                return(FHUGE);
854 >        }
855 >        if (op->xfac) {         /* put ray in local coordinates */
856 >                multp3(xorg, rorg, op->xfb.b.xfm);
857 >                multv3(xdir, rdir, op->xfb.b.xfm);
858 >                VCOPY(darr, xorg); VCOPY(darr+3, xdir);
859 >        } else {
860 >                VCOPY(darr, rorg); VCOPY(darr+3, rdir);
861 >        }
862 >                                /* trace it */
863 >        if (process(op->rtp, (char *)darr, (char *)darr, sizeof(double),
864 >                        6*sizeof(double)) != sizeof(double))
865 >                error(SYSTEM, "rtrace communication error");
866 >                                /* return distance */
867 >        if (darr[0] >= .99*FHUGE)
868 >                return(FHUGE);
869 >        return(darr[0]*op->xfb.f.sca);
870   }
871  
872  
873 + int
874   dobj_render()                   /* render our objects in OpenGL */
875   {
876 +        int     nrendered = 0;
877          GLboolean       normalizing;
878          GLfloat vec[4];
879 +        FVECT   v1;
880          register DOBJECT        *op;
881          register int    i;
882                                          /* anything to render? */
# Line 812 | Line 884 | dobj_render()                  /* render our objects in OpenGL */
884                  if (op->drawcode != DO_HIDE)
885                          break;
886          if (op == NULL)
887 <                return(1);
887 >                return(0);
888                                          /* set up general rendering params */
889          glGetBooleanv(GL_NORMALIZE, &normalizing);
890 <        glPushAttrib(GL_LIGHTING_BIT|GL_TRANSFORM_BIT|
891 <                        GL_DEPTH_BUFFER_BIT|GL_POLYGON_BIT);
890 >        glPushAttrib(GL_LIGHTING_BIT|GL_TRANSFORM_BIT|GL_ENABLE_BIT|
891 >                GL_COLOR_BUFFER_BIT|GL_DEPTH_BUFFER_BIT|GL_POLYGON_BIT);
892          glDepthFunc(GL_LESS);
893          glEnable(GL_DEPTH_TEST);
894          glLightModeli(GL_LIGHT_MODEL_TWO_SIDE, GL_TRUE);
# Line 836 | Line 908 | dobj_render()                  /* render our objects in OpenGL */
908                  if (op->drawcode == DO_LIGHT && op->ol != NULL) {
909                          BYTE    pval;
910                          double  expval, d;
911 <
911 >                                                /* use computed sources */
912                          if (tmMapPixels(&pval, &op->ol->larb, TM_NOCHROM, 1)
913                                          != TM_E_OK)
914                                  error(CONSISTENCY, "dobj_render w/o tone map");
# Line 862 | Line 934 | dobj_render()                  /* render our objects in OpenGL */
934                                  glLightfv(glightid[i], GL_AMBIENT, vec);
935                                  glEnable(glightid[i]);
936                          }
937 <                } else {
938 <                        vec[0] = vec[1] = vec[2] = 1.; vec[3] = 1.;
937 >                } else {                        /* fake lighting */
938 >                        vec[0] = vec[1] = vec[2] = 0.; vec[3] = 1.;
939                          glLightModelfv(GL_LIGHT_MODEL_AMBIENT, vec);
940 +                        getdcent(v1, op);
941 +                        VSUB(v1, odev.v.vp, v1);
942 +                        if (normalize(v1) <= getdrad(op)) {
943 +                                vec[0] = -odev.v.vdir[0];
944 +                                vec[1] = -odev.v.vdir[1];
945 +                                vec[2] = -odev.v.vdir[2];
946 +                        } else
947 +                                VCOPY(vec, v1);
948 +                        vec[3] = 0.;
949 +                        glLightfv(GL_LIGHT0, GL_POSITION, vec);
950 +                        vec[0] = vec[1] = vec[2] = .7; vec[3] = 1.;
951 +                        glLightfv(GL_LIGHT0, GL_SPECULAR, vec);
952 +                        glLightfv(GL_LIGHT0, GL_DIFFUSE, vec);
953 +                        vec[0] = vec[1] = vec[2] = .3; vec[3] = 1.;
954 +                        glLightfv(GL_LIGHT0, GL_AMBIENT, vec);
955 +                        glEnable(GL_LIGHT0);
956                  }
957                                          /* set up object transform */
958                  if (op->xfac) {
# Line 882 | Line 970 | dobj_render()                  /* render our objects in OpenGL */
970                  }
971                                          /* render the display list */
972                  glCallList(op->listid);
973 +                nrendered++;
974                                          /* restore matrix */
975                  if (op->xfac) {
976                          glMatrixMode(GL_MODELVIEW);
# Line 893 | Line 982 | dobj_render()                  /* render our objects in OpenGL */
982                  if (op->drawcode == DO_LIGHT && op->ol != NULL)
983                          for (i = op->ol->nl; i--; )
984                                  glDisable(glightid[i]);
985 +                else
986 +                        glDisable(GL_LIGHT0);
987                                          /* check errors */
897 –                rgl_checkerr("rendering object in dobj_render");
988          }
989          glPopAttrib();                  /* restore rendering params */
990 <        return(1);
990 >        rgl_checkerr("rendering objects in dobj_render");
991 >        return(nrendered);
992   }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines