Thanks!
On Thu, May 20, 2021 at 9:53 AM Robert Kern

On Thu, May 20, 2021 at 9:47 AM Neal Becker

wrote: This seems like something that can be done with indexing, but I haven't found the solution.

out is a 2D array is initialized to zeros. x is a 1D array whose values correspond to the columns of out. For each row in out, set out[row,x[row]] = 1. Here is working code: def orthogonal_mod (x, nbits): out = np.zeros ((len(x), 1<

Any idea to do this without an explicit python loop?

i = np.arange(len(x)) j = x[i] out[i, j] = 1

-- Robert Kern _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@python.org https://mail.python.org/mailman/listinfo/numpy-discussion

-- Those who don't understand recursion are doomed to repeat it