We implement conservative interval arithmetic in Lean, on top of a software implementation of
floating point (since Lean's Float
is untrusted). The key types are
Floating
: 64+64 bit software floating point, with 64 bits of mantissa (roughly, as the high bit is explicit) and 64 bits of exponent. All arithmetic is conservatively rounded up or down.Interval
: Lower and upperFloating
bounds representing a conservative interval.Box
: Real and imaginaryInterval
s approximating a complex number.Approx A R
: Typeclass that says thatA
is a conservative approximation ofR
. The key examples areApprox Interval ℝ
andApprox Box ℂ
. There are a variety of associated typeclasses showing that particular arithmetic operations are conservative (respect the approximation), such asApproxAdd
andApproxField
.
On top of Interval
, we implement conservative versions of field arithmetic and various special
functions (exp
, log
, powers).
- Install
elan
(brew install elan-init
on Mac) lake build
Here is a partial list of contributors, in alphabetical order:
Plus migration fixes from Yury Kudryashov, David Renshaw, and Scott Morrison. Thank you!
I'm going to keep a list here, since otherwise I will forget them.
-- Tracing options
set_option trace.aesop true in
-- Print compiler IR, such as to see how well inlining worked:
set_option trace.compiler.ir.result true in
def ...