kryst 4.0.3

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

//! Legacy AMG-focused distributed CSR (diag/off splitting).
//!
//! NOTE: `ParCsrMatrix` is maintained for backward compatibility but is
//! considered secondary to `matrix::DistCsrOp`. New distributed code should
//! prefer `DistCsrOp`, and ParCsrMatrix will be gradually reworked to build on
//! the same halo infrastructure.
//!
//! Migration plan (active):
//! 1. `ParCsrMatrix` matvec now delegates to an internal canonical
//!    `DistCsrOp` adapter.
//! 2. Halo exchange and ownership semantics are unified through
//!    `matrix::dist::halo::HaloPlan` used by `DistCsrOp`.
//! 3. Solver/PC entry points should consume `LinOp::dist_layout()` /
//!    `LinOp::comm()` metadata; `ParCsrOp` forwards these from the canonical
//!    operator.
//! 4. Legacy diag/off data structures remain only as compatibility adapters and
//!    are marked deprecated (since 1.1.0) with planned removal after
//!    2026-12-31 in favor of canonical conversion helpers.
//! 5. Test parity gates validate that legacy `ParCsrMatrix` SpMV behavior
//!    matches canonical distributed SpMV before removing legacy routes.

pub type Global = u64;
pub type Local = usize;

pub mod builder;
pub mod halo;
pub mod mat;

pub use halo::HaloPlan;
pub use mat::{ParCsrMatrix, ParCsrOp};