kryst 3.2.1

Krylov subspace and preconditioned iterative solvers for dense and sparse linear systems, with shared and distributed memory parallelism.
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

#![cfg(all(feature = "backend-faer", not(feature = "complex")))]
use std::sync::Arc;

use faer::Mat;
use kryst::algebra::prelude::*;
use kryst::matrix::backend::DefaultBackend;
use kryst::matrix::format::AsFormat;
use kryst::matrix::op::DenseOp;

#[test]
fn denseop_csc_cache_includes_values_id() {
    // 2x2 dense with trivial values
    let m = Mat::<R>::from_fn(2, 2, |i, j| {
        if i == j {
            S::one().real()
        } else {
            R::default()
        }
    });
    let dop = DenseOp::<f64>::new(Arc::new(m));
    // First conversion
    let c1 = <DenseOp<f64> as AsFormat<f64, DefaultBackend>>::to_csc_cached(&dop, R::default());
    // Mark numeric change (simulate in-place edits)
    dop.mark_values_changed();
    // Second conversion should not hit the same cache entry
    let c2 = <DenseOp<f64> as AsFormat<f64, DefaultBackend>>::to_csc_cached(&dop, R::default());

    // The arcs should be distinct instances because key includes ValuesId
    let p1 = Arc::as_ptr(&c1) as usize;
    let p2 = Arc::as_ptr(&c2) as usize;
    assert_ne!(p1, p2, "CSC cache entry not invalidated on values change");
}

#[cfg(not(feature = "mat-values-fingerprint"))]
#[test]
fn raw_mat_csc_cache_is_stable_without_values_tracking() {
    // Raw Mat has values_id == 0 -> same key if structure and drop_tol stay the same
    let m = Mat::<R>::from_fn(2, 2, |i, j| {
        if i == j {
            S::one().real()
        } else {
            R::default()
        }
    });
    let c1 = <Mat<R> as AsFormat<f64, DefaultBackend>>::to_csc_cached(&m, R::default());
    let c2 = <Mat<R> as AsFormat<f64, DefaultBackend>>::to_csc_cached(&m, R::default());
    let p1 = Arc::as_ptr(&c1) as usize;
    let p2 = Arc::as_ptr(&c2) as usize;
    assert_eq!(
        p1, p2,
        "Raw Mat should reuse the same CSC when values_id is unknown (0)"
    );
}

#[test]
fn denseop_csc_cache_includes_structure_id() {
    let m = Mat::<R>::from_fn(2, 2, |i, j| {
        if i == j {
            S::one().real()
        } else {
            R::default()
        }
    });
    let dop = DenseOp::<f64>::new(Arc::new(m));
    let c1 = <DenseOp<f64> as AsFormat<f64, DefaultBackend>>::to_csc_cached(&dop, R::default());
    dop.mark_structure_changed();
    let c2 = <DenseOp<f64> as AsFormat<f64, DefaultBackend>>::to_csc_cached(&dop, R::default());
    let p1 = Arc::as_ptr(&c1) as usize;
    let p2 = Arc::as_ptr(&c2) as usize;
    assert_ne!(p1, p2, "CSC cache entry not invalidated on structure change");
}