dspk-synapse 0.0.5

Electrophysiology for dspk: recording and sorting formats, probes, spike detection, waveforms, features, clustering, localization, quality metrics, and the shared parts of the sorters
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

//! Kernels of [`crate::sorter::components::detection::matched`].

use cubecl::prelude::*;
use dspk_core::compute::element_position;

/// One unit per element `(c, i)` of `out` (`[channels, samples]`): `Σₖ proto[k] · x[c, i + k]`
/// for `i + len ≤ samples` (the correlation's valid part), 0 past it.
#[cube(launch)]
pub fn correlate_prototype_kernel<F: Float>(x: &[F], proto: &[F], out: &mut [F], samples: u32, len: u32, total: u32) {
    let e = element_position() as u32;
    if e < total {
        let c = e / samples;
        let i = e - c * samples;
        let mut acc = F::new(0.0f32);
        if i + len <= samples {
            let base = c * samples + i;
            let mut k = 0u32;
            while k < len {
                acc += proto[k as usize] * x[(base + k) as usize];
                k += 1u32;
            }
        }
        out[e as usize] = acc;
    }
}

/// One unit per element `(r, i)` of `out` (`[rows, samples]`): `Σ w · corr[c, i]` over row `r`'s
/// channels (CSR: `offsets[r]..offsets[r + 1]` of `cols` / `weights`).
#[cube(launch)]
#[allow(clippy::too_many_arguments)]
pub fn spatial_sum_kernel<F: Float>(
    corr: &[F],
    offsets: &[u32],
    cols: &[u32],
    weights: &[F],
    out: &mut [F],
    samples: u32,
    total: u32,
) {
    let e = element_position() as u32;
    if e < total {
        let r = e / samples;
        let i = e - r * samples;
        let mut acc = F::new(0.0f32);
        let mut k = offsets[r as usize];
        let end = offsets[(r + 1u32) as usize];
        while k < end {
            acc += weights[k as usize] * corr[(cols[k as usize] * samples + i) as usize];
            k += 1u32;
        }
        out[e as usize] = acc;
    }
}

/// `out[i] = x[channels[i] · samples + at[i]]`: values of a `[channels, samples]` buffer at points.
#[cube(launch)]
pub fn gather_points_kernel<F: Float>(x: &[F], channels: &[u32], at: &[u32], out: &mut [F], samples: u32, n: u32) {
    let i = element_position() as u32;
    if i < n {
        out[i as usize] = x[(channels[i as usize] * samples + at[i as usize]) as usize];
    }
}