Struct ParallelAssembler 

pub struct ParallelAssembler {
    pub ndofs: usize,
}

Assembles per-element stiffness contributions into a global CSR matrix.

Thread safety

Elements are processed in parallel (Rayon). Each non-zero entry’s position in the CSR values array is pre-determined from the sparsity pattern, so scatter is lock-free at the element level. A per-row Mutex guards concurrent writes to overlapping rows (e.g. when elements share nodes on different threads).

Fields

ndofs: usize

Total number of global DOFs.

Implementations

impl ParallelAssembler

pub fn new(ndofs: usize) -> Self

Create an assembler for a problem with ndofs degrees of freedom.

pub fn assemble(&self, tasks: &[AssemblyTask]) -> CsrMatrix

Assemble all element stiffness matrices into a global CSR matrix.

The sparsity pattern is computed from the union of all element connectivity patterns. Values are accumulated atomically using per-row mutexes.

Panics

Panics if any global_dof index exceeds ndofs.

pub fn assemble_rhs( &self, element_dofs: &[Vec<usize>], element_forces: &[Vec<f64>], ) -> Vec<f64>

Assemble a load vector (right-hand side) from per-element force vectors.

element_forces[e] must have the same DOF count as element_dofs[e].

pub fn assemble_colored( &self, tasks: &[AssemblyTask], element_dofs: &[Vec<usize>], ) -> CsrMatrix

Assemble using element graph coloring for better cache locality.

Computes a graph coloring of the elements (no two elements of the same color share a DOF) and assembles each color group in parallel. The result is identical to ParallelAssembler::assemble.

Trait Implementations

impl Debug for ParallelAssembler

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.

impl Default for ParallelAssembler

fn default() -> ParallelAssembler

Returns the “default value” for a type.