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.8 by gwlarson, Tue Oct 6 18:16:54 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)
50 >
51 > mark_active_tris(qtptr,fptr)
52   QUADTREE *qtptr;
53 < char *arg;
53 > int *fptr;
54   {
55    QUADTREE qt = *qtptr;
56 <  OBJECT os[QT_MAX_SET+1],*optr;
56 >  OBJECT *os,*optr;
57    register int i,t_id;
58 +  TRI *tri;
59  
60 <  if (!QT_IS_LEAF(qt))
61 <    return(TRUE);
60 >
61 >  if(!QT_FLAG_FILL_TRI(*fptr))
62 >     (*fptr)++;
63 >
64 >  if(QT_IS_EMPTY(qt) || QT_LEAF_IS_FLAG(qt))
65 >    return;
66    /* For each triangle in the set, set the which flag*/
67 <  qtgetset(os,qt);
67 >  os = qtqueryset(qt);
68  
69    for (i = QT_SET_CNT(os), optr = QT_SET_PTR(os); i > 0; i--)
70    {
71      t_id = QT_SET_NEXT_ELEM(optr);
72      /* Set the render flag */
73 <    if(SM_IS_NTH_T_BASE(smMesh,t_id))
73 >    tri = SM_NTH_TRI(smMesh,t_id);
74 >    if(!T_IS_VALID(tri) || SM_IS_NTH_T_BASE(smMesh,t_id))
75          continue;
76      SM_SET_NTH_T_ACTIVE(smMesh,t_id);
77 <    /* FOR NOW:Also set the LRU clock bit: MAY WANT TO CHANGE: */      
78 <    SM_SET_NTH_T_LRU(smMesh,t_id);
77 >    /* Set the Active bits of the Vertices */
78 >    S_SET_FLAG(T_NTH_V(tri,0));
79 >    S_SET_FLAG(T_NTH_V(tri,1));
80 >    S_SET_FLAG(T_NTH_V(tri,2));
81 >
82    }
72 –  return(TRUE);
83   }
84  
85 + #define mark_active_interior mark_active_tris
86 +
87   mark_tris_in_frustum(view)
88   VIEW *view;
89   {
90      FVECT nr[4],far[4];
91 <
91 >    FPEQ peq;
92 >    int debug=0;
93      /* Mark triangles in approx. view frustum as being active:set
94         LRU counter: for use in discarding samples when out
95         of space
# Line 87 | Line 100 | VIEW *view;
100      /* First clear all the quadtree node and triangle active flags */
101      qtClearAllFlags();
102      smClear_flags(smMesh,T_ACTIVE_FLAG);
103 +    /* Clear all of the active sample flags*/
104 +    sClear_all_flags(smMesh->samples);
105  
106 +
107      /* calculate the world space coordinates of the view frustum */
108      calculate_view_frustum(view->vp,view->hvec,view->vvec,view->horiz,
109                             view->vert, dev_zmin,dev_zmax,nr,far);
110  
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 –
111   #ifdef TEST_DRIVER
112      VCOPY(FrustumFar[0],far[0]);
113      VCOPY(FrustumFar[1],far[1]);
# Line 126 | Line 118 | VIEW *view;
118      VCOPY(FrustumNear[2],nr[2]);
119      VCOPY(FrustumNear[3],nr[3]);
120   #endif
121 +    /* Project the view frustum onto the spherical quadtree */
122 +    /* For every cell intersected by the projection of the faces
123 +
124 +       of the frustum: mark all triangles in the cell as ACTIVE-
125 +       Also set the triangles LRU clock counter
126 +       */
127 +    
128 +    if(EQUAL_VEC3(view->vp,SM_VIEW_CENTER(smMesh)))
129 +    {/* Near face triangles */
130 +      smLocator_apply_func(smMesh,nr[3],nr[2],nr[0],mark_active_tris,
131 +                         mark_active_interior,NULL);
132 +      smLocator_apply_func(smMesh,nr[1],nr[0],nr[2],mark_active_tris,
133 +                         mark_active_interior,NULL);
134 +      return;
135 +    }
136 +
137 +    /* Test the view against the planes: and swap orientation if inside:*/
138 +    tri_plane_equation(nr[0],nr[2],nr[3], &peq,FALSE);
139 +    if(PT_ON_PLANE(SM_VIEW_CENTER(smMesh),peq) > 0.0)
140 +    {/* Near face triangles */
141 +      smLocator_apply_func(smMesh,nr[3],nr[2],nr[0],mark_active_tris,
142 +                           mark_active_interior,NULL);
143 +      smLocator_apply_func(smMesh,nr[1],nr[0],nr[2],mark_active_tris,
144 +                         mark_active_interior,NULL);
145 +    }
146 +    else
147 +    {/* Near face triangles */
148 +      smLocator_apply_func(smMesh,nr[0],nr[2],nr[3],mark_active_tris,
149 +                           mark_active_interior,NULL);
150 +      smLocator_apply_func(smMesh,nr[2],nr[0],nr[1],mark_active_tris,
151 +                         mark_active_interior,NULL);
152 +    }
153 +    tri_plane_equation(nr[0],far[3],far[0], &peq,FALSE);
154 +    if(PT_ON_PLANE(SM_VIEW_CENTER(smMesh),peq) > 0.0)
155 +    { /* Right face triangles */
156 +      smLocator_apply_func(smMesh,far[0],far[3],nr[0],mark_active_tris,
157 +                           mark_active_interior,NULL);
158 +      smLocator_apply_func(smMesh,nr[3],nr[0],far[3],mark_active_tris,
159 +                           mark_active_interior,NULL);
160 +    }
161 +    else
162 +    {/* Right face triangles */
163 +      smLocator_apply_func(smMesh,nr[0],far[3],far[0],mark_active_tris,
164 +                           mark_active_interior,NULL);
165 +      smLocator_apply_func(smMesh,far[3],nr[0],nr[3],mark_active_tris,
166 +                           mark_active_interior,NULL);
167 +    }
168 +
169 +    tri_plane_equation(nr[1],far[2],nr[2], &peq,FALSE);
170 +    if(PT_ON_PLANE(SM_VIEW_CENTER(smMesh),peq) > 0.0)
171 +    { /* Left face triangles */
172 +      smLocator_apply_func(smMesh,nr[2],far[2],nr[1],mark_active_tris,
173 +                           mark_active_interior,NULL);
174 +      smLocator_apply_func(smMesh,far[1],nr[1],far[2],mark_active_tris,
175 +                           mark_active_interior,NULL);  
176 +    }
177 +    else
178 +    { /* Left face triangles */
179 +      smLocator_apply_func(smMesh,nr[1],far[2],nr[2],mark_active_tris,
180 +                           mark_active_interior,NULL);
181 +      smLocator_apply_func(smMesh,far[2],nr[1],far[1],mark_active_tris,
182 +                   mark_active_interior,NULL);
183 +
184 +    }
185 +    tri_plane_equation(nr[0],far[0],nr[1], &peq,FALSE);
186 +    if(PT_ON_PLANE(SM_VIEW_CENTER(smMesh),peq) > 0.0)
187 +    {/* Top face triangles */
188 +      smLocator_apply_func(smMesh,nr[1],far[0],nr[0],mark_active_tris,
189 +                           mark_active_interior,NULL);
190 +      smLocator_apply_func(smMesh,far[1],far[0],nr[1],mark_active_tris,
191 +                           mark_active_interior,NULL);
192 +    }
193 +    else
194 +    {/* Top face triangles */
195 +      smLocator_apply_func(smMesh,nr[0],far[0],nr[1],mark_active_tris,
196 +                           mark_active_interior,NULL);
197 +      smLocator_apply_func(smMesh,nr[1],far[0],far[1],mark_active_tris,
198 +                           mark_active_interior,NULL);
199 +    }
200 +    tri_plane_equation(nr[3],nr[2],far[3], &peq,FALSE);
201 +    if(PT_ON_PLANE(SM_VIEW_CENTER(smMesh),peq) > 0.0)
202 +    {/* Bottom face triangles */
203 +      smLocator_apply_func(smMesh,far[3],nr[2],nr[3],mark_active_tris,
204 +                           mark_active_interior,NULL);
205 +      smLocator_apply_func(smMesh,far[3],far[2],nr[2],mark_active_tris,
206 +                           mark_active_interior,NULL);
207 +    }
208 +    else
209 +    { /* Bottom face triangles */
210 +      smLocator_apply_func(smMesh,nr[3],nr[2],far[3],mark_active_tris,
211 +                           mark_active_interior,NULL);
212 +      smLocator_apply_func(smMesh,nr[2],far[2],far[3],mark_active_tris,
213 +                           mark_active_interior,NULL);
214 +    }
215 +     tri_plane_equation(far[2],far[0],far[1], &peq,FALSE);
216 +    if(PT_ON_PLANE(SM_VIEW_CENTER(smMesh),peq) > 0.0)
217 +    {/* Far face triangles */
218 +      smLocator_apply_func(smMesh,far[0],far[2],far[1],mark_active_tris,
219 +                           mark_active_interior,NULL);
220 +      smLocator_apply_func(smMesh,far[2],far[0],far[3],mark_active_tris,
221 +                           mark_active_interior,NULL);
222 +    }
223 +    else
224 +    {/* Far face triangles */
225 +      smLocator_apply_func(smMesh,far[1],far[2],far[0],mark_active_tris,
226 +                           mark_active_interior,NULL);
227 +      smLocator_apply_func(smMesh,far[3],far[0],far[2],mark_active_tris,
228 +                           mark_active_interior,NULL);
229 +    }
230 +
231   }
232  
233   /*
# Line 140 | Line 242 | smClean()
242   }
243  
244   int
245 < qtHash_init(nel)                /* initialize for at least nel elements */
245 > qtCache_init(nel)               /* initialize for at least nel elements */
246   int     nel;
247   {
248          static int  hsiztab[] = {
# Line 171 | Line 273 | int    nel;
273   }
274  
275   QT_LUENT *
276 < qtHash_find(qt)                 /* find a quadtree table entry */
276 > qtCache_find(qt)                /* find a quadtree table entry */
277   QUADTREE qt;
278   {
279          int     i, n;
280          register int    ndx;
281          register QT_LUENT       *le;
282  
283 <        if (qt_hsiz == 0)
284 <                qtHash_init(1);
283 >        if (qt_hsiz == 0 && !qtCache_init(1))
284 >                return(NULL);
285   tryagain:                               /* hash table lookup */
286          ndx = (unsigned long)qt % qt_hsiz;
287          for (i = 0, n = 1; i < qt_hsiz; i++, n += 2) {
# Line 192 | Line 294 | tryagain:                              /* hash table lookup */
294                                          /* table is full, reallocate */
295          le = qt_htbl;
296          ndx = qt_hsiz;
297 <        if (!qtHash_init(ndx+1)) {      /* no more memory! */
297 >        if (!qtCache_init(ndx+1)) {     /* no more memory! */
298                  qt_htbl = le;
299                  qt_hsiz = ndx;
300                  return(NULL);
301          }
302 <        if (!ndx)
201 <                goto tryagain;
202 <                                        /* copy old table to new */
302 >                                        /* copy old table to new and free */
303          while (ndx--)
304                  if (!QT_IS_EMPTY(le[ndx].qt))
305 <                        copystruct(qtHash_find(le[ndx].qt), &le[ndx]);
305 >                        copystruct(qtCache_find(le[ndx].qt), &le[ndx]);
306          free((char *)le);
307          goto tryagain;                  /* should happen only once! */
308   }
# Line 236 | Line 336 | int lvl;
336    FVECT a,b,c;
337    register QT_LUENT *le;
338    QTRAVG *rc[4];
339 +  TRI *tri;
340    
341    if (QT_IS_EMPTY(qt))                          /* empty leaf node */
342      return(NULL);
343    if (QT_IS_TREE(qt) && !QT_IS_FLAG(qt))        /* not in our frustum */
344      return(NULL);
345                                          /* else look up node */
346 <  if ((le = qtHash_find(qt)) == NULL)
346 >  if ((le = qtCache_find(qt)) == NULL)
347      error(SYSTEM, "out of memory in qtRender_level");
348    if (QT_IS_TREE(qt) && (QT_IS_EMPTY(le->qt) || lvl > 0))
349    {                                     /* compute children */
# Line 271 | Line 372 | int lvl;
372      }
373      else
374      {                                   /* from triangle set */
375 <      OBJECT os[QT_MAX_SET+1];
375 >      OBJECT *os;
376        int s0, s1, s2;
377  
378 <      qtgetset(os,qt);
378 >      os = qtqueryset(qt);
379        for (n = os[0]; n; n--)
380        {
381 <        qtTri_from_id(os[n],a,b,c,NULL,NULL,NULL,&s0,&s1,&s2);
381 >        tri = SM_NTH_TRI(sm,os[n]);
382 >        if(!T_IS_VALID(tri))
383 >          continue;
384 >
385 >        s0 = T_NTH_V(tri,0);
386 >        s1 = T_NTH_V(tri,1);
387 >        s2 = T_NTH_V(tri,2);
388 >        VCOPY(a,SM_NTH_WV(sm,s0));
389 >        VCOPY(b,SM_NTH_WV(sm,s1));
390 >        VCOPY(c,SM_NTH_WV(sm,s2));            
391          distsum += SM_BG_SAMPLE(sm,s0) ? dev_zmax
392                                  : sqrt(dist2(a,SM_VIEW_CENTER(sm)));
393          distsum += SM_BG_SAMPLE(sm,s1) ? dev_zmax
# Line 325 | Line 435 | int lvl;
435    QUADTREE root;
436    int i;
437    FVECT t0,t1,t2;
438 +  STREE *st;
439  
440 +  
441    if (lvl < 1)
442      return;
443 +  st = SM_LOCATOR(sm);
444    glPushAttrib(GL_LIGHTING_BIT);
445    glShadeModel(GL_FLAT);
446    glBegin(GL_TRIANGLES);
447 <  for(i=0; i < 4; i++)
447 >  for(i=0; i < ST_NUM_ROOT_NODES; i++)
448    {
449 <    root = ST_NTH_ROOT(SM_LOCATOR(sm),i);
450 <    stNth_base_verts(SM_LOCATOR(sm),i,t0,t1,t2);
451 <    qtRender_level(root,t0,t1,t2,sm,lvl-1);
449 >    root = ST_NTH_ROOT(st,i);
450 >    qtRender_level(root,ST_NTH_V(st,i,0),ST_NTH_V(st,i,1),ST_NTH_V(st,i,2),
451 >                   sm,lvl-1);
452    }
453    glEnd();
454    glPopAttrib();
# Line 360 | Line 473 | int qual;
473      return;
474    for (i = QT_MAX_LEVELS; i--; )
475      lvlcnt[i] = 0;
476 <  stCount_level_leaves(lvlcnt, SM_LOCATOR(sm)->root);
476 >  stCount_level_leaves(lvlcnt, ST_TOP_ROOT(SM_LOCATOR(sm)));
477 >  stCount_level_leaves(lvlcnt, ST_BOTTOM_ROOT(SM_LOCATOR(sm)));  
478    for (ntarget -= lvlcnt[i=0]; i < QT_MAX_LEVELS-1; ntarget -= lvlcnt[++i])
479      if (ntarget < lvlcnt[i+1])
480        break;
# Line 380 | Line 494 | int clr;
494    int j;
495  
496    tri = SM_NTH_TRI(sm,i);
497 <  if (clr) SM_CLEAR_NTH_T_NEW(sm,i);
497 >  if (clr) SM_CLR_NTH_T_NEW(sm,i);
498  
499    /* NOTE:Triangles are defined clockwise:historical relative to spherical
500       tris: could change
# Line 411 | Line 525 | int clr;
525    int rgb[3];
526  
527    tri = SM_NTH_TRI(sm,i);
528 <  if (clr) SM_CLEAR_NTH_T_NEW(sm,i);
528 >  if (clr) SM_CLR_NTH_T_NEW(sm,i);
529  
530    /* NOTE:Triangles are defined clockwise:historical relative to spherical
531       tris: could change
# Line 469 | Line 583 | int clr;
583    TRI *tri;
584    
585    tri = SM_NTH_TRI(sm,i);
586 <  if (clr) SM_CLEAR_NTH_T_NEW(sm,i);
586 >  if (clr) SM_CLR_NTH_T_NEW(sm,i);
587  
588    /* NOTE:Triangles are defined clockwise:historical relative to spherical
589       tris: could change
# Line 505 | Line 619 | int clr;
619    glPushAttrib(GL_DEPTH_BUFFER_BIT);
620    
621    /* First draw background polygons */
508 –
622    glDisable(GL_DEPTH_TEST);
623    glBegin(GL_TRIANGLES);
624    SM_FOR_ALL_ACTIVE_BG_TRIS(sm,i)
625       smRender_bg_tri(sm,i,vp,d,clr);
626    glEnd();
514 –  
627    glEnable(GL_DEPTH_TEST);
628    glBegin(GL_TRIANGLES);
629    SM_FOR_ALL_ACTIVE_FG_TRIS(sm,i)
630    {
631 +      if(SM_BG_TRI(sm,i))
632 +         continue;
633        if(!SM_MIXED_TRI(sm,i))
634          smRender_tri(sm,i,vp,clr);
635 <     else
635 >   else
636          smRender_mixed_tri(sm,i,vp,clr);
637    }
638    glEnd();
# Line 584 | Line 698 | T_DEPTH *td;
698          continue;
699      }
700      tri = SM_NTH_TRI(sm,t_id);
701 + #ifdef DEBUG
702 +    if(i >= smNew_tri_cnt)
703 +    {
704 +      eputs("smDepth_sort_tris():More tris than reported by smNew_tri_cnt\n");
705 +      break;
706 +    }
707 + #endif
708      td[i].tri = t_id;
709      min_d = -1;
710      for(j=0;j < 3;j++)
# Line 646 | Line 767 | int clr;
767       smRender_bg_tri(sm,pop_list(&bglist),vp,d,clr);
768    glEnd();
769  
770 <  glEnable(GL_DEPTH_TEST);
770 >
771    glBegin(GL_TRIANGLES);
772    i=0;
773    while(td[i].tri != -1)
# Line 684 | Line 805 | smUpdate(view,qual)
805       epsilon is calculated as running avg of distance of sample points
806       from canonical view: m = 1/(AVG(1/r)): some fraction of this
807     */
808 +
809 +  if(!smMesh)
810 +    return;
811    d = DIST(view->vp,SM_VIEW_CENTER(smMesh));
812    if(qual >= 100 && d > SM_ALLOWED_VIEW_CHANGE(smMesh))
813    {
814        /* Re-build the mesh */
815   #ifdef TEST_DRIVER
816      odev.v = *view;
817 < #endif    
818 <      smRebuild_mesh(smMesh,view->vp);
817 > #endif  
818 >      mark_tris_in_frustum(view);
819 >      smRebuild_mesh(smMesh,view);
820    }
821    /* This is our final update iff qual==100 and view==&odev.v */
822    last_update = qual>=100 && view==&(odev.v);
823    /* Check if we should draw ALL triangles in current frustum */
824 <  if(smClean_notify || smNew_tri_cnt > SM_NUM_TRIS(smMesh)*SM_INC_PERCENT)
824 >  if(smClean_notify || smNew_tri_cnt > SM_SAMPLE_TRIS(smMesh)*SM_INC_PERCENT)
825    {
826   #ifdef TEST_DRIVER
827      glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
828   #else
829 <    if (SM_TONE_MAP(smMesh) < SM_NUM_SAMP(smMesh))
829 >    if ( SM_TONE_MAP(smMesh) < SM_NUM_SAMP(smMesh))
830      {
831         tmClearHisto();
832         tmAddHisto(SM_BRT(smMesh),SM_NUM_SAMP(smMesh),1);
# Line 723 | Line 848 | smUpdate(view,qual)
848    }
849    /* Do an incremental update instead */
850    else
851 <  {  
852 <    if(!smNew_tri_cnt)
851 >  {
852 >      if(!smNew_tri_cnt)
853        return;
854   #ifdef TEST_DRIVER
855      glDrawBuffer(GL_FRONT);
# Line 752 | Line 877 | smUpdate(view,qual)
877    {
878      smClean_notify = FALSE;
879      smNew_tri_cnt = 0;
880 <    qtHash_init(0);
880 >    smClear_flags(smMesh,T_NEW_FLAG);
881 >    qtCache_init(0);
882    }
883  
884   }

Diff Legend

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