code: def test_geodesic_inv_intermediate__del_s__numpy(del_s_fact, flag):
      geod = Geod(ellps="clrk66")
      point_count = 5
      lons = numpy.empty(point_count)
      lats = numpy.empty(point_count)
      azis = numpy.empty(point_count)
      dist = 4164192.7080994663
      del_s = dist / (point_count - 1)
      res = geod.inv_intermediate(
          out_lons=lons,
          out_lats=lats,
          out_azis=azis,
          lon1=_BOSTON_LON,
          lat1=_BOSTON_LAT,
          lon2=_PORTLAND_LON,
          lat2=_PORTLAND_LAT,
          del_s=del_s * del_s_fact,
          initial_idx=0,
          terminus_idx=0,
          flags=flag,
          return_back_azimuth=False,
      )
      assert res.npts == point_count
      assert_almost_equal(res.del_s, del_s)
      assert_almost_equal(res.dist, dist)
      assert_almost_equal(
          res.lons,
          [
              _BOSTON_LON,
              -83.34061499,
              -96.62666098,
              -110.34292364,
              _PORTLAND_LON,
          ],
      )
      assert_almost_equal(
          res.lats,
          [
              _BOSTON_LAT,
              45.35049848,
              47.01597626,
              47.07350417,
              _PORTLAND_LAT,
          ],
      )
      assert_almost_equal(
          res.azis[:-1],
          [-66.53059479, -75.01125433, -84.6195982, -94.68069764],
      )
      assert res.lons is lons
      assert res.lats is lats
      assert res.azis is azis
  
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 6 for calling a global function with example usage: test_geodesic_inv_intermediate__del_s__numpy(del_s_fact, flag) without return types