mesh-sieve 4.0.2

Modular, high-performance Rust library for mesh and data management, designed for scientific computing and PDE codes.
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
31
32
33
34
35
36

use mesh_sieve::algs::{communicator::RayonComm, completion::complete_section};
use mesh_sieve::data::{atlas::Atlas, section::Section, storage::VecStorage};
use mesh_sieve::overlap::delta::CopyDelta;
use mesh_sieve::overlap::overlap::Overlap;
use mesh_sieve::topology::point::PointId;

#[test]
fn ghost_update_self() {
    // This test simulates a “rank” that owns P0 and also ghosts P0:
    // i.e. it adds a self‐link so that complete_section does one send
    // and one receive, but both target the same rank (yourself).

    let mut ovlp = Overlap::default();
    let p0 = PointId::new(1).unwrap();

    // Add a self‐overlap: owner and ghost are the same rank (0).
    ovlp.add_link(p0, /*remote_rank=*/ 0, p0);

    // Build a section with one value = 42
    let mut atlas = Atlas::default();
    atlas
        .try_insert(p0, 1)
        .expect("Failed to insert point into atlas");
    let mut sec = Section::<u32, VecStorage<u32>>::new(atlas);
    sec.try_set(p0, &[42]).expect("Failed to set section value");

    let comm = RayonComm::new(0, 1);

    // Should complete without deadlock and leave the value intact.
    let _ = complete_section::<u32, VecStorage<u32>, CopyDelta, _>(&mut sec, &ovlp, &comm, 0);

    assert_eq!(
        sec.try_restrict(p0).expect("Failed to restrict section")[0],
        42
    );
}