Ticket #6047: ticket6047_round_trunc_01242013.diff

trunc.hpp

@@ -16 +16 @@
 
 namespace boost{ namespace math{
 
+template < class T, class TargetType,  class Policy >
+inline TargetType trunc_cast(T const& x, Policy const& pol)
+{
+BOOST_STATIC_ASSERT(std::numeric_limits<TargetType>::is_specialized);  
+T r = trunc(x);
+if ((r> (std::numeric_limits<TargetType>::max)()) || 
+    (r < (std::numeric_limits<TargetType>::min)()))
+  return policies::raise_rounding_error
+  ("boost::math::round_cast<%1%,TargetType>(%%1%%)>",
+   "Value %1% exceeds the numeric limits",
+   r,
+   0,
+   pol);
+  return static_cast<TargetType> (r);
+  
+}
+  
+template <class T, class TargetType>
+inline TargetType trunc_cast(T const& x)  
+{
+  typedef typename policies::policy<
+  boost::math::policies::rounding_error<policies::throw_on_error>
+  > Policy;
+  return trunc_cast<T, TargetType, Policy>(x, Policy());
+}
+
+  
 template <class T, class Policy>
 inline T trunc(const T& v, const Policy& pol)
 {

round.hpp

@@ -14 +14 @@
 #include <boost/math/policies/error_handling.hpp>
 #include <boost/math/special_functions/fpclassify.hpp>
 
+
 namespace boost{ namespace math{
 
+
+template < class T, class TargetType,  class Policy >
+inline TargetType round_cast(T const& x, Policy const& pol)
+{
+  BOOST_STATIC_ASSERT(std::numeric_limits<TargetType>::is_specialized);  
+  T r = round(x);
+  if ((r> (std::numeric_limits<TargetType>::max)()) || 
+      (r < (std::numeric_limits<TargetType>::min)()))
+    return policies::raise_rounding_error
+      ("boost::math::round_cast<%1%,TargetType>(%%1%%)>",
+        "Value %1% exceeds the numeric limits",
+        r,
+        0,
+        pol);
+  return static_cast<TargetType> (r);
+
+}
+
+template <class T, class TargetType>
+inline TargetType round_cast(T const& x)  
+{
+  typedef typename policies::policy<
+  boost::math::policies::rounding_error<policies::throw_on_error>
+  > Policy;
+  return round_cast<T, TargetType, Policy>(x, Policy());
+}
+  
 template <class T, class Policy>
 inline T round(const T& v, const Policy& pol)
 {