aeon-tk 0.3.3

Toolkit for running finite difference simulations with adaptive mesh refinement
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
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54

use std::cell::UnsafeCell;

use bumpalo::Bump;
use thread_local::ThreadLocal;

/// A per-thread store that contains memory pools and caches for multi-threaded
/// workloads.
#[derive(Debug, Default)]
pub struct MeshStore {
    /// Linear allocator for each thread.
    arena: Bump,
}

impl MeshStore {
    /// Allocates scratch data for use by the current thread.
    pub fn scratch<T: Default>(&self, len: usize) -> &mut [T] {
        self.arena.alloc_slice_fill_default(len)
    }

    /// Resets memory cached in mesh store.
    pub fn reset(&mut self) {
        self.arena.reset();
    }
}

/// A data type which stores a pool of thread local variables.
/// Thus allowing each thread to access one copy of `T` mutably.
#[derive(Debug)]
pub struct UnsafeThreadCache<T: Send> {
    pool: ThreadLocal<UnsafeCell<T>>,
}

impl<T: Send> UnsafeThreadCache<T> {
    /// Constructs an empty thread cache.
    pub fn new() -> Self {
        Self::default()
    }
}

impl<T: Send + Default> UnsafeThreadCache<T> {
    /// Retrieves the object `T` associated with this thread, initializing the
    /// default value in place if this has not already been done.
    pub unsafe fn get_or_default(&self) -> &mut T {
        unsafe { &mut *self.pool.get_or_default().get() }
    }
}

impl<T: Send> Default for UnsafeThreadCache<T> {
    fn default() -> Self {
        Self {
            pool: ThreadLocal::default(),
        }
    }
}