Join GitHub today
GitHub is home to over 31 million developers working together to host and review code, manage projects, and build software together.Sign up
[operators.cg] fix first order integration #431
1st order integration for boundary integrals in operators.cg was broken.
All dimensions of size 1 of the array of shape functions, evaluated at integration points, were removed using "np.squeeze". This works fine as long as there is more than one integration point. When you integrate first order, there is only one integration point and the dimension for the integration point number is removed from the array.