Discussion :

[gbardy]

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]

Je ne voit pas de pb chez moi...

j'obtiens 5.0000000001659828e-005

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

[gbardy]

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]

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

[gbardy]

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]

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

[thewho]

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...).