gmp unterstützt meines wissens leider kein sinus, cosinus usw.

Dein Problem wirst Du also mit einem genaueren double auch nicht lösen. Was bedeuten $\theta_A$, $\theta_B$, $\lambda$ und $\zeta$?

Gruß
Werner

Wie kommst Du da auf 2000m?

Gruß
Werner

// p1.lat, p2.lat, p1.lon, p2.lon sind floats 
// (ist aber ncih wichtig, wird ja dann eh in double umgewandelt)

double c = sin(p1.lat)*sin(p2.lat)+cos(p1.lat)*cos(p2.lat)*(cos(p1.lon-p2.lon));
return Radius*acos(c);

///
Radius = 6370000.0f; // so ähnlich, aber genauer in meinem 
// programm, hab ich grad nich zur hand

da kommen bei mir 2000meter raus. ja, und ich rechne im Bogenmaß.

float DegToRad(float d)
{
return d*pi/180.0f; // pi auf 15 Stellen aus dem Windows-taschenrechner
}

Dein Code sieht auf den ersten Blick völlig korrekt aus; ist er aber nicht. Man muss dazu wissen, dass es die trigonometrischen Funktionen sowohl für float als auch für double gibt und Du mischt beide Typen. Du berechnest den Wert für 'c' noch als float aber durch die Zuweisung nach double werden die letzen Binärstellen von 'c' nicht mit 0-en sondern mit dem Rundungsfehler von float belegt. Mit dem Ergebnis gehst Du jetzt in eine Funktion der Signatur acos( double ) hinein. die Funktion arccos hat für x=1 die Steigung unendlich. Numerisch wirkt sich also jede kleinste Änderung von x sehr stark beim Ergebnis aus.

Dein Code würde funktionieren, wenn man durchgängig float oder durchgängig double benutzt. Nur Mischen darst Du nicht.
Also entweder

float c = ..;

oder

double lat, lon;

Gruß
Werner

wenn ich vor den sin, cos... aufrufen die floats in double uwandle funktionierts komischerweise...

float sin(float x);
// und
double sin(double x);

D.h. wenn du float reintust kommt auch float raus

naj,a hat sich ja nun geklärt Rechnung funktioniert auch sehr gut