code: def test_geodesic_bvp_vec(self, n_reps, n_times, atol):
      initial_point, end_point = self.data_generator.random_point(2)
      time = get_random_times(n_times)
  
      expected = self.space.metric.geodesic(
          initial_point, end_point=end_point
      )(time)
  
      vec_data = generate_vectorization_data(
          data=[
              dict(
                  initial_point=initial_point,
                  end_point=end_point,
                  time=time,
                  expected=expected,
                  atol=atol,
              )
          ],
          arg_names=["initial_point", "end_point"],
          expected_name="expected",
          n_reps=n_reps,
      )
      self._test_vectorization(vec_data, test_fnc_name="test_geodesic")
  
imports: import math
  import pytest
  import geomstats.backend
  from geomstats.test.random import RandomDataGenerator, get_random_times
  from geomstats.test.test_case import TestCase
  from geomstats.test.utils import IdentityPointTransformer
  from geomstats.test.vectorization import generate_vectorization_data
  from geomstats.vectorization import get_batch_shape
label_desc: Use image node_id 16 for calling the ConnectionTestCase obj's underlying member method code with example usage: obj.test_geodesic_bvp_vec(n_reps, n_times, atol) without return types