Struct LevelSetTopOpt 

pub struct LevelSetTopOpt {
    pub grid: GridSize,
    pub phi: Vec<f64>,
    pub volume_fraction: f64,
    pub dt: f64,
    pub material: MaterialParams,
    pub compliance_history: Vec<f64>,
}

Level-set topology optimization using the Hamilton-Jacobi equation.

The structural boundary is implicitly represented as the zero contour of the level-set function φ. The Hamilton-Jacobi PDE governs its evolution: ∂φ/∂t + V_n |∇φ| = 0.

Fields

grid: GridSize

Grid size.

phi: Vec<f64>

Level-set function values (one per grid node, length = (nx+1)*(ny+1)).

volume_fraction: f64

Volume fraction target.

dt: f64

Time step size for Hamilton-Jacobi update.

material: MaterialParams

Material parameters.

compliance_history: Vec<f64>

Compliance history.

Implementations

impl LevelSetTopOpt

pub fn new( grid: GridSize, volume_fraction: f64, dt: f64, material: MaterialParams, ) -> Self

Creates a new LevelSetTopOpt with a circular initial inclusion.

pub fn element_density(&self) -> Vec<f64>

Extracts element density from level-set (Heaviside of nodal averages).

pub fn mock_compliance(&self, velocity: &[f64]) -> f64

Computes mock compliance given a velocity field.

pub fn hj_update(&mut self, velocity: &[f64])

Performs one Hamilton-Jacobi update step.

φ_new = φ - dt * V_n * |∇φ| (using upwind difference for |∇φ|).

pub fn reinitialize(&mut self, n_iter: usize)

Re-initializes phi as a signed distance function using fast-marching approximation.

pub fn run_mock(&mut self, n_iter: usize)

Runs mock level-set optimization.

Trait Implementations

impl Clone for LevelSetTopOpt

fn clone(&self) -> LevelSetTopOpt

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl Debug for LevelSetTopOpt

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

Formats the value using the given formatter.