PointIntegralSolver

Equations resulting from a point source integral are fully decoupled in space. Instead of assemble a full matrix it is more efficient to iterate over the points and solve a local problem.

For linear equations it means solving a single dense linear system at each node, which we can use Armadillo for.
For nonlinear systems we need to use Newtons method together with linear solves.