A parallel implementation of a high-order accurate solution technique for variable coefficient Helmholtz problems