PBD implementation of Matthias Müller paper