Distanzberechnung zwischen zwei Koordinaten gibt manchmal eine Fehlermeldung aus

[fox67]

Hallo
ich bin gerade Dabei ein Programm das eine Datei aus einem Fluglogger(aus einem Segelflieger) ausließt und dann eine Strecke raus optiemieren soll(wie z.B hier). Im Internet findet man nicht viel wie so ein Algorithmus funktionieren könnte, aber ich hab ein Quellcode gefunden auf github. Allerdings in C jetzt bin ich gerade dabei das Schritt für Schritt in Delphi umzusetzen. Aber beim berechnen der Distanzmatrix kommt bei einem bestimmten Punkt die Fehlermeldung daß es eine ungültige Gleitkommaoperation sei, obwohl es davor schon einige 1000 distanzen fehlerfrei ausgerechnet hat. Um den Fehler nachzuvollziehen gebe ich den Index der Zwei Punkte aus zwischen denen die Distanz berechnet wird in einer Memo aus. Der Fehler tritt zum erstenmal auf wenn er den 12 mit dem 51 Punkt vergleicht. Leider kann ich absolut nicht nachvollziehen warum der Fehler auftritt schließlich hat er ja davor schon einige tausen distanzen fehlerfrei berechnet. Im Anhang habe ich das komplette Projekt mit quellcode und der Logger datei(.IGC).
Das ist die Funktion in der der Fehler auftritt.

Delphi-Quellcode:

function TFlug.dpointtopoint(index1: Integer; index2: Integer) : integer;
var
dist,x : double;
sli,cli,lri : double;

begin
  sli := Punkte[index1].sinlat;
  cli := Punkte[index1].coslat;
  lri := Punkte[index1].lonrad;
  x := sli*Punkte[index2].sinlat;
  x := x + cli*Punkte[index2].coslat * cos(lri - Punkte[index2].lonrad);
  form10.Memo1.Text := 'i: '+ inttostr(index1) +'|j: '+inttostr(index2);
  Application.ProcessMessages;
  if (index1=12) and (index2=51) then
  begin
    dist := 0;

  end
  else dist := FAK*arccos(x);

  result := round(dist);
end;

[Sir Rufo]

Und jetzt?

Wir sollen uns also dein Projekt herunterladen und ausführen und selber herausfinden, welcher Fehler auftritt und auch selber die Werte der Koordinaten herausfinden um dir zu helfen?

Ist bei dir heute Ostern und Weihnachten zugleich?
Oder klemmt die Zwischenablage?

Also:

Fehlermeldung per CopyPaste hier posten und die Werte der beiden Koordinaten, bei denen dieser Fehler kommt, dann können wir uns den Download komplett schenken.

[fox67]

Nein natürlich nicht, aber ich kann den Fehler absolut nicht nachvollziehen und hab gehofft des jemand hier im Forum villeicht so ein ähnlich es problem schon mal hatte? Der Fehler ensteht erst wenn der arcuscosinus berechnet wird. Was ich nicht verstehe davor hat das programm zig mal für den selben wert den arcuscosinus errechnet ohne probleme und da kommt aufeinnmal der Fehler. Die Distanz sollte 0 sein da x = 1 ist. Die Koordinaten sind: 51°18.014'N 014°07.87'E und der Punkt2 hat die selben Koordinaten, deshalb auch der Abstand 0.

[himitsu]

Zitat von arccos:

Die Werte von X müssen zwischen -1 und 1 liegen.

Und X ist da 1.

Aber wer benutzt schon den Debugger?
Haltepunkt auf den Anfang der Funktion und als Bedingung (index1 = 12) and (index2 = 55) angeben. (das sind übrigens die Werte, wo es immer knallt)

[fox67]

Die Fehler Nachricht ist im Anhang

[fox67]

Den Debugger habe ich benutzt, das x den wert 1 hat ist ganz normal das ist ca. bei den ersten 10 logs so und jedesmal wenn irgentein ereigniswar z.B: wenn ein neuer Flieger vom Kollisionswarngerät erfasst wurde. Da das Programm die ersten 10 werte geschluckt hat und beiden den war der wert x öfters 1 und auch bis zur fehlermeldung ist es bestimmt mehr als hundert mal aufgetreten das x= 1 ist daran kann der fehler eigentlich nicht liegen.

[fox67]

So ganz schlüßig ist mir das Problem nicht anscheinend manchmal passiert nichts wenn x = 1 ist und manchmal kommt die fehlermeldung wenn x= 1 ist, dann hilft aber auch eine if(if x=1 then) abfrage nicht den er sieht die bedingung x=1 als false an(warum auch immer) Was alerdings hilft ist das x mithilfe von trunc in eine Integervariable umzuwandeln und dann zu überprüfen ob x = 1 ist dann klappt es einwandfrei. Gibt es eine erklärung warum 1 als double nicht immer gleich 1 ist?