Fast covariance estimation for innovations computed from a spatial Gibbs point process