Struct LevelSetField 

pub struct LevelSetField {
    pub nx: usize,
    pub ny: usize,
    pub nz: usize,
    pub dx: f64,
    pub origin: [f64; 3],
    pub data: Vec<f64>,
}

Signed distance function sampled on a uniform 3-D Cartesian grid.

The grid has nx × ny × nz cells. Values are stored in row-major order as data[iz * ny * nx + iy * nx + ix]. Positive values are outside the surface; negative values are inside.

Fields

nx: usize

Number of grid points along x.

ny: usize

Number of grid points along y.

nz: usize

Number of grid points along z.

dx: f64

Grid spacing (uniform in all directions).

origin: [f64; 3]

Origin of the grid (coordinates of node [0,0,0]).

data: Vec<f64>

Flat array of signed distance values, length nx * ny * nz.

Implementations

impl LevelSetField

pub fn new(nx: usize, ny: usize, nz: usize, dx: f64, origin: [f64; 3]) -> Self

Create a new LevelSetField filled with f64::MAX (uninitialized).

Arguments

• nx, ny, nz – grid dimensions
• dx – uniform grid spacing
• origin – world-space coordinates of node (0,0,0)

pub fn idx(&self, ix: usize, iy: usize, iz: usize) -> usize

Linear index for grid node (ix, iy, iz).

pub fn node_pos(&self, ix: usize, iy: usize, iz: usize) -> [f64; 3]

World-space coordinates of grid node (ix, iy, iz).

pub fn get(&self, ix: usize, iy: usize, iz: usize) -> f64

Get the level-set value at grid node (ix, iy, iz).

Returns f64::MAX for out-of-bounds nodes.

pub fn set(&mut self, ix: usize, iy: usize, iz: usize, value: f64)

Set the level-set value at grid node (ix, iy, iz).

Silently ignores out-of-bounds writes.

pub fn interpolate(&self, p: [f64; 3]) -> f64

Tri-linear interpolation of the level-set field at world position p.

Returns f64::MAX for positions outside the grid.

pub fn gradient(&self, ix: usize, iy: usize, iz: usize) -> [f64; 3]

Numerical gradient of the level-set field at grid node (ix, iy, iz).

Uses central differences where possible, one-sided differences at boundaries.

pub fn mean_curvature(&self, ix: usize, iy: usize, iz: usize) -> f64

Mean curvature κ = div(∇φ / |∇φ|) at grid node (ix, iy, iz).

Computed via second-order finite differences.

pub fn init_sphere(&mut self, center: [f64; 3], radius: f64)

Initialize the level-set as a sphere SDF.

Sets φ(x) = |x − center| − radius for all nodes.

pub fn init_box(&mut self, lo: [f64; 3], hi: [f64; 3])

Initialize the level-set as an axis-aligned box SDF.

Sets φ(x) to the signed distance to the box [lo, hi].

pub fn count_inside(&self) -> usize

Count nodes with φ ≤ 0 (inside or on the interface).

pub fn min_value(&self) -> f64

Minimum value of φ over all nodes.

pub fn max_value(&self) -> f64

Maximum value of φ over all nodes.