code: def test_mean_euclid(rndstate, complex_valued):
      """Test the Euclidean mean for generic matrices"""
      n_matrices, n_dim0, n_dim1 = 10, 3, 4
      mats = rndstate.randn(n_matrices, n_dim0, n_dim1)
      if complex_valued:
          mats = mats + 1j * rndstate.randn(n_matrices, n_dim0, n_dim1)
      assert mean_euclid(mats) == approx(mats.mean(axis=0))
  
imports: import numpy
  from scipy.stats import gmean
  import pytest
  from pytest import approx
  from pyriemann.utils.geodesic import geodesic_riemann
  from pyriemann.utils.mean import mean_covariance, mean_ale, mean_alm, mean_euclid, mean_harmonic, mean_identity, mean_kullback_sym, mean_logdet, mean_logeuclid, mean_power, mean_riemann, mean_wasserstein, maskedmean_riemann, nanmean_riemann
label_desc: Use image node_id 9 for calling a global function with example usage: test_mean_euclid(rndstate, complex_valued) without return types