scpgen.hpp

/*
 *    This file is part of CasADi.
 *
 *    CasADi -- A symbolic framework for dynamic optimization.
 *    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl,
 *                            KU Leuven. All rights reserved.
 *    Copyright (C) 2011-2014 Greg Horn
 *
 *    CasADi is free software; you can redistribute it and/or
 *    modify it under the terms of the GNU Lesser General Public
 *    License as published by the Free Software Foundation; either
 *    version 3 of the License, or (at your option) any later version.
 *
 *    CasADi is distributed in the hope that it will be useful,
 *    but WITHOUT ANY WARRANTY; without even the implied warranty of
 *    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 *    Lesser General Public License for more details.
 *
 *    You should have received a copy of the GNU Lesser General Public
 *    License along with CasADi; if not, write to the Free Software
 *    Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
 *
 */
 
 
#ifndef CASADI_SCPGEN_HPP
#define CASADI_SCPGEN_HPP
 
#include "casadi/core/nlpsol_impl.hpp"
#include <casadi/solvers/casadi_nlpsol_scpgen_export.h>
 
 
namespace casadi {
 
  struct  ScpgenMemory : public NlpsolMemory {
    // Work vectors, nonlifted problem
    double *dxk, *dlam, *gfk, *gL, *b_gn;
    // Memory for lifted variables
    struct VarMem {
      casadi_int n;
      double *dx, *x0, *x, *lam, *dlam;
      double *res, *resL;
    };
    std::vector<VarMem> lifted_mem;
    // Penalty parameter of merit function
    double sigma;
    // 1-norm of last primal step
    double pr_step;
    // 1-norm of last dual step
    double du_step;
    // Regularization
    double reg;
    // Message applying to a particular iteration
    const char* iteration_note;
    // QP
    double *qpH, *qpA, *qpB, *qpL, *qpG, *qpH_times_du;
    // QP solver
    double *lbdz, *ubdz;
    // Linesearch parameters
    double* merit_mem;
    casadi_int merit_ind;
    // Timers
    double t_eval_mat, t_eval_res, t_eval_vec, t_eval_exp, t_solve_qp, t_mainloop;
    // Current iteration
    casadi_int iter_count;
  };
 
  class  Scpgen : public Nlpsol {
  public:
    explicit Scpgen(const std::string& name, const Function& nlp);
    ~Scpgen() override;
 
    // Get name of the plugin
    const char* plugin_name() const override { return "scpgen";}
 
    // Name of the class
    std::string class_name() const override { return "Scpgen";}
 
    static Nlpsol* creator(const std::string& name, const Function& nlp) {
      return new Scpgen(name, nlp);
    }
 
 
    static const Options options_;
    const Options& get_options() const override { return options_;}
 
    // Initialize the solver
    void init(const Dict& opts) override;
 
    void* alloc_mem() const override { return new ScpgenMemory();}
 
    int init_mem(void* mem) const override;
 
    void free_mem(void *mem) const override { delete static_cast<ScpgenMemory*>(mem);}
 
    Dict get_stats(void* mem) const override;
 
    void set_work(void* mem, const double**& arg, double**& res,
                          casadi_int*& iw, double*& w) const override;
 
    // Solve the NLP
    int solve(void* mem) const override;
 
    // Calculate the L1-norm of the primal infeasibility
    double primalInfeasibility(ScpgenMemory* m) const;
 
    // Calculate the L1-norm of the dual infeasibility
    double dualInfeasibility(ScpgenMemory* m) const;
 
    // Print iteration header
    void printIteration(ScpgenMemory* m, std::ostream &stream) const;
 
    // Print iteration
    void printIteration(ScpgenMemory* m, std::ostream &stream, casadi_int iter, double obj,
                        double pr_inf, double du_inf,
                        double rg, casadi_int ls_trials, bool ls_success) const;
 
    // Evaluate the matrices in the condensed QP
    void eval_mat(ScpgenMemory* m) const;
 
    // Evaluate the vectors in the condensed QP
    void eval_vec(ScpgenMemory* m) const;
 
    // Evaluate the residual function
    void eval_res(ScpgenMemory* m) const;
 
    // Regularize the condensed QP
    void regularize(ScpgenMemory* m) const;
 
    // Solve the QP to get the (full) step
    void solve_qp(ScpgenMemory* m) const;
 
    // Perform the line-search to take the step
    void line_search(ScpgenMemory* m, casadi_int& ls_iter, bool& ls_success) const;
 
    // Evaluate the step expansion
    void eval_exp(ScpgenMemory* m) const;
 
    Function qpsol_;
 
    bool gauss_newton_;
 
    casadi_int max_iter_;
 
    casadi_int lbfgs_memory_;
 
    double tol_pr_;
 
    double tol_du_;
 
    double tol_reg_;
 
    double tol_pr_step_;
 
    double tol_gl_;
 
    double c1_;
    double beta_;
    casadi_int max_iter_ls_;
    casadi_int merit_memsize_;
    double merit_start_;
 
    bool codegen_;
 
    const Function getConic() const { return qpsol_;}
 
    bool regularize_;
 
    // Number of gauss_newton equations
    casadi_int ngn_;
 
    // Options
    double reg_threshold_;
 
    bool print_time_;
 
    Function vinit_fcn_;
 
    Function res_fcn_;
 
    // Function to calculate the matrices in the reduced QP
    Function mat_fcn_;
    casadi_int mat_jac_, mat_hes_;
 
    Function vec_fcn_;
    casadi_int vec_gf_, vec_g_;
 
    Function exp_fcn_;
 
    // Residual function io indices
    casadi_int res_x_, res_p_, res_g_lam_, res_p_lam_, res_p_d_;
    casadi_int res_f_, res_gl_, res_g_;
 
    // Modifier function io indices
    casadi_int mod_x_, mod_p_, mod_g_lam_;
    casadi_int mod_f_, mod_gl_, mod_g_;
    casadi_int mod_du_, mod_dlam_g_;
 
    struct Var {
      casadi_int n;
      MX v, v_def, v_lam, v_defL;
      MX d, d_def, d_lam, d_defL;
 
      // Indices of function inputs and outputs
      casadi_int res_var, res_lam, res_d, res_lam_d;
      casadi_int mod_var, mod_lam, mod_def, mod_defL;
      casadi_int exp_def, exp_defL;
    };
 
    std::vector<Var> v_;
 
    // Names of the components
    std::vector<std::string> name_x_;
 
    // Components to print
    std::vector<casadi_int> print_x_;
 
    // QP sparsity
    Sparsity spH_, spA_, spL_;
 
    // Print options
    bool print_header_;
 
    static const std::string meta_doc;
  };
 
} // namespace casadi
 
 
#endif // CASADI_SCPGEN_HPP