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

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

//! The Gegenbauer polynomials are defined in Abramowitz & Stegun, Chapter 22, where they are known as Ultraspherical polynomials.

use crate::Value;
use std::mem::MaybeUninit;

/// This function evaluates the Gegenbauer polynomials C^{(\lambda)}_n(x) using explicit representations for n =1, 2, 3.
#[doc(alias = "gsl_sf_gegenpoly_1")]
pub fn gegenpoly_1(lambda: f64, x: f64) -> f64 {
    unsafe { sys::gsl_sf_gegenpoly_1(lambda, x) }
}

/// This function evaluates the Gegenbauer polynomials C^{(\lambda)}_n(x) using explicit representations for n =1, 2, 3.
#[doc(alias = "gsl_sf_gegenpoly_2")]
pub fn gegenpoly_2(lambda: f64, x: f64) -> f64 {
    unsafe { sys::gsl_sf_gegenpoly_2(lambda, x) }
}

/// This function evaluates the Gegenbauer polynomials C^{(\lambda)}_n(x) using explicit representations for n =1, 2, 3.
#[doc(alias = "gsl_sf_gegenpoly_3")]
pub fn gegenpoly_3(lambda: f64, x: f64) -> f64 {
    unsafe { sys::gsl_sf_gegenpoly_3(lambda, x) }
}

/// This function evaluates the Gegenbauer polynomials C^{(\lambda)}_n(x) using explicit representations for n =1, 2, 3.
#[doc(alias = "gsl_sf_gegenpoly_1_e")]
pub fn gegenpoly_1_e(lambda: f64, x: f64) -> (Value, ::types::Result) {
    let mut result = MaybeUninit::<sys::gsl_sf_result>::uninit();
    let ret = unsafe { sys::gsl_sf_gegenpoly_1_e(lambda, x, result.as_mut_ptr()) };

    (::Value::from(ret), unsafe { result.assume_init() }.into())
}

/// This function evaluates the Gegenbauer polynomials C^{(\lambda)}_n(x) using explicit representations for n =1, 2, 3.
#[doc(alias = "gsl_sf_gegenpoly_2_e")]
pub fn gegenpoly_2_e(lambda: f64, x: f64) -> (Value, ::types::Result) {
    let mut result = MaybeUninit::<sys::gsl_sf_result>::uninit();
    let ret = unsafe { sys::gsl_sf_gegenpoly_2_e(lambda, x, result.as_mut_ptr()) };

    (::Value::from(ret), unsafe { result.assume_init() }.into())
}

/// This function evaluates the Gegenbauer polynomials C^{(\lambda)}_n(x) using explicit representations for n =1, 2, 3.
#[doc(alias = "gsl_sf_gegenpoly_3_e")]
pub fn gegenpoly_3_e(lambda: f64, x: f64) -> (Value, ::types::Result) {
    let mut result = MaybeUninit::<sys::gsl_sf_result>::uninit();
    let ret = unsafe { sys::gsl_sf_gegenpoly_3_e(lambda, x, result.as_mut_ptr()) };

    (::Value::from(ret), unsafe { result.assume_init() }.into())
}

/// This function evaluates the Gegenbauer polynomial C^{(\lambda)}_n(x) for a specific value of n, lambda, x subject to \lambda > -1/2, n >= 0.
#[doc(alias = "gsl_sf_gegenpoly_n")]
pub fn gegenpoly_n(n: i32, lambda: f64, x: f64) -> f64 {
    unsafe { sys::gsl_sf_gegenpoly_n(n, lambda, x) }
}

/// This function evaluates the Gegenbauer polynomial C^{(\lambda)}_n(x) for a specific value of n, lambda, x subject to \lambda > -1/2, n >= 0.
#[doc(alias = "gsl_sf_gegenpoly_n_e")]
pub fn gegenpoly_n_e(n: i32, lambda: f64, x: f64) -> (Value, ::types::Result) {
    let mut result = MaybeUninit::<sys::gsl_sf_result>::uninit();
    let ret = unsafe { sys::gsl_sf_gegenpoly_n_e(n, lambda, x, result.as_mut_ptr()) };

    (::Value::from(ret), unsafe { result.assume_init() }.into())
}

/// This function computes an array of Gegenbauer polynomials C^{(\lambda)}_n(x) for n = 0, 1, 2, \dots, nmax, subject to \lambda > -1/2, nmax >= 0.
#[doc(alias = "gsl_sf_gegenpoly_array")]
pub fn gegenpoly_array(lambda: f64, x: f64, result_array: &mut [f64]) -> Value {
    Value::from(unsafe {
        sys::gsl_sf_gegenpoly_array(
            result_array.len() as _,
            lambda,
            x,
            result_array.as_mut_ptr(),
        )
    })
}