Fixup for DomainTuples

nifty5/operators/linear_interpolation.py

@@ -37,7 +37,6 @@ class LinearInterpolator(LinearOperator):
         Multilinear interpolation for points in an RGSpace
 
         :param domain:
-            RGSpace
         :param positions:
             positions at which to interpolate
             Field with UnstructuredDomain, shape (dim, ndata)
@@ -66,7 +65,7 @@ class LinearInterpolator(LinearOperator):
         pos = positions/dist
         excess = pos-pos.astype(int64)
         pos = pos.astype(int64)
-        max_index = array(self.domain[0].shape).reshape((-1,)+(1,))
+        max_index = array(self.domain.shape).reshape(-1, 1)
         data = zeros((len(mg[0]), N_points))
         ii = zeros((len(mg[0]), N_points), dtype=int64)
         jj = zeros((len(mg[0]), N_points), dtype=int64)