#13591 closed Bugs (obsolete)
maybe a bug in math::constants library
| Reported by: | Owned by: | John Maddock | |
|---|---|---|---|
| Milestone: | Boost 1.67.0 | Component: | math | 
| Version: | Boost 1.67.0 | Severity: | Problem | 
| Keywords: | Cc: | 
Description
The following templated code will generate incorrect result in visual studio community 2017:
#include<iostream>
#include<boost/math/constants/constants.hpp>
// convienent defination of math constants
template<typename T> const T pi = boost::math::constants::pi<T>(); // PI
template<typename T> const T two_pi = boost::math::constants::two_pi<T>(); // 2*PI
// physical constants
template<typename T> const T mu_0 = 4.0*pi<T>*1.0e-4; // permeability of free space in H / km
int main()
{
	using T = double;
	auto w = two_pi<T> * 0.01;
	//auto t1 = pi<T>;
	auto t2 = mu_0<T>;
	std::cout<<"omega="<<w<<std::endl;
	std::cout<<"mu_0="<<t2<<std::endl;
	system("pause");
	return 0;
}
In VS2017, mu_0 will be 0.000. However, if I uncomment
auto t1 = pi<T>
in main, the result is correct. This problem does not happen when compiled with gcc 8.1 in Manjaro Linux. The boost library version is 1.67.0
Change History (2)
comment:1 by , 4 years ago
| Resolution: | → obsolete | 
|---|---|
| Status: | new → closed | 
comment:2 by , 4 years ago
Not sure, but I think this is an order of initialization error, with mu_0 being initialized before pi.
Whatever I can reduce the issue down to:
template <class T>
T get_pi()
{
   return 3.14;
}
template <class T>
T get_2_pi()
{
   return 6.28;
}
template<typename T> const T pi = get_pi<T>(); // PI
template<typename T> const T two_pi = get_2_pi<T>(); // 2*PI
                                                                           // physical constants
template<typename T> const T mu_0 = 4.0*pi<T>*1.0e-4;
int main()
{
   using T = double;
   auto w = two_pi<T> * 0.01;
   //auto t1 = pi<T>;   // uncomment this line and everything works OK.
   auto t2 = mu_0<T>;
   std::cout << "omega=" << w << std::endl;
   std::cout << "mu_0=" << t2 << std::endl;
   system("pause");
   return 0;
}
So it's not our bug if it's a bug at all?
Note: please reply on the linked github issue as this Trac is now closed.
  Note:
 See   TracTickets
 for help on using tickets.
    
Moved to https://github.com/boostorg/math/issues/144