In direct and large eddy simulations, very small space steps are used close to the solid walls in order to resolve the boundary-layer structures. Due to the restrictive CFL stability criteria of explicit time-stepping schemes, the maximum allowable time step is also very small, leading to high computational costs, notably for converging flow statistics. The use of an implicit integration scheme may overcome this limitation at the price of an increased computational cost per step. Moreover, the most commonly used fully implicit schemes induce higher errors due to the necessary approximations and poor dispersion and dissipation properties. As a compromise, a fourth-order implicit residual smoothing scheme (IRS4), successfully validated for a finite volume solver in Cinnella and Content (2016); Hoarau and Cinnella (2020), has been introduced in a multiblock high-order finite-difference solver. Several improvements are proposed to ensure better dissipation properties, a more efficient treatment of physical boundaries and an accurate and stable parallel multiblock implementation. For moderate CFL numbers, a similar accuracy as the explicit method is obtained with substantial savings in terms of computational time.