code: def pinvh(a, rtol=None):
      """Compute a pseudo-inverse of a Hermitian matrix."""
      s, u = np.linalg.eigh(a)
      del a
      if rtol is None:
          rtol = s.size * np.finfo(s.dtype).eps
      maxS = np.max(np.abs(s))
      above_cutoff = abs(s) > maxS * rtol
      psigma_diag = 1.0 / s[above_cutoff]
      u = u[:, above_cutoff]
      return (u * psigma_diag) @ u.conj().T
  
imports: import inspect
  import operator
  import os
  import warnings
  from contextlib import contextmanager
  from io import StringIO
  from math import log
  from pprint import pprint
  import numpy
label_desc: Use image node_id 20 for calling a global function with example usage: pinvh(a, rtol) and returns: unknown