Use reduced integration in flux elements as well

(have to be consistent otherwise Newton method won't converge)