Extremely slow but does not overflow

Additive approximation now works, but would take about a week on a
cluster.

Co-occurrence did not work properly as the quartet_to_index was
producing too large numbers. Now the function outputs a structure with a
list of encountered quartets and related co-occurrence bins. Furter
optimization urgently needed.