src/glu/sgi/libnurbs/internals/reader.cc

   1 /*
   2 ** License Applicability. Except to the extent portions of this file are
   3 ** made subject to an alternative license as permitted in the SGI Free
   4 ** Software License B, Version 1.1 (the "License"), the contents of this
   5 ** file are subject only to the provisions of the License. You may not use
   6 ** this file except in compliance with the License. You may obtain a copy
   7 ** of the License at Silicon Graphics, Inc., attn: Legal Services, 1600
   8 ** Amphitheatre Parkway, Mountain View, CA 94043-1351, or at:
   9 **
  10 ** http://oss.sgi.com/projects/FreeB
  11 **
  12 ** Note that, as provided in the License, the Software is distributed on an
  13 ** "AS IS" basis, with ALL EXPRESS AND IMPLIED WARRANTIES AND CONDITIONS
  14 ** DISCLAIMED, INCLUDING, WITHOUT LIMITATION, ANY IMPLIED WARRANTIES AND
  15 ** CONDITIONS OF MERCHANTABILITY, SATISFACTORY QUALITY, FITNESS FOR A
  16 ** PARTICULAR PURPOSE, AND NON-INFRINGEMENT.
  17 **
  18 ** Original Code. The Original Code is: OpenGL Sample Implementation,
  19 ** Version 1.2.1, released January 26, 2000, developed by Silicon Graphics,
  20 ** Inc. The Original Code is Copyright (c) 1991-2000 Silicon Graphics, Inc.
  21 ** Copyright in any portions created by third parties is as indicated
  22 ** elsewhere herein. All Rights Reserved.
  23 **
  24 ** Additional Notice Provisions: The application programming interfaces
  25 ** established by SGI in conjunction with the Original Code are The
  26 ** OpenGL(R) Graphics System: A Specification (Version 1.2.1), released
  27 ** April 1, 1999; The OpenGL(R) Graphics System Utility Library (Version
  28 ** 1.3), released November 4, 1998; and OpenGL(R) Graphics with the X
  29 ** Window System(R) (Version 1.3), released October 19, 1998. This software
  30 ** was created using the OpenGL(R) version 1.2.1 Sample Implementation
  31 ** published by SGI, but has not been independently verified as being
  32 ** compliant with the OpenGL(R) version 1.2.1 Specification.
  33 */
  34
  35 /*
  36  *  reader.c++
  37  *
  38  */
  39
  40 #include <stdio.h>
  41 #include "glimports.h"
  42 #include "nurbsconsts.h"
  43 #include "reader.h"
  44 #include "trimvertex.h"
  45 #include "simplemath.h"
  46
  47 //when read a pwlCurve, if two consecutive points are the same, then
  48 //eliminate one of them. This makes the tessellator more robust. The spec
  49 //assumes the application makes sure there are no redundant points.
  50 //but in Inspector, the trim curves seem to have redundant points a lot.
  51 //I guess other similar users may have the same problem.
  52
  53 #define ELIMINATE_REDUNDANT_POINTS
  54
  55 #ifdef  ELIMINATE_REDUNDANT_POINTS
  56 #define equal(x,y) ( glu_abs(x-y) <= 0.00001)
  57 #endif
  58
  59 #ifdef ELIMINATE_REDUNDANT_POINTS
  60 O_pwlcurve::O_pwlcurve( long _type, long count, INREAL *array, long byte_stride, TrimVertex *trimpts )
  61 {
  62     next = 0;
  63     used = 0;
  64     owner = 0;
  65     pts = trimpts;
  66     npts = (int) count;
  67     save = 0;
  68     int i;
  69
  70     /* copy user data into internal trimming data structures */
  71     switch( _type ) {
  72         case N_P2D: {
  73             TrimVertex *v = pts;
  74             TrimVertex *prev = NULL;
  75             int num = 0;
  76             int doit;
  77             for(i=0; i<count; i++) {
  78                 doit = 1;
  79                 if(prev != NULL)
  80                   {
  81                     if(equal(prev->param[0], array[0]) && equal(prev->param[1], array[1]))
  82                       {
  83                         doit = 0;
  84                       }
  85                   }
  86
  87                 if(doit)
  88                   {
  89                     v->param[0] = (REAL) array[0];
  90                     v->param[1] = (REAL) array[1];
  91                     prev = v;
  92                     v++;
  93                     num++;
  94                   }
  95                 array = (INREAL *) (((char *) array) + byte_stride);
  96               }
  97             npts = num;
  98             break;
  99         }
 100         case N_P2DR: {
 101             TrimVertex *v = pts;
 102             for( TrimVertex *lastv = v + count; v != lastv; v++ ) {
 103                 v->param[0] = (REAL) array[0] / (REAL) array[2];
 104                 v->param[1] = (REAL) array[1] / (REAL) array[2];
 105                 array = (INREAL *) (((char *) array) + byte_stride);
 106             }
 107             break;
 108         }
 109     }
 110 }
 111 #else
 112 O_pwlcurve::O_pwlcurve( long _type, long count, INREAL *array, long byte_stride, TrimVertex *trimpts )
 113 {
 114     next = 0;
 115     used = 0;
 116     owner = 0;
 117     pts = trimpts;
 118     npts = (int) count;
 119     save = 0;
 120
 121     /* copy user data into internal trimming data structures */
 122     switch( _type ) {
 123         case N_P2D: {
 124             TrimVertex *v = pts;
 125             for( TrimVertex *lastv = v + count; v != lastv; v++ ) {
 126                 v->param[0] = (REAL) array[0];
 127                 v->param[1] = (REAL) array[1];
 128                 array = (INREAL *) (((char *) array) + byte_stride);
 129             }
 130             break;
 131         }
 132         case N_P2DR: {
 133             TrimVertex *v = pts;
 134             for( TrimVertex *lastv = v + count; v != lastv; v++ ) {
 135                 v->param[0] = (REAL) array[0] / (REAL) array[2];
 136                 v->param[1] = (REAL) array[1] / (REAL) array[2];
 137                 array = (INREAL *) (((char *) array) + byte_stride);
 138             }
 139             break;
 140         }
 141     }
 142 }
 143 #endif
 144
 145
 146
 147
 148