made sure the arcCos si always computed for a value between -1 and 1

pharo-GIS · Feb 28, 2022 · 774a4b3 · 774a4b3
1 parent 22dabc4
commit 774a4b3
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/OGC-Spheric/OGCSphericSystem.class.st b/OGC-Spheric/OGCSphericSystem.class.st
@@ -35,5 +35,5 @@ OGCSphericSystem class >> distanceFromPoint: pointA toPoint: pointB [
 	sinValue := (latA sin ) * (latB sin).
 	cosValue := (latA cos ) * (latB cos) * ((lonA - lonB) cos).
 
-	^ ray * ((sinValue + cosValue) arcCos)
+	^ ray * (((sinValue + cosValue) max: -1.0) min: 1.0) arcCos
 ]