QPSolverObject.hxx

Go to the documentation of this file.

 
#ifndef CONICBUNDLE_QPSOLVEROBJECT_HXX
#define CONICBUNDLE_QPSOLVEROBJECT_HXX

#include "QPModelDataObject.hxx"

namespace ConicBundle {



class QPSolverParametersObject
{
public:
  virtual ~QPSolverParametersObject();
};

class QPSolverProxObject
{
public:
  virtual ~QPSolverProxObject();

  virtual bool is_DLR() const =0;

  virtual int add_H(CH_Matrix_Classes::Symmatrix& big_sym,CH_Matrix_Classes::Integer start_index=0) const=0;

  virtual CH_Matrix_Classes::Matrix& add_Hx(const CH_Matrix_Classes::Matrix& x,
                                            CH_Matrix_Classes::Matrix& outplusHx,
                                            CH_Matrix_Classes::Real alpha=1.) const =0;

  virtual CH_Matrix_Classes::Matrix& apply_Hinv(CH_Matrix_Classes::Matrix& x) const =0;

  virtual void get_precond(CH_Matrix_Classes::Matrix& D, 
                           const CH_Matrix_Classes::Matrix*& Vp) const=0;
  

  virtual void set_short_QPsteps(CH_Matrix_Classes::Integer shortQPst)=0;
  
  virtual CH_Matrix_Classes::Integer get_short_QPsteps()=0;

};
  
  
  class QPSolverObject: public virtual QPModelDataPointer
{
public:
  virtual void QPclear()=0;

  QPSolverObject(CBout* cb=0,int cbinc=-1):QPModelDataPointer(cb,cbinc){}
 
  virtual ~QPSolverObject();

  virtual int QPset_parameters(QPSolverParametersObject* params)=0;

  virtual bool QPsupports_yfixing()=0;

  virtual bool QPsupports_updates()=0;

  virtual int QPsolve(const CH_Matrix_Classes::Matrix& center_y,
                      CH_Matrix_Classes::Real lower_bound,
                      CH_Matrix_Classes::Real upper_bound,
                      CH_Matrix_Classes::Real relprec,
                      QPSolverProxObject* Hp,
                      const MinorantPointer& gs_aggr,
                      CH_Matrix_Classes::Indexmatrix* yfixed)=0;
  
  virtual int QPupdate(const CH_Matrix_Classes::Matrix& center_y,
                       CH_Matrix_Classes::Real lower_bound,
                       CH_Matrix_Classes::Real upper_bound,
                       CH_Matrix_Classes::Real relprec,
                       QPSolverProxObject* Hp,
                       const MinorantPointer& gs_aggr,
                       CH_Matrix_Classes::Indexmatrix* yfixed,
                       const MinorantPointer& delta_gs_aggr,
                       const CH_Matrix_Classes::Indexmatrix& delta_index)=0;

  virtual int QPresolve(CH_Matrix_Classes::Real lower_bound,
                      CH_Matrix_Classes::Real upper_bound,
                      CH_Matrix_Classes::Real relprec)=0;

  virtual int QPget_solution(CH_Matrix_Classes::Real& augval_lb,
                             CH_Matrix_Classes::Real& augval_ub,
                             CH_Matrix_Classes::Matrix& new_point,
                             CH_Matrix_Classes::Real& gs_aggr_offset,
                             CH_Matrix_Classes::Matrix& gs_aggr_gradient)=0;

  virtual CH_Matrix_Classes::Real QPget_lower_bound()=0;

  virtual GroundsetModification* QPstart_modification() =0;

  virtual int QPapply_modification(const GroundsetModification& mdf)=0;

  virtual bool QPis_feasible(const CH_Matrix_Classes::Matrix& y,
                             CH_Matrix_Classes::Real relprec=1e-10)=0;

  virtual int QPensure_feasibility(CH_Matrix_Classes::Matrix& y,
                                   bool& ychanged,
                                   QPSolverProxObject* Hp,
                                   CH_Matrix_Classes::Real relprec=1e-10)=0;
  
  virtual bool QPprefer_UQPSolver(QPSolverProxObject*) const=0;

  virtual bool QPconstrained() const=0;

  virtual bool QPboxconstrained(const CH_Matrix_Classes::Matrix*& lb,
                                const CH_Matrix_Classes::Matrix*& ub,
                                const CH_Matrix_Classes::Indexmatrix*& lbind,
                                const CH_Matrix_Classes::Indexmatrix*& ubind) const =0;

  virtual std::ostream& QPprint_statistics(std::ostream& out,int printlevel=0)=0;
  
 };
  
 


}

#endif