Struct RayRestoration 

pub struct RayRestoration {
    pub block: usize,
    pub rho_first: usize,
    pub rho_count: usize,
    pub log_strength_ratio: f64,
    pub likelihood_slope: f64,
    pub penalty_slope: f64,
    pub block_step_inf: f64,
}

The block-level reading of a ray the joint Newton was descending when it stopped: along the last accepted step δ, the likelihood term slopes down by likelihood_slope = ∇(−ℓ)·δ < 0 while block block’s penalty slopes up by only penalty_slope = (λ_b S_b β)·δ > 0. The penalized objective is stationary along δ at the strength ratio r = −likelihood_slope / penalty_slope > 1, so raising every log strength of that block by log_strength_ratio = ln r closes the ray at the iterate the solve stopped on. The ratio is read off the two slopes the solve already had; nothing here is a step size.

Fields

block: usize

Index of the parameter block carrying the ray.

rho_first: usize

The outer ρ coordinates of that block’s penalties: rho_count of them, contiguous from rho_first (ρ is laid out block by block).

rho_count: usize

log_strength_ratio: f64

ln r, the amount every one of rho_indices has to rise.

likelihood_slope: f64

∇(−ℓ)·δ along the accepted step (negative: the likelihood descends).

penalty_slope: f64

(λ_b S_b β)·δ along the accepted step (positive: the penalty resists).

block_step_inf: f64

‖δ_b‖∞, the block’s share of the accepted step.

Implementations

impl RayRestoration

pub fn rho_indices(&self) -> Range<usize>

The outer ρ coordinates this restoration raises.

Trait Implementations

impl Clone for RayRestoration

fn clone(&self) -> RayRestoration

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl Copy for RayRestoration

impl Debug for RayRestoration

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

Formats the value using the given formatter.

impl Display for RayRestoration

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

Formats the value using the given formatter.

impl PartialEq for RayRestoration

fn eq(&self, other: &RayRestoration) -> bool

Equality operator ==.

fn ne(&self, other: &Rhs) -> bool

Inequality operator !=.

impl StructuralPartialEq for RayRestoration