Remove submodules and use numpydoc for documentation
Implement kf-like solver for block-banded symmetric positive-definite matrices
Add unittests for block-banded matrices
First version of Gaussian Markov chain VB node with a demo and a test