208{
209 double result = 0;
210#if POSTGIS_DEBUG_LEVEL >= 4
211 double result2 = 0;
212#endif
213
214 if ( (lat1==lat2) && (long1 == long2) )
215 {
216 return 0.0;
217 }
218
220
221#if POSTGIS_DEBUG_LEVEL >= 4
223
224 POSTGIS_DEBUGF(4, "delta = %lf, skae says: %.15lf,2 circle says: %.15lf",
225 (result2-result),result,result2);
226 POSTGIS_DEBUGF(4, "2 circle says: %.15lf",result2);
227#endif
228
229 if (result != result)
230
231
232 {
234 lat2,long2, sphere);
235 }
236
237 return result;
238}
double distance_ellipse_calculation(double lat1, double long1, double lat2, double long2, SPHEROID *sphere)
double distance_sphere_method(double lat1, double long1, double lat2, double long2, SPHEROID *sphere)