Getting largest indices of array less than or equal to an array of numbers

[schniefen]

Current solution:

y_i=np.linspace(np.zeros(d),np.ones(d),m+1) 
y=np.random.uniform(0,1,d)                      
ind=np.zeros([len(y),d]) 
for k in range(d):
  k_i=y_i[:,k]
  i=np.argmax(np.repeat(k_i[:, np.newaxis], len(y), axis=1) >= y, axis=0)
  ind[:,k]=i