code: def hilbert_schmidt_inner_product(m1: np.ndarray, m2: np.ndarray) -> complex:
      
      return np.einsum('ij,ij', m1.conj(), m2)
docs: Computes Hilbert-Schmidt inner product of two matrices.
  
      Linear in second argument.