Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

boost::numeric::interval constructor seems not creating an interval enclosing its arguments #981

Open
Shengyis opened this issue Nov 16, 2024 · 0 comments

Comments

@Shengyis
Copy link

In the document of boost interval https://www.boost.org/doc/libs/1_86_0/libs/numeric/interval/doc/interval.htm#interval
If the interval constructor has only one argument, "The constructors create an interval enclosing their arguments."
But in the following code, seems not working like this, it just sets interval's upper and lower bounds as the same value of my input. This is NOT what we expect for interval arithmetic, since the input may not be precisely stored in a computer.

#include <iostream>
#include <boost/numeric/interval.hpp>

using namespace std;

int main() {
    cout << setprecision(20);
    boost::numeric::interval<double> I = boost::numeric::interval<double>(0.1);
    cout << I.upper() << endl; // print: 0.10000000000000000555
    cout << I.lower() << endl; // print: 0.10000000000000000555
}

another strange thing is when I add a double number, sometimes it gives the desired bound, but sometimes NOT! See below, adding 0.1 to the interval and adding 1.1 to the interval gives different type of the result. adding 1.1 gives the desired result.

#include <iostream>
#include <boost/numeric/interval.hpp>

using namespace std;

int main() {
    cout << setprecision(20);
    boost::numeric::interval<double> I = boost::numeric::interval<double>(0.1);
    I += 0.1 // add 0.1
    cout << I.upper() << endl; // print: 0.2000000000000000111
    cout << I.lower() << endl; // print: 0.2000000000000000111
}
#include <iostream>
#include <boost/numeric/interval.hpp>

using namespace std;

int main() {
    cout << setprecision(20);
    boost::numeric::interval<double> I = boost::numeric::interval<double>(0.1);
    I += 1.1 // add 1.1
    cout << I.upper() << endl; // print: 1.2000000000000001776
    cout << I.lower() << endl; // print: 1.1999999999999999556
}

When compiling I add -frounding-math, and the boost version is v1.85.0.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant