code: def test_mean_alm(kind, get_mats):
      """Test the ALM mean"""
      n_matrices, n_channels = 3, 3
      mats = get_mats(n_matrices, n_channels, kind)
      C_alm = mean_alm(mats)
      assert C_alm.shape == (n_channels, n_channels)
      C_riem = mean_riemann(mats)
      assert C_alm == approx(C_riem, abs=1e-6, rel=1e-3)
  
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 7 for calling a global function with example usage: test_mean_alm(kind, get_mats) without return types