Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition