sarkka branch: getting the dispersion matrix right
nothing
simpler code: Q computation following Sarkka (removed abasset confusing matrix)
nothing
merge sarkka