evc4:division des doubles

**gbardy** · 07/08/2006, 10h41

un petit problème pour les calculs des doubles.
j'utilise evc4, voici mon code:
double lat1=45.8541;
double lat2=45.854;
double deltaLat=lat1-lat2; //cette ligne donne bien 0.0001
deltaLat=deltaLat/2; //et là 5.00000000016
vous savez qu'est-ce qu'il s'est passé ??

**superspag** · 07/08/2006, 11h13

Je ne voit pas de pb chez moi...

j'obtiens 5.0000000001659828e-005

Tu n'aurrais pas oublié de regarder le e-005 ?

**gbardy** · 07/08/2006, 11h42

oui, effectivement

mais problème n'est pas résolu, en fait, j'ai fait une fonction qui calcul la distance entre 2 points latitude/longitude, mais le résultat est totalement faux.

#include <math.h>
double Distance(double Lat1,double Lng1,double Lat2,double Lng2)
{
double rayonTerre=6371;
double deltaLat=Lat1-Lat2;
double deltaLng=Lng1-Lng2; //jusqu'à ici c'est juste

double a=sin(deltaLat/2)*sin(deltaLat/2)+
cos(Lat1)*cos(Lat2)*sin(deltaLng/2)*sin(deltaLng/2);
double c=2*atan2(sqrt(a),sqrt(1-a));
double distance=rayonTerre*c;
return distance; //le résultat est complètement faux, du gerne 0.7 au lieu de 0.02
}

**superspag** · 07/08/2006, 12h09

Quelles sont tes données d'entrées ?

**gbardy** · 07/08/2006, 12h21

Les données entrées:
Lat1=45.854 Lng1=4.7319
Lat2=45.8541 Lng2=4.7321

J'attend comme résultat 0.01907 (en théorie) mais j'ai 0.740652288866354

L'incertitude (740 mètres par rapport à 20 mètres) est quand même inacceptable...

**superspag** · 07/08/2006, 12h26

Je pense que tu devrais convertir tes angles en radians...

**thewho** · 07/08/2006, 12h29

Bonjour,

Je ne prendrai pas la peine de lire ton code.

Une telle erreur vient forcément de l'algorithme, car si c'était le compilateur ou/et le processeur, ça se saurait depuis longtemps, n'est-ce pas ?

Donc revois ton calcul, en faisant quand même attention aux approximations inévitables quand on calcule avec des réels (par exemple, faire la différence entre deux valeurs très proches peut donner un résultat loin de la valeur attendue, idem diviser par une valeur très faible...).

**superspag** · 07/08/2006, 12h33

Je suis d'accord avec thewho

2 choses me choquent dans ton algorithme :

- Tu utilises apparemment des degrés plutôt que des radians
- Tu renvois des km plutôt que des m

Place toi toujours dans un système universel

**gbardy** · 07/08/2006, 13h14

j'utilise le formulaire d'Haversine, les lat/lng sont sous forme de degré décimale.
R = earth’s radius (6371km)
Δlat = lat2− lat1
Δlong = long2− long1
a = sin²(Δlat/2) + cos(lat1).cos(lat2).sin²(Δlong/2)
c = 2.atan2(√a, √(1−a))
d = R.c
donc l'algorithme n'a aucun erreur.

**Matthieu Brucher** · 07/08/2006, 13h36

sauf que les fonction sin, atan2, ... prennent des radians et non des degrés comme la formule le laisse à penser. Fais le test au lieu de dire que ce n'est pas ça.