id,summary,reporter,owner,description,type,status,milestone,component,version,severity,resolution,keywords,cc
448,interval evaluation bug,nobody,Guillaume Melquiond,"{{{
evaluation of sin of [0,1] with interval template 
specification

typedef boost::numeric::interval<
		double, 
		boost::numeric::interval_lib::policies<		
			boost::numeric::interval_lib::
save_state<
				boost::numeric::interval_lib::
rounded_transc_exact<double, 
					boost::numeric::
interval_lib::rounded_arith_opp<double> 
				> 
			>,
			boost::numeric::interval_lib::
checking_strict<double> 
			//boost::numeric::interval_lib::
checking_base<double>
		>
	> Interval;


and rint definition in ms Windows

static double inline rint( double x) 
// Copyright (C) 2001 Tor M. Aamodt, University of 
Toronto 
// Permisssion to use for all purposes commercial and 
otherwise granted. 
// THIS MATERIAL IS PROVIDED ""AS IS"" WITHOUT 
WARRANTY, OR ANY CONDITION OR 
// OTHER TERM OF ANY KIND INCLUDING, WITHOUT 
LIMITATION, ANY WARRANTY 
// OF MERCHANTABILITY, SATISFACTORY QUALITY, 
OR FITNESS FOR A PARTICULAR 
// PURPOSE. 
{ 
    if( x > 0 ) { 
        __int64 xint = (__int64) (x+0.5); 
        if( xint % 2 ) { 
            // then we might have an even number... 
            double diff = x - (double)xint; 
            if( diff == -0.5 ) 
                return double(xint-1); 
        } 
        return double(xint); 
    } else { 
        __int64 xint = (__int64) (x-0.5); 
        if( xint % 2 ) { 
            // then we might have an even number... 
            double diff = x - (double)xint; 
            if( diff == 0.5 ) 
                return double(xint+1); 
        } 
        return double(xint); 
    } 
}


gives incorrect result:
so sin([0,1]) = [ 0.841471 , -1.60814e-016 ]

Thats looks incredible!!! How can i fix it out?
}}}",Bugs,closed,,None,None,,Invalid,,