code: def test_geodesic_euclid(rndstate, complex_valued):
      """Test Euclidean geodesic for generic matrices"""
      n_matrices, n_dim0, n_dim1 = 2, 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)
      A, B = mats[0], mats[1]
      assert geodesic_euclid(A, B).shape == A.shape
  
file_path: pyRiemann/tests/test_utils_geodesic.py
imports: import numpy
  import pytest
  from pytest import approx
  from pyriemann.utils.geodesic import geodesic, geodesic_euclid, geodesic_logeuclid, geodesic_riemann
  from pyriemann.utils.mean import mean_euclid, mean_logeuclid, mean_riemann
label_desc: Use image node_id 8 for calling a global function with example usage: test_geodesic_euclid(rndstate, complex_valued) without return types