This will be an awesome capability.  At the moment I can’t use mesh adaptivity when I have a lot of dirac’s because it takes far too long to search for the elements - I mean the simulations run perhaps 10^3 or 10^4 times slower because of the search.  But with this capability they’ll be zooming along!

a