  The stochastic Laplacian Heaviside (LapH) method has proven to be successful
in hadronic calculations. In this work, with charm-light spectroscopy in mind,
we set up and optimise the LapH procedure limiting ourselves to the evaluation
of two-point mesonic correlation functions. The calculations are performed on
CLS ensembles with $N_f=2+1$ Wilson-Clover fermions on a $32^3\times64$ lattice
with open boundary conditions. We analyse the interplay between the LapH
parameters and the boundary effects, and implement a fitting procedure to
isolate excitations coming from the border.
