mirror of
https://github.com/AdaCore/why3.git
synced 2026-02-12 12:34:55 -08:00
ae5ddf9151
This adds an optional dependency on MPFR to run why3execute for floats. It is also used for reals (represented as intervals of floats). This commit does the following changes: - update the configure/Makefile to allow MPFR dependency - adds a MPFR wrapper so that why3execute can still be compiled if MPFR is not installed. In this case, an exception is raised when executing on real/float. - some updates are made to the standard library so that real/float functions are part of the program world (and can be executed). - pinterp changes to make elementary functions from real and float executable. - add some tests under bench/interp for manual testing of this feature Note that the correctness of the results given for reals depends on the precision. A too low precision may give unexpected results.
75 lines
1.8 KiB
Plaintext
75 lines
1.8 KiB
Plaintext
module N
|
|
|
|
use ieee_float.Float64
|
|
|
|
let f (x: t)
|
|
=
|
|
add RNE x (1.0: t)
|
|
|
|
let g () =
|
|
add RNE (24.0: t) (1.0: t)
|
|
|
|
exception BenchFailure
|
|
|
|
let bench1 () raises { BenchFailure -> false } =
|
|
let x = f (1.0: t) in
|
|
if not (eq x (2.0: t)) then raise BenchFailure;
|
|
x
|
|
|
|
let bench2 () raises { BenchFailure -> false } =
|
|
let o = (0x1p-52: t) in
|
|
let x = sqrt RNE ((1.0: t) .+ o .- (1.0: t) .- o) in (* Expected result: 0.0 *)
|
|
if not (eq x (0.0: t)) then raise BenchFailure;
|
|
let o = (0x1p-53: t) in
|
|
let x = sqrt RNE ((1.0: t) .+ o .- (1.0: t) .- o) in (* Expected result: Nan *)
|
|
if eq x x then raise BenchFailure;
|
|
x
|
|
|
|
let bench3 ()
|
|
raises { BenchFailure -> false }
|
|
ensures { eq result (0.0: t)} =
|
|
let o = (0x1p-1023: t) in
|
|
let v = o .* o in
|
|
if not (eq v (0.0:t)) then raise BenchFailure;
|
|
v
|
|
|
|
let bench4 ()
|
|
raises { BenchFailure -> false }
|
|
ensures { is_infinite result } =
|
|
let o = (0x1p+1023: t) in
|
|
let res = o .* o in
|
|
if not (is_infinite res) then raise BenchFailure;
|
|
res
|
|
|
|
let bench5 ()
|
|
raises { BenchFailure -> false }
|
|
ensures { eq result (1.0:t) } =
|
|
let o = add RNE (1.0:t) (0x1p-53:t) in
|
|
if not (eq o (1.0:t)) then raise BenchFailure;
|
|
o
|
|
|
|
let bench6 ()
|
|
raises { BenchFailure -> false }
|
|
ensures { not (eq result (1.0:t)) } =
|
|
let o = add RNE (1.0:t) (0x1p-52:t) in
|
|
if eq o (1.0:t) then raise BenchFailure;
|
|
o
|
|
|
|
let bench7 ()
|
|
raises { BenchFailure -> false }
|
|
ensures { eq result (0.0: t)} =
|
|
let o = (0x1p-1074: t) in
|
|
let o2 = o ./ (2.0: t) in
|
|
if not (eq o2 (0.0:t)) then raise BenchFailure;
|
|
o2
|
|
|
|
let bench8 ()
|
|
raises { BenchFailure -> false }
|
|
ensures { not (eq result (0.0: t))} =
|
|
let o = (0x1p-1073: t) in
|
|
let o2 = o ./ (2.0: t) in
|
|
if (eq o2 (0.0:t)) then raise BenchFailure;
|
|
o2
|
|
|
|
end
|