Struct NonLocalRegularization 

pub struct NonLocalRegularization {
    pub lc: f64,
}

Non-local integral averaging for damage regularization.

Computes a weighted average of a local field (e.g., equivalent strain) over a neighborhood defined by a characteristic length l_c.

ε̄(x) = Σ w(|x - xᵢ|) * ε(xᵢ) / Σ w(|x - xᵢ|)

where w(r) = exp(-r² / (2*l_c²)) is the Gaussian weight.

Fields

lc: f64

Characteristic (internal) length l_c (m).

Implementations

impl NonLocalRegularization

pub fn new(lc: f64) -> Self

Create a new non-local regularization with characteristic length l_c.

pub fn weight(&self, distance: f64) -> f64

Gaussian weight function: w(r) = exp(-r² / (2*l_c²)).

pub fn weight_bell(&self, distance: f64) -> f64

Bell-shaped (polynomial) weight function (compact support at r = l_c).

w(r) = (1 - (r/l_c)²)² for r < l_c, 0 otherwise

pub fn nonlocal_average( &self, point: &[f64; 3], positions: &[[f64; 3]], values: &[f64], ) -> f64

Compute non-local average of a field at a given point.

Arguments

• point - Position at which to evaluate the non-local average
• positions - Positions of all integration points
• values - Local field values at each integration point

pub fn nonlocal_average_all( &self, positions: &[[f64; 3]], values: &[f64], ) -> Vec<f64>

Compute non-local averages for all points simultaneously.

Returns a vector of non-local averages, one per point.

Trait Implementations

impl Clone for NonLocalRegularization

fn clone(&self) -> NonLocalRegularization

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Debug for NonLocalRegularization

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

Formats the value using the given formatter.