switched to inline heaps
[swftools.git] / lib / art / art_svp_point.c
1 /* Libart_LGPL - library of basic graphic primitives
2  * Copyright (C) 1999 Raph Levien
3  *
4  * This library is free software; you can redistribute it and/or
5  * modify it under the terms of the GNU Library General Public
6  * License as published by the Free Software Foundation; either
7  * version 2 of the License, or (at your option) any later version.
8  *
9  * This library is distributed in the hope that it will be useful,
10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
12  * Library General Public License for more details.
13  *
14  * You should have received a copy of the GNU Library General Public
15  * License along with this library; if not, write to the
16  * Free Software Foundation, Inc., 59 Temple Place - Suite 330,
17  * Boston, MA 02111-1307, USA.
18  */
19
20 #include "config.h"
21 #include "art_svp_point.h"
22
23 #include <math.h>
24 #include "art_misc.h"
25
26 #include "art_svp.h"
27
28 /* Determine whether a point is inside, or near, an svp. */
29
30 /* return winding number of point wrt svp */
31 /**
32  * art_svp_point_wind: Determine winding number of a point with respect to svp.
33  * @svp: The svp.
34  * @x: The X coordinate of the point.
35  * @y: The Y coordinate of the point.
36  *
37  * Determine the winding number of the point @x, @y with respect to @svp.
38  *
39  * Return value: the winding number.
40  **/
41 int
42 art_svp_point_wind (ArtSVP *svp, double x, double y)
43 {
44   int i, j;
45   int wind = 0;
46
47   for (i = 0; i < svp->n_segs; i++)
48     {
49       ArtSVPSeg *seg = &svp->segs[i];
50
51       if (seg->bbox.y0 > y)
52         break;
53
54       if (seg->bbox.y1 > y)
55         {
56           if (seg->bbox.x1 < x)
57             wind += seg->dir ? 1 : -1;
58           else if (seg->bbox.x0 <= x)
59             {
60               double x0, y0, x1, y1, dx, dy;
61
62               for (j = 0; j < seg->n_points - 1; j++)
63                 {
64                   if (seg->points[j + 1].y > y)
65                     break;
66                 }
67               x0 = seg->points[j].x;
68               y0 = seg->points[j].y;
69               x1 = seg->points[j + 1].x;
70               y1 = seg->points[j + 1].y;
71
72               dx = x1 - x0;
73               dy = y1 - y0;
74               if ((x - x0) * dy > (y - y0) * dx)
75                 wind += seg->dir ? 1 : -1;
76             }
77         }
78     }
79
80   return wind;
81 }
82
83 /**
84  * art_svp_point_dist: Determine distance between point and svp.
85  * @svp: The svp.
86  * @x: The X coordinate of the point.
87  * @y: The Y coordinate of the point.
88  *
89  * Determines the distance of the point @x, @y to the closest edge in
90  * @svp. A large number is returned if @svp is empty.
91  *
92  * Return value: the distance.
93  **/
94 double
95 art_svp_point_dist (ArtSVP *svp, double x, double y)
96 {
97   int i, j;
98   double dist_sq;
99   double best_sq = -1;
100
101   for (i = 0; i < svp->n_segs; i++)
102     {
103       ArtSVPSeg *seg = &svp->segs[i];
104       for (j = 0; j < seg->n_points - 1; j++)
105         {
106           double x0 = seg->points[j].x;
107           double y0 = seg->points[j].y;
108           double x1 = seg->points[j + 1].x;
109           double y1 = seg->points[j + 1].y;
110
111           double dx = x1 - x0;
112           double dy = y1 - y0;
113
114           double dxx0 = x - x0;
115           double dyy0 = y - y0;
116
117           double dot = dxx0 * dx + dyy0 * dy;
118
119           if (dot < 0)
120             dist_sq = dxx0 * dxx0 + dyy0 * dyy0;
121           else
122             {
123               double rr = dx * dx + dy * dy;
124
125               if (dot > rr)
126                 dist_sq = (x - x1) * (x - x1) + (y - y1) * (y - y1);
127               else
128                 {
129                   double perp = (y - y0) * dx - (x - x0) * dy;
130
131                   dist_sq = perp * perp / rr;
132                 }
133             }
134           if (best_sq < 0 || dist_sq < best_sq)
135             best_sq = dist_sq;
136         }
137     }
138
139   if (best_sq >= 0)
140     return sqrt (best_sq);
141   else
142     return 1e12;
143 }
144