docstr: Check geodesic belongs to manifold.
  
  This is for geodesics defined by the initial value problem (ivp).
  
  Parameters
  ----------
  n_points : int
      Number of random points to generate.
  atol : float
      Absolute tolerance.
docstr_tok: ["Check","geodesic","belongs","to","manifold",".","This","is","for","geodesics","defined","by","the","initial","value","problem","(","ivp",")",".","Parameters","--","--","--","--","--","n_points",":","int","Number","of","random","points","to","generate",".","atol",":","float","Absolute","tolerance","."]
code: def test_geodesic_ivp_belongs(self, n_points, n_times, atol):
      """Check geodesic belongs to manifold.
  
      This is for geodesics defined by the initial value problem (ivp).
  
      Parameters
      ----------
      n_points : int
          Number of random points to generate.
      atol : float
          Absolute tolerance.
      """
      initial_point = self.data_generator.random_point(n_points)
      initial_tangent_vec = self.data_generator.random_tangent_vec(
          initial_point
      )
  
      time = get_random_times(n_times)
  
      geod_func = self.space.metric.geodesic(
          initial_point, initial_tangent_vec=initial_tangent_vec
      )
  
      points = geod_func(time)
  
      res = self.space.belongs(
          gs.reshape(points, (-1, *self.space.shape)), atol=atol
      )
  
      expected_shape = (
          math.prod(
              get_batch_shape(self.space.point_ndim, initial_point)
          )
          * n_times,
      )
  
      expected = gs.ones(expected_shape, dtype=bool)
      self.assertAllEqual(res, expected)
  
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 21 for calling the ConnectionTestCase obj's underlying member method code with example usage: obj.test_geodesic_ivp_belongs(n_points, n_times, atol) without return types