code: def _value_and_jacobian(fun, point_ndim=1):
      def _value_and_jacobian_vec(x):
          if x.ndim == point_ndim:
              return _value_and_jacobian_op(fun)(x)
  
          ans = []
          jac = []
          for one_x in x:
              ans_, jac_ = _value_and_jacobian_op(fun)(one_x)
              ans.append(ans_)
              jac.append(jac_)
  
          return _np.stack(ans), _np.stack(jac)
  
      return _value_and_jacobian_vec
  
file_path: geomstats/geomstats/_backend/autograd/autodiff.py
imports: import autograd
  import autograd.numpy
  from autograd import jacobian
label_desc: Use image node_id 6 for calling a global function with example usage: _value_and_jacobian(fun, point_ndim) and returns: _value_and_jacobian_vec, _np,  _np, _value_and_jacobian_op