The computational approach for solving the Faddeev-Merkuriev equations in total orbital momentum representation is presented. These equations describe a system of three quantum charged particles and are widely used in bound state and scattering calculations. The approach is based on the spline collocation method and exploits intensively the tensor product form of discretized operators and preconditioner, which leads to a drastic economy in both computer resources and time.