class_docstr: Geodesic initial value problem solver.
  
  Parameters
  ----------
  integrator : ODEIVPSolver
      Instance of ODEIVP integrator.
class_docstr_tok: ["Geodesic","initial","value","problem","solver",".","Parameters","--","--","--","--","--","integrator",":","ODEIVPSolver","Instance","of","ODEIVP","integrator","."]
docstr: Geodesic curve for initial value problem.
  
  Parameters
  ----------
  space : Manifold
      Equipped manifold.
  tangent_vec : array-like, shape=[..., *space.shape]
      Tangent vector at the base point.
  base_point : array-like, shape=[..., *space.shape]
      Point on the manifold.
  
  Returns
  -------
  path : callable
      Time parametrized geodesic curve. `f(t)`.
code: def geodesic_ivp(self, space, tangent_vec, base_point):
      """Geodesic curve for initial value problem.
  
      Parameters
      ----------
      space : Manifold
          Equipped manifold.
      tangent_vec : array-like, shape=[..., *space.shape]
          Tangent vector at the base point.
      base_point : array-like, shape=[..., *space.shape]
          Point on the manifold.
  
      Returns
      -------
      path : callable
          Time parametrized geodesic curve. `f(t)`.
      """
      base_point = gs.broadcast_to(base_point, tangent_vec.shape)
  
      def path(t):
          """Time parametrized geodesic curve.
  
          Parameters
          ----------
          t : float or array-like, shape=[n_times,]
  
          Returns
          -------
          geodesic_points : array-like, shape=[..., n_times, *space.shape]
              Geodesic points evaluated at t.
          """
          if not gs.is_array(t):
              t = gs.array([t])
  
          if gs.ndim(t) == 0:
              t = gs.expand_dims(t, axis=0)
  
          result = self._solve(space, tangent_vec, base_point, t_eval=t)
          return self._simplify_result_t(result, space)
  
      return path
  
imports: import math
  from abc import ABC, abstractmethod
  import geomstats.backend
  from geomstats.numerics.bvp import ScipySolveBVP
  from geomstats.numerics.ivp import GSIVPIntegrator
  from geomstats.numerics.optimizers import ScipyMinimize
  from geomstats.vectorization import get_batch_shape
label_desc: Use image node_id 4 for calling the ExpODESolver obj's underlying member method code with example usage: obj.geodesic_ivp(space, tangent_vec, base_point) and returns: path, self