Extending gaussian_dlm_obs

I need support for time-varying system matrices when fitting dynamic linear models, a feature not currently supported by gaussian_dlm_obs(), so I have filed an issue #2854 in the stan repository. I plan to implement this myself, and I figured I should let people know I’m working on this.