code: def test_circular_mean_scipy(mean):
      """Test our `_circular_mean()` function gives same result than Scipy version."""
      rvs = st.vonmises.rvs(loc=mean, kappa=1, size=1000)
      mean_az = _circular_mean(rvs)
      mean_sp = st.circmean(rvs, low=-np.pi, high=np.pi)
      np.testing.assert_almost_equal(mean_az, mean_sp)
  
imports: from unittest.mock import Mock
  import numpy
  import pytest
  import scipy.stats
  from ...data import dict_to_dataset, from_dict, load_arviz_data
  from ...stats.density_utils import _circular_mean, _normalize_angle, _find_hdi_contours
  from ...utils import _stack, _subset_list, _var_names, expand_dims, flatten_inference_data_to_dict, one_de, two_de
  from ..helpers import TestRandomVariable
label_desc: Use image node_id 22 for calling a global function with example usage: test_circular_mean_scipy(mean) without return types