code: def test_geod__reverse_azimuth(radians):
      f = math.pi / 180 if radians else 1
      xy = numpy.array(
          [
              [0, 0 + 180],
              [180, 180 - 180],
              [-180, -180 + 180],
              [10, 10 - 180],
              [20, 20 - 180],
              [-10, -10 + 180],
          ]
      )
      for x, y in xy:
          assert_almost_equal(
              reverse_azimuth(x * f, radians=radians), y * f
          )
  
      xx = xy.T[0]
      yy = xy.T[1]
      assert_almost_equal(
          reverse_azimuth(xx * f, radians=radians), yy * f
      )
  
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 46 for calling a global function with example usage: test_geod__reverse_azimuth(radians) without return types