Added changelog
More efficient and better behaving magnitude computation

Previously 'magnitude e' had to compute all magnitudes smaller than e
O(e). Now it has the same complexity as '^' O(log e) but caches
intermediate results for later re-use.