BUG: Correct invocation of expand_dims

eric-wieser · web-flow · commit 5cea9ce4d4ee · 2017-05-11T15:00:54.000+01:00
`alpha` is 2d, yet here it was requested to insert a fourth dimension.

This rightly produces an error in numpy 1.13.0rc1, although we might soften this to a deprecation warning
diff --git a/lib/matplotlib/tri/triinterpolate.py b/lib/matplotlib/tri/triinterpolate.py
@@ -898,7 +898,7 @@ def get_bending_matrices(self, J, ecc):
         pts = self.gauss_pts
         for igauss in range(self.n_gauss):
             alpha = np.tile(pts[igauss, :], n).reshape(n, 3)
-            alpha = np.expand_dims(alpha, 3)
+            alpha = np.expand_dims(alpha, 2)
             weight = weights[igauss]
             d2Skdksi2 = self.get_d2Sidksij2(alpha, ecc)
             d2Skdx2 = _prod_vectorized(d2Skdksi2, H_rot)
