(BUGFIX) missing None when evaluating gradient

Author: lachlangrose

LoopStructural/interpolators/structured_tetra.py

@@ -130,7 +130,7 @@ def evaluate_gradient(self, pos, prop):
         vertices, element_gradients, tetras, inside = self.get_tetra_gradient_for_location(pos)
         vertex_vals = self.properties[prop][tetras]
         #grads = np.zeros(tetras.shape)
-        values[inside,:] = (element_gradients[inside,:,:]*self.properties[prop][tetras[inside,:]]).sum(2)
+        values[inside,:] = (element_gradients[inside,:,:]*self.properties[prop][tetras[inside,None,:]]).sum(2)
         length = np.sum(values[inside,:],axis=1)
         # values[inside,:] /= length[:,None]
         return values