interval-arithmetic
An implementation of an algebraically closed interval system of the extended real number set
Description
An interval
is a pair of numbers which represents all the numbers between them, closed
means that the bounds are also included in the representation, extended real
because the
real number system
is extended with two elements: -∞
and +∞
representing negative infinity
and positive infinity respectively.
The implementation is a modified port of the Boost's interval arithmetic library, the modifications are based on some guidelines from the following papers/presentations:
- Interval Arithmetic: from Principles to Implementation - T. Hickey, Q. Ju, M.H. van Emden
- Interval Arithmetic: Python Implementation and Applications - Stefano Taschini
- The Boost interval arithmetic library - Hervé Brönnimann, Guillaume Melquiond, Sylvain Pion
- Graphing equations with generalized interval arithmetic - Jeffrey Allen Tupper
floating point operations
Floating point is a way to represent a real number in an approximate way (due to the finite space existing on a computer), most calculations with real numbers will produce quantities that cannot be exactly represented with the space allocated and therefore this operation needs to be rounded in order to fit back into its finite representation, such errors are described in more detail here
Interval arithmetic
Instead of using a single floating point number as an approximation of a real number, interval
arithmetic represents the approximated value as a set of possible values (considering the numbers
that floating point cannot represent), let's say we want to represent the number 1 / 3
, as a single
floating point number it's approximated as 0.333333333333...
, in the end there will be some 333...
decimals that will be lost due to the nature of floating point, instead we can represent this
number with the interval [0.2, 0.4]
, with this interval we're completely sure that 1 / 3
is within
the interval (although the interval is also representing many more numbers), to improve the scope
of the interval we have to understand that numbers in JavaScript are represented with 64 bits,
therefore to get the next floating point number of a single precision number the last bit
needs to be incremented to get the upper bound, and the last bit also needs to be decremented
to get the lower point
Notable modifications
- next/previous IEEE754 floating point number implementation based on Typed Arrays
-
division
when both intervals contain zero creates a whole interval -
cosine, tangent
works with positive/negative values out of the box
Interval arithmetic evaluator
Due to the expressive nature of the way the methods interact with intervals it's sad that even the simplest
operation needs a lot of characters to be typed, let's consider evaluating the result of 1 + 2
expressed
with intervals
Interval.add(new Interval(1, 1), new Interval(2, 2))
This gets worse when the expression to be evaluated becomes complex like sin(exp(x)) + tan(x) - 1/cos(PI) * [1, 3]^2
:
const x = Interval(0, 1);
Interval.add(
Interval.sin(Interval.exp(x)),
Interval.sub(
Interval.tan(x),
Interval.mul(
Interval.div(Interval.ONE, Interval.cos(Interval.PI)),
Interval.pow(Interval(1, 3), 2)
)
)
);
To avoid this 'expressiveness' mess there's an interval arithmetic evaluator module which I've created to deal with all the work of parsing/evaluating expressions like the one above
const compile = require('interval-arithmetic-eval');
compile('sin(exp(x)) + tan(x) - 1/cos(PI) * [1, 3]^2').eval({ x: [0, 1] })
Installation & Usage
$ npm install --save interval-arithmetic
import Interval from 'interval-arithmetic'
Interval.add(Interval(1), Interval(2))
API
See the homepage
Development
npm test
Deployment steps
// after the working directory is clean
npm test
npm run build
npm version (major|minor|patch)
git push origin master
// if everything went well
npm publish
npm run deploy
2015-2022 © Mauricio Poppe