adding, but slower than p451.py