Add cpp_double_fp and exercise arithmetic_tests

performance/Jamfile.v2

@@ -71,7 +71,8 @@ test-suite performance :
             performance_test_files/test14.cpp  performance_test_files/test31.cpp  performance_test_files/test48.cpp
             performance_test_files/test15.cpp  performance_test_files/test32.cpp  performance_test_files/test49.cpp
             performance_test_files/test16.cpp  performance_test_files/test33.cpp  performance_test_files/test50.cpp
-            performance_test_files/test17.cpp  performance_test_files/test34.cpp  performance_test_files/test51.cpp            
+            performance_test_files/test17.cpp  performance_test_files/test34.cpp  performance_test_files/test51.cpp
+            performance_test_files/test52.cpp
             /boost/system//boost_system
           : release
           [ check-target-builds ../config//has_gmp : <define>TEST_MPF <define>TEST_MPZ <define>TEST_MPQ <source>gmp : ]

performance/performance_test.cpp

@@ -23,6 +23,7 @@
 // TEST_CPP_INT
 // TEST_CPP_INT_RATIONAL
 // TEST_CPP_BIN_FLOAT
+// TEST_CPP_DOUBLE_FLOAT
 //
 
 #include "performance_test.hpp"
@@ -214,6 +215,7 @@ int main()
    test49();
    test50();
    test51();
+   test52();
 
    quickbook_results();
    return 0;

performance/performance_test.hpp

@@ -12,7 +12,7 @@
 #if !defined(TEST_MPF) && !defined(TEST_MPZ) && \
     !defined(TEST_CPP_DEC_FLOAT) && !defined(TEST_MPFR) && !defined(TEST_MPQ) && !defined(TEST_TOMMATH) && \
     !defined(TEST_TOMMATH_BOOST_RATIONAL) && !defined(TEST_MPZ_BOOST_RATIONAL) && !defined(TEST_CPP_INT) && \
-    !defined(TEST_CPP_INT_RATIONAL) && !defined(TEST_CPP_BIN_FLOAT)
+    !defined(TEST_CPP_INT_RATIONAL) && !defined(TEST_CPP_BIN_FLOAT) && !defined(TEST_CPP_DOUBLE_FLOAT)
 #define TEST_MPF
 #define TEST_MPZ
 #define TEST_MPQ
@@ -23,6 +23,7 @@
 #define TEST_CPP_INT
 #define TEST_CPP_INT_RATIONAL
 #define TEST_CPP_BIN_FLOAT
+#define TEST_CPP_DOUBLE_FLOAT
 
 #ifdef _MSC_VER
 #pragma message("CAUTION!!: No backend type specified so testing everything.... this will take some time!!")
@@ -767,4 +768,5 @@ void test48();
 void test49();
 void test50();
 void test51();
+void test52();
 

performance/performance_test_df.cpp

@@ -0,0 +1,194 @@
+///////////////////////////////////////////////////////////////
+//  Copyright John Maddock 2019.
+//  Copyright Christopher Kormanyos 2021 - 2023.
+//  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
+
+#if !defined(TEST_CPP_DOUBLE_FLOAT)
+#define TEST_CPP_DOUBLE_FLOAT
+#endif
+
+#if !defined(TEST_CPP_BIN_FLOAT)
+#define TEST_CPP_BIN_FLOAT
+#endif
+
+//
+// This is the main entry point for our operator performance test suite.
+// In order to build this program, you must compile and link this file against
+// all the libs/multiprecision/performance/performance_test_files/*.cpp files.
+// 
+
+#include <boost/version.hpp>
+
+#if defined(TEST_CPP_DOUBLE_FLOAT)
+#include <boost/multiprecision/cpp_double_fp.hpp>
+#endif
+
+#if defined(TEST_CPP_BIN_FLOAT)
+#include <boost/multiprecision/cpp_bin_float.hpp>
+#endif
+
+#include "performance_test.hpp"
+
+//
+// Keys in order are:
+// Category
+// Operator
+// Type
+// Precision
+// Time
+//
+std::map<std::string, std::map<std::string, std::map<std::string, std::map<int, double> > > > result_table;
+
+unsigned bits_wanted; // for integer types
+
+void quickbook_results()
+{
+   //
+   // Keys in order are:
+   // Category
+   // Operator
+   // Type
+   // Precision
+   // Time
+   //
+   typedef std::map<std::string, std::map<std::string, std::map<std::string, std::map<int, double> > > >::const_iterator category_iterator;
+   typedef std::map<std::string, std::map<std::string, std::map<int, double> > >::const_iterator                         operator_iterator;
+   typedef std::map<std::string, std::map<int, double> >::const_iterator                                                 type_iterator;
+   typedef std::map<int, double>::const_iterator                                                                         precision_iterator;
+
+   for (category_iterator i = result_table.begin(); i != result_table.end(); ++i)
+   {
+      std::string cat = i->first;
+      cat[0]          = (char)std::toupper((char)cat[0]);
+      std::cout << "[section:" << i->first << "_performance " << cat << " Type Perfomance]" << std::endl;
+
+      for (operator_iterator j = i->second.begin(); j != i->second.end(); ++j)
+      {
+         std::string op = j->first;
+         std::cout << "[table Operator " << op << std::endl;
+         std::cout << "[[Backend]";
+
+         for (precision_iterator k = j->second.begin()->second.begin(); k != j->second.begin()->second.end(); ++k)
+         {
+            std::cout << "[" << k->first << " Bits]";
+         }
+         std::cout << "]\n";
+
+         std::vector<double> best_times(j->second.begin()->second.size(), (std::numeric_limits<double>::max)());
+         for (unsigned m = 0; m < j->second.begin()->second.size(); ++m)
+         {
+            for (type_iterator k = j->second.begin(); k != j->second.end(); ++k)
+            {
+               if (m < k->second.size())
+               {
+                  precision_iterator l = k->second.begin();
+                  std::advance(l, m);
+                  if (best_times[m] > l->second)
+                     best_times[m] = l->second ? l->second : best_times[m];
+               }
+            }
+         }
+
+         for (type_iterator k = j->second.begin(); k != j->second.end(); ++k)
+         {
+            std::cout << "[[" << k->first << "]";
+
+            unsigned m = 0;
+            for (precision_iterator l = k->second.begin(); l != k->second.end(); ++l)
+            {
+               double rel_time = l->second / best_times[m];
+               if (rel_time == 1)
+                  std::cout << "[[*" << rel_time << "]";
+               else
+                  std::cout << "[" << rel_time;
+               std::cout << " (" << l->second << "s)]";
+               ++m;
+            }
+
+            std::cout << "]\n";
+         }
+
+         std::cout << "]\n";
+      }
+
+      std::cout << "[endsect]" << std::endl;
+   }
+}
+
+#if defined(__HAS_INCLUDE)
+#if __has_include(<sys/utsname.h>)
+#define HAS_UTSNAME
+#include <sys/utsname.h>
+#endif
+#endif
+#ifdef _WIN32
+#include <windows.h>
+#endif
+
+void quickbook_platform_details()
+{
+   std::cout << "[table:platform Platform Details\n[[Platform][";
+#ifdef HAS_UTSNAME
+   utsname name;
+   uname(&name);
+   std::cout << name.sysname << " " << name.release << ", version " << name.version << ", " << name.machine << "]]\n";
+#elif defined(_WIN32)
+   std::cout << "Windows ";
+#ifdef _M_AMD64
+   std::cout << "x64";
+#elif defined(_M_IX86)
+   std::cout << "x86";
+#endif
+   std::cout << "]]\n";
+#endif
+   std::cout << "[[Compiler][" << BOOST_COMPILER << "]]\n";
+#ifdef TEST_MPZ
+   std::cout << "[[GMP][" << gmp_version << "]]\n";
+#endif
+#ifdef TEST_MPFR
+   std::cout << "[[MPFR][" << MPFR_VERSION << "]]\n";
+#endif
+   std::cout << "[[Boost][" << BOOST_VERSION << "]]\n";
+   std::cout << "[[Run date][" << __DATE__ << "]]\n";
+   std::cout << "]\n\n";
+}
+
+int main()
+{
+   quickbook_platform_details();
+
+   test32();
+   test52();
+
+   quickbook_results();
+}
+
+#ifdef TEST_CPP_DOUBLE_FLOAT
+using double_float_of_double_type = boost::multiprecision::cpp_double_double;
+
+constexpr auto digits10_for_performance_test = std::numeric_limits<double_float_of_double_type>::digits10;
+#else
+constexpr auto digits10_for_performance_test = 31;
+#endif
+
+void test32()
+{
+#ifdef TEST_CPP_BIN_FLOAT
+   static_assert(digits10_for_performance_test >= 31, "Error: Too few digits for performance comparison");
+
+   using cpp_bin_float_type =
+      boost::multiprecision::number<boost::multiprecision::backends::cpp_bin_float<digits10_for_performance_test>,
+                                    boost::multiprecision::et_off>;
+
+   test<cpp_bin_float_type>("cpp_bin_float", digits10_for_performance_test);
+#endif
+}
+
+void test52()
+{
+#ifdef TEST_CPP_DOUBLE_FLOAT
+   test<double_float_of_double_type>("cpp_double_fp_backend<double>", digits10_for_performance_test);
+#endif
+}

performance/performance_test_files/test52.cpp

@@ -0,0 +1,20 @@
+///////////////////////////////////////////////////////////////
+//  Copyright John Maddock 2019.
+//  Copyright Christopher Kormanyos 2021 - 2023.
+//  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
+
+#include "../performance_test.hpp"
+#if defined(TEST_CPP_DOUBLE_FLOAT)
+#include <boost/multiprecision/cpp_double_fp.hpp>
+#endif
+
+void test52()
+{
+#ifdef TEST_CPP_DOUBLE_FLOAT
+   using boost::multiprecision::cpp_double_double;
+
+   test<cpp_double_double>("cpp_double_fp_backend<double>", std::numeric_limits<cpp_double_double>::digits10);
+#endif
+}

test/Jamfile.v2

@@ -116,6 +116,7 @@ test-suite arithmetic_tests :
    [ run test_arithmetic_backend_concept.cpp no_eh_support ]
    [ compile test_arithmetic_skeleton.cpp ]
 
+   [ run test_arithmetic_df.cpp              no_eh_support  : : : <toolset>msvc:<cxxflags>-bigobj [ check-target-builds ../config//has_float128 : <source>quadmath ] ]
    [ run test_arithmetic_cpp_dec_float_1.cpp no_eh_support  : : : <toolset>msvc:<cxxflags>-bigobj [ check-target-builds ../config//has_float128 : <source>quadmath ] ]
    [ run test_arithmetic_cpp_dec_float_2.cpp no_eh_support  : : : <toolset>msvc:<cxxflags>-bigobj [ check-target-builds ../config//has_float128 : <source>quadmath ] ]
    [ run test_arithmetic_cpp_dec_float_3.cpp no_eh_support  : : : <toolset>msvc:<cxxflags>-bigobj [ check-target-builds ../config//has_float128 : <source>quadmath ] ]
@@ -260,7 +261,7 @@ alias precision_tests : test_preserve_source_precision_gmp test_preserve_source_
 rule get_function_tests
 {
    local result ;
-   for local source in test_exp.cpp test_log.cpp test_pow.cpp test_sinh.cpp test_sqrt.cpp test_cosh.cpp test_tanh.cpp test_sin.cpp test_cos.cpp test_tan.cpp test_asin.cpp test_acos.cpp test_atan.cpp test_round.cpp test_fpclassify.cpp test_sf_import_c99.cpp
+   for local source in test_exp.cpp test_log.cpp test_pow.cpp test_pow_df.cpp test_sinh.cpp test_sqrt.cpp test_cosh.cpp test_tanh.cpp test_sin.cpp test_cos.cpp test_tan.cpp test_asin.cpp test_acos.cpp test_atan.cpp test_round.cpp test_fpclassify.cpp test_sf_import_c99.cpp
    {
          result += [ run $(source) gmp no_eh_support
               : # command line
@@ -295,6 +296,13 @@ rule get_function_tests
               : # requirements
                <define>TEST_CPP_BIN_FLOAT
               : $(source:B)_cpp_bin_float ] ;
+         result += [ run $(source) quadmath no_eh_support
+              : # command line
+              : # input files
+              : # requirements
+               [ check-target-builds ../config//has_float128 : : <build>no ]
+               <define>TEST_CPP_DOUBLE_FLOAT
+              : $(source:B)_cpp_double_fp ] ;
          result += [ run $(source) quadmath no_eh_support
               : # command line
               : # input files
@@ -384,6 +392,13 @@ test-suite functions_and_limits :
               <define>TEST_CPP_BIN_FLOAT <define>PRINT_MAX_DIGITS10
               : test_numeric_limits_cpp_bin_float ]
 
+      [ run test_numeric_limits.cpp no_eh_support
+              : # command line
+              : # input files
+              : # requirements
+              <define>TEST_CPP_DOUBLE_FLOAT <define>PRINT_MAX_DIGITS10
+              : test_numeric_limits_cpp_double_float ]
+
       [ run test_numeric_limits.cpp $(TOMMATH) no_eh_support
               : # command line
               : # input files

test/concepts/sf_concept_check_basic.cpp

@@ -1,4 +1,5 @@
 //  Copyright John Maddock 2012.
+//  Copyright Christopher Kormanyos 2021 - 2023.
 //  Use, modification and distribution are subject to 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)
@@ -19,7 +20,7 @@
 #include <boost/container_hash/hash.hpp>
 #include <libs/math/test/compile_test/poison.hpp>
 
-#if !defined(TEST_MPF_50) && !defined(TEST_BACKEND) && !defined(TEST_MPZ) && !defined(TEST_CPP_DEC_FLOAT) && !defined(TEST_MPFR_50) && !defined(TEST_MPFR_6) && !defined(TEST_MPFR_15) && !defined(TEST_MPFR_17) && !defined(TEST_MPFR_30) && !defined(TEST_CPP_DEC_FLOAT_NO_ET) && !defined(TEST_LOGGED_ADAPTER) && !defined(TEST_CPP_BIN_FLOAT)
+#if !defined(TEST_MPF_50) && !defined(TEST_BACKEND) && !defined(TEST_MPZ) && !defined(TEST_CPP_DEC_FLOAT) && !defined(TEST_MPFR_50) && !defined(TEST_MPFR_6) && !defined(TEST_MPFR_15) && !defined(TEST_MPFR_17) && !defined(TEST_MPFR_30) && !defined(TEST_CPP_DEC_FLOAT_NO_ET) && !defined(TEST_LOGGED_ADAPTER) && !defined(TEST_CPP_BIN_FLOAT) && !defined(TEST_CPP_DOUBLE_FLOAT)
 #define TEST_MPF_50
 #define TEST_BACKEND
 #define TEST_MPZ
@@ -32,6 +33,7 @@
 #define TEST_CPP_DEC_FLOAT_NO_ET
 #define TEST_LOGGED_ADAPTER
 #define TEST_CPP_BIN_FLOAT
+#define TEST_CPP_DOUBLE_FLOAT
 
 #ifdef _MSC_VER
 #pragma message("CAUTION!!: No backend type specified so testing everything.... this will take some time!!")
@@ -60,6 +62,9 @@
 #ifdef TEST_LOGGED_ADAPTER
 #include <boost/multiprecision/logged_adaptor.hpp>
 #endif
+#ifdef TEST_CPP_DOUBLE_FLOAT
+#include <boost/multiprecision/cpp_double_fp.hpp>
+#endif
 
 #include <boost/math/special_functions.hpp>
 
@@ -157,6 +162,11 @@ void foo()
    typedef boost::multiprecision::number<boost::multiprecision::logged_adaptor<boost::multiprecision::cpp_dec_float<50> > > num_t;
    test_extra(num_t());
 #endif
+#ifdef TEST_CPP_DOUBLE_FLOAT
+   using cpp_double_float_of_double_type = boost::multiprecision::cpp_double_double;
+
+   test_extra(cpp_double_float_of_double_type());
+#endif
 }
 
 int main()

test/concepts/sf_concept_check_bessel.cpp

@@ -1,4 +1,5 @@
 //  Copyright John Maddock 2012.
+//  Copyright Christopher Kormanyos 2021 - 2023.
 //  Use, modification and distribution are subject to 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)
@@ -19,7 +20,7 @@
 #include <boost/container_hash/hash.hpp>
 #include <libs/math/test/compile_test/poison.hpp>
 
-#if !defined(TEST_MPF_50) && !defined(TEST_BACKEND) && !defined(TEST_MPZ) && !defined(TEST_CPP_DEC_FLOAT) && !defined(TEST_MPFR_50) && !defined(TEST_MPFR_6) && !defined(TEST_MPFR_15) && !defined(TEST_MPFR_17) && !defined(TEST_MPFR_30) && !defined(TEST_CPP_DEC_FLOAT_NO_ET) && !defined(TEST_LOGGED_ADAPTER) && !defined(TEST_CPP_BIN_FLOAT)
+#if !defined(TEST_MPF_50) && !defined(TEST_BACKEND) && !defined(TEST_MPZ) && !defined(TEST_CPP_DEC_FLOAT) && !defined(TEST_MPFR_50) && !defined(TEST_MPFR_6) && !defined(TEST_MPFR_15) && !defined(TEST_MPFR_17) && !defined(TEST_MPFR_30) && !defined(TEST_CPP_DEC_FLOAT_NO_ET) && !defined(TEST_LOGGED_ADAPTER) && !defined(TEST_CPP_BIN_FLOAT) && !defined(TEST_CPP_DOUBLE_FLOAT)
 #define TEST_MPF_50
 #define TEST_BACKEND
 #define TEST_MPZ
@@ -32,6 +33,7 @@
 #define TEST_CPP_DEC_FLOAT_NO_ET
 #define TEST_LOGGED_ADAPTER
 #define TEST_CPP_BIN_FLOAT
+#define TEST_CPP_DOUBLE_FLOAT
 
 #ifdef _MSC_VER
 #pragma message("CAUTION!!: No backend type specified so testing everything.... this will take some time!!")
@@ -60,6 +62,9 @@
 #ifdef TEST_LOGGED_ADAPTER
 #include <boost/multiprecision/logged_adaptor.hpp>
 #endif
+#ifdef TEST_CPP_DOUBLE_FLOAT
+#include <boost/multiprecision/cpp_double_fp.hpp>
+#endif
 
 #include <boost/math/special_functions.hpp>
 
@@ -124,6 +129,11 @@ void foo()
    typedef boost::multiprecision::number<boost::multiprecision::logged_adaptor<boost::multiprecision::cpp_dec_float<50> > > num_t;
    test_extra(num_t());
 #endif
+#ifdef TEST_CPP_DOUBLE_FLOAT
+   using cpp_double_float_of_double_type = boost::multiprecision::cpp_double_double;
+
+   test_extra(cpp_double_float_of_double_type());
+#endif
 }
 
 int main()