Struct MorseComplex 

pub struct MorseComplex {
    pub critical_points: Vec<(Vec<f64>, f64)>,
}

A Morse complex built from a collection of critical points of a Morse function.

Each critical point is stored as (coordinates, function_value).

Fields

critical_points: Vec<(Vec<f64>, f64)>

Critical points: (position_coords, function_value).

Implementations

impl MorseComplex

pub fn new(critical_points: Vec<(Vec<f64>, f64)>) -> Self

Create a new MorseComplex with the given critical points.

pub fn morse_index(&self, pt: &[f64]) -> usize

Compute the Morse index of a critical point by finite-difference Hessian.

The Morse index is the number of negative eigenvalues of the Hessian. This implementation approximates the Hessian via the ordering of surrounding critical points.

pub fn gradient_flow(&self) -> Vec<&(Vec<f64>, f64)>

Simulate gradient flow: sort critical points by function value.

Returns the critical points in ascending order of function value, representing the flow from minima to maxima.

pub fn pair_critical_points(&self) -> Vec<(usize, usize)>

Pair critical points using Morse cancellation heuristic.

Returns pairs (min_idx, max_idx) of critical points that can be cancelled (adjacent in function value with index difference 1).

Trait Implementations

impl Clone for MorseComplex

fn clone(&self) -> MorseComplex

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl Debug for MorseComplex

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

Formats the value using the given formatter.