poincare-lib 0.1.0

GPU-accelerated 3D plotting library for mathematical functions and scientific visualisation
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

use crate::coordinate::CoordinateSystem;
use crate::domain::{DataBounds, Domain};
use crate::plot_object::{PlotComponent, PlotGeometry, PlotObject};
use crate::resolution::Resolution;
use crate::style::PlotStyle;

struct Piece {
    domain: Domain,
    plot: Box<dyn PlotObject>,
}

/// A piecewise composition of independently-defined plot objects.
pub struct PiecewisePlot {
    pieces: Vec<Piece>,
    style: PlotStyle,
}

impl PiecewisePlot {
    pub fn new() -> Self {
        Self {
            pieces: Vec::new(),
            style: PlotStyle::default(),
        }
    }

    pub fn add_piece(&mut self, domain: Domain, plot: impl PlotObject + 'static) {
        self.pieces.push(Piece {
            domain,
            plot: Box::new(plot),
        });
    }

    pub fn with_piece(mut self, domain: Domain, plot: impl PlotObject + 'static) -> Self {
        self.add_piece(domain, plot);
        self
    }
}

impl Default for PiecewisePlot {
    fn default() -> Self {
        Self::new()
    }
}

impl PlotObject for PiecewisePlot {
    fn coordinate_system(&self) -> CoordinateSystem {
        CoordinateSystem::Cartesian
    }

    fn natural_bounds(&self) -> Option<DataBounds> {
        let mut iter = self.pieces.iter();
        let first = iter.next()?;
        let mut x_min = *first.domain.x.start();
        let mut x_max = *first.domain.x.end();
        let mut y_min = *first.domain.y.start();
        let mut y_max = *first.domain.y.end();
        let mut z_min = *first.domain.z.start();
        let mut z_max = *first.domain.z.end();

        for piece in iter {
            x_min = x_min.min(*piece.domain.x.start());
            x_max = x_max.max(*piece.domain.x.end());
            y_min = y_min.min(*piece.domain.y.start());
            y_max = y_max.max(*piece.domain.y.end());
            z_min = z_min.min(*piece.domain.z.start());
            z_max = z_max.max(*piece.domain.z.end());
        }

        Some(DataBounds {
            x: x_min..=x_max,
            y: y_min..=y_max,
            z: z_min..=z_max,
        })
    }

    fn generate(&self, _domain: &Domain, resolution: Resolution) -> PlotGeometry {
        let components = self
            .pieces
            .iter()
            .map(|piece| PlotComponent {
                geometry: piece.plot.generate(&piece.domain, resolution),
                style: piece.plot.style().clone(),
            })
            .collect();
        PlotGeometry::Composite(components)
    }

    fn style(&self) -> &PlotStyle {
        &self.style
    }

    fn resolution(&self) -> Resolution {
        Resolution::default()
    }
}