X-Git-Url: http://dxcluster.net/gitweb/gitweb.cgi?a=blobdiff_plain;f=perl%2FDXBearing.pm;h=db6f4686d0a75ed89c32e0c32e28509ae233bb5a;hb=8106275008346b0cceb090f7317ff4240a3ba44f;hp=b7101a53ce627ad1f535a413784ff98e6021a2a4;hpb=74803d3f4937220f536cd67b515faa15da6e0520;p=spider.git

diff --git a/perl/DXBearing.pm b/perl/DXBearing.pm
index b7101a53..db6f4686 100644
--- a/perl/DXBearing.pm
+++ b/perl/DXBearing.pm
@@ -8,25 +8,19 @@
 #
 # Copyright (c) 1998 - Dirk Koopman G1TLH
 #
-# $Id$
+#
 #
 
 package DXBearing;
 
-use POSIX;
 use DXUtil;
+use POSIX qw(:math_h);
 
 use strict;
 use vars qw($pi);
 
 $pi = 3.14159265358979;
 
-use vars qw($VERSION $BRANCH);
-$VERSION = sprintf( "%d.%03d", q$Revision$ =~ /(\d+)\.(\d+)/ );
-$BRANCH = sprintf( "%d.%03d", q$Revision$ =~ /\d+\.\d+\.(\d+)\.(\d+)/  || (0,0));
-$main::build += $VERSION;
-$main::branch += $BRANCH;
-
 # convert a qra locator into lat/long in DEGREES
 sub qratoll
 {
@@ -95,22 +89,36 @@ sub bdist
 	my $e = dr(shift);
 	return (0, 0) if $hn == $n && $he == $e;
 	my $co = cos($he-$e)*cos($hn)*cos($n)+sin($hn)*sin($n);
-	my $ca = atan(abs(sqrt(1-$co*$co)/$co));
+	my $ca = $co ? atan(abs(sqrt(1-$co*$co)/$co)) : $pi;
 	$ca = $pi-$ca if $co < 0;
 	my $dx = 6367*$ca;
 	my $si = sin($e-$he)*cos($n)*cos($hn);
 	$co = sin($n)-sin($hn)*cos($ca);
-	my $az = atan(abs($si/$co));
+	my $az = $co ? atan(abs($si/$co)) : $pi;
 	$az = $pi - $az if $co < 0;
 	$az = -$az if $si < 0;
 	$az = $az+2*$pi if $az < 0;
 	return (rd($az), $dx);
 }
 
+# just the distance - parameters as above
+sub distance
+{
+	my $hn = dr(shift);
+	my $he = dr(shift);
+	my $n = dr(shift);
+	my $e = dr(shift);
+	return (0, 0) if $hn == $n && $he == $e;
+	my $co = cos($he-$e)*cos($hn)*cos($n)+sin($hn)*sin($n);
+	my $ca = $co ? atan(abs(sqrt(1-$co*$co)/$co)) : $pi;
+	$ca = $pi-$ca if $co < 0;
+	my $dx = 6367*$ca;
+}
+
 # turn a lat long string into floating point lat and long
 sub stoll
 {
-	my ($latd, $latm, $latl, $longd, $longm, $longl) = $_[0] =~ /(\d{1,2})\s+(\d{1,2})\s*([NnSs])\s+(\d{1,3})\s+(\d{1,2})\s*([EeWw])/;
+	my ($latd, $latm, $latl, $longd, $longm, $longl) = $_[0] =~ /(\d{1,2})\s+(\d{1,2})\s*([NnSs])\s+(1?\d{1,2})\s+(\d{1,2})\s*([EeWw])/;
 	
 	$longd += ($longm/60);
 	$longd = 0-$longd if (uc $longl) eq 'W';