Solving Stochastic PDEs with Approximate Gaussian Markov Random Fields using Different Programming Environments