Add sparse triangular solver. Make A_ldiv_B a polyalgorithm, some reorganization and some tests.