use numpy.broadcast_arrays instead of broadcast_to for simplicity

Author: scivision

setup.cfg

@@ -1,6 +1,6 @@
 [metadata]
 name = pymap3d
-version = 2.7.1
+version = 2.7.2
 author = Michael Hirsch, Ph.D.
 author_email = [email protected]
 description = pure Python (no prereqs) coordinate conversions, following convention of several popular Matlab routines.

src/pymap3d/lox.py

@@ -4,18 +4,7 @@
 import typing
 
 try:
-    from numpy import (
-        radians,
-        degrees,
-        cos,
-        arctan2 as atan2,
-        tan,
-        pi,
-        array,
-        atleast_1d,
-        broadcast_arrays,
-        broadcast_to,
-    )
+    from numpy import radians, degrees, cos, arctan2 as atan2, tan, pi, array, broadcast_arrays
 except ImportError:
     from math import radians, degrees, cos, atan2, tan, pi  # type: ignore
 
@@ -227,16 +216,7 @@ def loxodrome_direct(
         lat1, lon1, a12 = radians(lat1), radians(lon1), radians(a12)
 
     try:
-        lat1 = atleast_1d(lat1)
-        rng = atleast_1d(rng)
-
-        if lat1.shape != rng.shape:
-            if rng.size == 1:
-                rng = broadcast_to(rng, lat1.shape)
-            elif lat1.size == 1:
-                lat1 = broadcast_to(lat1, rng.shape)
-            else:
-                raise ValueError("lat1 and rng must each be scalars or same shape")
+        lat1, rng = broadcast_arrays(lat1, rng)
         if (abs(lat1) > pi / 2).any():  # type: ignore
             raise ValueError("-90 <= latitude <= 90")
         if (rng < 0).any():  # type: ignore

src/pymap3d/tests/test_rhumb.py

@@ -78,9 +78,7 @@ def test_numpy_loxodrome_inverse():
 
 def test_numpy_2d_loxodrome_inverse():
     pytest.importorskip("numpy")
-    d, a = lox.loxodrome_inverse(
-        [[40, 40], [40, 40]], [[-80, -80], [-80, -80]], 65, -148
-    )
+    d, a = lox.loxodrome_inverse([[40, 40], [40, 40]], [[-80, -80], [-80, -80]], 65, -148)
     assert d == approx(5248666.209)
     assert a == approx(302.00567)
 
@@ -90,9 +88,7 @@ def test_numpy_2d_loxodrome_inverse():
     d, a = lox.loxodrome_inverse(
         [[40, 40], [40, 40]], [[-80, -80], [-80, -80]], 65, [[-148, -148], [-148, -148]]
     )
-    d, a = lox.loxodrome_inverse(
-        40, -80, [[65, 65], [65, 65]], [[-148, -148], [-148, -148]]
-    )
+    d, a = lox.loxodrome_inverse(40, -80, [[65, 65], [65, 65]], [[-148, -148], [-148, -148]])
 
 
 @pytest.mark.parametrize(