Skip to content

Instantly share code, notes, and snippets.

@smanek
Created January 28, 2010 15:45
Show Gist options
  • Save smanek/288845 to your computer and use it in GitHub Desktop.
Save smanek/288845 to your computer and use it in GitHub Desktop.
//per David Beazley's suggestion modified the import statements
//yielding a ~25% improvement
from math import *
import time
radius = 6371
def distance(latA, lngA, latB, lngB):
latAr = radians(latA)
lngAr = radians(lngA)
latBr = radians(latB)
lngBr = radians(lngB)
deltaLat = latBr - latAr
deltaLng = lngBr - lngAr
return radius * 2 * asin(math.sqrt(
sin(deltaLat/2)**2 + cos(latAr)
* cos(latBr) * (sin(deltaLng/2)**2)))
def bench():
increment = 2.5
start = time.time()
latA = -90
while(latA <= 90):
lngA = -180
while(lngA <= 180):
latB = -90
while(latB <= 90):
lngB = -180
while(lngB <= 180):
distance(latA, lngA, latB, lngB)
lngB = lngB + increment
latB = latB + increment
lngA = lngA + increment
latA = latA + increment
end = time.time()
print end-start
if __name__ == '__main__':
bench()
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment