first 'working' version with scipy lin algebra routines - still not accurate (on test data)
commented out scipy routines and replaced with numpy - code now works again: should be easy to debug the cholesky solve usage