code: def test_l1_basis_orientations():
      """Test that angles between the basis components matches orientations."""
      binname = (
          fil_path
          / "sub-noise_ses-001_task-noise220622_run-001_meg.bin"
      )
      raw = read_raw_fil(binname)
      raw.info["bads"].extend([b for b in bads])
      projs = compute_proj_hfc(raw.info, accuracy="point")
      basis = np.hstack([p["data"]["data"].T for p in projs])
      picks = pick_types(raw.info, meg="mag")
      assert len(picks) == 68
      assert basis.shape == (len(picks), 3)
      ang_model = _angle_between_each(basis)
      n_ang = len(picks) ** 2
      assert ang_model.shape == (n_ang,)
  
      chs = pick_info(raw.info, picks)["chs"]
      ori_sens = np.array([ch["loc"][-3:] for ch in chs])
      # match the normalization that our projectors get
      ori_sens /= np.linalg.norm(ori_sens, axis=0, keepdims=True)
      assert ori_sens.shape == (len(picks), 3)
      ang_sens = _angle_between_each(ori_sens)
      assert ang_sens.shape == (n_ang,)
  
      assert_allclose(ang_sens, ang_model, atol=1e-7)
  
file_path: mne-python/mne/preprocessing/tests/test_hfc.py
imports: from pathlib import Path
  import numpy
  import pytest
  from numpy.testing import assert_allclose
  from scipy.io import loadmat
  from mne._fiff.pick import pick_channels, pick_info, pick_types
  from mne.datasets import testing
  from mne.io import read_info, read_raw_fil
  from mne.preprocessing.hfc import compute_proj_hfc
label_desc: Use image node_id 4 for calling a global function with example usage: test_l1_basis_orientations() without return types