The authors consider the solution of delay differential equations using the boundary value (BVM) methods that have been introduced in recent years to solve initial value problems. The authors begin by explaining the background and introduce the BVM scheme, showing that its application leads to the need to solve a large linear system. They introduce Strang-type preconditioners for solving the linear system using the generalized minimal residual algorithm and perform a spectral analysis to assess the convergence rate of the resulting method. Finally they give numerical examples using Matlab to demonstrate the efficacy of their approach.