code: def test_wrap_ufunc_out_shape_multi_output_same():
      func = lambda x: (np.random.rand(1, 2), np.random.rand(1, 2))
      ary = np.ones((4, 100))
      res1, res2 = wrap_xarray_ufunc(
          func,
          ary,
          func_kwargs={"out_shape": ((1, 2), (1, 2))},
          ufunc_kwargs={"n_dims": 1, "n_output": 2},
      )
      assert res1.shape == (*ary.shape[:-1], 1, 2)
      assert res2.shape == (*ary.shape[:-1], 1, 2)
  
imports: import numpy
  import pytest
  from numpy.testing import assert_array_almost_equal
  from scipy.special import logsumexp
  from scipy.stats import circstd
  from ...data import from_dict, load_arviz_data
  from ...stats.density_utils import histogram
  from ...stats.stats_utils import ELPDData, _angle, _circfunc, _circular_standard_deviation, _sqrt, get_log_likelihood
  from ...stats.stats_utils import logsumexp
  from ...stats.stats_utils import make_ufunc, not_valid, stats_variance_2d, wrap_xarray_ufunc
label_desc: Use image node_id 6 for calling a global function with example usage: test_wrap_ufunc_out_shape_multi_output_same() without return types