Added Triangular expression to extract upper or lower (strictly or not) part of a matrix.
Triangular also provide an optimised method for forward and backward substitution. Further optimizations regarding assignments and products might come later.

Updated determinant() to take into account triangular matrices.

Started the QR module with a QR decompostion algorithm.

Help needed to build a QR algorithm (eigen solver) based on it.