Benchmark: Haversine Distance tests

Script Preparation code:

var EARTH_RADIUS_M = 6371e3;
var EARTH_DIAMETER_M = EARTH_RADIUS_M * 2;
var DEG_TO_RAD = Math.PI / 180;
var p1 = { lat: 43.63625, lon: -79.40454 };
var p2 = { lat: 43.69851, lon: -79.39699 };

function haversineClassic(lat1, lon1, lat2, lon2) {
  var φ1 = lat1 * DEG_TO_RAD;
  var φ2 = lat2 * DEG_TO_RAD;
  var Δφ = (lat2-lat1) * DEG_TO_RAD;
  var Δλ = (lon2-lon1) * DEG_TO_RAD;

  var a = Math.sin(Δφ/2) * Math.sin(Δφ/2) +
          Math.cos(φ1) * Math.cos(φ2) *
          Math.sin(Δλ/2) * Math.sin(Δλ/2);
  var c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a));

  var d = EARTH_RADIUS_M * c;
  
  return d;
}

function haversineOptimized(lat1, lon1, lat2, lon2) {
  var p = 0.017453292519943295;    // Math.PI / 180
  var c = Math.cos;
  var a = 0.5 - c((lat2 - lat1) * p)/2 + 
      c(lat1 * p) * c(lat2 * p) * 
      (1 - c((lon2 - lon1) * p))/2;

  return 12742 * Math.asin(Math.sqrt(a)); // 2 * R; R = 6371 km  
}

function haversineMineATan(lat1, lon1, lat2, lon2) {
  var fDP = Math.sin( (lat2 - lat1) * DEG_TO_RAD / 2 );
  var fDL = Math.sin( (lon2 - lon1) * DEG_TO_RAD / 2 );
  var fP1 = Math.cos( lat1 * DEG_TO_RAD );
  var fP2 = Math.cos( lat2 * DEG_TO_RAD );
  var a = fDP * fDP + fP1 * fP2 * fDL * fDL;
  return EARTH_DIAMETER_M * Math.atan2(Math.sqrt(a), Math.sqrt(1-a));
}
function haversineMineASin(lat1, lon1, lat2, lon2) {
  var fDP = Math.sin( (lat2 - lat1) * DEG_TO_RAD / 2 );
  var fDL = Math.sin( (lon2 - lon1) * DEG_TO_RAD / 2 );
  var fP1 = Math.cos( lat1 * DEG_TO_RAD );
  var fP2 = Math.cos( lat2 * DEG_TO_RAD );
  var a = fDP * fDP + fP1 * fP2 * fDL * fDL;
  return EARTH_DIAMETER_M * Math.asin(Math.sqrt(a));
}

function spherical(lat1, lon1, lat2, lon2) {
  var φ1 = lat1 * DEG_TO_RAD, φ2 = lat2 * DEG_TO_RAD, Δλ = (lon2-lon1) * DEG_TO_RAD;
  return Math.acos( Math.sin(φ1) * Math.sin(φ2) + Math.cos(φ1) * Math.cos(φ2) * Math.cos(Δλ) ) * EARTH_RADIUS_M;
}

function equirectangular(lat1, lon1, lat2, lon2) {
  var x = (lon2-lon1)*DEG_TO_RAD * Math.cos( (lat1+lat2)*DEG_TO_RAD/2 );
  var y = (lat2-lat1)*DEG_TO_RAD;
  return EARTH_RADIUS_M * Math.sqrt( x*x + y*y )
}

Tests:

Haversine Classic
haversineClassic(p1.lat, p1.lon, p2.lat, p2.lon);
Haversine Optimized
haversineOptimized(p1.lat, p1.lon, p2.lat, p2.lon);
Haversine Mine ATan
haversineMineATan(p1.lat, p1.lon, p2.lat, p2.lon);
Haversine Mine ASin
haversineMineASin(p1.lat, p1.lon, p2.lat, p2.lon);
Spherical Law of Cosines
spherical(p1.lat, p1.lon, p2.lat, p2.lon);
Equirectangular (a little less accurate)
equirectangular(p1.lat, p1.lon, p2.lat, p2.lon);

Rendered benchmark preparation results:

Suite status: <idle, ready to run>

Previous results

Test case name	Result
Haversine Classic
Haversine Optimized
Haversine Mine ATan
Haversine Mine ASin
Spherical Law of Cosines
Equirectangular (a little less accurate)

Fastest: N/A

Slowest: N/A

ⓘ Get performance results from real users on your website with Request Metrics.

Website Performance Monitoring by Request Metrics

Latest run results:

Run details: (Test run date: 4 months ago)

User agent: Mozilla/5.0 (Windows NT 10.0; Win64; x64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/120.0.0.0 Safari/537.36

Browser/OS: Chrome 120 on Windows

View result in a separate tab

Test name	Executions per second
Haversine Classic	426100.7 Ops/sec
Haversine Optimized	1009763.0 Ops/sec
Haversine Mine ATan	381929.8 Ops/sec
Haversine Mine ASin	435878.5 Ops/sec
Spherical Law of Cosines	485645.0 Ops/sec
Equirectangular (a little less accurate)	807031.6 Ops/sec

Haversine Distance tests (version: 4)

Comparing performance of: Haversine Classic vs Haversine Optimized vs Haversine Mine ATan vs Haversine Mine ASin vs Spherical Law of Cosines vs Equirectangular (a little less accurate)

Created: 6 years ago by: Registered User

Jump to the latest result

Haversine Classic

Haversine Optimized

Haversine Mine ATan

Haversine Mine ASin

Spherical Law of Cosines

Equirectangular (a little less accurate)

Suite status: <idle, ready to run>

Experimental features:

Fastest: N/A

Slowest: N/A