code: def test_fast_cross_3d():
      """Test cross product with lots of elements."""
      x = rng.rand(100000, 3)
      y = rng.rand(1, 3)
      z = np.cross(x, y)
      zz = fast_cross_3d(x, y)
      assert_array_equal(z, zz)
      # broadcasting and non-2D
      zz = fast_cross_3d(x[:, np.newaxis], y[0])
      assert_array_equal(z, zz[:, 0])
  
imports: from pathlib import Path
  import numpy
  import pytest
  from numpy.testing import assert_allclose, assert_array_equal, assert_equal
  from mne import decimate_surface, dig_mri_distances, get_montage_volume_labels, pick_types, read_surface, write_surface
  from mne._fiff.constants import FIFF
  from mne.channels import make_dig_montage
  from mne.datasets import testing
  from mne.io import read_info
  from mne.surface import _compute_nearest, _get_ico_surface, _marching_cubes, _normal_orth, _project_onto_surface, _read_patch, _tessellate_sphere, _voxel_neighbors, fast_cross_3d, get_head_surf, get_meg_helmet_surf, read_curvature
  from mne.transforms import _get_trans
  from mne.utils import _record_warnings, catch_logging, object_diff, requires_freesurfer
label_desc: Use image node_id 3 for calling a global function with example usage: test_fast_cross_3d() without return types