working helmholtz solver in Julia