gem.Sum with multiple arguments

Author: pbrubeck

finat/physically_mapped.py

@@ -268,7 +268,7 @@ def matvec(self, table):
         ii = gem.indices(len(table.shape)-1)
         phi = [gem.Indexed(table, (j, *ii)) for j in range(self.M.shape[1])]
         # the sum approach is faster than calling numpy.dot or gem.IndexSum
-        exprs = [gem.ComponentTensor(reduce(gem.Sum, (self.M.array[i, j] * phi[j] for j in js)), ii)
+        exprs = [gem.ComponentTensor(gem.Sum(self.M.array[i, j] * phi[j] for j in js)), ii)
                  for i, js in enumerate(self.csr)]
 
         val = gem.ListTensor(exprs)