code: def test_mvoxel_gqi():
      data, gtab = dsi_voxels()
      sphere = get_sphere("symmetric724")
  
      gq = GeneralizedQSamplingModel(gtab, "standard")
      gqfit = gq.fit(data)
      all_odfs = gqfit.odf(sphere)
  
      # Check that the first and last voxels each have 2 peaks
      odf = all_odfs[0, 0, 0]
      directions, values, indices = peak_directions(
          odf, sphere, 0.35, 25
      )
      assert_equal(directions.shape[0], 2)
      odf = all_odfs[-1, -1, -1]
      directions, values, indices = peak_directions(
          odf, sphere, 0.35, 25
      )
      assert_equal(directions.shape[0], 2)
  
file_path: dipy/dipy/reconst/tests/test_gqi.py
imports: import numpy
  from dipy.data import get_fnames, dsi_voxels, get_sphere, default_sphere
  from dipy.core.gradients import gradient_table
  from dipy.core.subdivide_octahedron import create_unit_sphere
  from dipy.core.sphere_stats import angular_similarity
  from dipy.reconst.gqi import GeneralizedQSamplingModel
  from dipy.reconst.tests.test_dsi import sticks_and_ball_dummies
  from dipy.sims.voxel import sticks_and_ball
  from dipy.reconst.odf import gfa
  from dipy.direction.peaks import peak_directions
  from numpy.testing import assert_equal, assert_almost_equal
label_desc: Use image node_id 2 for calling a global function with example usage: test_mvoxel_gqi() without return types