Run results for: geosailing2

Test name	Executions per second
✓ case1	146,516,928 Ops/sec
case6	145,745,600 Ops/sec
case7	145,584,896 Ops/sec
case2	5,831,424 Ops/sec
case3	786,914 Ops/sec
case5	627,956 Ops/sec
case4	580,408 Ops/sec

Script Preparation code:

function geosailing(latA, lngA, latB, lngB) {
  try {
    return (Math.acos(Math.sin(Math.atan(Math.tan(latA * (Math.PI / 180)))) * Math.sin(Math.atan(Math.tan(latB * (Math.PI / 180)))) + Math.cos(Math.atan(Math.tan(latA * (Math.PI / 180)))) * Math.cos(Math.atan(Math.tan(latB * (Math.PI / 180)))) * Math.cos(lngA * (Math.PI / 180) - lngB * (Math.PI / 180))) * 6371.008);
  } catch {
    return 0;
  }
}
var dist_cache = {};

function geosailing2(latA, lngA, latB, lngB) {
  if (dist_cache[latA] && dist_cache[latA][lngA] && dist_cache[latA][lngA][latB] && dist_cache[latA][lngA][latB][lngB]) {
    return dist_cache[latA][lngA][latB][lngB];
  } else {
    var dist = Math.acos(Math.sin(Math.atan(Math.tan(latA * (Math.PI / 180)))) * Math.sin(Math.atan(Math.tan(latB * (Math.PI / 180)))) + Math.cos(Math.atan(Math.tan(latA * (Math.PI / 180)))) * Math.cos(Math.atan(Math.tan(latB * (Math.PI / 180)))) * Math.cos(lngA * (Math.PI / 180) - lngB * (Math.PI / 180))) * 6371.008;
    if (!dist_cache[latA]) dist_cache[latA] = {};
    if (!dist_cache[latA][lngA]) dist_cache[latA][lngA] = {};
    if (!dist_cache[latA][lngA][latB]) dist_cache[latA][lngA][latB] = {};
    dist_cache[latA][lngA][latB][lngB] = dist;
    return dist;
  }
}
function geosailing3(a,b,c,d){
  with(Math)return acos(sin(a*PI/180)*sin(c*PI/180)+cos(a*PI/180)*cos(c*PI/180)*cos(b*PI/180-d*PI/180))*6371.008
}
function geosailing4(a,b,c,d){
  with(Math)return acos(sin(a*(i=PI/180))*sin(c*i)+cos(a*i)*cos(c*i)*cos(b*i-d*i))*6371.008
}
function geosailing5(a,b,c,d){
  with(Math)return acos(sin(a*=(i=PI/180))*sin(c*=i)+cos(a)*cos(c)*cos(b*i-d*i))*6371.008
}
  const R = Math.PI / 180;
  function geosailing6(lat1, lng1, lat2, lng2) {
  lat1 *= R;
  lng1 *= R;
  lat2 *= R;
  lng2 *= R;
  return 6371 * Math.acos(Math.cos(lat1) * Math.cos(lat2) * Math.cos(lng2 - lng1) + Math.sin(lat1) * Math.sin(lat2));
}
   function geosailing7(lat1,lon1,lat2,lon2) {

  var dLat = toRad(lon2 - lon1);
  var dLon = toRad(lat2 - lat1);
  var lat1 = toRad(lon1);
  var lat2 = toRad(lon2);

  var a = Math.pow(Math.sin(dLat/2), 2) +
          Math.pow(Math.sin(dLon/2), 2) * Math.cos(lat1) * Math.cos(lat2);
  var c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a));

  var R = 6373;
 
  var distance = R * c;
  return distance;
};

function toRad(degree) {
  return degree * Math.PI / 180;
}

Run results for: geosailing2

case1

case2

case3

case4

case5

case6

case7