Vous utilisez la formule de la loi sphérique des cosinus , pas la Formule Haversine (ce qui est légèrement plus lent).
Faire le calcul dans MySQL sera probablement beaucoup plus rapide, je vous conseille de lire cette question J'ai demandé il y a quelque temps, comme pour les optimisations de vitesse, vous devriez certainement lire l'excellent Geo Proximity Rechercher avec MySQL papier, portez une attention particulière aux pages 8-14 et 19 .