code: def _generate_icosahedron():
      x = np.array(
          [[0, -1, -phi], [0, -1, +phi], [0, +1, -phi], [0, +1, +phi]]
      )
      return np.concatenate([np.roll(x, i, axis=1) for i in range(3)])
  
imports: import numpy
  import itertools
  from numpy.testing import assert_equal, assert_almost_equal, assert_array_equal, assert_array_almost_equal
  import pytest
  from pytest import raises
  from scipy.spatial import SphericalVoronoi, distance
  from scipy.optimize import linear_sum_assignment
  from scipy.constants import golden
  from scipy.special import gamma
label_desc: Use image node_id 5 for calling a global function with example usage: _generate_icosahedron() and returns: np