reachability code is now rectangular-capable
Sparse triangular solver with trapezoidal matrices
LU passe pour n <= m (mais pas m > n)
Fixing bug in test
partially fixing problem with rectangular LU