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.5 by gwlarson, Fri Sep 11 11:52:26 1998 UTC vs.
Revision 3.9 by gwlarson, Wed Nov 11 12:05:39 1998 UTC

# Line 16 | Line 16 | static char SCCSid[] = "$SunId$ SGI";
16   #include <GL/glu.h>
17   #include <glut.h>
18   #endif
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 44 | Line 47 | 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 < int *arg;
53 > int *fptr;
54   {
55    QUADTREE qt = *qtptr;
56    OBJECT *os,*optr;
57    register int i,t_id;
58 +  TRI *tri;
59  
60 <  if (!QT_IS_LEAF(qt))
61 <    return(TRUE);
60 >  if(!QT_FLAG_FILL_TRI(*fptr))
61 >     (*fptr)++;
62 >  if(QT_IS_EMPTY(qt) || QT_LEAF_IS_FLAG(qt))
63 >    return;
64 >
65    /* For each triangle in the set, set the which flag*/
66    os = qtqueryset(qt);
67  
# Line 62 | Line 69 | int *arg;
69    {
70      t_id = QT_SET_NEXT_ELEM(optr);
71      /* Set the render flag */
72 <    if(SM_IS_NTH_T_BASE(smMesh,t_id))
72 >    tri = SM_NTH_TRI(smMesh,t_id);
73 >    if(!T_IS_VALID(tri) || SM_IS_NTH_T_BASE(smMesh,t_id))
74          continue;
75      SM_SET_NTH_T_ACTIVE(smMesh,t_id);
76 <    /* FOR NOW:Also set the LRU clock bit: MAY WANT TO CHANGE: */      
77 <    SM_SET_NTH_T_LRU(smMesh,t_id);
76 >    /* Set the Active bits of the Vertices */
77 >    S_SET_FLAG(T_NTH_V(tri,0));
78 >    S_SET_FLAG(T_NTH_V(tri,1));
79 >    S_SET_FLAG(T_NTH_V(tri,2));
80 >
81    }
71 –  return(TRUE);
82   }
83  
84 + #define mark_active_interior mark_active_tris
85 +
86   mark_tris_in_frustum(view)
87   VIEW *view;
88   {
89      FVECT nr[4],far[4];
90 <
90 >    FPEQ peq;
91 >    int debug=0;
92      /* Mark triangles in approx. view frustum as being active:set
93         LRU counter: for use in discarding samples when out
94         of space
# Line 86 | Line 99 | VIEW *view;
99      /* First clear all the quadtree node and triangle active flags */
100      qtClearAllFlags();
101      smClear_flags(smMesh,T_ACTIVE_FLAG);
102 +    /* Clear all of the active sample flags*/
103 +    sClear_all_flags(SM_SAMP(smMesh));
104  
105 +
106      /* calculate the world space coordinates of the view frustum */
107      calculate_view_frustum(view->vp,view->hvec,view->vvec,view->horiz,
108                             view->vert, dev_zmin,dev_zmax,nr,far);
109  
94 –    /* Project the view frustum onto the spherical quadtree */
95 –    /* For every cell intersected by the projection of the faces
96 –       of the frustum: mark all triangles in the cell as ACTIVE-
97 –       Also set the triangles LRU clock counter
98 –       */
99 –    /* Near face triangles */
100 –    smLocator_apply_func(smMesh,nr[0],nr[2],nr[3],mark_active_tris,NULL);
101 –    smLocator_apply_func(smMesh,nr[2],nr[0],nr[1],mark_active_tris,NULL);
102 –    /* Right face triangles */
103 –    smLocator_apply_func(smMesh,nr[0],far[3],far[0],mark_active_tris,NULL);
104 –    smLocator_apply_func(smMesh,far[3],nr[0],nr[3],mark_active_tris,NULL);
105 –    /* Left face triangles */
106 –    smLocator_apply_func(smMesh,nr[1],far[2],nr[2],mark_active_tris,NULL);
107 –    smLocator_apply_func(smMesh,far[2],nr[1],far[1],mark_active_tris,NULL);
108 –    /* Top face triangles */
109 –    smLocator_apply_func(smMesh,nr[0],far[0],nr[1],mark_active_tris,NULL);
110 –    smLocator_apply_func(smMesh,nr[1],far[0],far[1],mark_active_tris,NULL);
111 –    /* Bottom face triangles */
112 –    smLocator_apply_func(smMesh,nr[3],nr[2],far[3],mark_active_tris,NULL);
113 –    smLocator_apply_func(smMesh,nr[2],far[2],far[3],mark_active_tris,NULL);
114 –    /* Far face triangles */
115 –    smLocator_apply_func(smMesh,far[0],far[2],far[1],mark_active_tris,NULL);
116 –    smLocator_apply_func(smMesh,far[2],far[0],far[3],mark_active_tris,NULL);
117 –
110   #ifdef TEST_DRIVER
111      VCOPY(FrustumFar[0],far[0]);
112      VCOPY(FrustumFar[1],far[1]);
# Line 125 | Line 117 | VIEW *view;
117      VCOPY(FrustumNear[2],nr[2]);
118      VCOPY(FrustumNear[3],nr[3]);
119   #endif
120 +    /* Project the view frustum onto the spherical quadtree */
121 +    /* For every cell intersected by the projection of the faces
122 +
123 +       of the frustum: mark all triangles in the cell as ACTIVE-
124 +       Also set the triangles LRU clock counter
125 +       */
126 +    
127 +    if(EQUAL_VEC3(view->vp,SM_VIEW_CENTER(smMesh)))
128 +    {/* Near face triangles */
129 +
130 +      smLocator_apply_func(smMesh,nr[0],nr[2],nr[3],mark_active_tris,
131 +                         mark_active_interior,NULL);
132 +      smLocator_apply_func(smMesh,nr[2],nr[0],nr[1],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 233 | 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);
# Line 274 | Line 378 | int lvl;
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 >        if(SM_IS_NTH_T_BASE(sm,os[n]))
382 >           continue;
383 >        tri = SM_NTH_TRI(sm,os[n]);
384 >        if(!T_IS_VALID(tri))
385 >          continue;
386 >        
387 >        s0 = T_NTH_V(tri,0);
388 >        s1 = T_NTH_V(tri,1);
389 >        s2 = T_NTH_V(tri,2);
390 >        VCOPY(a,SM_NTH_WV(sm,s0));
391 >        VCOPY(b,SM_NTH_WV(sm,s1));
392 >        VCOPY(c,SM_NTH_WV(sm,s2));            
393          distsum += SM_BG_SAMPLE(sm,s0) ? dev_zmax
394                                  : sqrt(dist2(a,SM_VIEW_CENTER(sm)));
395          distsum += SM_BG_SAMPLE(sm,s1) ? dev_zmax
# Line 306 | Line 421 | int lvl;
421      VSUM(c,SM_VIEW_CENTER(sm),c,le->av.dist);
422                                          /* draw triangle */
423      glColor3ub(le->av.rgb[0],le->av.rgb[1],le->av.rgb[2]);
309 –    /* NOTE: Triangle vertex order may change */
310 –    glVertex3d(c[0],c[1],c[2]);
311 –    glVertex3d(b[0],b[1],b[2]);
424      glVertex3d(a[0],a[1],a[2]);
425 +    glVertex3d(b[0],b[1],b[2]);
426 +    glVertex3d(c[0],c[1],c[2]);
427 +
428    }
429    return(&le->av);
430   }
# Line 322 | Line 437 | int lvl;
437    QUADTREE root;
438    int i;
439    FVECT t0,t1,t2;
440 +  STREE *st;
441  
442 +  
443    if (lvl < 1)
444      return;
445 +  st = SM_LOCATOR(sm);
446    glPushAttrib(GL_LIGHTING_BIT);
447    glShadeModel(GL_FLAT);
448    glBegin(GL_TRIANGLES);
449 <  for(i=0; i < 4; i++)
449 >  for(i=0; i < ST_NUM_ROOT_NODES; i++)
450    {
451 <    root = ST_NTH_ROOT(SM_LOCATOR(sm),i);
452 <    stNth_base_verts(SM_LOCATOR(sm),i,t0,t1,t2);
453 <    qtRender_level(root,t0,t1,t2,sm,lvl-1);
451 >    root = ST_NTH_ROOT(st,i);
452 >    qtRender_level(root,ST_NTH_V(st,i,0),ST_NTH_V(st,i,1),ST_NTH_V(st,i,2),
453 >                   sm,lvl-1);
454    }
455    glEnd();
456    glPopAttrib();
# Line 357 | Line 475 | int qual;
475      return;
476    for (i = QT_MAX_LEVELS; i--; )
477      lvlcnt[i] = 0;
478 <  stCount_level_leaves(lvlcnt, SM_LOCATOR(sm)->root);
478 >  stCount_level_leaves(lvlcnt, ST_TOP_ROOT(SM_LOCATOR(sm)));
479 >  stCount_level_leaves(lvlcnt, ST_BOTTOM_ROOT(SM_LOCATOR(sm)));  
480    for (ntarget -= lvlcnt[i=0]; i < QT_MAX_LEVELS-1; ntarget -= lvlcnt[++i])
481      if (ntarget < lvlcnt[i+1])
482        break;
# Line 377 | Line 496 | int clr;
496    int j;
497  
498    tri = SM_NTH_TRI(sm,i);
499 <  if (clr) SM_CLEAR_NTH_T_NEW(sm,i);
500 <
382 <  /* NOTE:Triangles are defined clockwise:historical relative to spherical
383 <     tris: could change
384 <     */
385 <  for(j=2; j>= 0; j--)
499 >  if (clr) SM_CLR_NTH_T_NEW(sm,i);
500 >  for(j=0; j <= 2; j++)
501    {
502   #ifdef DEBUG
503      if(SM_BG_SAMPLE(sm,T_NTH_V(tri,j)))
# Line 408 | Line 523 | int clr;
523    int rgb[3];
524  
525    tri = SM_NTH_TRI(sm,i);
526 <  if (clr) SM_CLEAR_NTH_T_NEW(sm,i);
526 >  if (clr) SM_CLR_NTH_T_NEW(sm,i);
527  
528    /* NOTE:Triangles are defined clockwise:historical relative to spherical
529       tris: could change
# Line 433 | Line 548 | int clr;
548      rgb[0]/=cnt; rgb[1]/=cnt; rgb[2]/=cnt;
549      d /= (double)cnt;
550    }
551 <  for(j=2; j>= 0; j--)
551 >  for(j=0; j <= 2; j++)
552    {
553      if(SM_BG_SAMPLE(sm,ids[j]))
554      {
# Line 466 | Line 581 | int clr;
581    TRI *tri;
582    
583    tri = SM_NTH_TRI(sm,i);
584 <  if (clr) SM_CLEAR_NTH_T_NEW(sm,i);
584 >  if (clr) SM_CLR_NTH_T_NEW(sm,i);
585  
586    /* NOTE:Triangles are defined clockwise:historical relative to spherical
587       tris: could change
588       */
589 <  for(j=2; j>= 0; j--)
589 >  for(j=0; j <= 2; j++)
590    {
591        id = T_NTH_V(tri,j);
592        glColor3ub(SM_NTH_RGB(sm,id)[0],SM_NTH_RGB(sm,id)[1],
# Line 502 | Line 617 | int clr;
617    glPushAttrib(GL_DEPTH_BUFFER_BIT);
618    
619    /* First draw background polygons */
505 –
620    glDisable(GL_DEPTH_TEST);
621    glBegin(GL_TRIANGLES);
622    SM_FOR_ALL_ACTIVE_BG_TRIS(sm,i)
623       smRender_bg_tri(sm,i,vp,d,clr);
624    glEnd();
511 –  
625    glEnable(GL_DEPTH_TEST);
626    glBegin(GL_TRIANGLES);
627    SM_FOR_ALL_ACTIVE_FG_TRIS(sm,i)
628    {
629        if(!SM_MIXED_TRI(sm,i))
630          smRender_tri(sm,i,vp,clr);
631 <     else
631 >   else
632          smRender_mixed_tri(sm,i,vp,clr);
633    }
634    glEnd();
# Line 534 | Line 647 | int i;
647  
648    tri = SM_NTH_TRI(sm,i);
649  
650 <  /* Triangles are defined clockwise:historical relative to spherical
538 <     tris: could change
539 <     */
540 <  for(j=2; j >=0; j--)
650 >  for(j=0; j <= 2; j++)
651    {
652      VCOPY(ptr,SM_NTH_WV(sm,T_NTH_V(tri,j)));
653      glVertex3d(ptr[0],ptr[1],ptr[2]);
# Line 581 | Line 691 | T_DEPTH *td;
691          continue;
692      }
693      tri = SM_NTH_TRI(sm,t_id);
694 + #ifdef DEBUG
695 +    if(i >= smNew_tri_cnt)
696 +    {
697 +      eputs("smDepth_sort_tris():More tris than reported by smNew_tri_cnt\n");
698 +      break;
699 +    }
700 + #endif
701      td[i].tri = t_id;
702      min_d = -1;
703      for(j=0;j < 3;j++)
# Line 634 | Line 751 | int clr;
751  
752    /* Turn Depth Test off -- using Painter's algorithm */
753    glPushAttrib(GL_DEPTH_BUFFER_BIT);
754 <  glDisable(GL_DEPTH_TEST);
754 >  glDepthFunc(GL_ALWAYS);
755    d = (dev_zmin+dev_zmax)/2.0;
756    /* Now render back-to front */
757    /* First render bg triangles */
# Line 643 | Line 760 | int clr;
760       smRender_bg_tri(sm,pop_list(&bglist),vp,d,clr);
761    glEnd();
762  
763 <  glEnable(GL_DEPTH_TEST);
763 >
764    glBegin(GL_TRIANGLES);
765    i=0;
766    while(td[i].tri != -1)
# Line 681 | Line 798 | smUpdate(view,qual)
798       epsilon is calculated as running avg of distance of sample points
799       from canonical view: m = 1/(AVG(1/r)): some fraction of this
800     */
801 +
802 +  if(!smMesh)
803 +    return;
804    d = DIST(view->vp,SM_VIEW_CENTER(smMesh));
805    if(qual >= 100 && d > SM_ALLOWED_VIEW_CHANGE(smMesh))
806    {
807        /* Re-build the mesh */
808   #ifdef TEST_DRIVER
809      odev.v = *view;
810 < #endif    
811 <      smRebuild_mesh(smMesh,view->vp);
810 > #endif  
811 >      mark_tris_in_frustum(view);
812 >      smRebuild_mesh(smMesh,view);
813    }
814    /* This is our final update iff qual==100 and view==&odev.v */
815    last_update = qual>=100 && view==&(odev.v);
816    /* Check if we should draw ALL triangles in current frustum */
817 <  if(smClean_notify || smNew_tri_cnt > SM_NUM_TRIS(smMesh)*SM_INC_PERCENT)
817 >  if(smClean_notify || smNew_tri_cnt > SM_SAMPLE_TRIS(smMesh)*SM_INC_PERCENT)
818    {
819   #ifdef TEST_DRIVER
820      glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
821   #else
822 <    if (SM_TONE_MAP(smMesh) < SM_NUM_SAMP(smMesh))
822 >    if ( SM_TONE_MAP(smMesh) < SM_NUM_SAMP(smMesh))
823      {
824         tmClearHisto();
825         tmAddHisto(SM_BRT(smMesh),SM_NUM_SAMP(smMesh),1);
# Line 720 | Line 841 | smUpdate(view,qual)
841    }
842    /* Do an incremental update instead */
843    else
844 <  {  
845 <    if(!smNew_tri_cnt)
844 >  {
845 >      if(!smNew_tri_cnt)
846        return;
847   #ifdef TEST_DRIVER
848      glDrawBuffer(GL_FRONT);
# Line 744 | Line 865 | smUpdate(view,qual)
865      glDrawBuffer(GL_BACK);
866   #endif
867    }
868 +
869    SM_TONE_MAP(smMesh) = SM_NUM_SAMP(smMesh);
870 +
871    if (last_update)
872    {
873      smClean_notify = FALSE;
874      smNew_tri_cnt = 0;
875 +    smClear_flags(smMesh,T_NEW_FLAG);
876      qtCache_init(0);
877    }
878  

Diff Legend

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