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

#![cfg(not(feature = "complex"))]
use kryst::algebra::prelude::*;
use kryst::parallel::{Comm, NoComm, UniverseComm};

#[test]
fn serial_nonblocking_request_lifecycle() {
    // Use serial backend and exercise irecv/isend/wait_all with the shared Request type
    let comm = UniverseComm::NoComm(NoComm);

    let zero = S::zero().real();
    let one = S::one().real();
    let mut recv = vec![zero; 4];
    let send = vec![one; 4];

    let mut reqs: Vec<<UniverseComm as Comm>::Request<'_>> = Vec::new();
    reqs.push(comm.irecv_from(&mut recv[..], 0));
    reqs.push(comm.isend_to(&send[..], 0));
    comm.wait_all(&mut reqs);
    // Ensure requests are dropped before reading recv buffer
    drop(reqs);

    // In serial, no data moves; the important part is type/lifetime integration.
    assert_eq!(recv, vec![zero; 4]);
}