2D acceptance using flood fill

[oscarxblanco]

Dear all,
E. Serra and I have written an implementation of the floodfill algorithm used by B.~Riemman et al. in this article.
https://journals.aps.org/prab/abstract/10.1103/PhysRevAccelBeams.27.094002

The matlab version tracks one particle at the time, while the pyat version is parallel. I attach the two files.
atfloodfill.m
floodfill_aceptance.py
floodfill_acceptance.py
atfloodfill.m.txt

This will allow you to look for the acceptance of a ring in a quick scan going from unstable to the stable boundary. The number of turns the particle does before being lost is part of the outputs, allowing you to have a look outside the acceptance.

I wonder if this could be added to AT matlab and pyat, and if so, where should I put it. The python implementation needs non-standard parallel tracking settings based on a multiprocessing Queue.