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","."]
code: def _simplify_exp_result(self, result, space):
      y = result.get_last_y()
  
      if self.integrator.state_is_raveled:
          dim_vec = math.prod(space.shape)
          exp = y[..., :dim_vec]
          if space.point_ndim > 1:
              return gs.reshape(exp, result.batch_shape + space.shape)
  
          return exp
  
      return y[0]
  
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 8 for calling the ExpODESolver obj's underlying member method code with example usage: obj._simplify_exp_result(result, space) and returns: y, exp, gs