polygon intersector: finished active list splay tree optimization
[swftools.git] / lib / gfxpoly / poly.h
1 #ifndef __poly_h__
2 #define __poly_h__
3
4 #include <stdint.h>
5 #include "../q.h"
6
7 //#define DEBUG
8 //#define CHECKS
9 #define SPLAY
10
11 typedef enum {DIR_UP, DIR_DOWN} segment_dir_t;
12 typedef enum {EVENT_CROSS, EVENT_END, EVENT_CORNER, EVENT_START, EVENT_HORIZONTAL} eventtype_t;
13 typedef enum {SLOPE_POSITIVE, SLOPE_NEGATIVE} slope_t;
14
15 typedef struct _point {
16     int32_t x;
17     int32_t y;
18 } point_t;
19
20 typedef struct _fillstyle {
21     char is_filled;
22 } fillstyle_t;
23
24 typedef struct _edge {
25     point_t a;
26     point_t b;
27     fillstyle_t*style;
28     int tmp;
29     struct _edge *next;
30 } edge_t;
31
32 typedef struct _windstate
33 {
34     char is_filled;
35     int wind_nr;
36     int num_polygons;
37 } windstate_t;
38
39 typedef struct _windrule
40 {
41     windstate_t (*start)(int num_polygons);
42     windstate_t (*add)(windstate_t left, fillstyle_t*edge, segment_dir_t dir, int polygon_nr);
43     fillstyle_t* (*diff)(windstate_t*left, windstate_t*right);
44 } windrule_t;
45
46 typedef struct _segment {
47     point_t a;
48     point_t b;
49     point_t delta;
50     double k; //k = a.x*b.y-a.y*b.x = delta.y*a.x - delta.x*a.y (=0 for points on the segment)
51     int minx, maxx;
52     
53     segment_dir_t dir;
54     fillstyle_t*fs;
55     fillstyle_t*fs_out;
56     char fs_out_ok;
57     
58     int polygon_nr;
59     windstate_t wind;
60     int nr;
61
62 #ifdef SPLAY
63     struct _segment*parent;
64     struct _segment*leftchild;
65     struct _segment*rightchild;
66 #endif
67     struct _segment*left;
68     struct _segment*right;
69     char changed;
70
71     point_t pos;
72
73     dict_t scheduled_crossings;
74 } segment_t;
75
76 #define LINE_EQ(p,s) ((double)(s)->delta.y*(p).x - (double)(s)->delta.x*(p).y - (s)->k)
77
78 /* x1 + ((x2-x1)*(y-y1)) / dy = 
79    (x1*(y2-y1) + (x2-x1)*(y-y1)) / dy =
80    (x1*(y2-y)  +  x2    *(y-y1)) / dy =
81    (x1*y2 - x2*y1 + x2*y - y*x1) / dy =
82    (k + x2*y - x1*y) / dy
83    (k + dx*y) / dy
84 */
85 //#define XPOS(s,ypos) ((s)->a.x + ((s)->delta.x * (double)((ypos) - (s)->a.y)) / (s)->delta.y)
86 #define XPOS(s,ypos) (((s)->k + (double)(s)->delta.x*ypos) / (s)->delta.y)
87
88 #define XPOS_INT(s,ypos) ((int)ceil(XPOS((s),ypos)))
89 #define XDIFF(s1,s2,ypos) (((s1)->k + (double)(s1)->delta.x*ypos)*(s2)->delta.y - \
90                            ((s2)->k + (double)(s2)->delta.x*ypos)*(s1)->delta.y)
91
92 typedef struct _gfxpoly {
93     double gridsize;
94     edge_t*edges;
95 } gfxpoly_t;
96
97 void gfxpoly_fail(char*expr, char*file, int line, const char*function);
98
99 gfxpoly_t* gfxpoly_new(double gridsize);
100 char gfxpoly_check(gfxpoly_t*poly);
101 int gfxpoly_size(gfxpoly_t*poly);
102 void gfxpoly_dump(gfxpoly_t*poly);
103 gfxpoly_t* gfxpoly_save(gfxpoly_t*poly, const char*filename);
104 gfxpoly_t* gfxpoly_process(gfxpoly_t*poly, windrule_t*windrule);
105
106 typedef struct _event {
107     eventtype_t type;
108     point_t p;
109     segment_t*s1;
110     segment_t*s2;
111 } event_t;
112
113 #ifndef CHECKS
114 #ifdef assert
115 #undef assert
116 #endif
117 #define assert(x)
118 #else
119 #define assert(x) ((x)?0:gfxpoly_fail(__STRING(x), __FILE__, __LINE__, __PRETTY_FUNCTION__))
120 #endif
121
122
123 #endif