code: def test_multivox_shore(rng):
      gtab = get_3shell_gtab()
  
      data = rng.random([20, 30, 1, gtab.gradients.shape[0]])
      radial_order = 4
      zeta = 700
      asm = ShoreModel(
          gtab,
          radial_order=radial_order,
          zeta=zeta,
          lambdaN=1e-8,
          lambdaL=1e-8,
      )
      with warnings.catch_warnings():
          warnings.filterwarnings(
              "ignore",
              message=descoteaux07_legacy_msg,
              category=PendingDeprecationWarning,
          )
          asmfit = asm.fit(data)
      c_shore = asmfit.shore_coeff
      npt.assert_equal(c_shore.shape[0:3], data.shape[0:3])
      npt.assert_equal(np.all(np.isreal(c_shore)), True)
  
file_path: dipy/dipy/reconst/tests/test_shore_odf.py
imports: import warnings
  import numpy
  import numpy.testing
  from dipy.data import default_sphere, get_isbi2013_2shell_gtab, get_3shell_gtab
  from dipy.reconst.shore import ShoreModel, shore_matrix
  from dipy.reconst.shm import sh_to_sf, descoteaux07_legacy_msg
  from dipy.direction.peaks import peak_directions
  from dipy.reconst.odf import gfa
  from dipy.sims.voxel import sticks_and_ball
  from dipy.core.subdivide_octahedron import create_unit_sphere
  from dipy.core.sphere_stats import angular_similarity
  from dipy.reconst.tests.test_dsi import sticks_and_ball_dummies
  from dipy.testing.decorators import set_random_number_generator
label_desc: Use image node_id 2 for calling a global function with example usage: test_multivox_shore(rng) without return types