code: def test_geod_inv_radians():
      g = Geod(ellps="clrk66")
      lon1 = 0
      lat1 = 0
      lon2 = 1
      lat2 = 1
      # the third output is in distance, so we don't want to change from deg-rad there
      out_rad = list(g.inv(lon1, lat1, lon2, lat2, radians=True))
      out_rad[0] *= 180 / numpy.pi
      out_rad[1] *= 180 / numpy.pi
      assert_almost_equal(
          out_rad,
          g.inv(
              lon1 * 180 / numpy.pi,
              lat1 * 180 / numpy.pi,
              lon2 * 180 / numpy.pi,
              lat2 * 180 / numpy.pi,
          ),
      )
  
imports: import math
  import pickle
  from contextlib import nullcontext
  from itertools import permutations
  import numpy
  import pytest
  from numpy.testing import assert_almost_equal, assert_array_equal
  from pyproj import Geod
  from pyproj.geod import GeodIntermediateFlag, reverse_azimuth
label_desc: Use image node_id 41 for calling a global function with example usage: test_geod_inv_radians() without return types