test float32 geodetic2ecef

Author: scivision

src/pymap3d/ecef.py

@@ -88,8 +88,7 @@ def geodetic2ecef(
     N = ell.semimajor_axis ** 2 / sqrt(
         ell.semimajor_axis ** 2 * cos(lat) ** 2 + ell.semiminor_axis ** 2 * sin(lat) ** 2
     )
-    # Compute cartesian (geocentric) coordinates given  (curvilinear) geodetic
-    # coordinates.
+    # Compute cartesian (geocentric) coordinates given (curvilinear) geodetic coordinates.
     x = (N + alt) * cos(lat) * cos(lon)
     y = (N + alt) * cos(lat) * sin(lon)
     z = (N * (ell.semiminor_axis / ell.semimajor_axis) ** 2 + alt) * sin(lat)

src/pymap3d/tests/test_geodetic.py

@@ -240,3 +240,21 @@ def test_numpy_ecef2geodetic(xyz, lla):
     assert lat[0] == approx(lla[0])
     assert lon[0] == approx(lla[1])
     assert alt[0] == approx(lla[2])
+
+
+@pytest.mark.parametrize("lla, xyz", llaxyz)
+def test_numpy_geodetic2ecef(lla, xyz):
+    np = pytest.importorskip("numpy")
+    x, y, z = pm.geodetic2ecef(
+        *np.array(
+            [
+                [lla],
+            ],
+            dtype=np.float32,
+        ).T
+    )
+
+    atol_dist = 1  # meters
+    assert x[0] == approx(xyz[0], abs=atol_dist)
+    assert y[0] == approx(xyz[1], abs=atol_dist)
+    assert z[0] == approx(xyz[2], abs=atol_dist)