sva-samples 0.7.13

Buffers, frames, the f64 machine, the collapse of a law onto a grid, and the measured observations over one
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

// Concern: runs each addend of a sum on its own row and adds the planes | Non-concern: choosing the rows (plan.rs), running one (collapse.rs) | IO: (Vec<Plan>) -> (Buffer, Label)

use super::{AliasScore, Extent, plan::Plan};
use crate::buffer::Buffer;
use crate::error::CollapseError;
use crate::label::{Detail, Label, Source};
use crate::profile::Profile;

/// Each addend takes the row its own form names; the planes add.
pub fn added(
    parts: Vec<Plan>,
    rate: u32,
    extent: Extent,
    profile: &Profile,
    len: usize,
    score: AliasScore,
) -> Result<(Buffer, Label), CollapseError> {
    let mut collapsed = Vec::with_capacity(parts.len());
    for part in parts {
        collapsed.push(super::run(part, rate, extent, profile, len, score)?);
    }
    let width = collapsed
        .iter()
        .map(|(buffer, _)| buffer.width)
        .max()
        .expect("a sum row is built from its addends, and a sum has at least one");
    let mut planes = vec![vec![0.0; len]; width];
    for (buffer, _) in &collapsed {
        for (c, plane) in planes.iter_mut().enumerate() {
            // One component broadcasts into every lane of the sum.
            let lane = match buffer.width {
                1 => 0,
                _ => c,
            };
            if lane >= buffer.width {
                continue;
            }
            for (out, held) in plane.iter_mut().zip(buffer.plane(lane)) {
                *out += held;
            }
        }
    }
    let labels: Vec<Label> = collapsed.into_iter().map(|(_, label)| label).collect();
    Ok((
        Buffer {
            rate,
            start: extent.start,
            width,
            planes,
        },
        label(&labels, profile, rate),
    ))
}

/// One measured addend measures the whole.
fn label(parts: &[Label], profile: &Profile, rate: u32) -> Label {
    let source = match parts.iter().all(|p| p.source == Source::Exact) {
        true => Source::Exact,
        false => Source::Measured,
    };
    Label::new(
        source,
        profile.name,
        rate,
        Detail::Added {
            parts: parts.iter().map(|p| p.detail.clone()).collect(),
        },
    )
}