converting to float128 from double
Hi, I have a function that takes as an argument a vector of doubles and tries to convert it to a vector of float128. When I accumulate the values in the float128 vector it looks like it does double additions and only casts up at the last step. The initial vector is defined as const std::vector<double> The float128 vector is defined as std::vector<float128> epoint(num) And I try to do the casting (inside a loop) via: epoint.at(i)=(float128(e8.at(i))); If I do the accumulation on the original vector with no casting into a float128, I get the same as accumulating the vector that should be of float128 type. How should I be doing the casting? Thanks,
On 03/05/2016 13:44, Cooper, Bridgette R D wrote:
Hi,
I have a function that takes as an argument a vector of doubles and tries to convert it to a vector of float128.
When I accumulate the values in the float128 vector it looks like it does double additions and only casts up at the last step.
The initial vector is defined as conststd::vector<double>
The float128 vector is defined as std::vector<float128> epoint(num)
And I try to do the casting (inside a loop) via: epoint.at(i)=(float128(e8.at(i)));
If I do the accumulation on the original vector with no casting into a float128, I get the same as accumulating the vector that should be of float128 type.
How should I be doing the casting?
Sorry but you'll have to post way more code than that to see where the mistake is. John.
Thanks,
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
Hi John, Thanks for the response. Here are some more snippets from my code. #include <boost/multiprecision/float128.hpp> #include <boost/math/cstdfloat/cstdfloat_types.hpp> #include <boost/math/cstdfloat/cstdfloat_limits.hpp> #include <boost/math/cstdfloat/cstdfloat_cmath.hpp> #include <boost/math/cstdfloat/cstdfloat_iostream.hpp> #include <boost/multiprecision/detail/float_string_cvt.hpp> #include <boost/multiprecision/detail/generic_interconvert.hpp> extern "C"{ #include <quadmath.h> } #include <iostream> #include <sstream> #include <vector> #include <iterator> using namespace boost::multi precision; double stieltjes_gamma(int num, const std::vector<double> & e8, const std::vector<double> & g8, int printlevel ) { std::vector<float128> epoint(num), gpoint(num), bcoef(num); for (int i=0; i < num; i++) { epoint.at(i)=(float128(e8.at(i))); gpoint.at(i)=(float128(g8.at(i))); } bcoef.at(0)=0.0Q; for (int i = 0; i < num; i++) { bcoef.at(0)+=gpoint.at(i) } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0)<< std::endl; return(bcoef.at(0)) } When I compare this to fortran code that does the same I don't see the same value for bcoef.at(0). 0.0345217724606016853085032291117978629 <- from above c++ code 3.452175962187407162907246287483402E-0002 <- Fortran correctly utilising quad precision on same input vector The input vectors in this case are 7976 so errors in precision are really quickly apparrant. If I print the gpoint vector, it gets filled with the same garbage for additional precision from double as is for the fortran. Thanks, ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Tuesday, May 3, 2016 5:06:47 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double On 03/05/2016 13:44, Cooper, Bridgette R D wrote:
Hi,
I have a function that takes as an argument a vector of doubles and tries to convert it to a vector of float128.
When I accumulate the values in the float128 vector it looks like it does double additions and only casts up at the last step.
The initial vector is defined as conststd::vector<double>
The float128 vector is defined as std::vector<float128> epoint(num)
And I try to do the casting (inside a loop) via: epoint.at(i)=(float128(e8.at(i)));
If I do the accumulation on the original vector with no casting into a float128, I get the same as accumulating the vector that should be of float128 type.
How should I be doing the casting?
Sorry but you'll have to post way more code than that to see where the mistake is. John.
Thanks,
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
On 03/05/2016 17:35, Cooper, Bridgette R D wrote:
Hi John,
Thanks for the response.
Here are some more snippets from my code.
I don't see any obvious errors, but in any case you must have extracted that from something else but it doesn't compile! Not least you can't implicitly return from a float128 to a double in return(bcoef.at(0)) So if that's compiling for you, then I suspect your declaration of bcoef. If I re-write and simplify in a slightly more C++ idiom how does it compare to your fortran code now: double stieltjes_gamma(const std::vector<double> & e8, const std::vector<double> & g8, int printlevel) { assert(e8.size() == g8.size()); // precondition? std::vector<float128> epoint(e8.begin(), e8.end()), gpoint(g8.begin(), g8.end()), bcoef(e8.size()); bcoef[0] = 0.0Q; for(int i = 0; i < e8.size(); i++) { bcoef[0] += gpoint[i]; } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0) << std::endl; return static_cast<double>(bcoef[0]); } And finally... note that summing at higher precision can not actually protect you from cancellation error (in any language) since the error is inherent in the input values. HTH, John.
#include <boost/multiprecision/float128.hpp> #include <boost/math/cstdfloat/cstdfloat_types.hpp> #include <boost/math/cstdfloat/cstdfloat_limits.hpp> #include <boost/math/cstdfloat/cstdfloat_cmath.hpp> #include <boost/math/cstdfloat/cstdfloat_iostream.hpp> #include <boost/multiprecision/detail/float_string_cvt.hpp> #include <boost/multiprecision/detail/generic_interconvert.hpp> extern "C"{ #include <quadmath.h> } #include <iostream> #include <sstream> #include <vector> #include <iterator>
using namespace boost::multi precision;
double stieltjes_gamma(int num, const std::vector<double> & e8, const std::vector<double> & g8, int printlevel ) {
std::vector<float128> epoint(num), gpoint(num), bcoef(num);
for (int i=0; i < num; i++) { epoint.at(i)=(float128(e8.at(i))); gpoint.at(i)=(float128(g8.at(i))); }
bcoef.at(0)=0.0Q; for (int i = 0; i < num; i++) { bcoef.at(0)+=gpoint.at(i) } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0)<< std::endl; return(bcoef.at(0)) }
When I compare this to fortran code that does the same I don't see the same value for bcoef.at(0).
0.0345217724606016853085032291117978629 <- from above c++ code 3.452175962187407162907246287483402E-0002 <- Fortran correctly utilising quad precision on same input vector
The input vectors in this case are 7976 so errors in precision are really quickly apparrant. If I print the gpoint vector, it gets filled with the same garbage for additional precision from double as is for the fortran.
Thanks, ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Tuesday, May 3, 2016 5:06:47 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
On 03/05/2016 13:44, Cooper, Bridgette R D wrote:
Hi,
I have a function that takes as an argument a vector of doubles and tries to convert it to a vector of float128.
When I accumulate the values in the float128 vector it looks like it does double additions and only casts up at the last step.
The initial vector is defined as conststd::vector<double>
The float128 vector is defined as std::vector<float128> epoint(num)
And I try to do the casting (inside a loop) via: epoint.at(i)=(float128(e8.at(i)));
If I do the accumulation on the original vector with no casting into a float128, I get the same as accumulating the vector that should be of float128 type.
How should I be doing the casting?
Sorry but you'll have to post way more code than that to see where the mistake is.
John.
Thanks,
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
Hi John, No, my code is several hundreds of lines long. I extracted what I thought would be quite self contained, but yeah you are right about the return type mis-match. I am not sure that this helps. What I really want is to do a lot of computation in float128 precision, and for some reason when I cast up a double vector to a float128 vector all the operations on the float128 vector are happening a double precision, and the cast up only happens at the final stage. I've changed the initiation of epoint and gpoint vectors to the same as your code snippet and still the same problem. So for example, if I accumulate the original double vector without casting the values to float128 type and accumulate into a float128 variable, the cast up only happens to the result. This is the same number that I see as the result of accumulating the float128vector values. Hope this clarifies things a bit ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Tuesday, May 3, 2016 6:27:39 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double On 03/05/2016 17:35, Cooper, Bridgette R D wrote:
Hi John,
Thanks for the response.
Here are some more snippets from my code.
I don't see any obvious errors, but in any case you must have extracted that from something else but it doesn't compile! Not least you can't implicitly return from a float128 to a double in return(bcoef.at(0)) So if that's compiling for you, then I suspect your declaration of bcoef. If I re-write and simplify in a slightly more C++ idiom how does it compare to your fortran code now: double stieltjes_gamma(const std::vector<double> & e8, const std::vector<double> & g8, int printlevel) { assert(e8.size() == g8.size()); // precondition? std::vector<float128> epoint(e8.begin(), e8.end()), gpoint(g8.begin(), g8.end()), bcoef(e8.size()); bcoef[0] = 0.0Q; for(int i = 0; i < e8.size(); i++) { bcoef[0] += gpoint[i]; } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0) << std::endl; return static_cast<double>(bcoef[0]); } And finally... note that summing at higher precision can not actually protect you from cancellation error (in any language) since the error is inherent in the input values. HTH, John.
#include <boost/multiprecision/float128.hpp> #include <boost/math/cstdfloat/cstdfloat_types.hpp> #include <boost/math/cstdfloat/cstdfloat_limits.hpp> #include <boost/math/cstdfloat/cstdfloat_cmath.hpp> #include <boost/math/cstdfloat/cstdfloat_iostream.hpp> #include <boost/multiprecision/detail/float_string_cvt.hpp> #include <boost/multiprecision/detail/generic_interconvert.hpp> extern "C"{ #include <quadmath.h> } #include <iostream> #include <sstream> #include <vector> #include <iterator>
using namespace boost::multi precision;
double stieltjes_gamma(int num, const std::vector<double> & e8, const std::vector<double> & g8, int printlevel ) {
std::vector<float128> epoint(num), gpoint(num), bcoef(num);
for (int i=0; i < num; i++) { epoint.at(i)=(float128(e8.at(i))); gpoint.at(i)=(float128(g8.at(i))); }
bcoef.at(0)=0.0Q; for (int i = 0; i < num; i++) { bcoef.at(0)+=gpoint.at(i) } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0)<< std::endl; return(bcoef.at(0)) }
When I compare this to fortran code that does the same I don't see the same value for bcoef.at(0).
0.0345217724606016853085032291117978629 <- from above c++ code 3.452175962187407162907246287483402E-0002 <- Fortran correctly utilising quad precision on same input vector
The input vectors in this case are 7976 so errors in precision are really quickly apparrant. If I print the gpoint vector, it gets filled with the same garbage for additional precision from double as is for the fortran.
Thanks, ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Tuesday, May 3, 2016 5:06:47 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
On 03/05/2016 13:44, Cooper, Bridgette R D wrote:
Hi,
I have a function that takes as an argument a vector of doubles and tries to convert it to a vector of float128.
When I accumulate the values in the float128 vector it looks like it does double additions and only casts up at the last step.
The initial vector is defined as conststd::vector<double>
The float128 vector is defined as std::vector<float128> epoint(num)
And I try to do the casting (inside a loop) via: epoint.at(i)=(float128(e8.at(i)));
If I do the accumulation on the original vector with no casting into a float128, I get the same as accumulating the vector that should be of float128 type.
How should I be doing the casting?
Sorry but you'll have to post way more code than that to see where the mistake is.
John.
Thanks,
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
On 03/05/2016 20:58, Cooper, Bridgette R D wrote:
Hi John,
No, my code is several hundreds of lines long. I extracted what I thought would be quite self contained, but yeah you are right about the return type mis-match.
I am not sure that this helps. What I really want is to do a lot of computation in float128 precision, and for some reason when I cast up a double vector to a float128 vector all the operations on the float128 vector are happening a double precision, and the cast up only happens at the final stage.
I've changed the initiation of epoint and gpoint vectors to the same as your code snippet and still the same problem.
So for example, if I accumulate the original double vector without casting the values to float128 type and accumulate into a float128 variable, the cast up only happens to the result. This is the same number that I see as the result of accumulating the float128vector values.
Hope this clarifies things a bit
Not really :( I think you're going to have to try to produce a self contained test case - my guess is you'll spot the error in the process of doing that? Best, John.
________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Tuesday, May 3, 2016 6:27:39 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
Hi John,
Thanks for the response.
Here are some more snippets from my code. I don't see any obvious errors, but in any case you must have extracted
On 03/05/2016 17:35, Cooper, Bridgette R D wrote: that from something else but it doesn't compile!
Not least you can't implicitly return from a float128 to a double in
return(bcoef.at(0))
So if that's compiling for you, then I suspect your declaration of bcoef.
If I re-write and simplify in a slightly more C++ idiom how does it compare to your fortran code now:
double stieltjes_gamma(const std::vector<double> & e8, const std::vector<double> & g8, int printlevel) { assert(e8.size() == g8.size()); // precondition? std::vector<float128> epoint(e8.begin(), e8.end()), gpoint(g8.begin(), g8.end()), bcoef(e8.size());
bcoef[0] = 0.0Q; for(int i = 0; i < e8.size(); i++) { bcoef[0] += gpoint[i]; } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0) << std::endl; return static_cast<double>(bcoef[0]); }
And finally... note that summing at higher precision can not actually protect you from cancellation error (in any language) since the error is inherent in the input values.
HTH, John.
#include <boost/multiprecision/float128.hpp> #include <boost/math/cstdfloat/cstdfloat_types.hpp> #include <boost/math/cstdfloat/cstdfloat_limits.hpp> #include <boost/math/cstdfloat/cstdfloat_cmath.hpp> #include <boost/math/cstdfloat/cstdfloat_iostream.hpp> #include <boost/multiprecision/detail/float_string_cvt.hpp> #include <boost/multiprecision/detail/generic_interconvert.hpp> extern "C"{ #include <quadmath.h> } #include <iostream> #include <sstream> #include <vector> #include <iterator>
using namespace boost::multi precision;
double stieltjes_gamma(int num, const std::vector<double> & e8, const std::vector<double> & g8, int printlevel ) {
std::vector<float128> epoint(num), gpoint(num), bcoef(num);
for (int i=0; i < num; i++) { epoint.at(i)=(float128(e8.at(i))); gpoint.at(i)=(float128(g8.at(i))); }
bcoef.at(0)=0.0Q; for (int i = 0; i < num; i++) { bcoef.at(0)+=gpoint.at(i) } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0)<< std::endl; return(bcoef.at(0)) }
When I compare this to fortran code that does the same I don't see the same value for bcoef.at(0).
0.0345217724606016853085032291117978629 <- from above c++ code 3.452175962187407162907246287483402E-0002 <- Fortran correctly utilising quad precision on same input vector
The input vectors in this case are 7976 so errors in precision are really quickly apparrant. If I print the gpoint vector, it gets filled with the same garbage for additional precision from double as is for the fortran.
Thanks, ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Tuesday, May 3, 2016 5:06:47 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
On 03/05/2016 13:44, Cooper, Bridgette R D wrote:
Hi,
I have a function that takes as an argument a vector of doubles and tries to convert it to a vector of float128.
When I accumulate the values in the float128 vector it looks like it does double additions and only casts up at the last step.
The initial vector is defined as conststd::vector<double>
The float128 vector is defined as std::vector<float128> epoint(num)
And I try to do the casting (inside a loop) via: epoint.at(i)=(float128(e8.at(i)));
If I do the accumulation on the original vector with no casting into a float128, I get the same as accumulating the vector that should be of float128 type.
How should I be doing the casting?
Sorry but you'll have to post way more code than that to see where the mistake is.
John.
Thanks,
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
Hi John, I tried again. I compile pretty much exactly the function you wrote: #include <boost/multiprecision/float128.hpp> #include <boost/math/cstdfloat/cstdfloat_types.hpp> #include <boost/math/cstdfloat/cstdfloat_limits.hpp> #include <boost/math/cstdfloat/cstdfloat_cmath.hpp> #include <boost/math/cstdfloat/cstdfloat_iostream.hpp> extern "C"{ #include <quadmath.h> } #include <iostream> #include <fstream> #include <sstream> #include <vector> #include <iterator> using namespace std; using namespace boost::multiprecision; double stieltjes_gamma(const std::vector<double> & e8, const std::vector<double> & g8, int printlevel) { assert(e8.size() == g8.size()); // precondition? std::vector<float128> epoint(e8.begin(), e8.end()), gpoint(g8.begin(), g8.end()), bcoef(e8.size()); bcoef[0] = 0.0Q; for(int i = 0; i < e8.size(); i++) { bcoef[0] += gpoint[i]; } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0) << std::endl; return (0.0); } The answer is still: 0.0345217724606016853085032291117978629 Whereas the fortran gives the answer to be: 3.452175962187407162907246287483402E-0002 I really appreciate you taking the time to respond to my emails. I really don't understand the discrepancy. I have also tried things like: epoint.at(i)=static_cast<float128>(e8.at(i)); But still cannot reproduce the fortran value. I also have code that can reproduce the quad precision sqrt to match fortran code. Firstly, I have a c++ routine that calculates the sqrt (a^2+b^2) in quad precision without using the library. This reproduces the quad value from fortran, and the quad value of sqrt function using boost library. However, if the explicit type of a or b is not appended with Q I don't reproduce the correct answer, which I can understand. But what if I wanted to start with a double value of a and b, cast them up to new variables c d in float128, and want the square root of the quad numbers. How do I do this casting so that I don't get just the result casted up to a quad precision number, the sqrtq function is called? Thanks again for all your help. ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Wednesday, May 4, 2016 12:29:31 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double On 03/05/2016 20:58, Cooper, Bridgette R D wrote:
Hi John,
No, my code is several hundreds of lines long. I extracted what I thought would be quite self contained, but yeah you are right about the return type mis-match.
I am not sure that this helps. What I really want is to do a lot of computation in float128 precision, and for some reason when I cast up a double vector to a float128 vector all the operations on the float128 vector are happening a double precision, and the cast up only happens at the final stage.
I've changed the initiation of epoint and gpoint vectors to the same as your code snippet and still the same problem.
So for example, if I accumulate the original double vector without casting the values to float128 type and accumulate into a float128 variable, the cast up only happens to the result. This is the same number that I see as the result of accumulating the float128vector values.
Hope this clarifies things a bit
Not really :( I think you're going to have to try to produce a self contained test case - my guess is you'll spot the error in the process of doing that? Best, John.
________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Tuesday, May 3, 2016 6:27:39 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
Hi John,
Thanks for the response.
Here are some more snippets from my code. I don't see any obvious errors, but in any case you must have extracted
On 03/05/2016 17:35, Cooper, Bridgette R D wrote: that from something else but it doesn't compile!
Not least you can't implicitly return from a float128 to a double in
return(bcoef.at(0))
So if that's compiling for you, then I suspect your declaration of bcoef.
If I re-write and simplify in a slightly more C++ idiom how does it compare to your fortran code now:
double stieltjes_gamma(const std::vector<double> & e8, const std::vector<double> & g8, int printlevel) { assert(e8.size() == g8.size()); // precondition? std::vector<float128> epoint(e8.begin(), e8.end()), gpoint(g8.begin(), g8.end()), bcoef(e8.size());
bcoef[0] = 0.0Q; for(int i = 0; i < e8.size(); i++) { bcoef[0] += gpoint[i]; } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0) << std::endl; return static_cast<double>(bcoef[0]); }
And finally... note that summing at higher precision can not actually protect you from cancellation error (in any language) since the error is inherent in the input values.
HTH, John.
#include <boost/multiprecision/float128.hpp> #include <boost/math/cstdfloat/cstdfloat_types.hpp> #include <boost/math/cstdfloat/cstdfloat_limits.hpp> #include <boost/math/cstdfloat/cstdfloat_cmath.hpp> #include <boost/math/cstdfloat/cstdfloat_iostream.hpp> #include <boost/multiprecision/detail/float_string_cvt.hpp> #include <boost/multiprecision/detail/generic_interconvert.hpp> extern "C"{ #include <quadmath.h> } #include <iostream> #include <sstream> #include <vector> #include <iterator>
using namespace boost::multi precision;
double stieltjes_gamma(int num, const std::vector<double> & e8, const std::vector<double> & g8, int printlevel ) {
std::vector<float128> epoint(num), gpoint(num), bcoef(num);
for (int i=0; i < num; i++) { epoint.at(i)=(float128(e8.at(i))); gpoint.at(i)=(float128(g8.at(i))); }
bcoef.at(0)=0.0Q; for (int i = 0; i < num; i++) { bcoef.at(0)+=gpoint.at(i) } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0)<< std::endl; return(bcoef.at(0)) }
When I compare this to fortran code that does the same I don't see the same value for bcoef.at(0).
0.0345217724606016853085032291117978629 <- from above c++ code 3.452175962187407162907246287483402E-0002 <- Fortran correctly utilising quad precision on same input vector
The input vectors in this case are 7976 so errors in precision are really quickly apparrant. If I print the gpoint vector, it gets filled with the same garbage for additional precision from double as is for the fortran.
Thanks, ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Tuesday, May 3, 2016 5:06:47 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
On 03/05/2016 13:44, Cooper, Bridgette R D wrote:
Hi,
I have a function that takes as an argument a vector of doubles and tries to convert it to a vector of float128.
When I accumulate the values in the float128 vector it looks like it does double additions and only casts up at the last step.
The initial vector is defined as conststd::vector<double>
The float128 vector is defined as std::vector<float128> epoint(num)
And I try to do the casting (inside a loop) via: epoint.at(i)=(float128(e8.at(i)));
If I do the accumulation on the original vector with no casting into a float128, I get the same as accumulating the vector that should be of float128 type.
How should I be doing the casting?
Sorry but you'll have to post way more code than that to see where the mistake is.
John.
Thanks,
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
On 04/05/2016 13:47, Cooper, Bridgette R D wrote:
Hi John,
I tried again. I compile pretty much exactly the function you wrote:
Bridgette, without the input values I still can't go about reproducing here - I need a *complete* test case plus the matching FORTRAN code. John.
#include <boost/multiprecision/float128.hpp> #include <boost/math/cstdfloat/cstdfloat_types.hpp> #include <boost/math/cstdfloat/cstdfloat_limits.hpp> #include <boost/math/cstdfloat/cstdfloat_cmath.hpp> #include <boost/math/cstdfloat/cstdfloat_iostream.hpp> extern "C"{ #include <quadmath.h> } #include <iostream> #include <fstream> #include <sstream> #include <vector> #include <iterator>
using namespace std; using namespace boost::multiprecision; double stieltjes_gamma(const std::vector<double> & e8, const std::vector<double> & g8, int printlevel) { assert(e8.size() == g8.size()); // precondition? std::vector<float128> epoint(e8.begin(), e8.end()), gpoint(g8.begin(), g8.end()), bcoef(e8.size());
bcoef[0] = 0.0Q; for(int i = 0; i < e8.size(); i++) { bcoef[0] += gpoint[i]; } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0) << std::endl; return (0.0); }
The answer is still:
0.0345217724606016853085032291117978629
Whereas the fortran gives the answer to be:
3.452175962187407162907246287483402E-0002
I really appreciate you taking the time to respond to my emails. I really don't understand the discrepancy.
I have also tried things like:
epoint.at(i)=static_cast<float128>(e8.at(i));
But still cannot reproduce the fortran value.
I also have code that can reproduce the quad precision sqrt to match fortran code. Firstly, I have a c++ routine that calculates the sqrt (a^2+b^2) in quad precision without using the library. This reproduces the quad value from fortran, and the quad value of sqrt function using boost library. However, if the explicit type of a or b is not appended with Q I don't reproduce the correct answer, which I can understand. But what if I wanted to start with a double value of a and b, cast them up to new variables c d in float128, and want the square root of the quad numbers. How do I do this casting so that I don't get just the result casted up to a quad precision number, the sqrtq function is called?
Thanks again for all your help. ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Wednesday, May 4, 2016 12:29:31 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
On 03/05/2016 20:58, Cooper, Bridgette R D wrote:
Hi John,
No, my code is several hundreds of lines long. I extracted what I thought would be quite self contained, but yeah you are right about the return type mis-match.
I am not sure that this helps. What I really want is to do a lot of computation in float128 precision, and for some reason when I cast up a double vector to a float128 vector all the operations on the float128 vector are happening a double precision, and the cast up only happens at the final stage.
I've changed the initiation of epoint and gpoint vectors to the same as your code snippet and still the same problem.
So for example, if I accumulate the original double vector without casting the values to float128 type and accumulate into a float128 variable, the cast up only happens to the result. This is the same number that I see as the result of accumulating the float128vector values.
Hope this clarifies things a bit Not really :(
I think you're going to have to try to produce a self contained test case - my guess is you'll spot the error in the process of doing that?
Best, John.
________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Tuesday, May 3, 2016 6:27:39 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
Hi John,
Thanks for the response.
Here are some more snippets from my code. I don't see any obvious errors, but in any case you must have extracted
On 03/05/2016 17:35, Cooper, Bridgette R D wrote: that from something else but it doesn't compile!
Not least you can't implicitly return from a float128 to a double in
return(bcoef.at(0))
So if that's compiling for you, then I suspect your declaration of bcoef.
If I re-write and simplify in a slightly more C++ idiom how does it compare to your fortran code now:
double stieltjes_gamma(const std::vector<double> & e8, const std::vector<double> & g8, int printlevel) { assert(e8.size() == g8.size()); // precondition? std::vector<float128> epoint(e8.begin(), e8.end()), gpoint(g8.begin(), g8.end()), bcoef(e8.size());
bcoef[0] = 0.0Q; for(int i = 0; i < e8.size(); i++) { bcoef[0] += gpoint[i]; } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0) << std::endl; return static_cast<double>(bcoef[0]); }
And finally... note that summing at higher precision can not actually protect you from cancellation error (in any language) since the error is inherent in the input values.
HTH, John.
#include <boost/multiprecision/float128.hpp> #include <boost/math/cstdfloat/cstdfloat_types.hpp> #include <boost/math/cstdfloat/cstdfloat_limits.hpp> #include <boost/math/cstdfloat/cstdfloat_cmath.hpp> #include <boost/math/cstdfloat/cstdfloat_iostream.hpp> #include <boost/multiprecision/detail/float_string_cvt.hpp> #include <boost/multiprecision/detail/generic_interconvert.hpp> extern "C"{ #include <quadmath.h> } #include <iostream> #include <sstream> #include <vector> #include <iterator>
using namespace boost::multi precision;
double stieltjes_gamma(int num, const std::vector<double> & e8, const std::vector<double> & g8, int printlevel ) {
std::vector<float128> epoint(num), gpoint(num), bcoef(num);
for (int i=0; i < num; i++) { epoint.at(i)=(float128(e8.at(i))); gpoint.at(i)=(float128(g8.at(i))); }
bcoef.at(0)=0.0Q; for (int i = 0; i < num; i++) { bcoef.at(0)+=gpoint.at(i) } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0)<< std::endl; return(bcoef.at(0)) }
When I compare this to fortran code that does the same I don't see the same value for bcoef.at(0).
0.0345217724606016853085032291117978629 <- from above c++ code 3.452175962187407162907246287483402E-0002 <- Fortran correctly utilising quad precision on same input vector
The input vectors in this case are 7976 so errors in precision are really quickly apparrant. If I print the gpoint vector, it gets filled with the same garbage for additional precision from double as is for the fortran.
Thanks, ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Tuesday, May 3, 2016 5:06:47 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
On 03/05/2016 13:44, Cooper, Bridgette R D wrote:
Hi,
I have a function that takes as an argument a vector of doubles and tries to convert it to a vector of float128.
When I accumulate the values in the float128 vector it looks like it does double additions and only casts up at the last step.
The initial vector is defined as conststd::vector<double>
The float128 vector is defined as std::vector<float128> epoint(num)
And I try to do the casting (inside a loop) via: epoint.at(i)=(float128(e8.at(i)));
If I do the accumulation on the original vector with no casting into a float128, I get the same as accumulating the vector that should be of float128 type.
How should I be doing the casting?
Sorry but you'll have to post way more code than that to see where the mistake is.
John.
Thanks,
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
program run_stieltjes implicit none integer num,printflag,unitout,io integer np,i parameter (np=15000) real*8 e8(np),g8(np),gamma0 io=0 num=0 open(unit=68,file='coup_VQZ+5jp_iv1.dat',status='old', readonly) do i=1,15000 read(68,*,iostat=io) e8(i),g8(i) if (io > 0) then write(6,*)'Check input file of correct format' exit else if (io < 0) then write(6,*)'Reach the end of file' write(6,*)'Number of energy points = ',num exit else num=num+1 end if end do close(68) open(unit=65,file='Stieltjes_VQZ+5jp_iv2.dat') unitout=65 printflag=1 call stieltjes(num,e8,g8,gamma0,65,printflag) write(6,*)'End of program' write(6,*)'Number of energy points = ',num close(65) end subroutine stieltjes (num,e8,g8,gamma0,unitout,printflag) implicit none integer :: num,nmax,np,npol,unitout,printflag,i parameter (np=15000,nmax=1000) real*16 :: e_point(np),g_point(np) real*16 :: bcoef(0:nmax) real*8 :: e8(np),g8(np),gamma0 do i=1,num e_point(i)=e8(i) g_point(i)=g8(i) if (printflag.ne.0) write(6,*) e_point(i),g_point(i),i end do write(6,*)' ' bcoef(0)=0.q0 do i=1,num bcoef(0)=bcoef(0)+g_point(i) end do write(6,*) bcoef(0) return end This compiles with the intel compiler, but not gfortran unfortunately. ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Wednesday, May 4, 2016 4:49:07 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double On 04/05/2016 13:47, Cooper, Bridgette R D wrote:
Hi John,
I tried again. I compile pretty much exactly the function you wrote:
Bridgette, without the input values I still can't go about reproducing here - I need a *complete* test case plus the matching FORTRAN code. John.
#include <boost/multiprecision/float128.hpp> #include <boost/math/cstdfloat/cstdfloat_types.hpp> #include <boost/math/cstdfloat/cstdfloat_limits.hpp> #include <boost/math/cstdfloat/cstdfloat_cmath.hpp> #include <boost/math/cstdfloat/cstdfloat_iostream.hpp> extern "C"{ #include <quadmath.h> } #include <iostream> #include <fstream> #include <sstream> #include <vector> #include <iterator>
using namespace std; using namespace boost::multiprecision; double stieltjes_gamma(const std::vector<double> & e8, const std::vector<double> & g8, int printlevel) { assert(e8.size() == g8.size()); // precondition? std::vector<float128> epoint(e8.begin(), e8.end()), gpoint(g8.begin(), g8.end()), bcoef(e8.size());
bcoef[0] = 0.0Q; for(int i = 0; i < e8.size(); i++) { bcoef[0] += gpoint[i]; } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0) << std::endl; return (0.0); }
The answer is still:
0.0345217724606016853085032291117978629
Whereas the fortran gives the answer to be:
3.452175962187407162907246287483402E-0002
I really appreciate you taking the time to respond to my emails. I really don't understand the discrepancy.
I have also tried things like:
epoint.at(i)=static_cast<float128>(e8.at(i));
But still cannot reproduce the fortran value.
I also have code that can reproduce the quad precision sqrt to match fortran code. Firstly, I have a c++ routine that calculates the sqrt (a^2+b^2) in quad precision without using the library. This reproduces the quad value from fortran, and the quad value of sqrt function using boost library. However, if the explicit type of a or b is not appended with Q I don't reproduce the correct answer, which I can understand. But what if I wanted to start with a double value of a and b, cast them up to new variables c d in float128, and want the square root of the quad numbers. How do I do this casting so that I don't get just the result casted up to a quad precision number, the sqrtq function is called?
Thanks again for all your help. ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Wednesday, May 4, 2016 12:29:31 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
On 03/05/2016 20:58, Cooper, Bridgette R D wrote:
Hi John,
No, my code is several hundreds of lines long. I extracted what I thought would be quite self contained, but yeah you are right about the return type mis-match.
I am not sure that this helps. What I really want is to do a lot of computation in float128 precision, and for some reason when I cast up a double vector to a float128 vector all the operations on the float128 vector are happening a double precision, and the cast up only happens at the final stage.
I've changed the initiation of epoint and gpoint vectors to the same as your code snippet and still the same problem.
So for example, if I accumulate the original double vector without casting the values to float128 type and accumulate into a float128 variable, the cast up only happens to the result. This is the same number that I see as the result of accumulating the float128vector values.
Hope this clarifies things a bit Not really :(
I think you're going to have to try to produce a self contained test case - my guess is you'll spot the error in the process of doing that?
Best, John.
________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Tuesday, May 3, 2016 6:27:39 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
Hi John,
Thanks for the response.
Here are some more snippets from my code. I don't see any obvious errors, but in any case you must have extracted
On 03/05/2016 17:35, Cooper, Bridgette R D wrote: that from something else but it doesn't compile!
Not least you can't implicitly return from a float128 to a double in
return(bcoef.at(0))
So if that's compiling for you, then I suspect your declaration of bcoef.
If I re-write and simplify in a slightly more C++ idiom how does it compare to your fortran code now:
double stieltjes_gamma(const std::vector<double> & e8, const std::vector<double> & g8, int printlevel) { assert(e8.size() == g8.size()); // precondition? std::vector<float128> epoint(e8.begin(), e8.end()), gpoint(g8.begin(), g8.end()), bcoef(e8.size());
bcoef[0] = 0.0Q; for(int i = 0; i < e8.size(); i++) { bcoef[0] += gpoint[i]; } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0) << std::endl; return static_cast<double>(bcoef[0]); }
And finally... note that summing at higher precision can not actually protect you from cancellation error (in any language) since the error is inherent in the input values.
HTH, John.
#include <boost/multiprecision/float128.hpp> #include <boost/math/cstdfloat/cstdfloat_types.hpp> #include <boost/math/cstdfloat/cstdfloat_limits.hpp> #include <boost/math/cstdfloat/cstdfloat_cmath.hpp> #include <boost/math/cstdfloat/cstdfloat_iostream.hpp> #include <boost/multiprecision/detail/float_string_cvt.hpp> #include <boost/multiprecision/detail/generic_interconvert.hpp> extern "C"{ #include <quadmath.h> } #include <iostream> #include <sstream> #include <vector> #include <iterator>
using namespace boost::multi precision;
double stieltjes_gamma(int num, const std::vector<double> & e8, const std::vector<double> & g8, int printlevel ) {
std::vector<float128> epoint(num), gpoint(num), bcoef(num);
for (int i=0; i < num; i++) { epoint.at(i)=(float128(e8.at(i))); gpoint.at(i)=(float128(g8.at(i))); }
bcoef.at(0)=0.0Q; for (int i = 0; i < num; i++) { bcoef.at(0)+=gpoint.at(i) } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0)<< std::endl; return(bcoef.at(0)) }
When I compare this to fortran code that does the same I don't see the same value for bcoef.at(0).
0.0345217724606016853085032291117978629 <- from above c++ code 3.452175962187407162907246287483402E-0002 <- Fortran correctly utilising quad precision on same input vector
The input vectors in this case are 7976 so errors in precision are really quickly apparrant. If I print the gpoint vector, it gets filled with the same garbage for additional precision from double as is for the fortran.
Thanks, ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Tuesday, May 3, 2016 5:06:47 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
On 03/05/2016 13:44, Cooper, Bridgette R D wrote:
Hi,
I have a function that takes as an argument a vector of doubles and tries to convert it to a vector of float128.
When I accumulate the values in the float128 vector it looks like it does double additions and only casts up at the last step.
The initial vector is defined as conststd::vector<double>
The float128 vector is defined as std::vector<float128> epoint(num)
And I try to do the casting (inside a loop) via: epoint.at(i)=(float128(e8.at(i)));
If I do the accumulation on the original vector with no casting into a float128, I get the same as accumulating the vector that should be of float128 type.
How should I be doing the casting?
Sorry but you'll have to post way more code than that to see where the mistake is.
John.
Thanks,
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
Hello, changing line 10 (the 'open' statement) to open(unit=68,file='coup_VQZ+5jp_iv1.dat',status='old',action='read') the program compiles fine for me using gfortran4.9.2. An input file 'coup_VQZ+5jp_iv1.dat' of the form 1 1 1 1 1 1 ... ... will allow the program to complete. For 15'000 e/g pairs the output will be ... ... ... 1.00000000000000000000000000000000000 1.00000000000000000000000000000000000 14999 1.00000000000000000000000000000000000 1.00000000000000000000000000000000000 15000 15000.0000000000000000000000000000000 End of program Number of energy points = 15000 Best Regards, Max On Wed, May 04, 2016 at 04:29:11PM +0000, Cooper, Bridgette R D wrote:
program run_stieltjes implicit none integer num,printflag,unitout,io integer np,i parameter (np=15000) real*8 e8(np),g8(np),gamma0
io=0 num=0 open(unit=68,file='coup_VQZ+5jp_iv1.dat',status='old', readonly) do i=1,15000 read(68,*,iostat=io) e8(i),g8(i) if (io > 0) then write(6,*)'Check input file of correct format' exit else if (io < 0) then write(6,*)'Reach the end of file' write(6,*)'Number of energy points = ',num exit else num=num+1 end if end do close(68) open(unit=65,file='Stieltjes_VQZ+5jp_iv2.dat') unitout=65 printflag=1 call stieltjes(num,e8,g8,gamma0,65,printflag) write(6,*)'End of program' write(6,*)'Number of energy points = ',num close(65) end
subroutine stieltjes (num,e8,g8,gamma0,unitout,printflag) implicit none integer :: num,nmax,np,npol,unitout,printflag,i parameter (np=15000,nmax=1000) real*16 :: e_point(np),g_point(np) real*16 :: bcoef(0:nmax) real*8 :: e8(np),g8(np),gamma0
do i=1,num e_point(i)=e8(i) g_point(i)=g8(i) if (printflag.ne.0) write(6,*) e_point(i),g_point(i),i end do write(6,*)' ' bcoef(0)=0.q0 do i=1,num bcoef(0)=bcoef(0)+g_point(i) end do write(6,*) bcoef(0) return end
This compiles with the intel compiler, but not gfortran unfortunately. ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Wednesday, May 4, 2016 4:49:07 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
On 04/05/2016 13:47, Cooper, Bridgette R D wrote:
Hi John,
I tried again. I compile pretty much exactly the function you wrote:
Bridgette, without the input values I still can't go about reproducing here - I need a *complete* test case plus the matching FORTRAN code.
John.
#include <boost/multiprecision/float128.hpp> #include <boost/math/cstdfloat/cstdfloat_types.hpp> #include <boost/math/cstdfloat/cstdfloat_limits.hpp> #include <boost/math/cstdfloat/cstdfloat_cmath.hpp> #include <boost/math/cstdfloat/cstdfloat_iostream.hpp> extern "C"{ #include <quadmath.h> } #include <iostream> #include <fstream> #include <sstream> #include <vector> #include <iterator>
using namespace std; using namespace boost::multiprecision; double stieltjes_gamma(const std::vector<double> & e8, const std::vector<double> & g8, int printlevel) { assert(e8.size() == g8.size()); // precondition? std::vector<float128> epoint(e8.begin(), e8.end()), gpoint(g8.begin(), g8.end()), bcoef(e8.size());
bcoef[0] = 0.0Q; for(int i = 0; i < e8.size(); i++) { bcoef[0] += gpoint[i]; } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0) << std::endl; return (0.0); }
The answer is still:
0.0345217724606016853085032291117978629
Whereas the fortran gives the answer to be:
3.452175962187407162907246287483402E-0002
I really appreciate you taking the time to respond to my emails. I really don't understand the discrepancy.
I have also tried things like:
epoint.at(i)=static_cast<float128>(e8.at(i));
But still cannot reproduce the fortran value.
I also have code that can reproduce the quad precision sqrt to match fortran code. Firstly, I have a c++ routine that calculates the sqrt (a^2+b^2) in quad precision without using the library. This reproduces the quad value from fortran, and the quad value of sqrt function using boost library. However, if the explicit type of a or b is not appended with Q I don't reproduce the correct answer, which I can understand. But what if I wanted to start with a double value of a and b, cast them up to new variables c d in float128, and want the square root of the quad numbers. How do I do this casting so that I don't get just the result casted up to a quad precision number, the sqrtq function is called?
Thanks again for all your help. ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Wednesday, May 4, 2016 12:29:31 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
On 03/05/2016 20:58, Cooper, Bridgette R D wrote:
Hi John,
No, my code is several hundreds of lines long. I extracted what I thought would be quite self contained, but yeah you are right about the return type mis-match.
I am not sure that this helps. What I really want is to do a lot of computation in float128 precision, and for some reason when I cast up a double vector to a float128 vector all the operations on the float128 vector are happening a double precision, and the cast up only happens at the final stage.
I've changed the initiation of epoint and gpoint vectors to the same as your code snippet and still the same problem.
So for example, if I accumulate the original double vector without casting the values to float128 type and accumulate into a float128 variable, the cast up only happens to the result. This is the same number that I see as the result of accumulating the float128vector values.
Hope this clarifies things a bit Not really :(
I think you're going to have to try to produce a self contained test case - my guess is you'll spot the error in the process of doing that?
Best, John.
________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Tuesday, May 3, 2016 6:27:39 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
Hi John,
Thanks for the response.
Here are some more snippets from my code. I don't see any obvious errors, but in any case you must have extracted
On 03/05/2016 17:35, Cooper, Bridgette R D wrote: that from something else but it doesn't compile!
Not least you can't implicitly return from a float128 to a double in
return(bcoef.at(0))
So if that's compiling for you, then I suspect your declaration of bcoef.
If I re-write and simplify in a slightly more C++ idiom how does it compare to your fortran code now:
double stieltjes_gamma(const std::vector<double> & e8, const std::vector<double> & g8, int printlevel) { assert(e8.size() == g8.size()); // precondition? std::vector<float128> epoint(e8.begin(), e8.end()), gpoint(g8.begin(), g8.end()), bcoef(e8.size());
bcoef[0] = 0.0Q; for(int i = 0; i < e8.size(); i++) { bcoef[0] += gpoint[i]; } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0) << std::endl; return static_cast<double>(bcoef[0]); }
And finally... note that summing at higher precision can not actually protect you from cancellation error (in any language) since the error is inherent in the input values.
HTH, John.
#include <boost/multiprecision/float128.hpp> #include <boost/math/cstdfloat/cstdfloat_types.hpp> #include <boost/math/cstdfloat/cstdfloat_limits.hpp> #include <boost/math/cstdfloat/cstdfloat_cmath.hpp> #include <boost/math/cstdfloat/cstdfloat_iostream.hpp> #include <boost/multiprecision/detail/float_string_cvt.hpp> #include <boost/multiprecision/detail/generic_interconvert.hpp> extern "C"{ #include <quadmath.h> } #include <iostream> #include <sstream> #include <vector> #include <iterator>
using namespace boost::multi precision;
double stieltjes_gamma(int num, const std::vector<double> & e8, const std::vector<double> & g8, int printlevel ) {
std::vector<float128> epoint(num), gpoint(num), bcoef(num);
for (int i=0; i < num; i++) { epoint.at(i)=(float128(e8.at(i))); gpoint.at(i)=(float128(g8.at(i))); }
bcoef.at(0)=0.0Q; for (int i = 0; i < num; i++) { bcoef.at(0)+=gpoint.at(i) } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0)<< std::endl; return(bcoef.at(0)) }
When I compare this to fortran code that does the same I don't see the same value for bcoef.at(0).
0.0345217724606016853085032291117978629 <- from above c++ code 3.452175962187407162907246287483402E-0002 <- Fortran correctly utilising quad precision on same input vector
The input vectors in this case are 7976 so errors in precision are really quickly apparrant. If I print the gpoint vector, it gets filled with the same garbage for additional precision from double as is for the fortran.
Thanks, ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Tuesday, May 3, 2016 5:06:47 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
On 03/05/2016 13:44, Cooper, Bridgette R D wrote:
Hi,
I have a function that takes as an argument a vector of doubles and tries to convert it to a vector of float128.
When I accumulate the values in the float128 vector it looks like it does double additions and only casts up at the last step.
The initial vector is defined as conststd::vector<double>
The float128 vector is defined as std::vector<float128> epoint(num)
And I try to do the casting (inside a loop) via: epoint.at(i)=(float128(e8.at(i)));
If I do the accumulation on the original vector with no casting into a float128, I get the same as accumulating the vector that should be of float128 type.
How should I be doing the casting?
Sorry but you'll have to post way more code than that to see where the mistake is.
John.
Thanks,
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
On 04/05/2016 17:29, Cooper, Bridgette R D wrote:
program run_stieltjes implicit none integer num,printflag,unitout,io integer np,i parameter (np=15000) real*8 e8(np),g8(np),gamma0
io=0 num=0 open(unit=68,file='coup_VQZ+5jp_iv1.dat',status='old', readonly) do i=1,15000
I must be missing some data then, because I see 7976 pairs of entries in that file not 15000. If I sum all the entries in the second column (from 1.813373152916055E-005 down to 1.283872761367943E-008) using: template <class T> T sum_as_type(const std::vector<double>& v) { T val = 0; for(unsigned i = 0; i < v.size(); ++i) val += T(v[i]); return val; } int main() { using namespace boost::multiprecision; std::ifstream is("M:\\data\\boost\\boost\\coup_VQZ+5jp_iv1.dat"); std::vector<double> data; while(is.good()) { double a, b; is >> a; is >> b; is >> std::ws; data.push_back(b); } std::cout << data.back() << std::endl; std::cout << data.size() << std::endl; std::cout << std::setprecision(35); std::cout << sum_as_type<double>(data) << std::endl; std::cout << float128(sum_as_type<__float128>(data)) << std::endl; std::cout << sum_as_type<float128>(data) << std::endl; std::cout << sum_as_type<cpp_bin_float_quad>(data) << std::endl; std::cout << sum_as_type<mpfi_float_100>(data) << std::endl; return 0; } Then I see as output: 0.034521759621874245627104471623169957 0.034521759621874071629072462874834025 0.034521759621874071629072462874834025 0.034521759621874071629072462874834025 0.034521759621874071629072462874834025 So both cpp_bin_float, float128, __float128 and mpfi all agree, and the last line indicates the result is exact (no interval). HTH, John.
read(68,*,iostat=io) e8(i),g8(i) if (io > 0) then write(6,*)'Check input file of correct format' exit else if (io < 0) then write(6,*)'Reach the end of file' write(6,*)'Number of energy points = ',num exit else num=num+1 end if end do close(68) open(unit=65,file='Stieltjes_VQZ+5jp_iv2.dat') unitout=65 printflag=1 call stieltjes(num,e8,g8,gamma0,65,printflag) write(6,*)'End of program' write(6,*)'Number of energy points = ',num close(65) end
subroutine stieltjes (num,e8,g8,gamma0,unitout,printflag) implicit none integer :: num,nmax,np,npol,unitout,printflag,i parameter (np=15000,nmax=1000) real*16 :: e_point(np),g_point(np) real*16 :: bcoef(0:nmax) real*8 :: e8(np),g8(np),gamma0
do i=1,num e_point(i)=e8(i) g_point(i)=g8(i) if (printflag.ne.0) write(6,*) e_point(i),g_point(i),i end do write(6,*)' ' bcoef(0)=0.q0 do i=1,num bcoef(0)=bcoef(0)+g_point(i) end do write(6,*) bcoef(0) return end
This compiles with the intel compiler, but not gfortran unfortunately. ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Wednesday, May 4, 2016 4:49:07 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
On 04/05/2016 13:47, Cooper, Bridgette R D wrote:
Hi John,
I tried again. I compile pretty much exactly the function you wrote: Bridgette, without the input values I still can't go about reproducing here - I need a *complete* test case plus the matching FORTRAN code.
John.
#include <boost/multiprecision/float128.hpp> #include <boost/math/cstdfloat/cstdfloat_types.hpp> #include <boost/math/cstdfloat/cstdfloat_limits.hpp> #include <boost/math/cstdfloat/cstdfloat_cmath.hpp> #include <boost/math/cstdfloat/cstdfloat_iostream.hpp> extern "C"{ #include <quadmath.h> } #include <iostream> #include <fstream> #include <sstream> #include <vector> #include <iterator>
using namespace std; using namespace boost::multiprecision; double stieltjes_gamma(const std::vector<double> & e8, const std::vector<double> & g8, int printlevel) { assert(e8.size() == g8.size()); // precondition? std::vector<float128> epoint(e8.begin(), e8.end()), gpoint(g8.begin(), g8.end()), bcoef(e8.size());
bcoef[0] = 0.0Q; for(int i = 0; i < e8.size(); i++) { bcoef[0] += gpoint[i]; } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0) << std::endl; return (0.0); }
The answer is still:
0.0345217724606016853085032291117978629
Whereas the fortran gives the answer to be:
3.452175962187407162907246287483402E-0002
I really appreciate you taking the time to respond to my emails. I really don't understand the discrepancy.
I have also tried things like:
epoint.at(i)=static_cast<float128>(e8.at(i));
But still cannot reproduce the fortran value.
I also have code that can reproduce the quad precision sqrt to match fortran code. Firstly, I have a c++ routine that calculates the sqrt (a^2+b^2) in quad precision without using the library. This reproduces the quad value from fortran, and the quad value of sqrt function using boost library. However, if the explicit type of a or b is not appended with Q I don't reproduce the correct answer, which I can understand. But what if I wanted to start with a double value of a and b, cast them up to new variables c d in float128, and want the square root of the quad numbers. How do I do this casting so that I don't get just the result casted up to a quad precision number, the sqrtq function is called?
Thanks again for all your help. ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Wednesday, May 4, 2016 12:29:31 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
On 03/05/2016 20:58, Cooper, Bridgette R D wrote:
Hi John,
No, my code is several hundreds of lines long. I extracted what I thought would be quite self contained, but yeah you are right about the return type mis-match.
I am not sure that this helps. What I really want is to do a lot of computation in float128 precision, and for some reason when I cast up a double vector to a float128 vector all the operations on the float128 vector are happening a double precision, and the cast up only happens at the final stage.
I've changed the initiation of epoint and gpoint vectors to the same as your code snippet and still the same problem.
So for example, if I accumulate the original double vector without casting the values to float128 type and accumulate into a float128 variable, the cast up only happens to the result. This is the same number that I see as the result of accumulating the float128vector values.
Hope this clarifies things a bit Not really :(
I think you're going to have to try to produce a self contained test case - my guess is you'll spot the error in the process of doing that?
Best, John.
________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Tuesday, May 3, 2016 6:27:39 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
Hi John,
Thanks for the response.
Here are some more snippets from my code. I don't see any obvious errors, but in any case you must have extracted
On 03/05/2016 17:35, Cooper, Bridgette R D wrote: that from something else but it doesn't compile!
Not least you can't implicitly return from a float128 to a double in
return(bcoef.at(0))
So if that's compiling for you, then I suspect your declaration of bcoef.
If I re-write and simplify in a slightly more C++ idiom how does it compare to your fortran code now:
double stieltjes_gamma(const std::vector<double> & e8, const std::vector<double> & g8, int printlevel) { assert(e8.size() == g8.size()); // precondition? std::vector<float128> epoint(e8.begin(), e8.end()), gpoint(g8.begin(), g8.end()), bcoef(e8.size());
bcoef[0] = 0.0Q; for(int i = 0; i < e8.size(); i++) { bcoef[0] += gpoint[i]; } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0) << std::endl; return static_cast<double>(bcoef[0]); }
And finally... note that summing at higher precision can not actually protect you from cancellation error (in any language) since the error is inherent in the input values.
HTH, John.
#include <boost/multiprecision/float128.hpp> #include <boost/math/cstdfloat/cstdfloat_types.hpp> #include <boost/math/cstdfloat/cstdfloat_limits.hpp> #include <boost/math/cstdfloat/cstdfloat_cmath.hpp> #include <boost/math/cstdfloat/cstdfloat_iostream.hpp> #include <boost/multiprecision/detail/float_string_cvt.hpp> #include <boost/multiprecision/detail/generic_interconvert.hpp> extern "C"{ #include <quadmath.h> } #include <iostream> #include <sstream> #include <vector> #include <iterator>
using namespace boost::multi precision;
double stieltjes_gamma(int num, const std::vector<double> & e8, const std::vector<double> & g8, int printlevel ) {
std::vector<float128> epoint(num), gpoint(num), bcoef(num);
for (int i=0; i < num; i++) { epoint.at(i)=(float128(e8.at(i))); gpoint.at(i)=(float128(g8.at(i))); }
bcoef.at(0)=0.0Q; for (int i = 0; i < num; i++) { bcoef.at(0)+=gpoint.at(i) } std::cout << std::setprecision(std::numeric_limits<float128>::max_digits10) << bcoef.at(0)<< std::endl; return(bcoef.at(0)) }
When I compare this to fortran code that does the same I don't see the same value for bcoef.at(0).
0.0345217724606016853085032291117978629 <- from above c++ code 3.452175962187407162907246287483402E-0002 <- Fortran correctly utilising quad precision on same input vector
The input vectors in this case are 7976 so errors in precision are really quickly apparrant. If I print the gpoint vector, it gets filled with the same garbage for additional precision from double as is for the fortran.
Thanks, ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of John Maddock <jz.maddock@googlemail.com> Sent: Tuesday, May 3, 2016 5:06:47 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
On 03/05/2016 13:44, Cooper, Bridgette R D wrote:
Hi,
I have a function that takes as an argument a vector of doubles and tries to convert it to a vector of float128.
When I accumulate the values in the float128 vector it looks like it does double additions and only casts up at the last step.
The initial vector is defined as conststd::vector<double>
The float128 vector is defined as std::vector<float128> epoint(num)
And I try to do the casting (inside a loop) via: epoint.at(i)=(float128(e8.at(i)));
If I do the accumulation on the original vector with no casting into a float128, I get the same as accumulating the vector that should be of float128 type.
How should I be doing the casting?
Sorry but you'll have to post way more code than that to see where the mistake is.
John.
Thanks,
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
On Wed, May 04, 2016 at 07:18:46PM +0100, John Maddock wrote:
I must be missing some data then, because I see 7976 pairs of entries in that file not 15000.
I have to apologize, I didn't see there was an input file attached. My less than two cents' worth were only about the compilation anyway.
No worries Max thanks for your help. I still can't get gfortran to compile it. I used the following options and fixed the action=read gfortran -std=legacy -lquadmath *.f -o test, but it fails with stieltjes.f:(.text+0x15c): undefined reference to `_gfortran_transfer_real128_write' John, there are only 7976 pairs, you aren't missing data. The last four values you have also match my fortran code. The problem I have is that I don't even have the first number for my c++ code I have 0.0345217724606016853085032291117978629 for adding the same vector. I will have a play around with what you have sent and see if I can reproduce the numbers. Thank you again for all your help John. It is really appreciated. ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of Maximilian Winkler <mbw_bst@mailbox.org> Sent: Wednesday, May 4, 2016 7:35:41 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double On Wed, May 04, 2016 at 07:18:46PM +0100, John Maddock wrote:
I must be missing some data then, because I see 7976 pairs of entries in that file not 15000.
I have to apologize, I didn't see there was an input file attached. My less than two cents' worth were only about the compilation anyway. _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
Ok I didn't know you were forced to write "FORTRAN", but of course I know that's how things work in quantum chemistry. However, using the flags you provided, the program still compiles for me. Of course the open statement has to be something like this: open(unit=68,file='coup_VQZ+5jp_iv1.dat',status='old', +action='read') Regards, Max On Wed, May 04, 2016 at 07:45:38PM +0000, Cooper, Bridgette R D wrote:
No worries Max thanks for your help. I still can't get gfortran to compile it. I used the following options and fixed the action=read gfortran -std=legacy -lquadmath *.f -o test, but it fails with stieltjes.f:(.text+0x15c): undefined reference to `_gfortran_transfer_real128_write'
John, there are only 7976 pairs, you aren't missing data. The last four values you have also match my fortran code. The problem I have is that I don't even have the first number for my c++ code I have 0.0345217724606016853085032291117978629 for adding the same vector.
I will have a play around with what you have sent and see if I can reproduce the numbers. Thank you again for all your help John. It is really appreciated. ________________________________________ From: Boost-users <boost-users-bounces@lists.boost.org> on behalf of Maximilian Winkler <mbw_bst@mailbox.org> Sent: Wednesday, May 4, 2016 7:35:41 PM To: boost-users@lists.boost.org Subject: Re: [Boost-users] converting to float128 from double
On Wed, May 04, 2016 at 07:18:46PM +0100, John Maddock wrote:
I must be missing some data then, because I see 7976 pairs of entries in that file not 15000.
I have to apologize, I didn't see there was an input file attached. My less than two cents' worth were only about the compilation anyway. _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users _______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users
participants (1)
-
Cooper, Bridgette R D -
John Maddock -
Maximilian Winkler -
Maximilian Winkler