code: def powermh(mat, power):
      """Compute the matrix power for a Hermitian matrix.
  
      Parameters
      ----------
      mat : array_like, shape=[..., n, n]
          Symmetric matrix with non-negative eigenvalues.
      power : float, list
          Power at which mat will be raised. If a list of powers is passed,
          a list of results will be returned.
  
      Returns
      -------
      powerm : array_like or list of arrays, shape=[..., n, n]
          Matrix power of mat.
      """
      if isinstance(power, list):
          power_ = [lambda ev, p=p: gs.power(ev, p) for p in power]
      else:
  
          def power_(ev):
              return gs.power(ev, power)
  
      return apply_func_to_eigvalsh(mat, power_, check_positive=False)
  
file_path: geomstats/geomstats/geometry/hermitian_matrices.py
imports: import logging
  import geomstats.backend
  from geomstats import algebra_utils
  from geomstats.geometry.base import ComplexMatrixVectorSpace
  from geomstats.geometry.complex_matrices import ComplexMatrices, ComplexMatricesMetric
  from geomstats.geometry.matrices import Matrices
label_desc: Use image node_id 2 for calling a global function with example usage: powermh(mat, power) and returns: apply_func_to_eigvalsh, gs