On computing gradients of potential field data in the space domain