SNEA solver written but not yet tested. Currently variance matrices are all identities