[Numpy-discussion] Any easy way to do this?

[Neal Becker]

given: w[i,j,k], y[l, k]

find:
d[l,i,j] = norm(w[i,j] - y[l])

for each triple (l,i,j), w[i,j]-y[l] is a vector, of which I want to find the 
norm, and store into d[l,i,j]