docs: def scale_8(b):
      """Translations necessary for class-8 coefficients
      to go from a basis with only non-redundant WCxf
      operators to a basis where the Wilson coefficients are symmetrized like
      the operators."""
      a = np.array(b, copy=True, dtype=complex)
      for i in range(3):
          a[0, 0, 1, i] = 1/2 * b[0, 0, 1, i]
          a[0, 0, 2, i] = 1/2 * b[0, 0, 2, i]
          a[0, 1, 1, i] = 1/2 * b[0, 1, 1, i]
          a[0, 1, 2, i] = 2/3 * b[0, 1, 2, i] - 1/6 * b[0, 2, 1, i] - 1/6 * b[1, 0, 2, i] + 1/6 * b[1, 2, 0, i]
          a[0, 2, 1, i] = - (1/6) * b[0, 1, 2, i] + 2/3 * b[0, 2, 1, i] + 1/6 * b[1, 0, 2, i] + 1/3 * b[1, 2, 0, i]
          a[0, 2, 2, i] = 1/2 * b[0, 2, 2, i]
          a[1, 0, 2, i] = - (1/6) * b[0, 1, 2, i] + 1/6 * b[0, 2, 1, i] + 2/3 * b[1, 0, 2, i] - 1/6 * b[1, 2, 0, i]
          a[1, 1, 2, i] = 1/2 * b[1, 1, 2, i]
          a[1, 2, 0, i] = 1/6 * b[0, 1, 2, i] + 1/3 * b[0, 2, 1, i] - 1/6 * b[1, 0, 2, i] + 2/3 * b[1, 2, 0, i]
          a[1, 2, 2, i] = 1/2 * b[1, 2, 2, i]
      return a