BundleHKWeight.hxx

#ifndef CONICBUNDLE_BUNDLEHKWEIGHT_HXX
#define CONICBUNDLE_BUNDLEHKWEIGHT_HXX

#include "BundleWeight.hxx"

namespace ConicBundle {



class BundleHKWeight: public BundleWeight
{
  Groundset* groundset;  
  BundleModel* model;    
  CH_Matrix_Classes::Integer iweight; 
  CH_Matrix_Classes::Real weight; 
  CH_Matrix_Classes::Real epsweight; 
  CH_Matrix_Classes::Real minweight; 
  CH_Matrix_Classes::Real maxweight; 
  CH_Matrix_Classes::Real modelmax; 
  
  bool weightchanged;    
  bool next_weight_set;  
  int nullstep_updates; 

  CH_Matrix_Classes::Matrix valuelevels; 
  CH_Matrix_Classes::Matrix ratiolevels; 
  
  CH_Matrix_Classes::Real mR; 
    
public:
  BundleHKWeight(CH_Matrix_Classes::Real mRin=.5,BundleWeight* bwp=0,const CBout* cbo=0,int incr=-1);
  ~BundleHKWeight(){}

 
  virtual void set_defaults(); 

  virtual void set_nullstep_updates(int nu=0)
  {nullstep_updates=nu;} 

  virtual void clear();        
  
  int init(CH_Matrix_Classes::Real aggr_dnmormsqr,Groundset* groundset,BundleModel* model);
  
  void set_next_weight(CH_Matrix_Classes::Real u)
  { if (u<=0.) return; 
  weight=CH_Matrix_Classes::max(u,1e-10); weightchanged=true;next_weight_set=true;
  modelmax=CB_minus_infinity;iweight=0;}
  
  void set_minweight(CH_Matrix_Classes::Real mw)
  { 
    minweight=mw; 
    if (minweight>0){
      if ((weight>0)&&(weight<minweight)) 
        weight=minweight;
      if ((maxweight>0)&&(maxweight<minweight)) 
        maxweight=minweight;
    }
  }

  bool get_next_weight_set() const {return next_weight_set;}

  CH_Matrix_Classes::Real get_minweight() const {return minweight;}
  
  virtual void set_maxweight(CH_Matrix_Classes::Real mw)
  { 
    maxweight=mw; 
    if (maxweight>0){
      if (weight>maxweight){
        weight=maxweight;
      }
      if ((minweight>0)&&(minweight>maxweight)) 
        minweight=maxweight;
    }
  }
  
  CH_Matrix_Classes::Real get_maxweight() const {return maxweight;}

  CH_Matrix_Classes::Real get_weight() const;
  
  bool weight_changed() const;
  
  int descent_update(CH_Matrix_Classes::Real newval,
                     CH_Matrix_Classes::Real oldval,
                     CH_Matrix_Classes::Real modelval,
                     const CH_Matrix_Classes::Matrix& y, 
                     const CH_Matrix_Classes::Matrix& newy,
                     CH_Matrix_Classes::Real normsubg2,
                     BundleProxObject* Hp);
  
  int nullstep_update(CH_Matrix_Classes::Real newval,
                      CH_Matrix_Classes::Real oldval,
                      CH_Matrix_Classes::Real modelval,
                      const CH_Matrix_Classes::Matrix& y, 
                      const CH_Matrix_Classes::Matrix& newy,
                      MinorantPointer& new_minorant,
                      MinorantPointer& aggregate,
                      CH_Matrix_Classes::Real nullstep_bound,
                      CH_Matrix_Classes::Real normsubg2,
                      BundleProxObject* Hp);
      
  virtual int apply_modification(const GroundsetModification& gsmdf);
 
};

}

#endif