Again fix the precomputed default behaviour
Add precision arg to integral function
[WIP] Add some more unit-testing