Accurate computations with Riordan arrays associated with Schröder matrices

Large and little Schröder matrices as well as with their inverses are triangular matrices arising in Combinatorics. Through their total positivity and bidiagonal decomposition it is shown that their singular values, inverses and some associated linear systems can be solved with high relative accuracy. Numerical experiments confirm the theoretical results.