code: def test_cosp_covariances(rndstate):
      """Test fit CospCovariances"""
      n_matrices, n_channels, n_times = 2, 3, 1000
      x = rndstate.randn(n_matrices, n_channels, n_times)
      cov = CospCovariances()
      cov.fit(x)
      covmats = cov.transform(x)
      assert cov.get_params() == dict(
          window=128, overlap=0.75, fmin=None, fmax=None, fs=None
      )
      n_freqs = 65
      assert covmats.shape == (
          n_matrices,
          n_channels,
          n_channels,
          n_freqs,
      )
      assert is_spsd(covmats.transpose(0, 3, 1, 2))
  
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 14 for calling a global function with example usage: test_cosp_covariances(rndstate) without return types