code: def taylor_exp_even_func(
      point, taylor_function, order=5, tol=EPSILON
  ):
      """Taylor Approximation of an even function around zero.
  
      Parameters
      ----------
      point : array-like
          Argument of the function to approximate.
      taylor_function : dict with following keys
          function : callable
              Even function to approximate around zero.
          coefficients : list
              Taylor coefficients of even order at zero.
      order : int
          Order of the Taylor approximation.
          Optional, Default: 5.
      tol : float
          Threshold to use the approximation instead of the function's value.
          Where `abs(point) <= tol`, the approximation is returned.
  
      Returns
      -------
      function_value: array-like
          Value of the function at point.
      """
      approx = gs.einsum(
          "k,k...->...",
          gs.array(taylor_function["coefficients"][:order]),
          gs.array([point**k for k in range(order)]),
      )
      point_ = gs.where(gs.abs(point) <= tol, tol, point)
      exact = taylor_function["function"](gs.sqrt(point_))
      result = gs.where(gs.abs(point) < tol, approx, exact)
      return result
  
imports: import math
  import geomstats.backend
label_desc: Use image node_id 2 for calling a global function with example usage: taylor_exp_even_func(point, taylor_function, order, tol) and returns: result