code: def test_compare_geodetic():
      pyproj = pytest.importorskip("pyproj")
  
      xyz = pm.geodetic2ecef(*lla0)
  
      ecef = pyproj.Proj(proj="geocent", ellps="WGS84", datum="WGS84")
      lla = pyproj.Proj(proj="latlong", ellps="WGS84", datum="WGS84")
  
      assert pyproj.transform(
          lla, ecef, lla0[1], lla0[0], lla0[2]
      ) == approx(xyz)
      assert pyproj.transform(ecef, lla, *xyz) == approx(
          (lla0[1], lla0[0], lla0[2])
      )
  
imports: import pymap3d
  import pytest
  from pymap3d.vincenty import vreckon
  from pytest import approx
label_desc: Use image node_id 2 for calling a global function with example usage: test_compare_geodetic() without return types