Learning and inference for Gaussian linear dynamical systems with fast Cython and BLAS/LAPACK implementations