code: def test_geodesic_inv__string_init(scalar_and_array):
      geod = Geod("+ellps=clrk66")
      az12, az21, dist = geod.inv(
          scalar_and_array(_BOSTON_LON),
          scalar_and_array(_BOSTON_LAT),
          scalar_and_array(_PORTLAND_LON),
          scalar_and_array(_PORTLAND_LAT),
      )
      assert_almost_equal(
          (az12, az21, dist),
          (
              scalar_and_array(-66.531),
              scalar_and_array(75.654),
              scalar_and_array(4164192.708),
          ),
          decimal=3,
      )
  
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 12 for calling a global function with example usage: test_geodesic_inv__string_init(scalar_and_array) without return types