Three-dimensional h-adaptivity for the multigroup neutron diffusion equations