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