Approximate implementation in Pyrex.

Mostly just a carbon copy of the python version, having some weird
issues with the numbering of roots/ghosts. Which is why
the test existed in the first place :).