Add "projection" updates

[marvinpfoertner]

When projecting a Gaussian random variable x via Ax onto some (lower-dimensional) subspace, the covariance matrix of the transformed random variable is given by \Sigma' = A \Sigma A^T. The corresponding update to the Cholesky factors can then be realized by means of a QR decomposition.