code: def test_inside_ecef2geodetic():
      np = pytest.importorskip("numpy")
      # test values with no points inside ellipsoid
      lla0_array = (
          np.array([lla0[0], lla0[0]]),
          np.array([lla0[1], lla0[1]]),
          np.array([lla0[2], lla0[2]]),
      )
      xyz = pm.geodetic2ecef(*lla0_array)
      lats, lons, alts = pm.ecef2geodetic(*xyz)
  
      assert lats == approx(lla0_array[0])
      assert lons == approx(lla0_array[1])
      assert alts == approx(lla0_array[2])
  
      # test values with some (but not all) points inside ellipsoid
      lla0_array_inside = (
          np.array([lla0[0], lla0[0]]),
          np.array([lla0[1], lla0[1]]),
          np.array([lla0[2], -lla0[2]]),
      )
      xyz = pm.geodetic2ecef(*lla0_array_inside)
      lats, lons, alts = pm.ecef2geodetic(*xyz)
  
      assert lats == approx(lla0_array_inside[0])
      assert lons == approx(lla0_array_inside[1])
      assert alts == approx(lla0_array_inside[2])
  
file_path: pymap3d/src/pymap3d/tests/test_geodetic.py
imports: from math import isnan, nan, radians, sqrt
  import pymap3d
  import pytest
  from pytest import approx
label_desc: Use image node_id 5 for calling a global function with example usage: test_inside_ecef2geodetic() without return types