Added a method to do a direct Q calculation via the dtsqr function