code: def expmh(mat):
      """Compute the matrix exponential for a Hermitian matrix.
  
      Parameters
      ----------
      mat : array_like, shape=[..., n, n]
          Symmetric matrix.
  
      Returns
      -------
      exponential : array_like, shape=[..., n, n]
          Exponential of mat.
      """
      n = mat.shape[-1]
      dim_3_mat = gs.reshape(mat, [-1, n, n])
      expm = apply_func_to_eigvalsh(dim_3_mat, gs.exp)
      return gs.reshape(expm, mat.shape)
  
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 1 for calling a global function with example usage: expmh(mat) and returns: gs