Is there a simple way to create a rectangular grid, like this loopy code: (This is for a QAM256 constellation, btw): import numpy as np u = np.arange (-7.5, 8) const = np.empty ((16,16), dtype=complex) for ix in range (16): for iy in range (16): const[ix, iy] = u[ix] + 1j*u[iy]