Fiddling and thinking about how to build QR decomp out of

the current householder/reflect/matrixmult code I've got instead of
doing it in one massive kernel like the book