3 # Any copyright is dedicated to the Public Domain.
4 # https://creativecommons.org/publicdomain/zero/1.0/
6 # Written by Francois Fleuret <francois@fleuret.org>
13 def eingather(op, src, *indexes):
14 s_src, s_dst = re.search("^([^ ]*) *-> *(.*)", op).groups()
15 s_indexes = re.findall("\(([^)]*)\)", s_src)
16 s_src = re.sub("\([^)]*\)", "_", s_src)
18 all_sizes = tuple(d for s in (src,) + indexes for d in s.size())
19 s_all = "".join([s_src] + s_indexes)
20 shape = tuple(all_sizes[s_all.index(v)] for v in s_dst)
26 for i in range(x.dim()):
29 idx.append(do(indexes[n_index], s_indexes[n_index]))
34 torch.arange(x.size(i))
35 .reshape((1,) * j + (-1,) + (1,) * (len(s_dst) - j - 1))
45 def lambda_eingather(op, src_shape, *indexes_shape):
46 s_src, s_dst = re.search("^([^ ]*) *-> *(.*)", op).groups()
47 s_indexes = re.findall("\(([^)]*)\)", s_src)
48 s_src = re.sub("\([^)]*\)", "_", s_src)
50 all_sizes = tuple(d for s in (src_shape,) + indexes_shape for d in s)
51 s_all = "".join([s_src] + s_indexes)
52 shape = tuple(all_sizes[s_all.index(v)] for v in s_dst)
58 for i in range(len(x_shape)):
61 f = do(indexes_shape[n_index], s_indexes[n_index])
62 idx.append(lambda indexes: indexes[n_index][f(indexes)])
67 torch.arange(x_shape[i])
68 .reshape((1,) * j + (-1,) + (1,) * (len(s_dst) - j - 1))
71 idx.append(lambda indexes: a)
74 return lambda indexes: [f(indexes) for f in idx]
76 f = do(src_shape, s_src)
78 return lambda src, *indexes: src[f(indexes)]
81 ######################################################################
83 # src = torch.rand(3, 5, 3)
85 # print(eingather("aba -> ab", src))
87 # f = lambda_eingather("aba -> ab", src.shape)
93 ######################################################################
95 src = torch.rand(3, 5, 7, 11)
96 index1 = torch.randint(src.size(2), (src.size(3), src.size(1), src.size(3)))
97 index2 = torch.randint(src.size(3), (src.size(1),))
99 # f = lambda_eingather("ca(eae)(a) -> ace", src.shape, index1.shape, index2.shape)
101 # print(f(src, index1, index2))
103 # result[a, c, e] = src[c, a, index1[e, a, e], index2[a]]
105 # result = eingather("ca(eae)(a) -> ace", src, index1, index2)
107 from functorch.dim import dims
110 result = src[c, a, index1[e, a, e], index2[a]].order(a, c, e)
116 for a in range(result.size(0)):
117 for c in range(result.size(1)):
118 for e in range(result.size(2)):
119 error += (result[a, c, e] - src[c, a, index1[e, a, e], index2[a]]).abs()