Banded Matrix Operators for Gaussian Markov Models in the Automatic Differentiation Era