Fast estimation of matrix exponential spatial models