Is there any efficient way to add many dirac deltas on a variational form?

Hi Nate,

I also have the same problem with a large number of point sources in nonlinear problem. Could you please give a simple example about how to use PointSource class with nonlinear problem? Since as far as I remember, you need to do

A,b = assemble_system(a, L)

and

pointsources.apply(b)

which works for the linear problem?

Thanks a lot in advance.

Kind regards,
Alice