code: def test_hankel_covariances_delays(delays, rndstate):
      n_matrices, n_channels, n_times = 2, 3, 100
      x = rndstate.randn(n_matrices, n_channels, n_times)
      cov = HankelCovariances(delays=delays).fit(x)
      covmats = cov.fit_transform(x)
      assert cov.get_params() == dict(estimator="scm", delays=delays)
      if isinstance(delays, list):
          n_delays = 1 + len(delays)
      elif isinstance(delays, int):
          n_delays = delays
      assert covmats.shape == (
          n_matrices,
          n_delays * n_channels,
          n_delays * n_channels,
      )
      assert is_spd(covmats)
  
imports: import numpy
  from numpy.testing import assert_array_almost_equal, assert_array_equal, assert_raises
  import pytest
  from pyriemann.estimation import Covariances, ERPCovariances, XdawnCovariances, CospCovariances, HankelCovariances, Coherences, Shrinkage, BlockCovariances, Kernels
  from pyriemann.utils.test import is_sym_pos_def, is_sym_pos_semi_def
label_desc: Use image node_id 3 for calling a global function with example usage: test_hankel_covariances_delays(delays, rndstate) without return types