X-Git-Url: https://fleuret.org/cgi-bin/gitweb/gitweb.cgi?a=blobdiff_plain;f=eingather.py;h=e73e63ff384640a8f5edec5f3e2f347c4a59d0e2;hb=a5ee870f51cbe71921166fb6466a12f1382e670e;hp=b27110022263045413c87e0e4327ba2902aca283;hpb=7d526f2fa3849a4a049b3a1da4dddb92e3132428;p=pytorch.git diff --git a/eingather.py b/eingather.py index b271100..e73e63f 100755 --- a/eingather.py +++ b/eingather.py @@ -19,40 +19,27 @@ def eingather(op, src, *indexes): s_all = "".join([s_src] + s_indexes) shape = tuple(all_sizes[s_all.index(v)] for v in s_dst) - idx = [] - n_index = 0 - - for i in range(src.dim()): - v = s_src[i] - if v == "_": - index, s_index = indexes[n_index], s_indexes[n_index] - n_index += 1 - - sub_idx = [] - - for i in range(index.dim()): - v = s_index[i] + def do(x,s_x): + idx = [] + n_index = 0 + + for i in range(x.dim()): + v = s_x[i] + if v == "_": + idx.append(do(indexes[n_index], s_indexes[n_index])) + n_index += 1 + else: j = s_dst.index(v) a = ( - torch.arange(index.size(i)) + torch.arange(x.size(i)) .reshape((1,) * j + (-1,) + (1,) * (len(s_dst) - j - 1)) .expand(shape) ) - sub_idx.append(a) - - index = index[sub_idx] - idx.append(index) - else: - j = s_dst.index(v) - a = ( - torch.arange(src.size(i)) - .reshape((1,) * j + (-1,) + (1,) * (len(s_dst) - j - 1)) - .expand(shape) - ) - idx.append(a) + idx.append(a) - return src[idx] + return x[idx] + return do(src,s_src) ####################### @@ -60,7 +47,7 @@ src = torch.rand(3, 5, 7, 11) index1 = torch.randint(src.size(2), (src.size(3), src.size(1), src.size(3))) index2 = torch.randint(src.size(3), (src.size(1),)) -# I want result[a, c, e] = src[c, a, index1[e, a, e], index2[a], e] +# I want result[a, c, e] = src[c, a, index1[e, a, e], index2[a]] result = eingather("ca(eae)(a) -> ace", src, index1, index2)