conspire 0.7.0

The Rust interface to conspire.
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

#[cfg(test)]
mod test;

use crate::{
    geometry::mesh::{Mesh, differential::laplace::Weighting},
    math::Scalar,
};

impl<const D: usize> Mesh<D> {
    pub fn taubin_smooth(
        &mut self,
        iterations: usize,
        pass_band: Scalar,
        scale: Scalar,
        weighting: Weighting,
        preserve_boundary: bool,
        preserve_interfaces: bool,
    ) {
        let scale_deflate = scale;
        let scale_inflate = scale / (pass_band * scale - 1.0);
        let adjacency = (preserve_boundary || preserve_interfaces)
            .then(|| self.constrained_adjacency(preserve_boundary, preserve_interfaces));
        for iteration in 0..iterations {
            let laplacian = match &adjacency {
                Some(adjacency) => self.laplacian_over(adjacency, weighting),
                None => self.laplacian(weighting),
            };
            let scale = if scale_inflate < 0.0 && iteration % 2 == 1 {
                scale_inflate
            } else {
                scale_deflate
            };
            self.coordinates
                .iter_mut()
                .zip(laplacian)
                .for_each(|(x, u)| *x -= u * scale)
        }
    }
}