import qualified Algebra.Algebraic as Algebraic ( C )
import Data.Vector.Fixed ( ifoldl )
import Data.Vector.Fixed.Cont ( Arity )
-import Debug.Trace
import NumericPrelude hiding ( (*) )
import Linear.Matrix (
-- factorization. We keep the pair updated by multiplying @q@ and
-- @r@ by the new rotator (or its transpose).
--
+-- We do not require that the diagonal elements of R are positive,
+-- so our factorization is a little less unique than usual.
+--
-- Examples:
--
-- >>> import Linear.Matrix