The strain and stress inside of an element is bilinear. Therefore the bmatrix needs to be calculated for every MP in the element. This wasn't being done. Abdel fixed it on his branch and those changes have been pullled onto the Anura3D_QuadMoore branch.