[Numpy-discussion] simple rectangular grid

[Neal Becker]

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]