pg-binary-matrix-factorization Using PG augmentation and VI to fit the model Z ~ N(0,I) X|Z ~ Bernoulli(sigmoid(AX))