NNCBoxSupportFunction.hxx

Go to the documentation of this file.

#ifndef CONICBUNDLE_NNCBOXSUPPORTFUNCTION_HXX
#define CONICBUNDLE_NNCBOXSUPPORTFUNCTION_HXX

//------------------------------------------------------------

#include <map>
#include "MatrixCBSolver.hxx"
#include "NNCBoxSupportModification.hxx"

//------------------------------------------------------------

namespace ConicBundle {


  class NNCBoxSupportFunction;

  class NNCBoxSupportMinorantExtender: public CBout, public MinorantExtender
  {
  private:
    NNCBoxSupportFunction* fun;

  public:
    NNCBoxSupportMinorantExtender(NNCBoxSupportFunction* fun);
                        
    ~NNCBoxSupportMinorantExtender();

    /*@brief called by ConicBundle to update internal Minorant objects, has to return 0 on success

        for each relevant index i the minorant value is set to the projection 
        of 0 onto the interval [lower_bound(i),upper_bound(i)]

        @param[in,out] minorant  (Minorant&)
            it holds a (possibly aggregated) minorant that was generated
            from minorants returned by oracle calls, e.g. as in 
            FunctionOracle::evaluate() 

        @param[in] n_coords (int)
            the number of coordinate positions that have to be filled in

        @param[out] new_subgradient_values  (DVector &)
            the indices of these coordinate positions (sorted in
            strictly increasing order)
      
        @return
           -  0 on success,
           -  1 if extension/update is impossible 
      
     */ 
    int extend(Minorant& minorant,int n_coords,const int* indices);

  };



  class NNCBoxSupportFunction: public MatrixFunctionOracle, public CBout
  {
  private:
    CH_Matrix_Classes::Matrix lb; 
    CH_Matrix_Classes::Matrix ub; 
    
    int apply_modification(const NNCBoxSupportModification& mod);
    
  public:
    NNCBoxSupportFunction(const CH_Matrix_Classes::Matrix& lb,
                       const CH_Matrix_Classes::Matrix& ub,
                       const CBout* cb=0,int incr=-1);    

    ~NNCBoxSupportFunction(){}


    //----------- Oracle Implementation of MatrixFunctionOracle ----------

    int evaluate(const  CH_Matrix_Classes::Matrix& current_point,
                 CH_Matrix_Classes::Real relprec,
                 CH_Matrix_Classes::Real& objective_value,
                 std::vector<Minorant*>&  minorants,
                 PrimalExtender*& primal_extender);


    virtual 
    int 
    apply_modification
    (
     const OracleModification& oracle_modification ,
     const CH_Matrix_Classes::Matrix* new_center,
     const CH_Matrix_Classes::Matrix* old_center,
     bool& discard_objective_in_center,
     bool& discard_model, 
     bool& discard_aggregates,
     MinorantExtender*& minorant_extender
    );

    
     virtual
    bool
    check_correctness() const
    {return true;}
    

    //------------------  routines for querying the problem ----------------

    const CH_Matrix_Classes::Matrix& get_lower_bounds() { return lb;}
    
    const CH_Matrix_Classes::Matrix& get_upper_bounds() {return ub;}
            


   //------------------  routines for IO ----------------

    void  set_out(std::ostream* o=0,int pril=1);
    void  set_cbout(const CBout* cb,int incr=-1);

    std::ostream& print_problem_data(std::ostream& out) const;

    std::istream& read_problem_data(std::istream& in);

    std::ostream& print_problem_data_to_mfile(std::ostream& out, CH_Matrix_Classes::Integer blocknr) const;
    //std::istream& read_problem_data_from_mfile(std::istream& in);


  };


}


#endif