docs: def neighbors(self, distance=2.0):
          """Get all neighbors for all neurons."""
          dgrid = self.distance_grid.reshape(self.num_neurons, self.num_neurons)
          for x, y in zip(*np.nonzero(dgrid <= distance)):
              if x != y:
                  yield x, y