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

Comparing ray/src/hd/sm_ogl.c (file contents):
Revision 3.4 by gwlarson, Thu Sep 3 11:29:00 1998 UTC vs.
Revision 3.10 by gwlarson, Mon Dec 28 18:07:35 1998 UTC

# Line 16 | Line 16 | static char SCCSid[] = "$SunId$ SGI";
16   #include <GL/glu.h>
17   #include <glut.h>
18   #endif
19 < #include "object.h"
19 > #include "sm_flag.h"
20   #include "sm_list.h"
21   #include "sm_geom.h"
22 + #include "sm_qtree.h"
23 + #include "sm_stree.h"
24   #include "sm.h"
25  
26   #ifdef TEST_DRIVER
# Line 41 | Line 43 | typedef struct {
43          QTRAVG          av;     /* node average */
44   } QT_LUENT;             /* lookup table entry */
45  
46 < static QT_LUENT *qt_htbl = NULL;        /* quadtree hash table */
47 < static int      qt_hsiz = 0;            /* quadtree hash table size */
46 > static QT_LUENT *qt_htbl = NULL;        /* quadtree cache */
47 > static int      qt_hsiz = 0;            /* quadtree cache size */
48  
49  
50 < int
51 < mark_active_tris(qtptr,arg)
52 < QUADTREE *qtptr;
53 < char *arg;
54 < {
55 <  QUADTREE qt = *qtptr;
54 <  OBJECT os[QT_MAX_SET+1],*optr;
55 <  register int i,t_id;
56 <
57 <  if (!QT_IS_LEAF(qt))
58 <    return(TRUE);
59 <  /* For each triangle in the set, set the which flag*/
60 <  qtgetset(os,qt);
61 <
62 <  for (i = QT_SET_CNT(os), optr = QT_SET_PTR(os); i > 0; i--)
63 <  {
64 <    t_id = QT_SET_NEXT_ELEM(optr);
65 <    /* Set the render flag */
66 <    if(SM_IS_NTH_T_BASE(smMesh,t_id))
67 <        continue;
68 <    SM_SET_NTH_T_ACTIVE(smMesh,t_id);
69 <    /* FOR NOW:Also set the LRU clock bit: MAY WANT TO CHANGE: */      
70 <    SM_SET_NTH_T_LRU(smMesh,t_id);
71 <  }
72 <  return(TRUE);
73 < }
74 <
75 < mark_tris_in_frustum(view)
76 < VIEW *view;
77 < {
78 <    FVECT nr[4],far[4];
79 <
80 <    /* Mark triangles in approx. view frustum as being active:set
81 <       LRU counter: for use in discarding samples when out
82 <       of space
83 <       Radiance often has no far clipping plane: but driver will set
84 <       dev_zmin,dev_zmax to satisfy OGL
85 <    */
86 <
87 <    /* First clear all the quadtree node and triangle active flags */
88 <    qtClearAllFlags();
89 <    smClear_flags(smMesh,T_ACTIVE_FLAG);
90 <
91 <    /* calculate the world space coordinates of the view frustum */
92 <    calculate_view_frustum(view->vp,view->hvec,view->vvec,view->horiz,
93 <                           view->vert, dev_zmin,dev_zmax,nr,far);
94 <
95 <    /* Project the view frustum onto the spherical quadtree */
96 <    /* For every cell intersected by the projection of the faces
97 <       of the frustum: mark all triangles in the cell as ACTIVE-
98 <       Also set the triangles LRU clock counter
99 <       */
100 <    /* Near face triangles */
101 <    smLocator_apply_func(smMesh,nr[0],nr[2],nr[3],mark_active_tris,NULL);
102 <    smLocator_apply_func(smMesh,nr[2],nr[0],nr[1],mark_active_tris,NULL);
103 <    /* Right face triangles */
104 <    smLocator_apply_func(smMesh,nr[0],far[3],far[0],mark_active_tris,NULL);
105 <    smLocator_apply_func(smMesh,far[3],nr[0],nr[3],mark_active_tris,NULL);
106 <    /* Left face triangles */
107 <    smLocator_apply_func(smMesh,nr[1],far[2],nr[2],mark_active_tris,NULL);
108 <    smLocator_apply_func(smMesh,far[2],nr[1],far[1],mark_active_tris,NULL);
109 <    /* Top face triangles */
110 <    smLocator_apply_func(smMesh,nr[0],far[0],nr[1],mark_active_tris,NULL);
111 <    smLocator_apply_func(smMesh,nr[1],far[0],far[1],mark_active_tris,NULL);
112 <    /* Bottom face triangles */
113 <    smLocator_apply_func(smMesh,nr[3],nr[2],far[3],mark_active_tris,NULL);
114 <    smLocator_apply_func(smMesh,nr[2],far[2],far[3],mark_active_tris,NULL);
115 <    /* Far face triangles */
116 <    smLocator_apply_func(smMesh,far[0],far[2],far[1],mark_active_tris,NULL);
117 <    smLocator_apply_func(smMesh,far[2],far[0],far[3],mark_active_tris,NULL);
118 <
119 < #ifdef TEST_DRIVER
120 <    VCOPY(FrustumFar[0],far[0]);
121 <    VCOPY(FrustumFar[1],far[1]);
122 <    VCOPY(FrustumFar[2],far[2]);
123 <    VCOPY(FrustumFar[3],far[3]);
124 <    VCOPY(FrustumNear[0],nr[0]);
125 <    VCOPY(FrustumNear[1],nr[1]);
126 <    VCOPY(FrustumNear[2],nr[2]);
127 <    VCOPY(FrustumNear[3],nr[3]);
128 < #endif
129 < }
130 <
131 < /*
132 < * smClean()            : display has been wiped clean
133 < *
134 < * Called after display has been effectively cleared, meaning that all
135 < * geometry must be resent down the pipeline in the next call to smUpdate().
136 < */
50 > /*
51 >  * smClean()           : display has been wiped clean
52 >  *
53 >  * Called after display has been effectively cleared, meaning that all
54 >  * geometry must be resent down the pipeline in the next call to smUpdate().
55 >  */
56   smClean()
57   {
58      smClean_notify = TRUE;
59   }
60  
61   int
62 < qtHash_init(nel)                /* initialize for at least nel elements */
62 > qtCache_init(nel)               /* initialize for at least nel elements */
63   int     nel;
64   {
65          static int  hsiztab[] = {
# Line 171 | Line 90 | int    nel;
90   }
91  
92   QT_LUENT *
93 < qtHash_find(qt)                 /* find a quadtree table entry */
93 > qtCache_find(qt)                /* find a quadtree table entry */
94   QUADTREE qt;
95   {
96          int     i, n;
97          register int    ndx;
98          register QT_LUENT       *le;
99  
100 <        if (qt_hsiz == 0)
101 <                qtHash_init(1);
100 >        if (qt_hsiz == 0 && !qtCache_init(1))
101 >                return(NULL);
102   tryagain:                               /* hash table lookup */
103          ndx = (unsigned long)qt % qt_hsiz;
104          for (i = 0, n = 1; i < qt_hsiz; i++, n += 2) {
# Line 192 | Line 111 | tryagain:                              /* hash table lookup */
111                                          /* table is full, reallocate */
112          le = qt_htbl;
113          ndx = qt_hsiz;
114 <        if (!qtHash_init(ndx+1)) {      /* no more memory! */
114 >        if (!qtCache_init(ndx+1)) {     /* no more memory! */
115                  qt_htbl = le;
116                  qt_hsiz = ndx;
117                  return(NULL);
118          }
119 <        if (!ndx)
201 <                goto tryagain;
202 <                                        /* copy old table to new */
119 >                                        /* copy old table to new and free */
120          while (ndx--)
121                  if (!QT_IS_EMPTY(le[ndx].qt))
122 <                        copystruct(qtHash_find(le[ndx].qt), &le[ndx]);
122 >                        copystruct(qtCache_find(le[ndx].qt), &le[ndx]);
123          free((char *)le);
124          goto tryagain;                  /* should happen only once! */
125   }
# Line 222 | Line 139 | register QUADTREE qt;
139      stCount_level_leaves(lcnt+1, QT_NTH_CHILD(qt,3));
140    }
141    else
142 <    lcnt[0]++;
142 >    if(QT_LEAF_IS_FLAG(qt))
143 >      lcnt[0]++;
144   }
145  
146  
# Line 236 | Line 154 | int lvl;
154    FVECT a,b,c;
155    register QT_LUENT *le;
156    QTRAVG *rc[4];
157 +  TRI *tri;
158    
159    if (QT_IS_EMPTY(qt))                          /* empty leaf node */
160      return(NULL);
161    if (QT_IS_TREE(qt) && !QT_IS_FLAG(qt))        /* not in our frustum */
162      return(NULL);
163 +  if(QT_IS_LEAF(qt)  && !QT_LEAF_IS_FLAG(qt))   /* not in our frustum */
164 +    return(NULL);
165                                          /* else look up node */
166 <  if ((le = qtHash_find(qt)) == NULL)
166 >  if ((le = qtCache_find(qt)) == NULL)
167      error(SYSTEM, "out of memory in qtRender_level");
168    if (QT_IS_TREE(qt) && (QT_IS_EMPTY(le->qt) || lvl > 0))
169    {                                     /* compute children */
# Line 271 | Line 192 | int lvl;
192      }
193      else
194      {                                   /* from triangle set */
195 <      OBJECT os[QT_MAX_SET+1];
195 >      OBJECT *os;
196        int s0, s1, s2;
197  
198 <      qtgetset(os,qt);
198 >      os = qtqueryset(qt);
199        for (n = os[0]; n; n--)
200        {
201 <        qtTri_from_id(os[n],a,b,c,NULL,NULL,NULL,&s0,&s1,&s2);
201 >        if(SM_IS_NTH_T_BASE(sm,os[n]))
202 >           continue;
203 >        tri = SM_NTH_TRI(sm,os[n]);
204 >        if(!T_IS_VALID(tri))
205 >          continue;
206 >        
207 >        s0 = T_NTH_V(tri,0);
208 >        s1 = T_NTH_V(tri,1);
209 >        s2 = T_NTH_V(tri,2);
210 >        VCOPY(a,SM_NTH_WV(sm,s0));
211 >        VCOPY(b,SM_NTH_WV(sm,s1));
212 >        VCOPY(c,SM_NTH_WV(sm,s2));            
213          distsum += SM_BG_SAMPLE(sm,s0) ? dev_zmax
214                                  : sqrt(dist2(a,SM_VIEW_CENTER(sm)));
215          distsum += SM_BG_SAMPLE(sm,s1) ? dev_zmax
# Line 309 | Line 241 | int lvl;
241      VSUM(c,SM_VIEW_CENTER(sm),c,le->av.dist);
242                                          /* draw triangle */
243      glColor3ub(le->av.rgb[0],le->av.rgb[1],le->av.rgb[2]);
312 –    /* NOTE: Triangle vertex order may change */
313 –    glVertex3d(c[0],c[1],c[2]);
314 –    glVertex3d(b[0],b[1],b[2]);
244      glVertex3d(a[0],a[1],a[2]);
245 +    glVertex3d(b[0],b[1],b[2]);
246 +    glVertex3d(c[0],c[1],c[2]);
247 +
248    }
249    return(&le->av);
250   }
# Line 322 | Line 254 | smRender_stree_level(sm,lvl)
254   SM *sm;
255   int lvl;
256   {
257 <  QUADTREE root;
257 >  QUADTREE qt;
258    int i;
259    FVECT t0,t1,t2;
260 +  STREE *st;
261  
262 +  
263    if (lvl < 1)
264      return;
265 +  st = SM_LOCATOR(sm);
266    glPushAttrib(GL_LIGHTING_BIT);
267    glShadeModel(GL_FLAT);
268    glBegin(GL_TRIANGLES);
269 <  for(i=0; i < 4; i++)
269 >  for(i=0; i < ST_NUM_ROOT_NODES; i++)
270    {
271 <    root = ST_NTH_ROOT(SM_LOCATOR(sm),i);
272 <    stNth_base_verts(SM_LOCATOR(sm),i,t0,t1,t2);
273 <    qtRender_level(root,t0,t1,t2,sm,lvl-1);
271 >    qt = ST_ROOT_QT(st,i);
272 >    qtRender_level(qt,ST_NTH_V(st,i,0),ST_NTH_V(st,i,1),ST_NTH_V(st,i,2),
273 >                   sm,lvl-1);
274    }
275    glEnd();
276    glPopAttrib();
# Line 348 | Line 283 | int qual;
283   {
284    int i, ntarget;
285    int lvlcnt[QT_MAX_LEVELS];
286 +  STREE *st;
287  
288    if (qual <= 0)
289      return;
# Line 360 | Line 296 | int qual;
296      return;
297    for (i = QT_MAX_LEVELS; i--; )
298      lvlcnt[i] = 0;
299 <  stCount_level_leaves(lvlcnt, SM_LOCATOR(sm)->root);
299 >  
300 >  st = SM_LOCATOR(sm);
301 >  for(i=0; i < ST_NUM_ROOT_NODES;i++)
302 >    stCount_level_leaves(lvlcnt, ST_ROOT_QT(st,i));
303 >  
304    for (ntarget -= lvlcnt[i=0]; i < QT_MAX_LEVELS-1; ntarget -= lvlcnt[++i])
305      if (ntarget < lvlcnt[i+1])
306        break;
# Line 380 | Line 320 | int clr;
320    int j;
321  
322    tri = SM_NTH_TRI(sm,i);
323 <  if (clr) SM_CLEAR_NTH_T_NEW(sm,i);
324 <
385 <  /* NOTE:Triangles are defined clockwise:historical relative to spherical
386 <     tris: could change
387 <     */
388 <  for(j=2; j>= 0; j--)
323 >  if (clr) SM_CLR_NTH_T_NEW(sm,i);
324 >  for(j=0; j <= 2; j++)
325    {
326   #ifdef DEBUG
327      if(SM_BG_SAMPLE(sm,T_NTH_V(tri,j)))
# Line 411 | Line 347 | int clr;
347    int rgb[3];
348  
349    tri = SM_NTH_TRI(sm,i);
350 <  if (clr) SM_CLEAR_NTH_T_NEW(sm,i);
350 >  if (clr) SM_CLR_NTH_T_NEW(sm,i);
351  
352    /* NOTE:Triangles are defined clockwise:historical relative to spherical
353       tris: could change
# Line 436 | Line 372 | int clr;
372      rgb[0]/=cnt; rgb[1]/=cnt; rgb[2]/=cnt;
373      d /= (double)cnt;
374    }
375 <  for(j=2; j>= 0; j--)
375 >  for(j=0; j <= 2; j++)
376    {
377      if(SM_BG_SAMPLE(sm,ids[j]))
378      {
# Line 469 | Line 405 | int clr;
405    TRI *tri;
406    
407    tri = SM_NTH_TRI(sm,i);
408 <  if (clr) SM_CLEAR_NTH_T_NEW(sm,i);
408 >  if (clr) SM_CLR_NTH_T_NEW(sm,i);
409  
410    /* NOTE:Triangles are defined clockwise:historical relative to spherical
411       tris: could change
412       */
413 <  for(j=2; j>= 0; j--)
413 >  for(j=0; j <= 2; j++)
414    {
415        id = T_NTH_V(tri,j);
416        glColor3ub(SM_NTH_RGB(sm,id)[0],SM_NTH_RGB(sm,id)[1],
# Line 505 | Line 441 | int clr;
441    glPushAttrib(GL_DEPTH_BUFFER_BIT);
442    
443    /* First draw background polygons */
508 –
444    glDisable(GL_DEPTH_TEST);
445    glBegin(GL_TRIANGLES);
446    SM_FOR_ALL_ACTIVE_BG_TRIS(sm,i)
447       smRender_bg_tri(sm,i,vp,d,clr);
448    glEnd();
514 –  
449    glEnable(GL_DEPTH_TEST);
450    glBegin(GL_TRIANGLES);
451    SM_FOR_ALL_ACTIVE_FG_TRIS(sm,i)
452    {
453        if(!SM_MIXED_TRI(sm,i))
454          smRender_tri(sm,i,vp,clr);
455 <     else
455 >   else
456          smRender_mixed_tri(sm,i,vp,clr);
457    }
458    glEnd();
# Line 537 | Line 471 | int i;
471  
472    tri = SM_NTH_TRI(sm,i);
473  
474 <  /* Triangles are defined clockwise:historical relative to spherical
541 <     tris: could change
542 <     */
543 <  for(j=2; j >=0; j--)
474 >  for(j=0; j <= 2; j++)
475    {
476      VCOPY(ptr,SM_NTH_WV(sm,T_NTH_V(tri,j)));
477      glVertex3d(ptr[0],ptr[1],ptr[2]);
# Line 584 | Line 515 | T_DEPTH *td;
515          continue;
516      }
517      tri = SM_NTH_TRI(sm,t_id);
518 + #ifdef DEBUG
519 +    if(i >= smNew_tri_cnt)
520 +    {
521 +      eputs("smDepth_sort_tris():More tris than reported by smNew_tri_cnt\n");
522 +      break;
523 +    }
524 + #endif
525      td[i].tri = t_id;
526      min_d = -1;
527      for(j=0;j < 3;j++)
# Line 637 | Line 575 | int clr;
575  
576    /* Turn Depth Test off -- using Painter's algorithm */
577    glPushAttrib(GL_DEPTH_BUFFER_BIT);
578 <  glDisable(GL_DEPTH_TEST);
578 >  glDepthFunc(GL_ALWAYS);
579    d = (dev_zmin+dev_zmax)/2.0;
580    /* Now render back-to front */
581    /* First render bg triangles */
# Line 646 | Line 584 | int clr;
584       smRender_bg_tri(sm,pop_list(&bglist),vp,d,clr);
585    glEnd();
586  
587 <  glEnable(GL_DEPTH_TEST);
587 >
588    glBegin(GL_TRIANGLES);
589    i=0;
590    while(td[i].tri != -1)
# Line 684 | Line 622 | smUpdate(view,qual)
622       epsilon is calculated as running avg of distance of sample points
623       from canonical view: m = 1/(AVG(1/r)): some fraction of this
624     */
625 +
626 +  if(!smMesh)
627 +    return;
628 +
629 +
630 +
631    d = DIST(view->vp,SM_VIEW_CENTER(smMesh));
632    if(qual >= 100 && d > SM_ALLOWED_VIEW_CHANGE(smMesh))
633    {
634        /* Re-build the mesh */
635   #ifdef TEST_DRIVER
636      odev.v = *view;
637 < #endif    
638 <      smRebuild_mesh(smMesh,view->vp);
637 > #endif  
638 >      mark_tris_in_frustum(view);
639 >      smRebuild_mesh(smMesh,view);
640    }
641    /* This is our final update iff qual==100 and view==&odev.v */
642    last_update = qual>=100 && view==&(odev.v);
643    /* Check if we should draw ALL triangles in current frustum */
644 <  if(smClean_notify || smNew_tri_cnt > SM_NUM_TRIS(smMesh)*SM_INC_PERCENT)
644 >  if(smClean_notify || smNew_tri_cnt > SM_SAMPLE_TRIS(smMesh)*SM_INC_PERCENT)
645    {
646   #ifdef TEST_DRIVER
647      glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
648   #else
649 <    if (SM_TONE_MAP(smMesh) < SM_NUM_SAMP(smMesh))
649 >    if ( SM_TONE_MAP(smMesh) < SM_NUM_SAMP(smMesh))
650      {
651         tmClearHisto();
652         tmAddHisto(SM_BRT(smMesh),SM_NUM_SAMP(smMesh),1);
653 <       if(tmComputeMapping(0.,0.,0.) != TM_E_OK ||
654 <                   tmMapPixels(SM_RGB(smMesh),SM_BRT(smMesh),SM_CHR(smMesh),
655 <                                SM_NUM_SAMP(smMesh)) != TM_E_OK)
711 <            return;
653 >       tmComputeMapping(0.,0.,0.);
654 >       tmMapPixels(SM_RGB(smMesh),SM_BRT(smMesh),SM_CHR(smMesh),
655 >                   SM_NUM_SAMP(smMesh));
656      }
657   #endif
658      mark_tris_in_frustum(view);
# Line 723 | Line 667 | smUpdate(view,qual)
667    }
668    /* Do an incremental update instead */
669    else
670 <  {  
671 <    if(!smNew_tri_cnt)
670 >  {
671 >      if(!smNew_tri_cnt)
672        return;
673   #ifdef TEST_DRIVER
674      glDrawBuffer(GL_FRONT);
# Line 737 | Line 681 | smUpdate(view,qual)
681          if(tmComputeMapping(0.,0.,0.) != TM_E_OK)
682             return;
683      }
684 <    if(tmMapPixels(SM_NTH_RGB(smMesh,t),&SM_NTH_BRT(smMesh,t),
685 <                   SM_NTH_CHR(smMesh,t), SM_NUM_SAMP(smMesh)-t) != TM_E_OK)
686 <          return;
684 >    tmMapPixels(SM_NTH_RGB(smMesh,t),&SM_NTH_BRT(smMesh,t),
685 >                   SM_NTH_CHR(smMesh,t), SM_NUM_SAMP(smMesh)-t);
686 >        
687   #endif    
688      smUpdate_Rendered_mesh(smMesh,view->vp,last_update);
689      
# Line 747 | Line 691 | smUpdate(view,qual)
691      glDrawBuffer(GL_BACK);
692   #endif
693    }
694 +
695    SM_TONE_MAP(smMesh) = SM_NUM_SAMP(smMesh);
696 +
697    if (last_update)
698    {
699      smClean_notify = FALSE;
700      smNew_tri_cnt = 0;
701 <    qtHash_init(0);
701 >    smClear_flags(smMesh,T_NEW_FLAG);
702 >    qtCache_init(0);
703    }
704  
705   }

Diff Legend

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