Discrete constitutive equations over hexahedral grids for eddy-current problems