code: def test_sample_gaussian_spd_dimsup(n_dim, n_jobs, sampling_method):
      """Test Riemannian Gaussian sampling for dim>2."""
      n_matrices, sigma = 5, 1.0
      mean = np.eye(n_dim)
      X = sample_gaussian_spd(
          n_matrices,
          mean,
          sigma,
          random_state=42,
          n_jobs=n_jobs,
          sampling_method=sampling_method,
      )
      assert X.shape == (n_matrices, n_dim, n_dim)  # X shape mismatch
      assert is_spd(X)
  
file_path: pyRiemann/tests/test_sampling.py
imports: import pytest
  import numpy
  from pyriemann.datasets.sampling import sample_gaussian_spd, generate_random_spd_matrix
  from pyriemann.utils.distance import distance_riemann
  from pyriemann.utils.test import is_sym_pos_def
label_desc: Use image node_id 2 for calling a global function with example usage: test_sample_gaussian_spd_dimsup(n_dim, n_jobs, sampling_method) without return types