Struct Tower2 

pub struct Tower2<const K: usize> {
    pub v: f64,
    pub g: [f64; K],
    pub h: [[f64; K]; K],
}

Truncated SECOND-order multivariate Taylor scalar in K variables.

This is the value/gradient/Hessian-only sibling of Tower4. Every channel it carries (v, g, h) is computed by the SAME formulas Tower4 uses for those orders, so for any program written over both towers the order-≤2 outputs are bit-identical: the order-2 Leibniz and Faà-di-Bruno terms read only the order-≤2 channels of their inputs (see Tower4::mul / Tower4::compose_unary — out.h never touches t3 or t4), so dropping the third/fourth tensors cannot perturb the value, gradient, or Hessian.

It exists purely for performance: an inner Newton step (and the value-only ρ-homotopy pre-warm) needs at most curvature, never the outer-κ/ψ third/fourth derivatives. Evaluating a row likelihood over Tower2 skips the K⁴ fourth-tensor product/composition arithmetic that dominates the cold marginal-slope fit, while returning the exact same (v, g, h).

Fields

v: f64

Value ℓ.

g: [f64; K]

Gradient ∂ℓ/∂p_a.

h: [[f64; K]; K]

Hessian ∂²ℓ/∂p_a∂p_b (symmetric).

Implementations

impl<const K: usize> Tower2<K>

pub fn zero() -> Self

The additive identity.

pub fn constant(c: f64) -> Self

A constant: value c, all derivatives zero.

pub fn variable(value: f64, idx: usize) -> Self

The seeded variable p_idx with current value value: unit first derivative in slot idx, zero elsewhere and above.

pub fn mul(&self, o: &Self) -> Self

Exact truncated (order ≤ 2) Leibniz product. The v/g/h channels match Tower4::mul term-for-term.

pub fn compose_unary(&self, d: [f64; 3]) -> Self

Exact (order ≤ 2) multivariate Faà di Bruno composition f ∘ self.

d = [f(u), f′(u), f″(u)] evaluated at u = self.v. The v/g/h channels match Tower4::compose_unary term-for-term (which uses only d[0..=2] for those orders), so this is a strict truncation, not an approximation. The full-order [f64; 5] derivative stacks the families already produce can be passed by slicing their first three entries.

pub fn scale(&self, s: f64) -> Self

Multiply every channel by a plain scalar.

pub fn exp(&self) -> Self

e^self.

pub fn sqrt(&self) -> Self

√self. Caller guarantees positivity.

Trait Implementations

impl<const K: usize> Add for Tower2<K>

type Output = Tower2<K>

The resulting type after applying the + operator.

fn add(self, o: Self) -> Self

Performs the + operation.

impl<const K: usize> Add<f64> for Tower2<K>

type Output = Tower2<K>

The resulting type after applying the + operator.

fn add(self, c: f64) -> Self

Performs the + operation.

impl<const K: usize> Clone for Tower2<K>

fn clone(&self) -> Tower2<K>

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl<const K: usize> Copy for Tower2<K>

impl<const K: usize> Debug for Tower2<K>

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

Formats the value using the given formatter.

impl<const K: usize> Mul for Tower2<K>

type Output = Tower2<K>

The resulting type after applying the * operator.

fn mul(self, o: Self) -> Self

Performs the * operation.

impl<const K: usize> Mul<f64> for Tower2<K>

type Output = Tower2<K>

The resulting type after applying the * operator.

fn mul(self, c: f64) -> Self

Performs the * operation.