only fix posterior feature expectations according to p(z|x) starting from the third iteration + refactoring