code: def _compare_bem_surfaces(surfs_1, surfs_2):
      """Compare BEM surfaces."""
      names = [
          "id",
          "nn",
          "rr",
          "coord_frame",
          "tris",
          "sigma",
          "ntri",
          "np",
      ]
      ignores = ["tri_cent", "tri_nn", "tri_area", "neighbor_tri"]
      for s0, s1 in zip(surfs_1, surfs_2):
          assert_equal(set(names), set(s0.keys()) - set(ignores))
          assert_equal(set(names), set(s1.keys()) - set(ignores))
          for name in names:
              assert_allclose(
                  s0[name],
                  s1[name],
                  rtol=1e-3,
                  atol=1e-6,
                  err_msg='Mismatch: "%s"' % name,
              )
  
imports: import re
  from copy import deepcopy
  from os import makedirs
  from pathlib import Path
  from shutil import copy
  import numpy
  import pytest
  from numpy.testing import assert_allclose, assert_equal
  import mne
  from mne import Info, Transform, make_bem_model, make_bem_solution, make_sphere_model, read_bem_solution, read_bem_surfaces, write_bem_solution, write_bem_surfaces, write_head_bem, write_surface
  from mne._fiff.constants import FIFF
  from mne.bem import _assert_complete_surface, _assert_inside, _bem_find_surface, _check_surface_size, _get_ico_map, _ico_downsample, _order_surfaces, distance_to_bem, fit_sphere_to_headshape, make_scalp_surfaces
  from mne.datasets import testing
  from mne.io import read_info
  from mne.surface import _get_ico_surface, read_surface
  from mne.transforms import translation
  from mne.utils import catch_logging, check_version
label_desc: Use image node_id 1 for calling a global function with example usage: _compare_bem_surfaces(surfs_1, surfs_2) without return types