codac_SepFixPoint.h

//============================================================================
//                                  I B E X
// File        : ibex_SepFixPoint.h
// Author      : Benoit Desrochers
// Copyright   : Benoit Desrochers
// License     : See the LICENSE file
// Created     : May 04, 2015
//============================================================================

#ifndef __IBEX_SEP_FIXPOINT_H__
#define __IBEX_SEP_FIXPOINT_H__


#include <ibex_IntervalVector.h>
#include <ibex_LargestFirst.h>
#include <ibex_Sep.h>
#include <vector>

using ibex::Sep;
using ibex::IntervalVector;

// Compute the fixed point of a separator
namespace codac {
class SepFixPoint : public Sep {

public:
    SepFixPoint(Sep& sep, double ratio=default_ratio);

    ~SepFixPoint();

    void separate(IntervalVector &x_in, IntervalVector &x_out);


protected:

    Sep& sep;

    bool impact_cin;
    bool impact_cout;

    IntervalVector *first_cin_boxes;

    IntervalVector *first_cout_boxes;

    int n_in;
    
    int n_out;


    double ratio;

    static const double default_ratio;

private:

    void clearFlags();
    void setCoutFlags(IntervalVector &x_out, IntervalVector &x0);

    void setCinFlags(IntervalVector &x_in, IntervalVector &x0);

    bool reconstruct(IntervalVector &x_in, IntervalVector &x_out, IntervalVector &x_old);

};

} // namespace pyibex

#endif // __IBEX_SEP_FIXPOINT_H__