[Boost-commit] svn:boost r66887 - trunk/libs/math/example

Subject: [Boost-commit] svn:boost r66887 - trunk/libs/math/example
From: pbristow_at_[hidden]
Date: 2010-11-30 05:05:17

Author: pbristow
Date: 2010-11-30 05:05:14 EST (Tue, 30 Nov 2010)
New Revision: 66887
URL: http://svn.boost.org/trac/boost/changeset/66887

Log:
Check before using max_digits10
Text files modified: 
   trunk/libs/math/example/geometric_examples.cpp                 |    14 ++++++++++----                          
   trunk/libs/math/example/inverse_chi_squared_example.cpp        |     2 --                                      
   trunk/libs/math/example/inverse_gamma_distribution_example.cpp |    13 +++++++++----                           
   trunk/libs/math/example/inverse_gamma_example.cpp              |    12 ++++++++++--                            
   4 files changed, 29 insertions(+), 12 deletions(-)

Modified: trunk/libs/math/example/geometric_examples.cpp
==============================================================================
--- trunk/libs/math/example/geometric_examples.cpp	(original)
+++ trunk/libs/math/example/geometric_examples.cpp	2010-11-30 05:05:14 EST (Tue, 30 Nov 2010)
@@ -156,11 +156,17 @@
 */
     geometric g05(0.5);  // Probability of success = 0.5 or 50%
     // Output all potentially significant digits for the type, here double.
-#if (_MSC_VER >= 1600)
-    cout.precision(std::numeric_limits<double>::max_digits10);
-#else
-    cout.precision(2 + (boost::math::policies::digits<double, boost::math::policies::policy<> >() * 30103UL) / 100000UL);
+
+#ifdef BOOST_NO_NUMERIC_LIMITS_LOWEST
+  int max_digits10 = 2 + (boost::math::policies::digits<double, boost::math::policies::policy<> >() * 30103UL) / 100000UL;
+  cout << "BOOST_NO_NUMERIC_LIMITS_LOWEST is defined" << endl; 
+#else 
+  int max_digits10 = std::numeric_limits<double>::max_digits10;
 #endif
+  cout << "Show all potentially significant decimal digits std::numeric_limits<double>::max_digits10 = "
+    << max_digits10 << endl; 
+  cout.precision(max_digits10); // 
+
     cout << cdf(g05, 0.0001) << endl; // returns 0.5000346561579232, not exact 0.5.
 /*`To get the R discrete behaviour, you simply need to round with,
 for example, the `floor` function.

Modified: trunk/libs/math/example/inverse_chi_squared_example.cpp
==============================================================================
--- trunk/libs/math/example/inverse_chi_squared_example.cpp	(original)
+++ trunk/libs/math/example/inverse_chi_squared_example.cpp	2010-11-30 05:05:14 EST (Tue, 30 Nov 2010)
@@ -92,8 +92,6 @@
     << max_digits10 << endl; 
   cout.precision(max_digits10); // 
 
-  cout << _MSC_VER << ' ' << _MSC_FULL_VER << ' ' << _CPPLIB_VER << endl;
-
   inverse_chi_squared ichsqdef; // All defaults  - not very useful!
   cout << "default df = " << ichsqdef.degrees_of_freedom()
     << ", default scale = " <<  ichsqdef.scale() << endl; //  default df = 1, default scale = 0.5

Modified: trunk/libs/math/example/inverse_gamma_distribution_example.cpp
==============================================================================
--- trunk/libs/math/example/inverse_gamma_distribution_example.cpp	(original)
+++ trunk/libs/math/example/inverse_gamma_distribution_example.cpp	2010-11-30 05:05:14 EST (Tue, 30 Nov 2010)
@@ -32,10 +32,15 @@
   cout << "Example using Inverse Gamma distribution. " << endl;
   // TODO - awaiting a real example using Bayesian statistics.
 
-  cout.precision(std::numeric_limits<double>::max_digits10); // 
-
-  int i = std::numeric_limits<double>::max_digits10;
-  cout << "std::numeric_limits<double>::max_digits10 = " << i << endl; 
+#ifdef BOOST_NO_NUMERIC_LIMITS_LOWEST
+  int max_digits10 = 2 + (boost::math::policies::digits<double, boost::math::policies::policy<> >() * 30103UL) / 100000UL;
+  cout << "BOOST_NO_NUMERIC_LIMITS_LOWEST is defined" << endl; 
+#else 
+  int max_digits10 = std::numeric_limits<double>::max_digits10;
+#endif
+  cout << "Show all potentially significant decimal digits std::numeric_limits<double>::max_digits10 = "
+    << max_digits10 << endl; 
+  cout.precision(max_digits10); // 
 
   double shape = 1.;
   double scale = 1.;

Modified: trunk/libs/math/example/inverse_gamma_example.cpp
==============================================================================
--- trunk/libs/math/example/inverse_gamma_example.cpp	(original)
+++ trunk/libs/math/example/inverse_gamma_example.cpp	2010-11-30 05:05:14 EST (Tue, 30 Nov 2010)
@@ -28,8 +28,16 @@
 {
   cout << "Example 1 using Inverse Gamma function. " << endl;
 
-  cout.precision(std::numeric_limits<double>::max_digits10);
-  
+  #ifdef BOOST_NO_NUMERIC_LIMITS_LOWEST
+  int max_digits10 = 2 + (boost::math::policies::digits<double, boost::math::policies::policy<> >() * 30103UL) / 100000UL;
+  cout << "BOOST_NO_NUMERIC_LIMITS_LOWEST is defined" << endl; 
+#else 
+  int max_digits10 = std::numeric_limits<double>::max_digits10;
+#endif
+  cout << "Show all potentially significant decimal digits std::numeric_limits<double>::max_digits10 = "
+    << max_digits10 << endl; 
+  cout.precision(max_digits10); // 
+
   double x = 1.;
   double a = 10;