Much faster additive approximation

Sped up additive approximation by computing the co-occurrence only from
the pixels affected by change.

Discovered an error in additive_approximation, where intersect instead
of union is still used.

No wet costs yet, the algorithm allows values -1 and 256!!