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
#![cfg_attr(test, feature(test))]
#[cfg(test)]
extern crate assert;
#[cfg(test)]
extern crate test;
#[cfg(test)]
extern crate rand;
extern crate libc;
extern crate num;
pub use beta::{ln_beta, inc_beta, inv_inc_beta};
mod beta;
#[link_name = "m"]
mod m {
use libc::{c_double, c_int};
extern {
pub fn erf(x: c_double) -> c_double;
pub fn erfc(x: c_double) -> c_double;
pub fn lgamma_r(x: c_double, sign: &mut c_int) -> c_double;
}
}
#[inline]
pub fn erf(x: f64) -> f64 {
unsafe { m::erf(x) }
}
#[inline]
pub fn erfc(x: f64) -> f64 {
unsafe { m::erfc(x) }
}
#[inline]
pub fn ln_gamma(x: f64) -> (f64, i32) {
let mut sign: i32 = 0;
let value = unsafe { m::lgamma_r(x, &mut sign) };
(value, sign)
}