return r

def org_app(r, x, y, c, i_max, j_max):
    for i in range(i_max):
        for j in range(j_max):
            r[i, j, x[i, j], y[i, j]] = c[i, j]
return r
Оптимизация операции индексирования numpy ndarray