Interval.h

/********************************************************************/
/* Copyright (c) 2015 DG Technologies Inc. and Yuzi Mizuno          */
/* All rights reserved.                                             */
/********************************************************************/
#ifndef _MGInterval_HH_
#define _MGInterval_HH_
#include "mg/EReal.h"
//  MGInterval.h
//  header for class MGInterval
// Forward Declaration
class MGIfstream;
class MGOfstream;


class MG_DLL_DECLR MGInterval{

public:

MG_DLL_DECLR friend bool operator>(const MGEReal&, const MGInterval&);
MG_DLL_DECLR friend MGInterval operator+ (double, const MGInterval& );

MG_DLL_DECLR friend MGInterval operator* (double a, const MGInterval& );

MG_DLL_DECLR friend bool operator<(const MGEReal&, const MGInterval&);
MG_DLL_DECLR friend bool operator>=(const MGEReal&, const MGInterval&);
MG_DLL_DECLR friend bool operator<=(const MGEReal&, const MGInterval&);
MG_DLL_DECLR friend bool operator>(double, const MGInterval& i);
MG_DLL_DECLR friend bool operator<(double, const MGInterval& i);

MG_DLL_DECLR friend std::ostream& operator<< (std::ostream&, const MGInterval&);


MGInterval();

explicit MGInterval(MGINTERVAL_TYPE, double=0.0);

MGInterval(const MGEReal& t1, const MGEReal& t2); 

MGInterval(const MGInterval&, const MGEReal&); 

explicit MGInterval(double t);


const MGEReal& operator() (int i)const;
MGEReal& operator() (int i);
const MGEReal& operator[] (int i)const;
MGEReal& operator[] (int i);

MGInterval operator+ (const MGInterval& ) const;

MGInterval operator+ (double) const;

MGInterval& operator+= (const MGInterval& );

MGInterval& operator+= (double);

MGInterval operator- () const;

MGInterval operator- (const MGInterval& ) const ;

MGInterval operator- (double) const;

MGInterval& operator-= (const MGInterval& );

MGInterval& operator-= (double);

MGInterval operator* (double a) const;

MGInterval& operator*= (double a);

MGInterval operator/ (double a) const; 

MGInterval& operator/= (double a);

MGInterval operator| (const MGInterval&) const ;

MGInterval& operator|= (const MGInterval&);

MGInterval operator& (const MGInterval&) const ;

MGInterval& operator&= (const MGInterval&);


bool operator&& (const MGInterval&)const;

bool operator>> (const MGInterval& intrval2) const;

bool operator>> (const MGEReal& t) const;

bool operator<< (const MGInterval& intrval2) const;

bool operator<< (const MGEReal& t) const;

bool operator== (const MGInterval&) const;

bool operator!= (const MGInterval&) const;

bool operator< (const MGInterval&) const;
bool operator> (const MGInterval&) const;
bool operator<= (const MGInterval&) const;
bool operator>= (const MGInterval&) const;

bool operator< (const MGEReal& ) const;
bool operator> (const MGEReal&) const;
bool operator<= (const MGEReal&) const;
bool operator>= (const MGEReal&) const;

//Return true if m_high<t;
bool operator< (double t) const;

//Return true if m_low>t;
bool operator> (double t) const;


void change_range(
    double t0,  
    double t1   
);

bool empty() const;

void expand(const MGEReal& val);

bool finite() const;

bool finite_above() const;
 
bool finite_below() const;

const MGEReal& high() const {return m_high;};
MGEReal& high() {return m_high;};

double high_point() const;

bool includes(double t)const;

bool infinite() const;

bool infinite_above() const;

bool infinite_below() const;

double interpolate(double t) const;

bool is_null() const{return empty();};

MGEReal length() const;

const MGEReal& low() const {return m_low;};
MGEReal& low() {return m_low;};

double low_point() const;

double mid_point() const;

double relative_error() const;

double round_into_interval(double t)const;

void set_length(double len);

void set_null(double t=0.);
void set_empty(double t=0.){set_null(t);};

void set_high_point(const MGEReal&);
void set_high(const MGEReal& t){set_high_point(t);};

void set_low_point(const MGEReal&);
void set_low(const MGEReal& t){set_low_point(t);};

MGINTERVAL_TYPE type() const;

int dump_size() const;

int dump(MGOfstream& ) const;

int restore(MGIfstream& );

private:


    MGEReal m_high; 
    MGEReal m_low;  

};
 // end of BASE group
#endif