trap.h

Go to the documentation of this file.

 
/*************************************************************************
 * Copyright (c) 2011 AT&T Intellectual Property 
 * All rights reserved. This program and the accompanying materials
 * are made available under the terms of the Eclipse Public License v1.0
 * which accompanies this distribution, and is available at
 * https://www.eclipse.org/legal/epl-v10.html
 *
 * Contributors: Details at https://graphviz.org
 *************************************************************************/
 
#pragma once
 
#include <math.h>
#include <stdbool.h>
#include <stddef.h>
#include <stdint.h>
#include <util/list.h>
 
/* Segment attributes */
 
typedef struct {
  pointf v0, v1;       /* two endpoints */ 
  bool is_inserted;      /* inserted in trapezoidation yet ? */
  size_t root0, root1;   
  int next;         /* Next logical segment */
  int prev;         /* Previous segment */
} segment_t;
 
 
/* Trapezoid attributes */
 
typedef struct {
  int lseg, rseg;       /* two adjoining segments */
  pointf hi, lo;       /* max/min y-values */ 
  size_t u0, u1;
  size_t d0, d1;
  size_t sink;          
  size_t usave;         
  int uside;            
  int state;
} trap_t; 
 
static inline bool is_valid_trap(size_t index) {
  return index != 0 && index != SIZE_MAX;
}
 
DEFINE_LIST(traps, trap_t)
 
#define ST_VALID 1      /* for trapezium state */
#define ST_INVALID 2
 
#define C_EPS 1.0e-7        /* tolerance value: Used for making */
                /* all decisions about collinearity or */
                /* left/right of segment. Decrease */
                /* this value if the input points are */
                /* spaced very close together */
 
static inline bool fp_equal(double s, double t) { return fabs(s - t) <= C_EPS; }
 
static inline int dfp_cmp(double f1, double f2) {
  double d = f1 - f2;
  if (d < -C_EPS)
    return -1;
  if (d > C_EPS)
    return 1;
  return 0;
}
 
static inline bool equal_to(pointf v0, pointf v1) {
  return fp_equal(v0.y, v1.y) && fp_equal(v0.x, v1.x);
}
 
static inline bool greater_than(pointf v0, pointf v1) {
  return v0.y > v1.y + C_EPS ? true
                             : (v0.y < v1.y - C_EPS ? false : v0.x > v1.x);
}
 
extern traps_t construct_trapezoids(int, segment_t*, int*);