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

Comparing ray/src/hd/rhd_odraw.c (file contents):
Revision 3.3 by gwlarson, Mon Dec 21 13:17:00 1998 UTC vs.
Revision 3.10 by gwlarson, Thu Dec 24 16:18:45 1998 UTC

# Line 14 | Line 14 | static char SCCSid[] = "$SunId$ SGI";
14   #include <GL/glx.h>
15   #include <GL/glu.h>
16  
17 #include "random.h"
17   #include "rhd_odraw.h"
18  
19   #ifndef DEPTHEPS
# Line 131 | Line 130 | int    n;
130                  DCHECK(count<=0 | nextsamp>=n,
131                                  CONSISTENCY, "counter botch in odInit");
132                  if (!i) count = j;
133 +                odView[i].sfirst = nextsamp;
134                  while (j--) {           /* initialize blocks & free lists */
135                          odView[i].bmap[j].pthresh = FHUGE;
136                          odView[i].bmap[j].first = k = nextsamp;
# Line 142 | Line 142 | int    n;
142                          odS.nextfree(k-1) = ENDFREE;
143                          odView[i].bmap[j].nused = 0;
144                  }
145 +                odView[i].snext = nextsamp;
146          }
147          CLR4ALL(odS.redraw, odS.nsamp);         /* clear redraw flags */
148          for (i = odS.nsamp; i--; ) {            /* clear values */
# Line 316 | Line 317 | int    newhist;
317   }
318  
319  
320 + odRedrawAll()                           /* mark all samples for redraw */
321 + {
322 +        register int    i;
323 +
324 +        if ((needmapping&(NEWMAP|NEWRGB)) == (NEWMAP|NEWRGB))
325 +                return;                 /* will be called later, anyway */
326 +        for (i = odS.nsamp; i--; )
327 +                if (odS.ip[i][0] >= 0)
328 +                        SET4(odS.redraw, i);
329 + }
330 +
331 +
332   odRedraw(vn, hmin, vmin, hmax, vmax)    /* redraw view region */
333   int     vn, hmin, vmin, hmax, vmax;
334   {
# Line 424 | Line 437 | GLfloat        *dm;
437   odUpdate(vn)                            /* update this view */
438   int     vn;
439   {
440 <        int     i, j;
428 <        register struct ODblock *bp;
429 <        register int    k;
440 >        register int    i, j;
441  
442          DCHECK(vn<0 | vn>=odNViews, CONSISTENCY,
443                          "bad view number in odUpdate");
# Line 435 | Line 446 | int    vn;
446                  if (needmapping & NEWMAP) {
447                          if (needmapping & NEWHIST)
448                                  tmClearHisto();
449 +                        needmapping &= ~NEWHIST;
450                          if (tmAddHisto(odS.brt,odS.nsamp,1) != TM_E_OK)
451                                  return;
452                          if (tmComputeMapping(0.,0.,0.) != TM_E_OK)
453                                  return;
454 <                        for (k = odS.nsamp; k--; )      /* redraw all */
455 <                                if (odS.ip[k][0] >= 0)
444 <                                        SET4(odS.redraw, k);
454 >                        needmapping &= ~NEWMAP;
455 >                        odRedrawAll();                  /* redraw everything */
456                  }
457                  if (tmMapPixels(odS.rgb,odS.brt,odS.chr,odS.nsamp) != TM_E_OK)
458                          return;
459 <                needmapping = 0;                /* reset flag */
459 >                needmapping &= ~NEWRGB;
460          }
461 <                                        /* draw each block in view */
462 <        for (j = odView[vn].vlow; j--; )
463 <                for (i = 0; i < odView[vn].hlow; i++) {
464 <                                        /* get block */
465 <                        bp = odView[vn].bmap + j*odView[vn].hlow + i;
466 <                                        /* do quick, conservative flag check */
456 <                        for (k = (bp->first+bp->nsamp+31)>>5;
457 <                                        k-- > bp->first>>5; )
458 <                                if (odS.redraw[k])
459 <                                        break;          /* non-zero flag */
460 <                        if (k < bp->first>>5)
461 <                                continue;               /* no flags set */
462 <                        for (k = bp->nsamp; k--; )      /* sample by sample */
463 <                                if (CHK4(odS.redraw, bp->first+k)) {
464 <                                        odDrawBlockSamp(vn, i, j, bp->first+k);
465 <                                        CLR4(odS.redraw, bp->first+k);
466 <                                }
461 >                                        /* this code segment was too slow */
462 > #if 0
463 >        for (i = odView[vn].sfirst; i < odView[vn].snext; i++)
464 >                if (CHK4(odS.redraw, i)) {
465 >                        odDrawSamp(vn, i);
466 >                        CLR4(odS.redraw, i);
467                  }
468 + #else
469 +                                        /* redraw samples at end */
470 +        for (i = odView[vn].snext-31; i < odView[vn].snext; i++)
471 +                if (CHK4(odS.redraw, i)) {
472 +                        odDrawSamp(vn, i);
473 +                        CLR4(odS.redraw, i);
474 +                }
475 +                                        /* faster flag checks in middle */
476 +        for (j = odView[vn].snext>>5; j-- > (odView[vn].sfirst+0x1f)>>5; )
477 +                for (i = 0; odS.redraw[j]; i++)         /* skips faster */
478 +                        if (odS.redraw[j] & 1L<<i) {
479 +                                odDrawSamp(vn, (j<<5)+i);
480 +                                odS.redraw[j] &= ~(1L<<i);
481 +                        }
482 +                                        /* redraw samples at beginning */
483 +        for (i = odView[vn].sfirst; i < odView[vn].sfirst+31; i++)
484 +                if (CHK4(odS.redraw, i)) {
485 +                        odDrawSamp(vn, i);
486 +                        CLR4(odS.redraw, i);
487 +                }
488 + #endif
489   }
490  
491  
492 +                                        /* this turned out to be unnecessary */
493   #if 0
494   static
495   clip_end(p, o, vp)                      /* clip line segment to view */
# Line 505 | Line 527 | register struct ODview *vp;
527   double  sz;
528   {
529          int     na, dv;
530 <        double  hrad, vrad, phi0, phi;
530 >        double  hrad, vrad, phi;
531          register int    i;
532  
533          DCHECK(sz > 1, CONSISTENCY, "super-unary size in make_arms");
534 <        na = MAXFAN*sz*sz + 0.5;                /* keep area constant */
534 >        na = MAXFAN*sz + 0.5;                   /* keep arc length constant */
535          if (na < MINFAN) na = MINFAN;
536          hrad = FANSIZE*sz*vp->hhi/vp->hlow;
537          vrad = FANSIZE*sz*vp->vhi/vp->vlow;
538          if (hrad*vrad < 2.25)
539                  hrad = vrad = 1.5;
518        phi0 = (2.*PI) * frandom();
540          dv = OMAXDEPTH*sz + 0.5;
541          for (i = 0; i < na; i++) {
542 <                phi = phi0 + (2.*PI)*i/na;
542 >                phi = (2.*PI)*i/na;
543                  ar[i][0] = cp[0] + tcos(phi)*hrad + 0.5;
544                  ar[i][1] = cp[1] + tsin(phi)*vrad + 0.5;
545                  ar[i][2] = dv;
# Line 571 | Line 592 | register struct ODview *vp;
592          n = rise + run;
593          while (n--)                     /* run out arm, checking depth */
594                  if (run2 > rise2) {
574                        if (depthchange(vp, x, y, x+xstep, y))
575                                break;
595                          x += xstep;
596                          rise2 += rise;
597 <                } else {
579 <                        if (depthchange(vp, x, y, x, y+ystep))
597 >                        if (depthchange(vp, x-xstep, y, x, y))
598                                  break;
599 +                } else {
600                          y += ystep;
601                          run2 += run;
602 +                        if (depthchange(vp, x, y-ystep, x, y))
603 +                                break;
604                  }
605          if (n < 0)                      /* found something? */
606                  return;
# Line 603 | Line 624 | register int   h, v;
624   }
625  
626  
627 < odDrawBlockSamp(vn, h, v, id)           /* draw sample in view block */
628 < int     vn, h, v;
627 > static int
628 > blockedge(vp, bi0, bi1)                 /* check for edge between blocks? */
629 > register struct ODview  *vp;
630 > register int    bi0, bi1;
631 > {
632 >        if (bi1 == bi0)
633 >                return(0);              /* same block */
634 >        if (bi1 < 0)
635 >                return(1);              /* end off view */
636 >        if (CHK4(vp->emap, bi1))
637 >                return(1);              /* end block has edges */
638 >        if (bi1 == bi0+1 || bi1 == bi0-1 ||
639 >                        bi1 == bi0+vp->hlow || bi1 == bi0-vp->hlow)
640 >                return(0);              /* end in adjacent block -- no edges */
641 >        return(1);                      /* conservative for rarer case */
642 > }
643 >
644 >
645 > odDrawSamp(vn, id)                      /* draw view sample */
646 > int     vn;
647   register int    id;
648   {
649          GLshort arm[MAXFAN][3];
650 <        int     narms, blockindex, bi1;
650 >        int     narms, blockindex;
651          register struct ODview  *vp;
652          double  size;
653          int     home_edges;
654          register int    i;
655  
656          vp = odView + vn;
657 <        blockindex = v*vp->hlow + h;
658 <        DCHECK(odS.ip[id][0]*vp->hlow/vp->hhi != h |
620 <                        odS.ip[id][1]*vp->vlow/vp->vhi != v,
621 <                        CONSISTENCY, "bad sample position in odDrawBlockSamp");
657 >        blockindex = getblock(vp, odS.ip[id][0], odS.ip[id][1]);
658 >        DCHECK(blockindex<0, CONSISTENCY, "bad sample handed to odDrawSamp");
659          DCHECK(vp->bmap[blockindex].nused <= 0,
660 <                        CONSISTENCY, "bad in-use count in odDrawBlockSamp");
660 >                        CONSISTENCY, "bad in-use count in odDrawSamp");
661                                          /* create triangle fan */
662          size = 1./sqrt((double)vp->bmap[blockindex].nused);
663          narms = make_arms(arm, odS.ip[id], vp, size);
664          if (vp->emap != NULL) {         /* check for edge collisions */
665                  home_edges = CHK4(vp->emap, blockindex);
666                  for (i = 0; i < narms; i++)
667 <                        /* the following test is flawed, because we could
668 <                         * be passing through a block on a diagonal run */
632 <                        if (home_edges ||
633 <                                ( (bi1 = getblock(vp, arm[i][0], arm[i][1]))
634 <                                                != blockindex &&
635 <                                        (bi1 < 0 || CHK4(vp->emap, bi1)) ))
667 >                        if (home_edges || blockedge(vp, blockindex,
668 >                                        getblock(vp, arm[i][0], arm[i][1])))
669                                  clip_edge(arm[i], odS.ip[id], vp);
670          }
671                                          /* draw triangle fan */

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines