codac_QInterProjF.h

//============================================================================
//                                  I B E X
// File        : Q-intersection contracto
// Author      : Clément Carbonnel, Benoit Desrochers
// Copyright   : Ecole des Mines de Nantes (France)
// License     : See the LICENSE file
// Created     : Apr 30, 2012
// Last Update : Apr 30, 2012
//============================================================================
 
#ifndef __IBEX_CTC_Q_INTER_2_H__
#define __IBEX_CTC_Q_INTER_2_H__
 
#include "ibex_Ctc.h"
#include "ibex_Sep.h"
#include "ibex_Array.h"
#include "ibex_IntervalMatrix.h"
 
using ibex::IntervalVector;
using ibex::IntervalMatrix;
using ibex::Array;
using ibex::Ctc;
using ibex::Sep;
using ibex::Interval;
 
 
 
namespace codac {
 
 
IntervalVector qinter_projf(const Array<IntervalVector>& _boxes, int q);
 
class CtcQInterProjF : public Ctc {
public:
    CtcQInterProjF(const Array<Ctc>& list, int q);
 
    virtual void contract(IntervalVector& box);
 
    Array<Ctc> list;
 
    int q;
 
protected:
    IntervalMatrix boxes; 
};
 
 
class SepQInterProjF : public Sep {
public:
    SepQInterProjF(const Array<Sep>& list, int q = 0);
 
 
    virtual void separate(IntervalVector& xin, IntervalVector& xout);
 
 
    void set_q(int q);
 
    int get_q();
 
protected:
 
    Array<Sep> list;
 
    IntervalMatrix boxes_in;
    IntervalMatrix boxes_out;
 
    int q;
 
};
 
/* ============================================================================
                             inline implementation
  ============================================================================*/
 
inline CtcQInterProjF::CtcQInterProjF(const Array<Ctc>& list, int q) :
        Ctc(list), list(list), q(q), boxes(list.size(), nb_var) { }
 
inline SepQInterProjF::SepQInterProjF(const Array<Sep>& list, int q) :
        Sep(list[0].nb_var),
        list(list),
        boxes_in(list.size(), list[0].nb_var),
        boxes_out(list.size(), list[0].nb_var)
    { this->set_q(q); }
 
 
inline void SepQInterProjF::set_q(int q){
    assert (q >= 0 || q < list.size());
    this->q = q;
}
 
inline int SepQInterProjF::get_q(){ return q; }
 
 
 
} // end namespace pyibex
#endif // __IBEX_CTC_Q_INTER_2_H__