code: def test_geodetic_alt_geocentric(geodetic_lat, alt_m, geocentric_lat):
      assert latitude.geod2geoc(geodetic_lat, alt_m) == approx(
          geocentric_lat
      )
  
      r = rcurve.geocentric_radius(geodetic_lat)
      assert latitude.geoc2geod(geocentric_lat, r) == approx(
          geodetic_lat
      )
      assert latitude.geoc2geod(geocentric_lat, 1e5 + r) == approx(
          latitude.geocentric2geodetic(geocentric_lat, 1e5 + alt_m)
      )
  
      assert latitude.geod2geoc(geodetic_lat, 1e5 + alt_m) == approx(
          latitude.geodetic2geocentric(geodetic_lat, 1e5 + alt_m)
      )
  
file_path: pymap3d/src/pymap3d/tests/test_latitude.py
imports: from math import inf, radians
  import pymap3d.latitude
  import pymap3d.rcurve
  import pytest
  from pytest import approx
label_desc: Use image node_id 1 for calling a global function with example usage: test_geodetic_alt_geocentric(geodetic_lat, alt_m, geocentric_lat) without return types