sparsetools 0.2.4

2-D sparse matrix package for numeric data
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
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
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168

use pretty_dtoa::{dtoa, FmtFloatConfig};
use std::{fmt, ops};

use crate::util::DEFAULT_FLOAT_CONFIG;

pub trait Integer:
    num_traits::PrimInt + ops::AddAssign + ops::SubAssign + fmt::Display + fmt::Debug
{
}

impl Integer for usize {}
impl Integer for u8 {}
impl Integer for u16 {}
impl Integer for u32 {}
impl Integer for u64 {}
impl Integer for u128 {}

impl Integer for isize {}
impl Integer for i8 {}
impl Integer for i16 {}
impl Integer for i32 {}
impl Integer for i64 {}
impl Integer for i128 {}

pub trait Float: Scalar {}

impl Float for f64 {}
impl Float for f32 {}

pub trait Complex<F: Float>: Scalar {
    fn new(re: F, im: F) -> Self;
    fn conj(&self) -> Self;
    fn real(&self) -> F;
    fn imag(&self) -> F;
}

impl Complex<f64> for num_complex::Complex<f64> {
    fn new(re: f64, im: f64) -> Self {
        Self::new(re, im)
    }

    fn conj(&self) -> Self {
        self.conj()
    }

    fn real(&self) -> f64 {
        self.re
    }

    fn imag(&self) -> f64 {
        self.im
    }
}

impl Complex<f32> for num_complex::Complex<f32> {
    fn new(re: f32, im: f32) -> Self {
        Self::new(re, im)
    }

    fn conj(&self) -> Self {
        self.conj()
    }

    fn real(&self) -> f32 {
        self.re
    }

    fn imag(&self) -> f32 {
        self.im
    }
}

pub trait Scalar<Output = Self>:
    PartialEq
    + Copy
    + Clone
    + ops::Add<Output = Self>
    + ops::Sub<Output = Self>
    + ops::Mul<Output = Self>
    + ops::Div<Output = Self>
    + ops::AddAssign
    + ops::SubAssign
    + ops::MulAssign
    + num_traits::Zero
    + num_traits::One<Output = Self>
    + fmt::Display
    + Nonzero
{
    fn pretty_string(&self, _config: Option<FmtFloatConfig>) -> String {
        format!("{}", self)
    }
}

impl Scalar for f64 {
    fn pretty_string(&self, config: Option<FmtFloatConfig>) -> String {
        dtoa(*self, config.unwrap_or(DEFAULT_FLOAT_CONFIG))
    }
}
impl Scalar for f32 {
    fn pretty_string(&self, config: Option<FmtFloatConfig>) -> String {
        dtoa(*self as f64, config.unwrap_or(DEFAULT_FLOAT_CONFIG))
    }
}

impl Scalar for num_complex::Complex<f64> {
    fn pretty_string(&self, config: Option<FmtFloatConfig>) -> String {
        format!(
            "{}{}j{}",
            dtoa(self.re, config.unwrap_or(DEFAULT_FLOAT_CONFIG)),
            if self.im.signum() < 0.0 { "-" } else { "+" },
            dtoa(self.im, config.unwrap_or(DEFAULT_FLOAT_CONFIG))
        )
        .to_string()
    }
}
impl Scalar for num_complex::Complex<f32> {
    fn pretty_string(&self, config: Option<FmtFloatConfig>) -> String {
        format!(
            "{}{}j{}",
            dtoa(self.re as f64, config.unwrap_or(DEFAULT_FLOAT_CONFIG)),
            if self.im.signum() < 0.0 { "-" } else { "+" },
            dtoa(self.im as f64, config.unwrap_or(DEFAULT_FLOAT_CONFIG))
        )
        .to_string()
    }
}

impl Scalar for usize {}
impl Scalar for u8 {}
impl Scalar for u16 {}
impl Scalar for u32 {}
impl Scalar for u64 {}
impl Scalar for u128 {}

impl Scalar for isize {}
impl Scalar for i8 {}
impl Scalar for i16 {}
impl Scalar for i32 {}
impl Scalar for i64 {}
impl Scalar for i128 {}

pub trait Nonzero: PartialEq + Default {
    fn nonzero(&self) -> bool {
        *self != Default::default()
    }
}

impl Nonzero for f64 {}
impl Nonzero for f32 {}

impl Nonzero for num_complex::Complex<f64> {}
impl Nonzero for num_complex::Complex<f32> {}

impl Nonzero for bool {}

impl Nonzero for usize {}
impl Nonzero for u8 {}
impl Nonzero for u16 {}
impl Nonzero for u32 {}
impl Nonzero for u64 {}
impl Nonzero for u128 {}

impl Nonzero for isize {}
impl Nonzero for i8 {}
impl Nonzero for i16 {}
impl Nonzero for i32 {}
impl Nonzero for i64 {}
impl Nonzero for i128 {}