Laguerre–Galerkin methods with reduced sum-products