QPSolverParameters.hxx

Go to the documentation of this file.

#ifndef CONICBUNDLE_QPSOLVERPARAMETERS_HXX
#define CONICBUNDLE_QPSOLVERPARAMETERS_HXX

#include "QPKKTSolverObject.hxx"

namespace ConicBundle {



class QPSolverParameters: public virtual QPSolverParametersObject, public virtual CBout
{
private:
  //--- preset termination parameters
  CH_Matrix_Classes::Real min_objective_relprec; 
  CH_Matrix_Classes::Integer maxiter; 

  //--- selecting system and system solver
  bool allow_unconstrained; 

  //--- adaptive termination parameters set by the code
  CH_Matrix_Classes::Real objective_gap_eps;   
  CH_Matrix_Classes::Real primal_infeasibility_eps;  
  CH_Matrix_Classes::Real dual_infeasibility_eps; 
  CH_Matrix_Classes::Real lower_bound;  
  CH_Matrix_Classes::Real upper_bound;  
  CH_Matrix_Classes::Real lower_bound_gap_eps; 
  CH_Matrix_Classes::Real upper_bound_gap_eps; 

  QPKKTSolverObject* KKTsolver; 
  
  bool use_predictor_corrector; 
  bool use_neighborhood; 
  CH_Matrix_Classes::Real nbh_ub; 
  CH_Matrix_Classes::Real nbh_lb; 
  
  bool use_socqp; 

  QPSolverParameters(const QPSolverParameters& /*params*/);  

  QPSolverParameters& operator=(const QPSolverParameters&);

public:
  QPSolverParameters(CBout* cb=0,int incr=-1);    
  virtual ~QPSolverParameters(){delete KKTsolver;}

  CH_Matrix_Classes::Real QPget_min_objective_relprec() const {return min_objective_relprec;}
  CH_Matrix_Classes::Real QPget_objective_gap_eps() const {return objective_gap_eps;}
  CH_Matrix_Classes::Real QPget_primal_infeasibility_eps() const {return primal_infeasibility_eps;}
  CH_Matrix_Classes::Real QPget_dual_infeasibility_eps() const {return dual_infeasibility_eps;}
  CH_Matrix_Classes::Real QPget_lower_bound() const {return lower_bound;}
  CH_Matrix_Classes::Real QPget_upper_bound() const {return upper_bound;}
  CH_Matrix_Classes::Real QPget_lower_bound_gap_eps() const {return lower_bound_gap_eps;}
  CH_Matrix_Classes::Real QPget_upper_bound_gap_eps() const {return upper_bound_gap_eps;}
  CH_Matrix_Classes::Integer QPget_maxiter() const {return maxiter;}
  QPKKTSolverObject* QPget_KKTsolver() {return KKTsolver;};
  bool QPget_use_predictor_corrector() const {return use_predictor_corrector;}
  bool QPget_use_neighborhood() const {return use_neighborhood;}
  CH_Matrix_Classes::Real QPget_nbh_ub() const {return nbh_ub;}
  CH_Matrix_Classes::Real QPget_nbh_lb() const {return nbh_lb;}
  bool QPget_use_socqp() const {return use_socqp;}


  int QPset_min_objective_relprec(CH_Matrix_Classes::Real eps)
  {min_objective_relprec=(eps>0)?eps:min_objective_relprec; return (eps<=0);}
  int QPset_objective_gap_eps(CH_Matrix_Classes::Real eps)
  {objective_gap_eps=(eps>0)?eps:objective_gap_eps; return (eps<=0);}
  int QPset_primal_infeasibility_eps(CH_Matrix_Classes::Real eps)
  {primal_infeasibility_eps=(eps>0)?eps:primal_infeasibility_eps; return (eps<=0);}
  int QPset_dual_infeasibility_eps(CH_Matrix_Classes::Real eps)
  {dual_infeasibility_eps=(eps>0)?eps:dual_infeasibility_eps; return (eps<=0);}
  int QPset_lower_and_upper_bounds(CH_Matrix_Classes::Real lb,CH_Matrix_Classes::Real ub)
  {if (lb>ub) return 1;lower_bound=lb;upper_bound=ub;return 0;}
  int QPset_lower_bound_gap_eps(CH_Matrix_Classes::Real eps)
  {lower_bound_gap_eps=(eps>0)?eps:lower_bound_gap_eps; return (eps<=0);}
  int QPset_upper_bound_gap_eps(CH_Matrix_Classes::Real eps)
  {upper_bound_gap_eps=(eps>0)?eps:upper_bound_gap_eps; return (eps<=0);}
  int QPset_maxiter(CH_Matrix_Classes::Integer mi)
  {maxiter=mi;return 0;}

  int QPset_KKTsolver(QPKKTSolverObject* in_KKTsolver)
  {delete KKTsolver; KKTsolver=in_KKTsolver;return 0;}

  int QPset_use_predictor_corrector(bool upc)
  {use_predictor_corrector=upc; return 0;}

  int QPset_use_neighborhood(bool nbh)
  {use_neighborhood=nbh; return 0;}

  int QPset_nbh_bounds(CH_Matrix_Classes::Real nbhlb,CH_Matrix_Classes::Real nbhub)
  {
    nbh_lb=CH_Matrix_Classes::max(CH_Matrix_Classes::eps_Real,nbhlb);
    nbh_ub=CH_Matrix_Classes::max(nbh_lb,nbhub);
    return 0;
  }

  int QPset_use_socqp(bool s)
  {use_socqp=s; return 0;}

  int QPset_allow_UQPSolver(bool allow)
  {allow_unconstrained=allow; return 0;}
  bool QPallow_UQPSolver()
  {return allow_unconstrained;}
};






}

#endif