@@ -581,9 +581,9 @@ def _compute_tri_eccentricities(tris_pts):
581581 The so-called eccentricity parameters [1] needed for
582582 HCT triangular element.
583583 """
584- a = np .expand_dims (tris_pts [:, 2 , :]- tris_pts [:, 1 , :], axis = 2 )
585- b = np .expand_dims (tris_pts [:, 0 , :]- tris_pts [:, 2 , :], axis = 2 )
586- c = np .expand_dims (tris_pts [:, 1 , :]- tris_pts [:, 0 , :], axis = 2 )
584+ a = np .expand_dims (tris_pts [:, 2 , :] - tris_pts [:, 1 , :], axis = 2 )
585+ b = np .expand_dims (tris_pts [:, 0 , :] - tris_pts [:, 2 , :], axis = 2 )
586+ c = np .expand_dims (tris_pts [:, 1 , :] - tris_pts [:, 0 , :], axis = 2 )
587587 # Do not use np.squeeze, this is dangerous if only one triangle
588588 # in the triangulation...
589589 dot_a = _prod_vectorized (_transpose_vectorized (a ), a )[:, 0 , 0 ]
@@ -1071,9 +1071,9 @@ def get_dof_vec(tri_z, tri_dz, J):
10711071 J1 = _prod_vectorized (_ReducedHCT_Element .J0_to_J1 , J )
10721072 J2 = _prod_vectorized (_ReducedHCT_Element .J0_to_J2 , J )
10731073
1074- col0 = _prod_vectorized (J , np .expand_dims (tri_dz [:, 0 , :], axis = 3 ))
1075- col1 = _prod_vectorized (J1 , np .expand_dims (tri_dz [:, 1 , :], axis = 3 ))
1076- col2 = _prod_vectorized (J2 , np .expand_dims (tri_dz [:, 2 , :], axis = 3 ))
1074+ col0 = _prod_vectorized (J , np .expand_dims (tri_dz [:, 0 , :], axis = 2 ))
1075+ col1 = _prod_vectorized (J1 , np .expand_dims (tri_dz [:, 1 , :], axis = 2 ))
1076+ col2 = _prod_vectorized (J2 , np .expand_dims (tri_dz [:, 2 , :], axis = 2 ))
10771077
10781078 dfdksi = _to_matrix_vectorized ([
10791079 [col0 [:, 0 , 0 ], col1 [:, 0 , 0 ], col2 [:, 0 , 0 ]],
0 commit comments