  Localized point sources (monopoles) in an acoustical domain are implemented
to a three dimensional non-singular Helmholtz boundary element method in the
frequency domain. It allows for the straightforward use of higher order surface
elements on the boundaries of the problem. It will been shown that the effect
of the monopole sources ends up on the right hand side of the resulting matrix
system. Some carefully selected examples are studied, such as point sources
near and within a concentric spherical core-shell scatterer (with theoretical
verification), near a curved focusing surface and near a multi-scale and
multi-domain acoustic lens.
