Source code for the ab-splat algorithm