Automate upper_bound_log_pt

[hai-dang-dau]

This isn't a bug but rather a suggestion. Currently the user needs to code themselves the method upper_bound_log_pt(self, t) of the class StateSpaceModel in order to use smoothing algorithms. However most models in this package have Gaussian transitions, for which the upper bound could be calculated automatically.