radiansTo method
Implementation
double radiansTo(LCGeoPoint point) {
double d2r = pi / 180.0;
double lat1rad = latitude * d2r;
double long1rad = longitude * d2r;
double lat2rad = point.latitude * d2r;
double long2rad = point.longitude * d2r;
double deltaLat = lat1rad - lat2rad;
double deltaLong = long1rad - long2rad;
double sinDeltaLatDiv2 = sin(deltaLat / 2);
double sinDeltaLongDiv2 = sin(deltaLong / 2);
double a = sinDeltaLatDiv2 * sinDeltaLatDiv2 +
cos(lat1rad) * cos(lat2rad) * sinDeltaLongDiv2 * sinDeltaLongDiv2;
a = min(1.0, a);
return 2 * sin(sqrt(a));
}