code: def _objective(self, velocity, space, point, base_point):
      velocity = gs.reshape(velocity, base_point.shape)
      delta = space.metric.exp(velocity, base_point) - point
      return gs.sum(delta**2)
  
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 3 for calling the _LogShootingSolverFlatten obj's underlying member method code with example usage: obj._objective(velocity, space, point, base_point) and returns: gs