code: def test_geodesic_bvp(self, initial_point, end_point, time, atol):
      initial_point_ = self.point_transformer.transform_point(
          initial_point
      )
      end_point_ = self.point_transformer.transform_point(end_point)
  
      res = self.space.metric.geodesic(
          initial_point, end_point=end_point
      )(time)
      res_ = self.other_space.metric.geodesic(
          initial_point_, end_point=end_point_
      )(time)
  
      res_ = self.point_transformer.inverse_transform_point(res_)
      self.assertAllClose(res, res_, atol=atol)
  
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 20 for calling the ConnectionComparisonTestCase obj's underlying member method code with example usage: obj.test_geodesic_bvp(initial_point, end_point, time, atol) without return types