obtaining jacobian matrix of each element

Hello everyone.
I try to employ non-local averaging formulation for the stress field in a domain using subroutines. To do so I need the value of det [J] for each element. How can I obtain this value at each integration point? Is there any utility routine to be called for this purpose?