Relativistic calculations of electric field gradients using the Douglas–Kroll method