code: def _generate_dodecahedron():
      x1 = _generate_cube()
      x2 = np.array(
          [
              [0, -phi, -1 / phi],
              [0, -phi, +1 / phi],
              [0, +phi, -1 / phi],
              [0, +phi, +1 / phi],
          ]
      )
      x3 = np.array(
          [
              [-1 / phi, 0, -phi],
              [+1 / phi, 0, -phi],
              [-1 / phi, 0, +phi],
              [+1 / phi, 0, +phi],
          ]
      )
      x4 = np.array(
          [
              [-phi, -1 / phi, 0],
              [-phi, +1 / phi, 0],
              [+phi, -1 / phi, 0],
              [+phi, +1 / phi, 0],
          ]
      )
      return np.concatenate((x1, x2, x3, x4))
  
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 4 for calling a global function with example usage: _generate_dodecahedron() and returns: np