Some progress with triangular solver for compressed_matrix: unit-lower-solve...

Some progress with triangular solver for compressed_matrix: unit-lower-solve and upper-solve working, but performance is rather poor. This is, however, expected, because the data structure is not well suited for that.