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
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116

//
// A rust binding for the GSL library by Guillaume Gomez (guillaume1.gomez@gmail.com)
//

/*!
The polygamma functions of order n are defined by

\psi^{(n)}(x) = (d/dx)^n \psi(x) = (d/dx)^{n+1} \log(\Gamma(x))

where \psi(x) = \Gamma'(x)/\Gamma(x) is known as the digamma function.
!*/

pub mod diagamma {
    use crate::{types, Value};
    use std::mem::MaybeUninit;

    /// This routine computes the digamma function \psi(n) for positive integer n. The digamma function is also called the Psi function.
    #[doc(alias = "gsl_sf_psi_int")]
    pub fn psi_int(n: i32) -> f64 {
        unsafe { sys::gsl_sf_psi_int(n) }
    }

    /// This routine computes the digamma function \psi(n) for positive integer n. The digamma function is also called the Psi function.
    #[doc(alias = "gsl_sf_psi_int_e")]
    pub fn psi_int_e(n: i32) -> Result<types::Result, Value> {
        let mut result = MaybeUninit::<sys::gsl_sf_result>::uninit();
        let ret = unsafe { sys::gsl_sf_psi_int_e(n, result.as_mut_ptr()) };

        result_handler!(ret, unsafe { result.assume_init() }.into())
    }

    /// This routine computes the digamma function \psi(x) for general x, x \ne 0.
    #[doc(alias = "gsl_sf_psi")]
    pub fn psi(x: f64) -> f64 {
        unsafe { sys::gsl_sf_psi(x) }
    }

    /// This routine computes the digamma function \psi(x) for general x, x \ne 0.
    #[doc(alias = "gsl_sf_psi_e")]
    pub fn psi_e(x: f64) -> Result<types::Result, Value> {
        let mut result = MaybeUninit::<sys::gsl_sf_result>::uninit();
        let ret = unsafe { sys::gsl_sf_psi_e(x, result.as_mut_ptr()) };

        result_handler!(ret, unsafe { result.assume_init() }.into())
    }

    /// This routine computes the real part of the digamma function on the line 1+i y, \Re[\psi(1 + i y)].
    #[doc(alias = "gsl_sf_psi_1piy")]
    pub fn psi_1piy(x: f64) -> f64 {
        unsafe { sys::gsl_sf_psi_1piy(x) }
    }

    /// This routine computes the real part of the digamma function on the line 1+i y, \Re[\psi(1 + i y)].
    #[doc(alias = "gsl_sf_psi_1piy_e")]
    pub fn psi_1piy_e(x: f64) -> Result<types::Result, Value> {
        let mut result = MaybeUninit::<sys::gsl_sf_result>::uninit();
        let ret = unsafe { sys::gsl_sf_psi_1piy_e(x, result.as_mut_ptr()) };

        result_handler!(ret, unsafe { result.assume_init() }.into())
    }
}

pub mod trigamma {
    use crate::{types, Value};
    use std::mem::MaybeUninit;

    /// This routine computes the Trigamma function \psi'(n) for positive integer n.
    #[doc(alias = "gsl_sf_psi_1_int")]
    pub fn psi_1_int(n: i32) -> f64 {
        unsafe { sys::gsl_sf_psi_1_int(n) }
    }

    /// This routine computes the Trigamma function \psi'(n) for positive integer n.
    #[doc(alias = "gsl_sf_psi_1_int_e")]
    pub fn psi_1_int_e(n: i32) -> Result<types::Result, Value> {
        let mut result = MaybeUninit::<sys::gsl_sf_result>::uninit();
        let ret = unsafe { sys::gsl_sf_psi_1_int_e(n, result.as_mut_ptr()) };

        result_handler!(ret, unsafe { result.assume_init() }.into())
    }

    /// This routine computes the Trigamma function \psi'(x) for general x.
    #[doc(alias = "gsl_sf_psi_1")]
    pub fn psi_1(x: f64) -> f64 {
        unsafe { sys::gsl_sf_psi_1(x) }
    }

    /// This routine computes the Trigamma function \psi'(x) for general x.
    #[doc(alias = "gsl_sf_psi_1_e")]
    pub fn psi_1_e(x: f64) -> Result<types::Result, Value> {
        let mut result = MaybeUninit::<sys::gsl_sf_result>::uninit();
        let ret = unsafe { sys::gsl_sf_psi_1_e(x, result.as_mut_ptr()) };

        result_handler!(ret, unsafe { result.assume_init() }.into())
    }
}

pub mod polygamma {
    use crate::{types, Value};
    use std::mem::MaybeUninit;

    /// This routine computes the polygamma function \psi^{(n)}(x) for n >= 0, x > 0.
    #[doc(alias = "gsl_sf_psi_n")]
    pub fn psi_n(n: i32, x: f64) -> f64 {
        unsafe { sys::gsl_sf_psi_n(n, x) }
    }

    /// This routine computes the polygamma function \psi^{(n)}(x) for n >= 0, x > 0.
    #[doc(alias = "gsl_sf_psi_n_e")]
    pub fn psi_n_e(n: i32, x: f64) -> Result<types::Result, Value> {
        let mut result = MaybeUninit::<sys::gsl_sf_result>::uninit();
        let ret = unsafe { sys::gsl_sf_psi_n_e(n, x, result.as_mut_ptr()) };

        result_handler!(ret, unsafe { result.assume_init() }.into())
    }
}