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

//! Matrix storage schemes.

mod band;
mod compressed;
mod dense;
mod diagonal;
mod packed;

pub use band::BandMatrix;
pub use compressed::{CompressedMatrix, CompressedFormat};
pub use dense::DenseMatrix;
pub use diagonal::DiagonalMatrix;
pub use packed::{PackedMatrix, PackedFormat};

/// An element of a matrix.
pub trait Element: Copy {
    /// Return the zero element.
    fn zero() -> Self;
}

macro_rules! implement(
    ($kind:ty, $zero:expr) => (
        impl Element for $kind {
            #[inline(always)]
            fn zero() -> Self {
                $zero
            }
        }
    );
);

implement!(u8, 0);
implement!(u16, 0);
implement!(u32, 0);
implement!(u64, 0);

implement!(i8, 0);
implement!(i16, 0);
implement!(i32, 0);
implement!(i64, 0);

implement!(f32, 0.0);
implement!(f64, 0.0);

implement!(isize, 0);
implement!(usize, 0);