semiflow 0.15.0-beta

no_std Rust PDE solver: evolution equations via Chernoff approximation of operator semigroups (heat, diffusion, Schrodinger)
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

//! Gate: `apply_into` allocates 0 heap bytes per step in steady state.
//!
//! ADR-0041 Wave 1 acceptance criterion AC-4:
//! After the first step (warm-up), subsequent steps must allocate 0 bytes when
//! `DiffusionChernoff` or `TruncatedExp4thDiffusionChernoff` is the inner kernel.
//!
//! Strategy: run `apply_into` once (warm-up), then measure allocations for a
//! second identical call with the same pre-warmed `ScratchPool`.
//!
//! The `allocation-counter` crate overrides the global allocator (test-only).

use allocation_counter::{self, AllocationInfo};
use semiflow::{
    scratch::ScratchPool, ChernoffFunction, ChernoffSemigroup, DiffusionChernoff, Grid1D, GridFn1D,
    TruncatedExp4thDiffusionChernoff,
};

fn a_half(_: f64) -> f64 {
    0.5
}

fn a_zero(_: f64) -> f64 {
    0.0
}

// ---------------------------------------------------------------------------
// DiffusionChernoff: 0 allocs after warm-up
// ---------------------------------------------------------------------------

#[test]
fn diffusion_zero_alloc_steady() {
    let n = 64;
    let grid = Grid1D::new(-4.0, 4.0, n).unwrap();
    let dx = grid.dx();
    let a0 = 0.5;
    let tau = 0.35 * dx * dx / a0;
    let op = DiffusionChernoff::new(a_half, a_zero, a_zero, a0, grid);
    let f0 = GridFn1D::from_fn(grid, |x| (-x * x).exp());

    // Warm-up: populates the scratch pool's free-list.
    let mut pool = ScratchPool::new();
    let mut dst = f0.zeroed_like();
    op.apply_into(tau, &f0, &mut dst, &mut pool).unwrap();

    // Steady-state measurement: pool already holds all needed capacity.
    let info: AllocationInfo = allocation_counter::measure(|| {
        op.apply_into(tau, &f0, &mut dst, &mut pool).unwrap();
    });

    assert_eq!(
        info.count_total, 0,
        "DiffusionChernoff::apply_into allocated {} time(s) in steady state (expected 0)",
        info.count_total
    );
}

// ---------------------------------------------------------------------------
// TruncatedExp4thDiffusionChernoff: 0 allocs after warm-up
// ---------------------------------------------------------------------------

#[test]
fn texp4_zero_alloc_steady() {
    let n = 64;
    let grid = Grid1D::new(-4.0, 4.0, n).unwrap();
    let dx = grid.dx();
    let a0 = 0.5;
    let tau = 0.25 * dx * dx / a0;
    let op = TruncatedExp4thDiffusionChernoff::new(a_half, a_zero, a_zero, a0, grid);
    let f0 = GridFn1D::from_fn(grid, |x| (-x * x).exp());

    let mut pool = ScratchPool::new();
    let mut dst = f0.zeroed_like();
    // Warm-up: take_vec × 4 → return_vec × 4 → 4 items in free-list.
    op.apply_into(tau, &f0, &mut dst, &mut pool).unwrap();

    let info: AllocationInfo = allocation_counter::measure(|| {
        op.apply_into(tau, &f0, &mut dst, &mut pool).unwrap();
    });

    assert_eq!(
        info.count_total, 0,
        "TruncatedExp4th::apply_into allocated {} time(s) in steady state (expected 0)",
        info.count_total
    );
}

// ---------------------------------------------------------------------------
// ChernoffSemigroup::evolve ping-pong: ≤ 2 allocs (setup only, 0 per step)
// ---------------------------------------------------------------------------

#[test]
fn chernoff_semigroup_ping_pong_allocs() {
    let n = 64;
    let grid = Grid1D::new(-4.0, 4.0, n).unwrap();
    let a0 = 0.5;
    let t = 0.01;
    let steps = 4_usize;
    #[allow(clippy::cast_precision_loss)]
    let tau = t / steps as f64;
    let dx = grid.dx();
    assert!(tau < 0.4 * dx * dx / a0, "CFL violated in test setup");

    let op = DiffusionChernoff::new(a_half, a_zero, a_zero, a0, grid);
    let semi = ChernoffSemigroup::new(op, steps).unwrap();
    let f0 = GridFn1D::from_fn(grid, |x| (-x * x).exp());

    // Measure allocations for a complete evolve call: the two ping-pong states and
    // the τ-keyed plan buffer (ADR-0215; it replaces the per-step ghost table and
    // is dropped with the local pool, so the free list never grows).
    let info: AllocationInfo = allocation_counter::measure(|| {
        let _result = semi.evolve(t, &f0).unwrap();
    });

    assert!(
        info.count_total <= 3,
        "ChernoffSemigroup::evolve allocated {} time(s) (expected ≤ 3: two states + plan)",
        info.count_total
    );
}