pounce-algorithm 0.4.0

Algorithm-side core for POUNCE (port of Ipopt's src/Algorithm/): IteratesVector, IpoptData, CalculatedQuantities, KKT solvers, line search, mu update, conv check, initializer, IpoptAlg main loop, AlgBuilder.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39

//! SQP iterate state — `(x, λ_g, λ_x)` plus the working set
//! carried across QP subproblem solves for warm-starting (the
//! §5/§6 design-note contract).
//!
//! Distinct from [`crate::iterates_vector::IteratesVector`] (which
//! is IPM-shaped with slacks and barrier multipliers). The SQP
//! iterate is simpler: primal `x`, constraint multipliers
//! `λ_g`, packed bound multipliers `λ_x = z_l − z_u`, and the
//! `pounce_qp::WorkingSet` from the previous QP solve.

use pounce_common::Number;
use pounce_qp::WorkingSet;

#[derive(Debug, Clone)]
pub struct SqpIterates {
    pub x: Vec<Number>,
    pub lambda_g: Vec<Number>,
    pub lambda_x: Vec<Number>,
    pub working: Option<WorkingSet>,
}

impl SqpIterates {
    pub fn cold(n: usize, m: usize) -> Self {
        Self {
            x: vec![0.0; n],
            lambda_g: vec![0.0; m],
            lambda_x: vec![0.0; n],
            working: None,
        }
    }

    pub fn n(&self) -> usize {
        self.x.len()
    }

    pub fn m(&self) -> usize {
        self.lambda_g.len()
    }
}