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(); }