Add point source (delta function) in bilinear form

Hi all,

I am solving a non-linear Poisson problem with Newton’s method, and there is a point source. I checked our old forum and only found this Point source of electric current in radially symmetric cylindrical coordinates - Kunena (ngsolve.org).

However, now it seems that I need to put the point source into the a = BilinearForm(V) according to the documentation page 3.3 Nonlinear problems — NGS-Py 6.2.2404-2-ge6a630133 documentation (ngsolve.org).

Is there any faster way to do this other than approximating the delta function by some smooth exponential function? Thanks in advance for any possible help!