code: def test_array_ecef2geodetic():
      np = pytest.importorskip("numpy")
  
      xyz = (
          np.asarray(xyz0[0]),
          np.asarray(xyz0[1]),
          np.asarray(xyz0[2]),
      )
      lla = pm.ecef2geodetic(*xyz)
      assert np.isclose(pm.geodetic2ecef(*lla), xyz).all()
  
      xyz = (
          np.atleast_1d(xyz0[0]),
          np.atleast_1d(xyz0[1]),
          np.atleast_1d(xyz0[2]),
      )
      lla = pm.ecef2geodetic(*xyz)
      assert np.isclose(pm.geodetic2ecef(*lla), xyz).all()
  
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 4 for calling a global function with example usage: test_array_ecef2geodetic() without return types