[Dev] Expose the duals of the JuMP path, boundary and box constraints, aligned with signs.md

Ouverte
  • Quelque chose ne va pas ?
#430 0 commentaires 0 réactions 1 personne assignée Voir sur GitHub

Les mainteneurs répondent en général sous 1 jour

@ocots y travaille déjà.

Depuis le 8/11/2026.

Évaluation

Cette issue n'a pas encore été évaluée.

Description

Context

#326 aligns the sign of the JuMP costate with OptimalControl (ct-conventions rules/signs.md, SGN1–SGN6, v0.11.0): the normalized multiplier of a constraint row is -JuMP.dual(ref), for a Min and a Max model alike. Today only the dynamics constraints are exposed — their duals are read as the costate. The duals of the other constraints the JuMP models create exist but are never read:

  • path constraints: lower <= g(x, u, v, t) <= upper rows, built at every node;
  • boundary constraints: lower <= g(x0, xf, v) <= upper rows;
  • box bounds on the state, control and variable components (LowerBoundRef / UpperBoundRef).

All of them are concatenated in model[:constraints] (dynamics, path, boundary); build.jl would have to keep the refs per group to expose them.

What to do (later)

  • Expose the duals of the path and boundary constraints, and of the box bounds, normalized as the costate (-JuMP.dual).
  • Conventions to apply (rules/signs.md v0.11.0):
    • SGN3: a two-sided row l ≤ c ≤ u carries one signed multiplier y = μU − μL (≥ 0 upper active, ≤ 0 lower active) — a JuMP Interval constraint gives one signed dual, mapped by -dual;
    • SGN2: for a box bound or a single-bounded row, the dual is the multiplier of that side, ≥ 0 — i.e. -dual(LowerBoundRef) / -dual(UpperBoundRef) give the nonnegative per-side values;
    • SGN5: the same constraint declared as a path constraint or as a box must have the same dual — the comparison tests should check it;
    • the scale of a time-dependent dual (density vs raw multiplier) is free in SGN5 — CTDirect divides the path multipliers by h; choose the same here for consistency.
  • Follow the ecosystem design for reading duals by constraint label (dual(sol, model, label)): CTModels.jl#105, CTFlows.jl#314 and discussion #312 — this issue waits for it, as #326 did for the costate.

Tests

Extend the sign tests of #326 (two layers, SGN6): the normalized duals must satisfy ∇F + Jᵀy = 0 on hand-written JuMP problems with known multipliers — active lower and upper sides of a two-sided row, bounds of variables, min f vs max -f (SGN1) — and match the OptimalControl duals of the same problem where both are exposed.

Related

#326, #170, control-toolbox/CTDirect.jl#405, control-toolbox/CTModels.jl#105.

Langage dominant
Julia
Étoiles
13
Forks
1
Merge moyen
3 h 33 min
PR mergées (30 j)
55

Préparer son environnement

Par où commencer

  1. Lisez l'issue en entier, puis le guide de contribution du projet.
  2. Signalez en commentaire que vous la prenez — cela évite que deux personnes fassent le même travail.
  3. Forkez le dépôt et travaillez sur une branche.
  4. Ouvrez une pull request qui référence le numéro de l'issue.

Autres issues de control-toolbox/OptimalControlProblems.jl

Toutes les issues de control-toolbox/OptimalControlProblems.jl

Issues similaires

Plus d'issues Julia

Recevez les nouvelles issues par e-mail

Un résumé court des issues GitHub adaptées aux débutants.