Initialization: use QR factorization of D^T
Merge pull request #8 from veeranjaneyulus/master

Initialization: use QR factorization of D^T