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

//! Inidividual enums for all bases who transform either
//! real-to-real, real-to-complex or complex-to-complex plus
//! a collective enum of enums, which combines these three types
//! and which can be used to put all Bases in a single array.
//!
//! Real-to-complex transforms implement the differentiate and
//! from ortho trait only for complex numbers
use crate::chebyshev::Chebyshev;
use crate::chebyshev::CompositeChebyshev;
use crate::fourier::FourierC2c;
use crate::fourier::FourierR2c;
use crate::traits::Differentiate;
use crate::traits::DifferentiatePar;
use crate::traits::FromOrtho;
use crate::traits::FromOrthoPar;
use crate::traits::Transform;
use crate::traits::TransformPar;
use crate::types::FloatNum;
use ndarray::prelude::*;
use num_complex::Complex;

#[enum_dispatch(Basics<T>, LaplacianInverse<T>)]
#[derive(Clone)]
/// Enum of enums which binds all bases
pub enum BaseAll<T: FloatNum> {
    /// Real-to-real bases
    BaseR2r(BaseR2r<T>),
    /// Real-to-complex bases
    BaseR2c(BaseR2c<T>),
    /// Complex-to-complex bases
    BaseC2c(BaseC2c<T>),
}

#[allow(clippy::large_enum_variant)]
#[enum_dispatch(Basics<T>, LaplacianInverse<T>)]
#[derive(Clone)]
/// All bases who transform: real-to-real
pub enum BaseR2r<T: FloatNum> {
    /// Chebyshev polynomials (orthogonal)
    Chebyshev(Chebyshev<T>),
    /// Chebyshev polynomials (composite)
    CompositeChebyshev(CompositeChebyshev<T>),
}

#[enum_dispatch(Basics<T>, LaplacianInverse<T>)]
#[derive(Clone)]
/// All bases who transform: real-to-complex
pub enum BaseR2c<T: FloatNum> {
    /// Fourier polynomials
    FourierR2c(FourierR2c<T>),
}

#[enum_dispatch(Basics<T>, LaplacianInverse<T>)]
#[derive(Clone)]
/// All bases who transform: complex-to-complex
pub enum BaseC2c<T: FloatNum> {
    /// Fourier polynomials
    FourierC2c(FourierC2c<T>),
}

// Implement traits on real-to-real
impl_transform_trait_for_base!(BaseR2r, A, A, Chebyshev, CompositeChebyshev);
impl_differentiate_trait_for_base!(BaseR2r, A, Chebyshev, CompositeChebyshev);
impl_differentiate_trait_for_base!(BaseR2r, Complex<A>, Chebyshev, CompositeChebyshev);
impl_from_ortho_trait_for_base!(BaseR2r, A, Chebyshev, CompositeChebyshev);
impl_from_ortho_trait_for_base!(BaseR2r, Complex<A>, Chebyshev, CompositeChebyshev);

// Implement traits on real-to-complex
impl_transform_trait_for_base!(BaseR2c, A, Complex<A>, FourierR2c);
impl_differentiate_trait_for_base!(BaseR2c, Complex<A>, FourierR2c);
impl_from_ortho_trait_for_base!(BaseR2c, Complex<A>, FourierR2c);

// Implement traits on complex-to-complex
impl_transform_trait_for_base!(BaseC2c, Complex<A>, Complex<A>, FourierC2c);
impl_differentiate_trait_for_base!(BaseC2c, Complex<A>, FourierC2c);
impl_from_ortho_trait_for_base!(BaseC2c, Complex<A>, FourierC2c);