QPKKTSolverObject.hxx

Go to the documentation of this file.

#ifndef CONICBUNDLE_QPKKTSOLVEROBJECT_HXX
#define CONICBUNDLE_QPKKTSOLVEROBJECT_HXX

#include "QPSolverObject.hxx"
#include "QPModelBlockObject.hxx"

namespace ConicBundle {

  class QPSolverParameters;


class QPKKTSolverObject: public virtual CBout
{
protected:

  //--- data describing the KKT system
  QPSolverProxObject* Hp;  
  QPModelBlockObject* model; 
  const CH_Matrix_Classes::Sparsemat* A;  
  const CH_Matrix_Classes::Indexmatrix* eq_indices; 

  CH_Matrix_Classes::Real blockH_norm;  
  CH_Matrix_Classes::Real blockA_norm;  
  
public:
  virtual void clear()
  {Hp=0;model=0;A=0;eq_indices=0;blockH_norm=1.;blockA_norm=1.;}

  QPKKTSolverObject(CBout* cb=0,int cbinc=-1):
    CBout(cb,cbinc),Hp(0),model(0),A(0),eq_indices(0),blockH_norm(1.),blockA_norm(1.)
  {}

  virtual ~QPKKTSolverObject();
  
  
  virtual int QPinit_KKTdata(QPSolverProxObject* Hp, 
                             QPModelBlockObject* model, 
                             const CH_Matrix_Classes::Sparsemat* A, 
                             const CH_Matrix_Classes::Indexmatrix* eq_indices 
                             )=0;

  virtual int QPinit_KKTsystem(const CH_Matrix_Classes::Matrix& KKTdiagx,
                               const CH_Matrix_Classes::Matrix& KKTdiagy,
                               CH_Matrix_Classes::Real Hfactor,
                               CH_Matrix_Classes::Real prec,
                               QPSolverParameters* params) =0;

  virtual int QPsolve_KKTsystem(CH_Matrix_Classes::Matrix& solx,
                                 CH_Matrix_Classes::Matrix& soly,
                                 const CH_Matrix_Classes::Matrix& primalrhs,
                                 const CH_Matrix_Classes::Matrix& dualrhs,
                                 CH_Matrix_Classes::Real rhsmu,
                                 CH_Matrix_Classes::Real rhscorr,
                                 CH_Matrix_Classes::Real prec,
                                 QPSolverParameters* params) =0;
  
  virtual CH_Matrix_Classes::Real QPget_blockH_norm()
  {return blockH_norm;}

  virtual CH_Matrix_Classes::Real QPget_blockA_norm()
  {return blockA_norm;}

  virtual CH_Matrix_Classes::Integer QPget_nmult() const
  {return 0;}

  virtual CH_Matrix_Classes::Real QPget_condition_number() 
  {return -1.;}

  virtual CH_Matrix_Classes::Integer QPget_precond_rank() 
  {return -1;}

  virtual CH_Matrix_Classes::Integer QPget_system_size() 
  {return 0;}
};






}

#endif