Implementation of direct and iterative methods for a class of Helmholtz wave problems