40# define _glu_dprintf printf
47#include "subdivider.h"
50#include "trimvertex.h"
51#include "simplemath.h"
56 return bbox( a->param[p], b->param[p], c->param[p],
57 a->param[1-p], b->param[1-p], c->param[1-p] );
61Subdivider::ccwTurn_sr( Arc_ptr j1, Arc_ptr j2 )
63 register TrimVertex *v1 = &j1->pwlArc->pts[j1->pwlArc->npts-1];
64 register TrimVertex *v1last = &j1->pwlArc->pts[0];
65 register TrimVertex *v2 = &j2->pwlArc->pts[0];
66 register TrimVertex *v2last = &j2->pwlArc->pts[j2->pwlArc->npts-1];
67 register TrimVertex *v1next = v1-1;
68 register TrimVertex *v2next = v2+1;
71 assert( v1 != v1last );
72 assert( v2 != v2last );
75 _glu_dprintf(
"arc_ccw_turn, p = %d\n", 0 );
79 if( v1->param[0] == v1next->param[0] && v2->param[0] == v2next->param[0] )
85 if( v1->param[1] < v2->param[1] )
87 else if( v1->param[1] > v2->param[1] )
91 if( v1next->param[0] < v2next->param[0] ) {
93 _glu_dprintf(
"case a\n" );
95 assert( v1->param[0] <= v1next->param[0] );
96 assert( v2->param[0] <= v1next->param[0] );
97 switch( bbox( v2, v2next, v1next, 1 ) ) {
101 sgn = ccw( v1next, v2, v2next );
106 _glu_dprintf(
"decr\n" );
111 _glu_dprintf(
"no good results\n" );
120 }
else if( v1next->param[0] > v2next->param[0] ) {
122 _glu_dprintf(
"case b\n" );
124 assert( v1->param[0] <= v2next->param[0] );
125 assert( v2->param[0] <= v2next->param[0] );
126 switch( bbox( v1, v1next, v2next, 1 ) ) {
130 sgn = ccw( v1next, v1, v2next );
135 _glu_dprintf(
"incr\n" );
140 _glu_dprintf(
"no good results\n" );
151 _glu_dprintf(
"case ab\n" );
153 if( v1next->param[1] < v2next->param[1] )
155 else if( v1next->param[1] > v2next->param[1] )
159 _glu_dprintf(
"incr\n" );
164 _glu_dprintf(
"no good results\n" );
174Subdivider::ccwTurn_sl( Arc_ptr j1, Arc_ptr j2 )
176 register TrimVertex *v1 = &j1->pwlArc->pts[j1->pwlArc->npts-1];
177 register TrimVertex *v1last = &j1->pwlArc->pts[0];
178 register TrimVertex *v2 = &j2->pwlArc->pts[0];
179 register TrimVertex *v2last = &j2->pwlArc->pts[j2->pwlArc->npts-1];
180 register TrimVertex *v1next = v1-1;
181 register TrimVertex *v2next = v2+1;
184 assert( v1 != v1last );
185 assert( v2 != v2last );
188 _glu_dprintf(
"arc_ccw_turn, p = %d\n", 0 );
192 if( v1->param[0] == v1next->param[0] && v2->param[0] == v2next->param[0] )
195 if( v2next->param[0] > v2->param[0] || v1next->param[0] > v1->param[0] )
199 if( v1->param[1] < v2->param[1] )
201 else if( v1->param[1] > v2->param[1] )
205 if( v1next->param[0] > v2next->param[0] ) {
207 _glu_dprintf(
"case c\n" );
209 assert( v1->param[0] >= v1next->param[0] );
210 assert( v2->param[0] >= v1next->param[0] );
211 switch( bbox( v2next, v2, v1next, 1 ) ) {
215 sgn = ccw( v1next, v2, v2next );
221 _glu_dprintf(
"decr\n" );
225 _glu_dprintf(
"no good results\n" );
234 }
else if( v1next->param[0] < v2next->param[0] ) {
236 _glu_dprintf(
"case d\n" );
238 assert( v1->param[0] >= v2next->param[0] );
239 assert( v2->param[0] >= v2next->param[0] );
240 switch( bbox( v1next, v1, v2next, 1 ) ) {
244 sgn = ccw( v1next, v1, v2next );
250 _glu_dprintf(
"incr\n" );
254 _glu_dprintf(
"no good results\n" );
265 _glu_dprintf(
"case cd\n" );
267 if( v1next->param[1] < v2next->param[1] )
269 else if( v1next->param[1] > v2next->param[1] )
274 _glu_dprintf(
"incr\n" );
278 _glu_dprintf(
"no good results\n" );
288Subdivider::ccwTurn_tr( Arc_ptr j1, Arc_ptr j2 )
290 register TrimVertex *v1 = &j1->pwlArc->pts[j1->pwlArc->npts-1];
291 register TrimVertex *v1last = &j1->pwlArc->pts[0];
292 register TrimVertex *v2 = &j2->pwlArc->pts[0];
293 register TrimVertex *v2last = &j2->pwlArc->pts[j2->pwlArc->npts-1];
294 register TrimVertex *v1next = v1-1;
295 register TrimVertex *v2next = v2+1;
298 assert( v1 != v1last );
299 assert( v2 != v2last );
302 _glu_dprintf(
"arc_ccw_turn, p = %d\n", 1 );
306 if( v1->param[1] == v1next->param[1] && v2->param[1] == v2next->param[1] )
309 if( v2next->param[1] < v2->param[1] || v1next->param[1] < v1->param[1] )
313 if( v1->param[0] < v2->param[0] )
315 else if( v1->param[0] > v2->param[0] )
319 if( v1next->param[1] < v2next->param[1] ) {
321 _glu_dprintf(
"case a\n" );
323 assert( v1->param[1] <= v1next->param[1] );
324 assert( v2->param[1] <= v1next->param[1] );
325 switch( bbox( v2, v2next, v1next, 0 ) ) {
329 sgn = ccw( v1next, v2, v2next );
334 _glu_dprintf(
"decr\n" );
339 _glu_dprintf(
"no good results\n" );
348 }
else if( v1next->param[1] > v2next->param[1] ) {
350 _glu_dprintf(
"case b\n" );
352 assert( v1->param[1] <= v2next->param[1] );
353 assert( v2->param[1] <= v2next->param[1] );
354 switch( bbox( v1, v1next, v2next, 0 ) ) {
358 sgn = ccw( v1next, v1, v2next );
363 _glu_dprintf(
"incr\n" );
368 _glu_dprintf(
"no good results\n" );
379 _glu_dprintf(
"case ab\n" );
381 if( v1next->param[0] < v2next->param[0] )
383 else if( v1next->param[0] > v2next->param[0] )
387 _glu_dprintf(
"incr\n" );
392 _glu_dprintf(
"no good results\n" );
402Subdivider::ccwTurn_tl( Arc_ptr j1, Arc_ptr j2 )
404 register TrimVertex *v1 = &j1->pwlArc->pts[j1->pwlArc->npts-1];
405 register TrimVertex *v1last = &j1->pwlArc->pts[0];
406 register TrimVertex *v2 = &j2->pwlArc->pts[0];
407 register TrimVertex *v2last = &j2->pwlArc->pts[j2->pwlArc->npts-1];
408 register TrimVertex *v1next = v1-1;
409 register TrimVertex *v2next = v2+1;
412 assert( v1 != v1last );
413 assert( v2 != v2last );
416 _glu_dprintf(
"arc_ccw_turn, p = %d\n", 1 );
420 if( v1->param[1] == v1next->param[1] && v2->param[1] == v2next->param[1] )
423 if( v2next->param[1] > v2->param[1] || v1next->param[1] > v1->param[1] )
427 if( v1->param[0] < v2->param[0] )
429 else if( v1->param[0] > v2->param[0] )
433 if( v1next->param[1] > v2next->param[1] ) {
435 _glu_dprintf(
"case c\n" );
437 assert( v1->param[1] >= v1next->param[1] );
438 assert( v2->param[1] >= v1next->param[1] );
439 switch( bbox( v2next, v2, v1next, 0 ) ) {
443 sgn = ccw( v1next, v2, v2next );
449 _glu_dprintf(
"decr\n" );
453 _glu_dprintf(
"no good results\n" );
462 }
else if( v1next->param[1] < v2next->param[1] ) {
464 _glu_dprintf(
"case d\n" );
465 assert( v1->param[1] >= v2next->param[1] );
466 assert( v2->param[1] >= v2next->param[1] );
468 switch( bbox( v1next, v1, v2next, 0 ) ) {
472 sgn = ccw( v1next, v1, v2next );
478 _glu_dprintf(
"incr\n" );
482 _glu_dprintf(
"no good results\n" );
493 _glu_dprintf(
"case cd\n" );
495 if( v1next->param[0] < v2next->param[0] )
497 else if( v1next->param[0] > v2next->param[0] )
502 _glu_dprintf(
"incr\n" );
506 _glu_dprintf(
"no good results\n" );
518Subdivider::bbox(
register REAL sa,
register REAL sb,
register REAL sc,
519 register REAL ta,
register REAL tb,
register REAL tc )
522Subdivider::bbox(
register REAL sa,
register REAL sb,
register REAL sc,
523 register REAL ,
register REAL ,
register REAL )
534 }
else if( sb <= sc ) {
539 }
else if( sa > sb ) {
542 }
else if( sb >= sc ) {
550 }
else if( sb > sc ) {
569 REAL d = det3( a, b, c );
570 if( glu_abs(d) < 0.0001 )
return -1;
571 return (d < 0.0) ? 0 : 1;