1 files changed, 42 insertions, 0 deletions
diff --git a/src/boost/libs/multiprecision/example/mpfi_snips.cpp b/src/boost/libs/multiprecision/example/mpfi_snips.cpp
new file mode 100644
index 00000000..b36e22b5
--- /dev/null
+++ b/src/boost/libs/multiprecision/example/mpfi_snips.cpp
@@ -0,0 +1,42 @@
+///////////////////////////////////////////////////////////////
+//  Copyright 2011 John Maddock. Distributed under the Boost
+//  Software License, Version 1.0. (See accompanying file
+//  LICENSE_1_0.txt or copy at https://www.boost.org/LICENSE_1_0.txt
+
+//[mpfi_eg
+#include <boost/multiprecision/mpfi.hpp>
+#include <boost/math/special_functions/gamma.hpp>
+#include <iostream>
+
+int main()
+{
+   using namespace boost::multiprecision;
+
+   // Operations at variable precision and no numeric_limits support:
+   mpfi_float a = 2;
+   mpfi_float::default_precision(1000);
+   std::cout << mpfi_float::default_precision() << std::endl;
+   std::cout << sqrt(a) << std::endl; // print root-2
+
+   // Operations at fixed precision and full numeric_limits support:
+   mpfi_float_100 b = 2;
+   std::cout << std::numeric_limits<mpfi_float_100>::digits << std::endl;
+   // We can use any C++ std lib function:
+   std::cout << log(b) << std::endl; // print log(2)
+
+   // Access the underlying data:
+   mpfi_t r;
+   mpfi_init(r);
+   mpfi_set(r, b.backend().data());
+
+   // Construct some explicit intervals and perform set operations:
+   mpfi_float_50 i1(1, 2), i2(1.5, 2.5);
+   std::cout << intersect(i1, i2) << std::endl;
+   std::cout << hull(i1, i2) << std::endl;
+   std::cout << overlap(i1, i2) << std::endl;
+   std::cout << subset(i1, i2) << std::endl;
+   mpfi_clear(r);
+   return 0;
+}
+//]
+