code: def test_gaussian_blobs():
      """Test function for sampling Gaussian blobs."""
      n_matrices, n_dim = 5, 4
      X, y = make_gaussian_blobs(
          n_matrices=n_matrices,
          n_dim=n_dim,
          class_sep=2.0,
          class_disp=1.0,
          return_centers=False,
          random_state=None,
      )
      assert X.shape == (
          2 * n_matrices,
          n_dim,
          n_dim,
      )  # X shape mismatch
      assert is_spd(X)  # X is an array of SPD matrices
      assert y.shape == (2 * n_matrices,)  # y shape mismatch
      assert np.unique(y).shape == (2,)  # Unexpected number of classes
      assert (
          sum(y == 0) == n_matrices
      )  # Unexpected number of samples in class 0
      assert (
          sum(y == 1) == n_matrices
      )  # Unexpected number of samples in class 1
      _, _, centers = make_gaussian_blobs(
          n_matrices=1,
          n_dim=n_dim,
          class_sep=2.0,
          class_disp=1.0,
          return_centers=True,
          random_state=None,
      )
      assert centers.shape == (2, n_dim, n_dim)
  
file_path: pyRiemann/tests/test_simulated.py
imports: import pytest
  import numpy
  from numpy.testing import assert_array_almost_equal
  from pyriemann.datasets.sampling import generate_random_spd_matrix
  from pyriemann.datasets.simulated import make_covariances, make_matrices, make_masks, make_gaussian_blobs, make_outliers
  from pyriemann.utils.test import is_real, is_sym_pos_def, is_sym_pos_semi_def, is_herm_pos_def, is_herm_pos_semi_def
label_desc: Use image node_id 5 for calling a global function with example usage: test_gaussian_blobs() without return types