UQPSumModelBlock.hxx

Go to the documentation of this file.

#ifndef CONICBUNDLE_UQPSUMMODELBLOCK_HXX
#define CONICBUNDLE_UQPSUMMODELBLOCK_HXX

#include "UQPModelBlock.hxx"

namespace ConicBundle {

  
  class UQPSumModelBlock: public QPSumModelDataObject, public UQPModelBlock 
{
private:
  std::vector<UQPModelBlock*> blocks;

  CH_Matrix_Classes::Integer xstart; 
  CH_Matrix_Classes::Integer xend;  
  CH_Matrix_Classes::Integer ystart; 
  CH_Matrix_Classes::Integer yend; 


public:
  void clear()
  {blocks.clear();xstart=0;xend=0;ystart=0;yend=0;UQPModelBlock::clear();}

  UQPSumModelBlock(CBout* cb=0,int incr=-1):QPSumModelDataObject(cb,incr),UQPModelBlock(cb,incr){clear();}
  
  ~UQPSumModelBlock();

  int append(QPModelDataObject* inblock);
  
  //-----------  UQPModelBlock routines

  CH_Matrix_Classes::Integer xdim() const;

  CH_Matrix_Classes::Integer ydim() const;

  int set_qp_xstart(CH_Matrix_Classes::Integer x_start_index);
  
  int set_qp_ystart(CH_Matrix_Classes::Integer y_start_index);

  int starting_x(CH_Matrix_Classes::Matrix& qp_x);

  int starting_y(CH_Matrix_Classes::Matrix& qp_y,
                 const CH_Matrix_Classes::Matrix& qp_Qx,
                 const CH_Matrix_Classes::Matrix& qp_c);

  CH_Matrix_Classes::Real get_local_primalcost() const; 

  CH_Matrix_Classes::Real get_local_dualcost() const; 

  int get_Ab(CH_Matrix_Classes::Matrix& qp_A,CH_Matrix_Classes::Matrix &qp_b) const;
      
  int restart_x(CH_Matrix_Classes::Matrix& qp_x,const CH_Matrix_Classes::Matrix& qp_c,const CH_Matrix_Classes::Matrix& qp_dc);
    

  int restart_y(CH_Matrix_Classes::Matrix& qp_y,
                const CH_Matrix_Classes::Matrix& qp_Qx,
                const CH_Matrix_Classes::Matrix& qp_c,
                const CH_Matrix_Classes::Matrix& qp_dc);

  int add_xinv_kron_z(CH_Matrix_Classes::Symmatrix& barQ);
    

  int add_local_sys(CH_Matrix_Classes::Symmatrix& sysdy,CH_Matrix_Classes::Matrix& rhs);

  int suggest_mu(CH_Matrix_Classes::Real& ip_xz,
                 CH_Matrix_Classes::Integer& mu_dim,
                 CH_Matrix_Classes::Real& sigma,
                 const CH_Matrix_Classes::Matrix& qp_dx,
                 const CH_Matrix_Classes::Matrix& qp_dy,
                 const CH_Matrix_Classes::Matrix& rhs_residual);
    

  int get_corr(CH_Matrix_Classes::Matrix& xcorr,
               CH_Matrix_Classes::Matrix& rhs,
               CH_Matrix_Classes::Real mu);
       

  int line_search(CH_Matrix_Classes::Real& alpha,
                  const CH_Matrix_Classes::Matrix& qp_dx,
                  const CH_Matrix_Classes::Matrix& qp_dy,
                  const CH_Matrix_Classes::Matrix& rhs_residual);
  

  int set_point(const CH_Matrix_Classes::Matrix& qp_x,
                const CH_Matrix_Classes::Matrix& qp_y,
                CH_Matrix_Classes::Real alpha);

   int add_modelx_aggregate(CH_Matrix_Classes::Real& offset,
                            CH_Matrix_Classes::Matrix& gradient);

  void set_out(std::ostream* o=0,int pril=1);

  void set_cbout(const CBout* cb,int incr=-1);
  
  //---------------- for debugging purposes

  CH_Matrix_Classes::Matrix& add_Bs(CH_Matrix_Classes::Matrix& qp_vec) const;

  CH_Matrix_Classes::Matrix& subtract_z(CH_Matrix_Classes::Matrix& dual_residual,bool with_step=false) const;

};


}

#endif