Draft: Adding Carlson's sqrt filter
The idea is to add an alternative method to the parametrised Kalman Filter, using the more numerically stable Carlson's sqrt filter. The branch will explore the differences in line selection and plots first, and then, if results are attractive, it will implement a version that keeps the covariance matrix decomposed to save memory and computation (along the necessary methods and data structures).