code: def test_geodesic_inv__multiple_points(
      lon_start,
      lat_start,
      lon_end,
      lat_end,
      res12,
      res21,
      resdist,
      scalar_and_array,
  ):
      geod = Geod(ellps="clrk66")
      o_az12, o_az21, o_dist = geod.inv(
          scalar_and_array(lon_start),
          scalar_and_array(lat_start),
          scalar_and_array(lon_end),
          scalar_and_array(lat_end),
      )
      assert_almost_equal(
          (o_az12, o_az21, o_dist),
          (
              scalar_and_array(res12),
              scalar_and_array(res21),
              scalar_and_array(resdist),
          ),
      )
  
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 2 for calling a global function with example usage: test_geodesic_inv__multiple_points(lon_start, lat_start, lon_end, lat_end, res12, res21, resdist, scalar_and_array) without return types