Method write_solution_to_file

interfaces/AMPL/AMPLModel.cpp:294–329

   }

   void AMPLModel::write_solution_to_file(Result& result) const {
      // write the primal-dual solution and status into a *.sol file
      this->asl->p.solve_code_ = 400; // limit
      if (result.solution_status == SolutionStatus::FEASIBLE_KKT_POINT) {
         this->asl->p.solve_code_ = 0;
      }
      else if (result.solution_status == SolutionStatus::FEASIBLE_SMALL_STEP) {
         this->asl->p.solve_code_ = 100;
      }
      else if (result.solution_status == SolutionStatus::INFEASIBLE_STATIONARY_POINT) {
         this->asl->p.solve_code_ = 200;
      }
      else if (result.solution_status == SolutionStatus::UNBOUNDED) {
         this->asl->p.solve_code_ = 300;
      }
      else if (result.solution_status == SolutionStatus::INFEASIBLE_SMALL_STEP) {
         this->asl->p.solve_code_ = 500;
      }

      // include the bound duals in the .sol file, using suffixes
      SufDecl lower_bound_suffix{const_cast<char*>("lower_bound_duals"), nullptr, ASL_Sufkind_var | ASL_Sufkind_real, 0};
      SufDecl upper_bound_suffix{const_cast<char*>("upper_bound_duals"), nullptr, ASL_Sufkind_var | ASL_Sufkind_real, 0};
      std::array<SufDecl, 2> suffixes{lower_bound_suffix, upper_bound_suffix};
      suf_declare_ASL(this->asl, suffixes.data(), suffixes.size());
      suf_rput_ASL(this->asl, "lower_bound_duals", ASL_Sufkind_var, result.lower_bound_dual_solution.data());
      suf_rput_ASL(this->asl, "upper_bound_duals", ASL_Sufkind_var, result.upper_bound_dual_solution.data());

      Option_Info option_info{};
      option_info.wantsol = 9; // write the solution without printing the message to stdout
      std::string message = "Uno ";
      message.append(Uno::current_version()).append(": ").append(solution_status_to_message(result.solution_status));

      // flip the signs of the constraint multipliers (different Lagrangian convention between ASL and Uno)
      Vector<double> ampl_dual_solution = -result.constraint_dual_solution;
      write_sol_ASL(this->asl, message.data(), result.primal_solution.data(), ampl_dual_solution.data(), &option_info);
   }

   size_t AMPLModel::number_jacobian_nonzeros() const {
      return static_cast<size_t>(this->asl->i.nzc_);