+ runtime check of rounding mode support

For example, following program

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

    int main()
    {
        boost::numeric::interval<double>::traits_type::rounding _;
    }
</code>

will work fine if is built with gcc using '-frounding-math', but will fail otherwise.
Running on Valgrind is another practical case when this program will crash (even if it was correctly built).

Author: AndreyG

include/boost/numeric/interval/hw_rounding.hpp

@@ -11,6 +11,8 @@
 #ifndef BOOST_NUMERIC_INTERVAL_HW_ROUNDING_HPP
 #define BOOST_NUMERIC_INTERVAL_HW_ROUNDING_HPP
 
+#include <cassert>
+
 #include <boost/numeric/interval/rounding.hpp>
 #include <boost/numeric/interval/rounded_arith.hpp>
 
@@ -46,23 +48,50 @@ namespace boost {
 namespace numeric {
 namespace interval_lib {
 
+namespace detail
+{
+  template<class T>
+  struct runtime_checker
+  {
+  protected:
+    runtime_checker()
+    {
+      static helper _;
+    }
+
+  private:
+    struct helper
+    {
+      helper()
+      {
+        T a = -(T(-1.1) * T(10.1));
+        T b = +(T(+1.1) * T(10.1));
+        assert(a != b);
+      }
+    };
+  };
+}
+
 /*
  * Three specializations of rounded_math<T>
  */
 
 template<>
 struct rounded_math<float>
   : save_state<rounded_arith_opp<float> >
+  , private detail::runtime_checker<float>
 {};
 
 template<>
 struct rounded_math<double>
   : save_state<rounded_arith_opp<double> >
+  , private detail::runtime_checker<double>
 {};
 
 template<>
 struct rounded_math<long double>
   : save_state<rounded_arith_opp<long double> >
+  , private detail::runtime_checker<long double>
 {};
 
 } // namespace interval_lib