BundleProxObject.hxx

Go to the documentation of this file.

#ifndef CONICBUNDLE_BUNDLEPROXOBJECT_HXX
#define CONICBUNDLE_BUNDLEPROXOBJECT_HXX


#include "AffineFunctionTransformation.hxx"
#include "QPSolverObject.hxx"
#include "VariableMetric.hxx"

namespace ConicBundle {

class BundleModel;
class Groundset;




  class BundleProxObject: public VariableMetric, public QPSolverProxObject 
{
protected:
  static const CH_Matrix_Classes::Integer xdim_threshold = 100; 

  CH_Matrix_Classes::Real factor;

  CH_Matrix_Classes::Integer short_QPsteps;

public:
  BundleProxObject(VariableMetricSelection* vp=0,bool local_scaling=false,bool bounds_scaling=false,CBout* cb=0,int cbincr=-1);
  //:DynamicScaling(dynamic_scaling,local_scaling,bounds_scaling){}
  //somehow inline did not work here, so I moved it to BundleProxObject.cxx

  virtual ~BundleProxObject(); //implemented in BundleScaling.cxx

  virtual void set_weightu(CH_Matrix_Classes::Real weightu)=0;

  virtual void set_short_QPsteps(CH_Matrix_Classes::Integer shortQPst)
  {short_QPsteps=shortQPst;}
  
  virtual CH_Matrix_Classes::Integer get_short_QPsteps()
  {return short_QPsteps;}
  
  
  int apply_factor(CH_Matrix_Classes::Real f)
  {if (f<=0.) return 1; factor*=f; return 0; }

  CH_Matrix_Classes::Real get_factor()
  {return factor;}
  
  virtual CH_Matrix_Classes::Real get_weightu() const=0;

  virtual CH_Matrix_Classes::Real get_term_corr(void) const {return 1.;}
   
  virtual CH_Matrix_Classes::Real norm_sqr(const CH_Matrix_Classes::Matrix& B) const =0;
  virtual CH_Matrix_Classes::Real dnorm_sqr(const MinorantPointer& B) const =0;

  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 int compute_QP_costs(CH_Matrix_Classes::Symmatrix& Q,
                               CH_Matrix_Classes::Matrix& d,
                               CH_Matrix_Classes::Real& offset,
                               const MinorantPointer& constant_minorant,
                               const MinorantBundle& bundle,
                               const CH_Matrix_Classes::Matrix& center_y,
                               const MinorantPointer& groundset_minorant,
                               CH_Matrix_Classes::Indexmatrix* yfixed) =0;

  

  virtual int update_QP_costs(CH_Matrix_Classes::Symmatrix& delta_Q, 
                              CH_Matrix_Classes::Matrix& delta_d,  
                              CH_Matrix_Classes::Real& delta_offset,
                              const MinorantPointer& constant_minorant,
                              const MinorantBundle& bundle,
                              const CH_Matrix_Classes::Matrix& center_y,
                              const MinorantPointer& groundset_minorant,
                              const MinorantPointer& delta_groundset_minorant,
                              const CH_Matrix_Classes::Indexmatrix& delta_index,
                              CH_Matrix_Classes::Indexmatrix* yfixed) =0;

  virtual int apply_modification(const GroundsetModification& gsmdf)=0;

  virtual BundleProxObject* projected_clone(const CH_Matrix_Classes::Indexmatrix& indices) =0;

};



}

#endif