Initialization: use QR factorization of D^T