diff options
| author | Bruce Momjian | 2003-08-04 00:43:34 +0000 |
|---|---|---|
| committer | Bruce Momjian | 2003-08-04 00:43:34 +0000 |
| commit | 089003fb462fcce46c02bf47322b429f73c33c50 (patch) | |
| tree | 77d78bc3a149df06f5603f60200a6ab363336624 /contrib/earthdistance/earthdistance.c | |
| parent | 63354a0228a1dbc4a0d5ddc8ecdd8326349d2100 (diff) | |
pgindent run.
Diffstat (limited to 'contrib/earthdistance/earthdistance.c')
| -rw-r--r-- | contrib/earthdistance/earthdistance.c | 13 |
1 files changed, 7 insertions, 6 deletions
diff --git a/contrib/earthdistance/earthdistance.c b/contrib/earthdistance/earthdistance.c index 19c81a5783f..28ce40378cf 100644 --- a/contrib/earthdistance/earthdistance.c +++ b/contrib/earthdistance/earthdistance.c @@ -6,8 +6,8 @@ /* Earth's radius is in statute miles. */ -const double EARTH_RADIUS = 3958.747716; -const double TWO_PI = 2.0 * M_PI; +const double EARTH_RADIUS = 3958.747716; +const double TWO_PI = 2.0 * M_PI; double *geo_distance(Point *pt1, Point *pt2); @@ -66,10 +66,11 @@ geo_distance(Point *pt1, Point *pt2) if (longdiff > M_PI) longdiff = TWO_PI - longdiff; - sino = sqrt(sin(fabs(lat1-lat2)/2.)*sin(fabs(lat1-lat2)/2.) + - cos(lat1) * cos(lat2) * sin(longdiff/2.)*sin(longdiff/2.)); - if (sino > 1.) sino = 1.; - *resultp = 2. * EARTH_RADIUS * asin(sino); + sino = sqrt(sin(fabs(lat1 - lat2) / 2.) * sin(fabs(lat1 - lat2) / 2.) + + cos(lat1) * cos(lat2) * sin(longdiff / 2.) * sin(longdiff / 2.)); + if (sino > 1.) + sino = 1.; + *resultp = 2. * EARTH_RADIUS * asin(sino); return resultp; } |