Matrix.QR.Householder
leastSquares
decompose
solve
det
data Reflection i a
reflectMatrix
reflectVector
data Upper i j a
matrixFromUpper
solveUpper
detUpper