19 #include <grass/gis.h>
20 #include <grass/glocale.h>
22 static double min4(
double,
double,
double,
double);
23 static double min2(
double,
double);
54 if (
st->factor <= 0.0) {
75 double G_distance(
double e1,
double n1,
double e2,
double n2)
77 if (
st->projection == PROJECTION_LL)
80 return st->factor * hypot(e1 - e2, n1 - n2);
92 double ax2,
double ay2,
93 double bx1,
double by1,
94 double bx2,
double by2)
101 bx1, by1, bx2, by2, &ra, &rb, &
x, &y) > 0)
121 double x1,
double y1,
double x2,
126 double xq, yq, ra, rb;
133 if (dx == 0.0 && dy == 0.0)
136 if (fabs(dy) > fabs(dx)) {
138 yq = (dx / dy) * (xp - xq) + yp;
142 xq = (dy / dx) * (yp - yq) + xp;
155 "code=%d P=(%f,%f) S=(%f,%f)(%f,%f)"),
156 "G_distance_point_to_line_segment",
t, xp, yp, x1, y1, x2, y2);
161 if (rb >= 0 && rb <= 1.0)
170 static double min4(
double a,
double b,
double c,
double d)
172 return min2(min2(a,
b), min2(c, d));
175 static double min2(
double a,
double b)
177 return a <
b ? a :
b;