UnconstrainedGroundset.hxx

Go to the documentation of this file.

#ifndef CONICBUNDLE_UNCONSTRAINEDGROUNDSET_HXX
#define CONICBUNDLE_UNCONSTRAINEDGROUNDSET_HXX


#include "Groundset.hxx"
#include "LPGroundsetModification.hxx"

namespace ConicBundle {


  class UnconstrainedGroundset:virtual public Groundset
{
protected:
  
  CH_Matrix_Classes::Integer dim;  

  CH_Matrix_Classes::Integer groundset_id;

  CH_Matrix_Classes::Matrix starting_point;

  MinorantPointer gs_aggregate;

  bool use_yfixing;                       

  CH_Matrix_Classes::Indexmatrix yfixed;

  QPSolverObject* qp_solver;

  //data for describing the QP
  const BundleProxObject* Hp;  
  CH_Matrix_Classes::Matrix c; 
  CH_Matrix_Classes::Real gamma; 
    

public:
  virtual void clear( CH_Matrix_Classes::Integer indim=0, CH_Matrix_Classes::Integer in_groundset_id=0);

  UnconstrainedGroundset(CH_Matrix_Classes::Integer indim=0,
            const CH_Matrix_Classes::Matrix* start_val=0,
            const CH_Matrix_Classes::Matrix* costs=0,
            const CH_Matrix_Classes::Real offset=0., 
            CH_Matrix_Classes::Integer in_groundset_id=0);

  virtual ~UnconstrainedGroundset();

  CH_Matrix_Classes::Integer get_groundset_id() const {return groundset_id;}

  void set_groundset_id(CH_Matrix_Classes::Integer gsid){groundset_id=gsid;}

  virtual CH_Matrix_Classes::Integer get_dim() const 
  {return dim;}

  virtual bool constrained() const {return false;}

  virtual bool is_feasible(CH_Matrix_Classes::Integer& in_groundset_id,
                            const CH_Matrix_Classes::Matrix& y,
                             CH_Matrix_Classes::Real relprec=1e-10);

  virtual int ensure_feasibility(CH_Matrix_Classes::Integer& in_groundset_id,
                                 CH_Matrix_Classes::Matrix& y,
                                 bool& ychanged,
                                 BundleProxObject* Hp=0,
                                 CH_Matrix_Classes::Real relprec=1e-10);

  virtual QPSolverObject* get_qp_solver(bool& solves_model_without_gs,
                                        BundleProxObject* Hp);

  int set_qp_solver_parameters(QPSolverParametersObject* /* in_qpsp */)
  {
    return 0;
  }

  virtual const CH_Matrix_Classes::Matrix& get_starting_point() const 
  {return starting_point;}

  virtual int set_starting_point(const CH_Matrix_Classes::Matrix& vec)
  {starting_point=vec; int dummy=-1; return (is_feasible(dummy,vec)==false);}

  virtual int candidate(CH_Matrix_Classes::Integer& gs_id,
                        CH_Matrix_Classes::Matrix& newy,
                        CH_Matrix_Classes::Real& cand_gs_val,
                        CH_Matrix_Classes::Real& linval,
                        CH_Matrix_Classes::Real& augval_lb,
                        CH_Matrix_Classes::Real& augval_ub,
                        CH_Matrix_Classes::Real& subgnorm2,
                        const CH_Matrix_Classes::Matrix& center_y,
                        CH_Matrix_Classes::Real center_value,
                        const MinorantPointer& model_minorant,
                        BundleProxObject* Hp,
                        MinorantPointer* delta_groundset_minorant=0,
                        CH_Matrix_Classes::Indexmatrix* delta_index=0,
                        CH_Matrix_Classes::Real relprec=1e-2);

  virtual const MinorantPointer& get_gs_aggregate() const { return gs_aggregate;}

  virtual const MinorantPointer& get_gs_minorant() const { return gs_aggregate;}

  virtual const CH_Matrix_Classes::Indexmatrix* get_yfixed() const
  {return (use_yfixing? &yfixed : 0);}

  virtual CH_Matrix_Classes::Indexmatrix* set_yfixed()
  {return (use_yfixing? &yfixed : 0);}

  bool get_use_yfixing() const {return use_yfixing;}

  void set_use_yfixing(bool uyf) {use_yfixing=uyf;}

  GroundsetModification* start_modification()
  { return new LPGroundsetModification(dim,0,this);}

  int apply_modification(const GroundsetModification& mdf);

  int mfile_data(std::ostream& out) const;

  void set_cbout(const CBout* cb, int incr=-1)
  {
    VariableMetricModel::set_cbout(cb,incr);
    if (qp_solver) 
      qp_solver->set_cbout(this);
  }
};




}

#endif