double EllipsoidalIntegrator::scale( const Interval &E, const Interval &X ) const{

    double TOL,ATOL;
    get( INTEGRATOR_TOLERANCE, TOL  );
    get( ABSOLUTE_TOLERANCE  , ATOL );

    return 0.5*acadoMax( ::fabs(E.l()), ::fabs(E.u()) )/
           ( TOL*acadoMax( ::fabs(X.l()), ::fabs(X.u()) ) + ATOL );
}
bool operator==(Interval i1, Interval i2)
{
    return i1.l() == i2.l() && i1.r() == i2.r();
}
bool operator>(Interval i1, Interval i2)
{
    if(i1.l() > i2.l()) return true;
    if(i1.l() < i2.l()) return false;
    return i1.r() > i2.r();
}