1 files changed, 90 insertions, 0 deletions
diff --git a/src/boost/libs/numeric/interval/examples/transc.cpp b/src/boost/libs/numeric/interval/examples/transc.cpp
new file mode 100644
index 00000000..2036f983
--- /dev/null
+++ b/src/boost/libs/numeric/interval/examples/transc.cpp
@@ -0,0 +1,90 @@
+/* Boost example/transc.cpp
+ * how to use an external library (GMP/MPFR in this case) in order to
+ * get correct transcendental functions on intervals
+ *
+ * Copyright 2003 Guillaume Melquiond
+ *
+ * Distributed under the Boost Software License, Version 1.0.
+ * (See accompanying file LICENSE_1_0.txt or
+ * copy at http://www.boost.org/LICENSE_1_0.txt)
+ */
+
+#include <boost/numeric/interval.hpp>
+//extern "C" {
+#include <gmp.h>
+#include <mpfr.h>
+//}
+#include <iostream>
+
+struct full_rounding:
+  boost::numeric::interval_lib::rounded_arith_opp<double>
+{
+private:
+  typedef int mpfr_func(mpfr_t, const __mpfr_struct*, mp_rnd_t);
+  double invoke_mpfr(double x, mpfr_func f, mp_rnd_t r) {
+    mpfr_t xx;
+    mpfr_init_set_d(xx, x, GMP_RNDN);
+    f(xx, xx, r);
+    double res = mpfr_get_d(xx, r);
+    mpfr_clear(xx);
+    return res;
+  }
+public:
+# define GENR_FUNC(name) \
+  double name##_down(double x) { return invoke_mpfr(x, mpfr_##name, GMP_RNDD); } \
+  double name##_up  (double x) { return invoke_mpfr(x, mpfr_##name, GMP_RNDU); }
+  GENR_FUNC(exp)
+  GENR_FUNC(log)
+  GENR_FUNC(sin)
+  GENR_FUNC(cos)
+  GENR_FUNC(tan)
+  GENR_FUNC(asin)
+  GENR_FUNC(acos)
+  GENR_FUNC(atan)
+  GENR_FUNC(sinh)
+  GENR_FUNC(cosh)
+  GENR_FUNC(tanh)
+  GENR_FUNC(asinh)
+  GENR_FUNC(acosh)
+  GENR_FUNC(atanh)
+};
+
+namespace dummy {
+  using namespace boost;
+  using namespace numeric;
+  using namespace interval_lib;
+  typedef save_state<full_rounding> R;
+  typedef checking_strict<double> P;
+  typedef interval<double, policies<R, P> > I;
+};
+
+typedef dummy::I I;
+
+template<class os_t>
+os_t& operator<<(os_t &os, const I &a) {
+  os << '[' << a.lower() << ',' << a.upper() << ']';
+  return os;
+}
+
+int main() {
+  I x(0.5, 2.5);
+  std::cout << "x = " << x << std::endl;
+  std::cout.precision(16);
+# define GENR_TEST(name) \
+  std::cout << #name "(x) = " << name(x) << std::endl
+  GENR_TEST(exp);
+  GENR_TEST(log);
+  GENR_TEST(sin);
+  GENR_TEST(cos);
+  GENR_TEST(tan);
+  GENR_TEST(asin);
+  GENR_TEST(acos);
+  GENR_TEST(atan);
+  GENR_TEST(sinh);
+  GENR_TEST(cosh);
+  GENR_TEST(tanh);
+  GENR_TEST(asinh);
+  GENR_TEST(acosh);
+  GENR_TEST(atanh);
+  return 0;
+}