numeris 0.5.1

Pure-Rust numerical algorithms library — high performance with SIMD support while also supporting no-std for embedded and WASM targets.
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

//! Cholesky rank-1 update and downdate.
//!
//! Delegates to the generic implementations in [`crate::linalg`].

use crate::linalg::{cholesky_rank1_update, cholesky_rank1_downdate};
use crate::matrix::vector::Vector;
use crate::traits::FloatScalar;
use crate::Matrix;

use super::EstimateError;

/// In-place Cholesky rank-1 update (`sign = +1`) or downdate (`sign = -1`).
///
/// `l` must be lower-triangular on entry (upper triangle is ignored/untouched).
/// `v` is used as workspace and is modified in place.
///
/// O(N²), no allocation.
#[allow(dead_code)]
pub(crate) fn cholupdate<T: FloatScalar, const N: usize>(
    l: &mut Matrix<T, N, N>,
    v: &mut Vector<T, N>,
    sign: T,
) -> Result<(), EstimateError> {
    let result = if sign >= T::zero() {
        cholesky_rank1_update(l, v.as_mut_slice())
    } else {
        cholesky_rank1_downdate(l, v.as_mut_slice())
    };
    result.map_err(|_| EstimateError::CholdowndateFailed)
}