Nonsymmetric and per-state remainder estimation

[agerlach]

Currently the library only supports a fixed, symmetric remainder estimate, e.g.

remainder estimation 0.001

This is equivalent to

remainder estimation [-0.001, 0.001]

Non-symmetric is possible via

remainder estimation [-0.001, 0.005]

or state-dependent via

remainder estimation { x :[0.1 ,0.101] , y :[ -0.01 ,0.06] }

See https://home.cs.colorado.edu/~xich8622/manual/manual-2.0.0.pdf section 3