Efficient methods for Gaussian Markov random fields under sparse linear constraints