render.c File Reference

#include "gluos.h"
#include <assert.h>
#include "tess.h"

Include dependency graph for render.c:

Go to the source code of this file.

Classes
struct	FaceCount

Macros
#define	TRUE   1
#define	FALSE   0
#define	Marked(f)   (! (f)->inside || (f)->marked)
#define	AddToTrail(f, t)   ((f)->trail = (t), (t) = (f), (f)->marked = TRUE)
#define	FreeTrail(t)
#define	IsEven(n)   (((n) & 1) == 0)
#define	SIGN_INCONSISTENT   2

Functions
static struct FaceCount	MaximumFan (GLUhalfEdge *eOrig)
static struct FaceCount	MaximumStrip (GLUhalfEdge *eOrig)
static void	RenderFan (GLUtesselator *tess, GLUhalfEdge *eStart, long size)
static void	RenderStrip (GLUtesselator *tess, GLUhalfEdge *eStart, long size)
static void	RenderTriangle (GLUtesselator *tess, GLUhalfEdge *eStart, long size)
static void	RenderMaximumFaceGroup (GLUtesselator *tess, GLUface *fOrig)
static void	RenderLonelyTriangles (GLUtesselator *tess, GLUface *head)
void	__gl_renderMesh (GLUtesselator *tess, GLUmesh *mesh)
void	__gl_renderBoundary (GLUtesselator *tess, GLUmesh *mesh)
static int	ComputeNormal (GLUtesselator *tess, GLdouble norm[3], int check)
GLboolean	__gl_renderCache (GLUtesselator *tess)

Macro Definition Documentation

AddToTrail

#define AddToTrail	(		f,
			t
	)		((f)->trail = (t), (t) = (f), (f)->marked = TRUE)

Definition at line 148 of file render.c.

FALSE

#define FALSE   0

Definition at line 46 of file render.c.

FreeTrail

#define FreeTrail

(

t

)

Value:

              do { \
              while( (t) != NULL ) { \
                (t)->marked = FALSE; t = (t)->trail; \
              } \
            } while(0) /* absorb trailing semicolon */

Definition at line 150 of file render.c.

IsEven

#define IsEven

(

n

)

(((n) & 1) == 0)

Definition at line 183 of file render.c.

Marked

#define Marked

(

f

)

(! (f)->inside || (f)->marked)

Definition at line 146 of file render.c.

SIGN_INCONSISTENT

#define SIGN_INCONSISTENT   2

Definition at line 360 of file render.c.

TRUE

#define TRUE   1

Definition at line 43 of file render.c.

Function Documentation

__gl_renderBoundary()

void __gl_renderBoundary	(	GLUtesselator *	tess,
		GLUmesh *	mesh
	)

Definition at line 339 of file render.c.

{
  GLUface *f;
  GLUhalfEdge *e;
 
  for( f = mesh->fHead.next; f != &mesh->fHead; f = f->next ) {
    if( f->inside ) {
      CALL_BEGIN_OR_BEGIN_DATA( GL_LINE_LOOP );
      e = f->anEdge;
      do {
        CALL_VERTEX_OR_VERTEX_DATA( e->Org->data ); 
    e = e->Lnext;
      } while( e != f->anEdge );
      CALL_END_OR_END_DATA();
    }
  }
}

Referenced by gluTessEndPolygon().

__gl_renderCache()

GLboolean __gl_renderCache

(

GLUtesselator *

tess

)

Definition at line 441 of file render.c.

{
  CachedVertex *v0 = tess->cache;
  CachedVertex *vn = v0 + tess->cacheCount;
  CachedVertex *vc;
  GLdouble norm[3];
  int sign;
 
  if( tess->cacheCount < 3 ) {
    /* Degenerate contour -- no output */
    return TRUE;
  }
 
  norm[0] = tess->normal[0];
  norm[1] = tess->normal[1];
  norm[2] = tess->normal[2];
  if( norm[0] == 0 && norm[1] == 0 && norm[2] == 0 ) {
    ComputeNormal( tess, norm, FALSE );
  }
 
  sign = ComputeNormal( tess, norm, TRUE );
  if( sign == SIGN_INCONSISTENT ) {
    /* Fan triangles did not have a consistent orientation */
    return FALSE;
  }
  if( sign == 0 ) {
    /* All triangles were degenerate */
    return TRUE;
  }
 
  /* Make sure we do the right thing for each winding rule */
  switch( tess->windingRule ) {
  case GLU_TESS_WINDING_ODD:
  case GLU_TESS_WINDING_NONZERO:
    break;
  case GLU_TESS_WINDING_POSITIVE:
    if( sign < 0 ) return TRUE;
    break;
  case GLU_TESS_WINDING_NEGATIVE:
    if( sign > 0 ) return TRUE;
    break;
  case GLU_TESS_WINDING_ABS_GEQ_TWO:
    return TRUE;
  }
 
  CALL_BEGIN_OR_BEGIN_DATA( tess->boundaryOnly ? GL_LINE_LOOP
              : (tess->cacheCount > 3) ? GL_TRIANGLE_FAN
              : GL_TRIANGLES );
 
  CALL_VERTEX_OR_VERTEX_DATA( v0->data ); 
  if( sign > 0 ) {
    for( vc = v0+1; vc < vn; ++vc ) {
      CALL_VERTEX_OR_VERTEX_DATA( vc->data ); 
    }
  } else {
    for( vc = vn-1; vc > v0; --vc ) {
      CALL_VERTEX_OR_VERTEX_DATA( vc->data ); 
    }
  }
  CALL_END_OR_END_DATA();
  return TRUE;
}

Referenced by gluTessEndPolygon().

__gl_renderMesh()

void __gl_renderMesh	(	GLUtesselator *	tess,
		GLUmesh *	mesh
	)

Definition at line 82 of file render.c.

{
  GLUface *f;
 
  /* Make a list of separate triangles so we can render them all at once */
  tess->lonelyTriList = NULL;
 
  for( f = mesh->fHead.next; f != &mesh->fHead; f = f->next ) {
    f->marked = FALSE;
  }
  for( f = mesh->fHead.next; f != &mesh->fHead; f = f->next ) {
 
    /* We examine all faces in an arbitrary order.  Whenever we find
     * an unprocessed face F, we output a group of faces including F
     * whose size is maximum.
     */
    if( f->inside && ! f->marked ) {
      RenderMaximumFaceGroup( tess, f );
      assert( f->marked );
    }
  }
  if( tess->lonelyTriList != NULL ) {
    RenderLonelyTriangles( tess, tess->lonelyTriList );
    tess->lonelyTriList = NULL;
  }
}

Referenced by gluTessEndPolygon().

ComputeNormal()

static int ComputeNormal ( GLUtesselator *  tess, GLdouble  norm[3], int  check  )

static

Definition at line 362 of file render.c.

{
  CachedVertex *v0 = tess->cache;
  CachedVertex *vn = v0 + tess->cacheCount;
  CachedVertex *vc;
  GLdouble dot, xc, yc, zc, xp, yp, zp, n[3];
  int sign = 0;
 
  /* Find the polygon normal.  It is important to get a reasonable
   * normal even when the polygon is self-intersecting (eg. a bowtie).
   * Otherwise, the computed normal could be very tiny, but perpendicular
   * to the true plane of the polygon due to numerical noise.  Then all
   * the triangles would appear to be degenerate and we would incorrectly
   * decompose the polygon as a fan (or simply not render it at all).
   *
   * We use a sum-of-triangles normal algorithm rather than the more
   * efficient sum-of-trapezoids method (used in CheckOrientation()
   * in normal.c).  This lets us explicitly reverse the signed area
   * of some triangles to get a reasonable normal in the self-intersecting
   * case.
   */
  if( ! check ) {
    norm[0] = norm[1] = norm[2] = 0.0;
  }
 
  vc = v0 + 1;
  xc = vc->coords[0] - v0->coords[0];
  yc = vc->coords[1] - v0->coords[1];
  zc = vc->coords[2] - v0->coords[2];
  while( ++vc < vn ) {
    xp = xc; yp = yc; zp = zc;
    xc = vc->coords[0] - v0->coords[0];
    yc = vc->coords[1] - v0->coords[1];
    zc = vc->coords[2] - v0->coords[2];
 
    /* Compute (vp - v0) cross (vc - v0) */
    n[0] = yp*zc - zp*yc;
    n[1] = zp*xc - xp*zc;
    n[2] = xp*yc - yp*xc;
 
    dot = n[0]*norm[0] + n[1]*norm[1] + n[2]*norm[2];
    if( ! check ) {
      /* Reverse the contribution of back-facing triangles to get
       * a reasonable normal for self-intersecting polygons (see above)
       */
      if( dot >= 0 ) {
    norm[0] += n[0]; norm[1] += n[1]; norm[2] += n[2];
      } else {
    norm[0] -= n[0]; norm[1] -= n[1]; norm[2] -= n[2];
      }
    } else if( dot != 0 ) {
      /* Check the new orientation for consistency with previous triangles */
      if( dot > 0 ) {
    if( sign < 0 ) return SIGN_INCONSISTENT;
    sign = 1;
      } else {
    if( sign > 0 ) return SIGN_INCONSISTENT;
    sign = -1;
      }
    }
  }
  return sign;
}

Referenced by __gl_renderCache().

MaximumFan()

static struct FaceCount MaximumFan ( GLUhalfEdge *  eOrig )

static

Definition at line 158 of file render.c.

{
  /* eOrig->Lface is the face we want to render.  We want to find the size
   * of a maximal fan around eOrig->Org.  To do this we just walk around
   * the origin vertex as far as possible in both directions.
   */
  struct FaceCount newFace = { 0, NULL, &RenderFan };
  GLUface *trail = NULL;
  GLUhalfEdge *e;
 
  for( e = eOrig; ! Marked( e->Lface ); e = e->Onext ) {
    AddToTrail( e->Lface, trail );
    ++newFace.size;
  }
  for( e = eOrig; ! Marked( e->Rface ); e = e->Oprev ) {
    AddToTrail( e->Rface, trail );
    ++newFace.size;
  }
  newFace.eStart = e;
  /*LINTED*/
  FreeTrail( trail );
  return newFace;
}

Referenced by RenderMaximumFaceGroup().

MaximumStrip()

static struct FaceCount MaximumStrip ( GLUhalfEdge *  eOrig )

static

Definition at line 185 of file render.c.

{
  /* Here we are looking for a maximal strip that contains the vertices
   * eOrig->Org, eOrig->Dst, eOrig->Lnext->Dst (in that order or the
   * reverse, such that all triangles are oriented CCW).
   *
   * Again we walk forward and backward as far as possible.  However for
   * strips there is a twist: to get CCW orientations, there must be
   * an *even* number of triangles in the strip on one side of eOrig.
   * We walk the strip starting on a side with an even number of triangles;
   * if both side have an odd number, we are forced to shorten one side.
   */
  struct FaceCount newFace = { 0, NULL, &RenderStrip };
  long headSize = 0, tailSize = 0;
  GLUface *trail = NULL;
  GLUhalfEdge *e, *eTail, *eHead;
 
  for( e = eOrig; ! Marked( e->Lface ); ++tailSize, e = e->Onext ) {
    AddToTrail( e->Lface, trail );
    ++tailSize;
    e = e->Dprev;
    if( Marked( e->Lface )) break;
    AddToTrail( e->Lface, trail );
  }
  eTail = e;
 
  for( e = eOrig; ! Marked( e->Rface ); ++headSize, e = e->Dnext ) {
    AddToTrail( e->Rface, trail );
    ++headSize;
    e = e->Oprev;
    if( Marked( e->Rface )) break;
    AddToTrail( e->Rface, trail );
  }
  eHead = e;
 
  newFace.size = tailSize + headSize;
  if( IsEven( tailSize )) {
    newFace.eStart = eTail->Sym;
  } else if( IsEven( headSize )) {
    newFace.eStart = eHead;
  } else {
    /* Both sides have odd length, we must shorten one of them.  In fact,
     * we must start from eHead to guarantee inclusion of eOrig->Lface.
     */
    --newFace.size;
    newFace.eStart = eHead->Onext;
  }
  /*LINTED*/
  FreeTrail( trail );
  return newFace;
}

Referenced by RenderMaximumFaceGroup().

RenderFan()

static void RenderFan ( GLUtesselator *  tess, GLUhalfEdge *  eStart, long  size  )

static

Definition at line 283 of file render.c.

{
  /* Render as many CCW triangles as possible in a fan starting from
   * edge "e".  The fan *should* contain exactly "size" triangles
   * (otherwise we've goofed up somewhere).
   */
  CALL_BEGIN_OR_BEGIN_DATA( GL_TRIANGLE_FAN ); 
  CALL_VERTEX_OR_VERTEX_DATA( e->Org->data ); 
  CALL_VERTEX_OR_VERTEX_DATA( e->Dst->data ); 
 
  while( ! Marked( e->Lface )) {
    e->Lface->marked = TRUE;
    --size;
    e = e->Onext;
    CALL_VERTEX_OR_VERTEX_DATA( e->Dst->data ); 
  }
 
  assert( size == 0 );
  CALL_END_OR_END_DATA();
}

Referenced by MaximumFan().

RenderLonelyTriangles()

static void RenderLonelyTriangles ( GLUtesselator *  tess, GLUface *  head  )

static

Definition at line 248 of file render.c.

{
  /* Now we render all the separate triangles which could not be
   * grouped into a triangle fan or strip.
   */
  GLUhalfEdge *e;
  int newState;
  int edgeState = -1;   /* force edge state output for first vertex */
 
  CALL_BEGIN_OR_BEGIN_DATA( GL_TRIANGLES );
 
  for( ; f != NULL; f = f->trail ) {
    /* Loop once for each edge (there will always be 3 edges) */
 
    e = f->anEdge;
    do {
      if( tess->flagBoundary ) {
    /* Set the "edge state" to TRUE just before we output the
     * first vertex of each edge on the polygon boundary.
     */
    newState = ! e->Rface->inside;
    if( edgeState != newState ) {
      edgeState = newState;
          CALL_EDGE_FLAG_OR_EDGE_FLAG_DATA( edgeState );
    }
      }
      CALL_VERTEX_OR_VERTEX_DATA( e->Org->data );
 
      e = e->Lnext;
    } while( e != f->anEdge );
  }
  CALL_END_OR_END_DATA();
}

Referenced by __gl_renderMesh().

RenderMaximumFaceGroup()

static void RenderMaximumFaceGroup ( GLUtesselator *  tess, GLUface *  fOrig  )

static

Definition at line 110 of file render.c.

{
  /* We want to find the largest triangle fan or strip of unmarked faces
   * which includes the given face fOrig.  There are 3 possible fans
   * passing through fOrig (one centered at each vertex), and 3 possible
   * strips (one for each CCW permutation of the vertices).  Our strategy
   * is to try all of these, and take the primitive which uses the most
   * triangles (a greedy approach).
   */
  GLUhalfEdge *e = fOrig->anEdge;
  struct FaceCount max, newFace;
 
  max.size = 1;
  max.eStart = e;
  max.render = &RenderTriangle;
 
  if( ! tess->flagBoundary ) {
    newFace = MaximumFan( e ); if( newFace.size > max.size ) { max = newFace; }
    newFace = MaximumFan( e->Lnext ); if( newFace.size > max.size ) { max = newFace; }
    newFace = MaximumFan( e->Lprev ); if( newFace.size > max.size ) { max = newFace; }
 
    newFace = MaximumStrip( e ); if( newFace.size > max.size ) { max = newFace; }
    newFace = MaximumStrip( e->Lnext ); if( newFace.size > max.size ) { max = newFace; }
    newFace = MaximumStrip( e->Lprev ); if( newFace.size > max.size ) { max = newFace; }
  }
  (*(max.render))( tess, max.eStart, max.size );
}

Referenced by __gl_renderMesh().

RenderStrip()

static void RenderStrip ( GLUtesselator *  tess, GLUhalfEdge *  eStart, long  size  )

static

Definition at line 305 of file render.c.

{
  /* Render as many CCW triangles as possible in a strip starting from
   * edge "e".  The strip *should* contain exactly "size" triangles
   * (otherwise we've goofed up somewhere).
   */
  CALL_BEGIN_OR_BEGIN_DATA( GL_TRIANGLE_STRIP );
  CALL_VERTEX_OR_VERTEX_DATA( e->Org->data ); 
  CALL_VERTEX_OR_VERTEX_DATA( e->Dst->data ); 
 
  while( ! Marked( e->Lface )) {
    e->Lface->marked = TRUE;
    --size;
    e = e->Dprev;
    CALL_VERTEX_OR_VERTEX_DATA( e->Org->data ); 
    if( Marked( e->Lface )) break;
 
    e->Lface->marked = TRUE;
    --size;
    e = e->Onext;
    CALL_VERTEX_OR_VERTEX_DATA( e->Dst->data ); 
  }
 
  assert( size == 0 );
  CALL_END_OR_END_DATA();
}

Referenced by MaximumStrip().

RenderTriangle()

static void RenderTriangle ( GLUtesselator *  tess, GLUhalfEdge *  eStart, long  size  )

static

Definition at line 238 of file render.c.

{
  /* Just add the triangle to a triangle list, so we can render all
   * the separate triangles at once.
   */
  assert( size == 1 );
  AddToTrail( e->Lface, tess->lonelyTriList );
}

Referenced by RenderMaximumFaceGroup().