QPKKTPrecondObject.hxx

Go to the documentation of this file.

#ifndef CONICBUNDLE_QPKKTPRECONDOBJECT_HXX
#define CONICBUNDLE_QPKKTPRECONDOBJECT_HXX

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

namespace ConicBundle {

class QPSolverParameters;

  

class QPKKTPrecondObject: 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; 
  bool SchurComplAineq; 

  CH_Matrix_Classes::Real Hfactor; 
  
public:
  virtual void clear()
  {Hp=0;model=0;A=0;eq_indices=0;}
  
  QPKKTPrecondObject(CBout* cb=0,int cbinc=-1):
    CBout(cb,cbinc),Hp(0),model(0),A(0),eq_indices(0),SchurComplAineq(false)
    {}

  virtual ~QPKKTPrecondObject();
  
  virtual int init_data(QPSolverProxObject* Hp, 
                        QPModelBlockObject* model, 
                        const CH_Matrix_Classes::Sparsemat* A, 
                        const CH_Matrix_Classes::Indexmatrix* eq_indices, 
                        bool SchurComplAineq 
                        )=0;
  
  virtual int init_system(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 CH_Matrix_Classes::Real get_lmin_invM1() {return 1.;}
  
  virtual int precondM1(CH_Matrix_Classes::Matrix& /* vec */) {return 0;}
  
  virtual int precondM2(CH_Matrix_Classes::Matrix& /* vec */) {return 0;}

  virtual int precond_invG1(CH_Matrix_Classes::Matrix& /* vec */) {return 0;}

  virtual int precond_invG1tran(CH_Matrix_Classes::Matrix& /* vec */) {return 0;}

  virtual CH_Matrix_Classes::Integer precond_size() {return -1;}
  
  virtual int cond_number_mult(CH_Matrix_Classes::Matrix& /* vec */,
                               const CH_Matrix_Classes::Matrix& /* KKTdiagx */,
                               const CH_Matrix_Classes::Matrix& /* KKTdiagy */)
  {return 1;}

  virtual  CH_Matrix_Classes::Real get_condition_number(const CH_Matrix_Classes::Matrix& KKTdiagx,
                                                        const CH_Matrix_Classes::Matrix& KKTdiagy);
  
  virtual CH_Matrix_Classes::Integer get_precond_rank()
  {return -1;}

  virtual CH_Tools::Microseconds get_t_precond_mult()
  {return 0;}

  virtual void reset_t_precond_mult()
  {}
};




  
}

#endif