honestintervals 0.2.0

Comprehensive, correctly rounding interval arithmetic library.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91

use libc::{c_char, c_double, c_int, c_long, c_longlong};

/// MPFR rounding mode enum.
#[repr(C)]
pub enum MpfrRnd {
    /// Round to the nearest representable number. Ties goes to the even number.
    HalfToEven = 0,
    /// Round towards zero.
    TowardsZero,
    /// Round towards positive infinity.
    Up,
    /// Round towards negative infinity.
    Down,
    /// Round away from zero.
    AwayFromZero,
    /// Round to the nearest representable number. Ties goes away from zero.
    HalfAwayFromZero = -1,
}

/// MPFR precision type.
pub type MpfrPrec = c_long;
/// MPFR sign type.
pub type MpfrSign = c_int;
/// MPFR exponent type.
pub type MpfrExp = c_long;
/// MPFR limb type.
pub type MpLimb = c_longlong;

/// Low-level MPFR struct.
#[repr(C)]
#[derive(Debug)]
pub struct MpfrStruct {
    pub _mpfr_prec: MpfrPrec,
    pub _mpfr_sign: MpfrSign,
    pub _mpfr_exp: MpfrExp,
    pub _mpfr_d: *mut MpLimb,
}

type MpfrPtr = *mut MpfrStruct;
type MpfrConstPtr = *const MpfrStruct;

/// External calls to MPFR library. See http://www.mpfr.org/mpfr-current/mpfr.html for details.
#[allow(missing_docs)]
#[link(name = "mpfr")]
extern "C" {
    // Initialization functions.
    pub fn mpfr_init2(x: MpfrPtr, prec: MpfrPrec);
    pub fn mpfr_clear(x: MpfrPtr);

    // Assignment functions.
    pub fn mpfr_set(rop: MpfrPtr, op: MpfrConstPtr, rnd: MpfrRnd) -> c_int;
    pub fn mpfr_set_nan(x: MpfrPtr);
    pub fn mpfr_set_inf(x: MpfrPtr, sign: c_int);
    pub fn mpfr_set_zero(x: MpfrPtr, sign: c_int);
    pub fn mpfr_set_d(rop: MpfrPtr, op: c_double, rnd: MpfrRnd) -> c_int;
    pub fn mpfr_set_str(rop: MpfrPtr, s: *const c_char, base: c_int, rnd: MpfrRnd) -> c_int;

    // Conversion functions.
    pub fn mpfr_get_d(op: MpfrConstPtr, rnd: MpfrRnd) -> c_double;

    // Comparison functions.
    pub fn mpfr_cmp(op1: MpfrConstPtr, op2: MpfrConstPtr) -> c_int;
    pub fn mpfr_nan_p(op: MpfrConstPtr) -> c_int;
    pub fn mpfr_inf_p(op: MpfrConstPtr) -> c_int;
    pub fn mpfr_number_p(op: MpfrConstPtr) -> c_int;
    pub fn mpfr_zero_p(op: MpfrConstPtr) -> c_int;
    pub fn mpfr_sgn(op: MpfrConstPtr) -> c_int;
    pub fn mpfr_greater_p(op1: MpfrConstPtr, op2: MpfrConstPtr) -> c_int;
    pub fn mpfr_greaterequal_p(op1: MpfrConstPtr, op2: MpfrConstPtr) -> c_int;
    pub fn mpfr_less_p(op1: MpfrConstPtr, op2: MpfrConstPtr) -> c_int;
    pub fn mpfr_lessequal_p(op1: MpfrConstPtr, op2: MpfrConstPtr) -> c_int;
    pub fn mpfr_equal_p(op1: MpfrConstPtr, op2: MpfrConstPtr) -> c_int;
    pub fn mpfr_unordered_p(op1: MpfrConstPtr, op2: MpfrConstPtr) -> c_int;

    // Basic arithmetic functions.
    pub fn mpfr_neg(rop: MpfrPtr, op: MpfrConstPtr, rnd: MpfrRnd) -> c_int;
    pub fn mpfr_abs(rop: MpfrPtr, op: MpfrConstPtr, rnd: MpfrRnd) -> c_int;
    pub fn mpfr_add(rop: MpfrPtr, op1: MpfrConstPtr, op2: MpfrConstPtr, rnd: MpfrRnd) -> c_int;
    pub fn mpfr_sub(rop: MpfrPtr, op1: MpfrConstPtr, op2: MpfrConstPtr, rnd: MpfrRnd) -> c_int;
    pub fn mpfr_mul(rop: MpfrPtr, op1: MpfrConstPtr, op2: MpfrConstPtr, rnd: MpfrRnd) -> c_int;
    pub fn mpfr_div(rop: MpfrPtr, op1: MpfrConstPtr, op2: MpfrConstPtr, rnd: MpfrRnd) -> c_int;
    pub fn mpfr_pow(rop: MpfrPtr, op1: MpfrConstPtr, op2: MpfrConstPtr, rnd: MpfrRnd) -> c_int;

    // Special functions.
    pub fn mpfr_log(rop: MpfrPtr, op: MpfrConstPtr, rnd: MpfrRnd) -> c_int;
    pub fn mpfr_exp(rop: MpfrPtr, op: MpfrConstPtr, rnd: MpfrRnd) -> c_int;

    // Miscellaneous Functions
    pub fn mpfr_min(rop: MpfrPtr, op1: MpfrConstPtr, op2: MpfrConstPtr, rnd: MpfrRnd) -> c_int;
    pub fn mpfr_max(rop: MpfrPtr, op1: MpfrConstPtr, op2: MpfrConstPtr, rnd: MpfrRnd) -> c_int;
}