PSCAffineModification.hxx

Go to the documentation of this file.

#ifndef CONICBUNDLE_PSCAffineMODIFICATION_HXX
#define CONICBUNDLE_PSCAffineMODIFICATION_HXX


#include <map>
#include "CBout.hxx"
#include "CBSolver.hxx"
#include "SparseCoeffmatMatrix.hxx"
#include "ModificationBase.hxx"

namespace ConicBundle {

  class PSCPrimal;




  class PSCAffineModification: public ModificationBase, public OracleModification
{
private:

  bool append_only;


  //first execute the changes in the variables then in the constrains

  //------ changes in variable space
  //first carry out the append step, afterwards the delete/reassign step

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

  CH_Matrix_Classes::Integer var_olddim;
  CH_Matrix_Classes::Integer var_newdim;

  CH_Matrix_Classes::Integer var_append_dim;

  SparseCoeffmatMatrix var_append;

  CH_Matrix_Classes::Indexmatrix* var_del_ind;
  CH_Matrix_Classes::Indexmatrix* var_map_to_old;
  CH_Matrix_Classes::Indexmatrix* var_new_ind;


  //------ changes in the block (the rows)
  //first carry out the append step, afterwards the reassign/delete step

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

  CH_Matrix_Classes::Indexmatrix block_olddim;
  CH_Matrix_Classes::Indexmatrix block_newdim;

  CH_Matrix_Classes::Indexmatrix block_append_dim;
  SparseCoeffmatMatrix offset_append;
  SparseCoeffmatMatrix block_append;
  CH_Matrix_Classes::Indexmatrix* block_del_ind;
  CH_Matrix_Classes::Indexmatrix* block_map_to_old;
  CH_Matrix_Classes::Indexmatrix* block_new_ind;


  //-----------------------------------------------------------------------
  bool reset_primal;
  PSCPrimal* generating_primal;
  bool skip_extension;


public:

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

  virtual ~PSCAffineModification();


  PSCAffineModification(CH_Matrix_Classes::Integer var_olddim,
                  const CH_Matrix_Classes::Indexmatrix& block_olddim,
                  const CBout* cb=0,int incr=0);

  int clear(CH_Matrix_Classes::Integer var_olddim,
            const CH_Matrix_Classes::Indexmatrix& block_olddim);
  
  
  //-----------------------------------------------------------------------
  

  int add_append_vars(CH_Matrix_Classes::Integer append_dim,
                      const SparseCoeffmatMatrix* append_cols);

  int add_reassign_vars(const CH_Matrix_Classes::Indexmatrix& map_to_old);

  int add_delete_vars(const CH_Matrix_Classes::Indexmatrix& del_ind,
                      CH_Matrix_Classes::Indexmatrix& map_to_old);




  int add_append_blocks(const CH_Matrix_Classes::Indexmatrix& append_dim,
                      const SparseCoeffmatMatrix* append_offsets,
                      const SparseCoeffmatMatrix* append_blocks);

  int add_reassign_blocks(const CH_Matrix_Classes::Indexmatrix& map_to_old);

  int add_delete_blocks(const CH_Matrix_Classes::Indexmatrix& del_ind,
                      CH_Matrix_Classes::Indexmatrix& map_to_old);


  int add_reset_generating_primal(PSCPrimal* new_generating_primal);

  int set_skip_extension(bool skip);


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

  int apply_to_PSCAffine(SparseCoeffmatMatrix* offset,SparseCoeffmatMatrix* matrix) const;


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

  bool no_modification() const;
  int set_append_to_old(bool append_only);
  bool append_to_old() const
  {return append_only;}
  bool deleted_variables_are_zero(const CH_Matrix_Classes::Matrix& oldpoint,const SparseCoeffmatMatrix& oldmat) const;
  bool new_variables_are_zero(const CH_Matrix_Classes::Matrix& newpoint,const SparseCoeffmatMatrix& newmat) const;
  bool mapped_variables_are_equal(const CH_Matrix_Classes::Matrix& newpoint,
                                  const CH_Matrix_Classes::Matrix& oldpoint) const;
  bool variable_modifications() const
  {return ((var_append_dim>0)||(var_map_to_old));}
  bool block_modifications() const
  {return ((block_append_dim.dim()>0)||(block_map_to_old));}


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

  CH_Matrix_Classes::Integer old_vardim() const {return var_olddim;} 
  CH_Matrix_Classes::Integer new_vardim() const {return var_newdim;}
  CH_Matrix_Classes::Integer appended_vardim() const {return var_append_dim;} 

  const CH_Matrix_Classes::Indexmatrix& old_blockdim() const {return block_olddim;} 
  const CH_Matrix_Classes::Indexmatrix& new_blockdim() const {return block_newdim;} 
  const CH_Matrix_Classes::Indexmatrix& appended_blockdim() const {return block_append_dim;} 


  const CH_Matrix_Classes::Indexmatrix* map_to_old_variables() const
  {return var_map_to_old;}

  const CH_Matrix_Classes::Indexmatrix* deleted_var_indices() const
  {return var_del_ind;}

  const CH_Matrix_Classes::Indexmatrix* new_var_indices() const
  {return var_new_ind;}

  const CH_Matrix_Classes::Indexmatrix* map_to_old_blocks() const
  {return block_map_to_old;}

  const CH_Matrix_Classes::Indexmatrix* deleted_block_indices() const
  {return block_del_ind;}

  const CH_Matrix_Classes::Indexmatrix* new_block_indices() const
  {return block_new_ind;}

  const SparseCoeffmatMatrix& get_var_append() const
  {return var_append;}
  const SparseCoeffmatMatrix& get_block_append() const
  {return block_append;}
  const SparseCoeffmatMatrix& get_offset_append() const
  {return offset_append;}
  bool get_reset_primal() const
  {return reset_primal;}
  const PSCPrimal* get_generating_primal() const
  {return generating_primal;}

  bool get_skip_extension() const
  {return skip_extension;}



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

  int get_old_vardim() const {return old_vardim();}
  int get_new_vardim() const {return new_vardim();}
  int get_appended_vardim() const {return appended_vardim();}

  const int* get_map_to_old_variables() const
  {return (map_to_old_variables())==0?0:(const int*)(map_to_old_variables()->get_store());}

  int incorporate(const OracleModification& m);

  OracleModification* new_initial_oraclemodification(int old_var_dim) const 
  {return new PSCAffineModification(CH_Matrix_Classes::Integer(old_var_dim),block_olddim,this,0);}
       
  int add_append_variables(int append_dim)
  {return add_append_vars(CH_Matrix_Classes::Integer(append_dim),0);}

    virtual int add_reassign_variables(int new_dim,const int* map_to_old_indices)
  { CH_Matrix_Classes::Indexmatrix maptoold(CH_Matrix_Classes::Integer(new_dim),1,map_to_old_indices); 
    return add_reassign_vars(maptoold);}


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

  void set_cbout(const CBout* out,int incr=-1)
  {CBout::set_cbout(out,incr);var_append.set_cbout(this);block_append.set_cbout(this);}

  void set_out(std::ostream* out=0,int print_level=1)
  {CBout::set_out(out,print_level);set_cbout(this,0);}


};






}

#endif