nalgebra/src/sparse/cs_utils.rs
2021-04-11 13:57:54 +02:00

19 lines
375 B
Rust

use crate::allocator::Allocator;
use crate::{DefaultAllocator, Dim, OVector};
pub fn cumsum<D: Dim>(a: &mut OVector<usize, D>, b: &mut OVector<usize, D>) -> usize
where
DefaultAllocator: Allocator<usize, D>,
{
assert!(a.len() == b.len());
let mut sum = 0;
for i in 0..a.len() {
b[i] = sum;
sum += a[i];
a[i] = b[i];
}
sum
}