I get different signs on the Jacobian's determinant for different Gauss integration points in a square 20 node hexa element, can that be correct? MWE

Question

I\'m trying to calculate the stiffness- and mass matrix contributions for a 3D case, using 20 node hexahedron elements.

I get very strange results though. The Jacobia