Optimizers_PADM.h

//
// Created by Henri on 22/04/2026.
//
 
#ifndef IDOL_OPTIMIZERS_PADM_H
#define IDOL_OPTIMIZERS_PADM_H
 
#include "idol/bilevel/modeling/Description.h"
#include "idol/general/optimizers/Algorithm.h"
 
namespace idol::Optimizers::Bilevel {
    class PADM;
}
 
class idol::Optimizers::Bilevel::PADM : public Algorithm {
    const idol::Bilevel::Description& m_description;
    std::unique_ptr<OptimizerFactory> m_deterministic_optimizer;
    std::unique_ptr<Model> m_single_level_reformulation;
    Annotation<unsigned int> m_decomposition;
    Annotation<double> m_initial_penalty_parameters;
 
    void throw_if_no_deterministic_model() const;
 
    void build_decomposition();
public:
    PADM(const Model& t_parent,
        const idol::Bilevel::Description& t_description,
        const OptimizerFactory& t_deterministic_optimizer);
 
    [[nodiscard]] std::string name() const override;
 
    [[nodiscard]] SolutionStatus get_status() const override;
 
    [[nodiscard]] SolutionReason get_reason() const override;
 
    [[nodiscard]] double get_best_obj() const override;
 
    [[nodiscard]] double get_best_bound() const override;
 
    [[nodiscard]] double get_var_primal(const Var &t_var) const override;
 
    [[nodiscard]] double get_var_reduced_cost(const Var &t_var) const override;
 
    [[nodiscard]] double get_var_ray(const Var &t_var) const override;
 
    [[nodiscard]] double get_ctr_dual(const Ctr &t_ctr) const override;
 
    [[nodiscard]] double get_ctr_farkas(const Ctr &t_ctr) const override;
 
    [[nodiscard]] unsigned int get_n_solutions() const override;
 
    [[nodiscard]] unsigned int get_solution_index() const override;
 
protected:
    void add(const Var &t_var) override;
 
    void add(const Ctr &t_ctr) override;
 
    void add(const QCtr &t_ctr) override;
 
    void remove(const Var &t_var) override;
 
    void remove(const Ctr &t_ctr) override;
 
    void remove(const QCtr &t_ctr) override;
 
    void update() override;
 
    void write(const std::string &t_name) override;
 
    void hook_optimize() override;
 
    void set_solution_index(unsigned int t_index) override;
 
    void update_obj_sense() override;
 
    void update_obj() override;
 
    void update_rhs() override;
 
    void update_obj_constant() override;
 
    void update_mat_coeff(const Ctr &t_ctr, const Var &t_var) override;
 
    void update_ctr_type(const Ctr &t_ctr) override;
 
    void update_ctr_rhs(const Ctr &t_ctr) override;
 
    void update_var_type(const Var &t_var) override;
 
    void update_var_lb(const Var &t_var) override;
 
    void update_var_ub(const Var &t_var) override;
 
    void update_var_obj(const Var &t_var) override;
};
 
#endif //IDOL_OPTIMIZERS_PADM_H