libs/math/test/test_lambert_w.cpp

   1 // Copyright Paul A. Bristow 2016, 2017, 2018.
   2 // Copyright John Maddock 2016.
   3
   4 // Use, modification and distribution are subject to the
   5 // Boost Software License, Version 1.0.
   6 // (See accompanying file LICENSE_1_0.txt
   7 // or copy at http://www.boost.org/LICENSE_1_0.txt)
   8
   9 // test_lambert_w.cpp
  10 //! \brief Basic sanity tests for Lambert W function using algorithms
  11 // informed by Thomas Luu, Darko Veberic and Tosio Fukushima for W0
  12 // and rational polynomials by John Maddock.
  13
  14 // #define BOOST_MATH_TEST_MULTIPRECISION  // Add tests for several multiprecision types (not just built-in).
  15 // #define BOOST_MATH_TEST_FLOAT128 // Add test using float128 type (GCC only, needing gnu++17 and quadmath library).
  16
  17 #ifdef BOOST_MATH_TEST_FLOAT128
  18 #include <boost/cstdfloat.hpp> // For float_64_t, float128_t. Must be first include!
  19 #endif // #ifdef #ifdef BOOST_MATH_TEST_FLOAT128
  20 // Needs gnu++17 for BOOST_HAS_FLOAT128
  21 #include <boost/config.hpp>   // for BOOST_MSVC definition etc.
  22 #include <boost/version.hpp>   // for BOOST_MSVC versions.
  23
  24 // Boost macros
  25 #define BOOST_TEST_MAIN
  26 #define BOOST_LIB_DIAGNOSTIC "on" // Report library file details.
  27 #include <boost/test/included/unit_test.hpp> // Boost.Test
  28 // #include <boost/test/unit_test.hpp> // Boost.Test
  29 #include <boost/test/tools/floating_point_comparison.hpp>
  30
  31 #include <boost/array.hpp>
  32 #include <boost/lexical_cast.hpp>
  33 #include <boost/type_traits/is_constructible.hpp>
  34
  35 #ifdef BOOST_MATH_TEST_MULTIPRECISION
  36 #include <boost/multiprecision/cpp_dec_float.hpp> // boost::multiprecision::cpp_dec_float_50
  37 using boost::multiprecision::cpp_dec_float_50;
  38
  39 #include <boost/multiprecision/cpp_bin_float.hpp>
  40 using boost::multiprecision::cpp_bin_float_quad;
  41
  42 #include <boost/math/concepts/real_concept.hpp>
  43
  44 #ifdef BOOST_MATH_TEST_FLOAT128
  45
  46 #ifdef BOOST_HAS_FLOAT128
  47 // Including this header below without float128 triggers:
  48 // fatal error C1189: #error:  "Sorry compiler is neither GCC, not Intel, don't know how to configure this header."
  49 #include <boost/multiprecision/float128.hpp>
  50 using boost::multiprecision::float128;
  51 #endif // ifdef BOOST_HAS_FLOAT128
  52 #endif // #ifdef #ifdef BOOST_MATH_TEST_FLOAT128
  53
  54 #endif //   #ifdef BOOST_MATH_TEST_MULTIPRECISION
  55
  56 //#include <boost/fixed_point/fixed_point.hpp> // If available.
  57
  58 #include <boost/math/concepts/real_concept.hpp> // for real_concept tests.
  59 #include <boost/math/special_functions/fpclassify.hpp> // isnan, ifinite.
  60 #include <boost/math/special_functions/next.hpp> // float_next, float_prior
  61 using boost::math::float_next;
  62 using boost::math::float_prior;
  63 #include <boost/math/special_functions/ulp.hpp>  // ulp
  64
  65 #include <boost/math/tools/test_value.hpp>  // for create_test_value and macro BOOST_MATH_TEST_VALUE.
  66 #include <boost/math/policies/policy.hpp>
  67 using boost::math::policies::digits2;
  68 using boost::math::policies::digits10;
  69 #include <boost/math/special_functions/lambert_w.hpp> // For Lambert W lambert_w function.
  70 using boost::math::lambert_wm1;
  71 using boost::math::lambert_w0;
  72
  73 #include "table_type.hpp"
  74
  75 #ifndef SC_
  76 #  define SC_(x) boost::lexical_cast<typename table_type<T>::type>(BOOST_STRINGIZE(x))
  77 #endif
  78
  79
  80 #include <limits>
  81 #include <cmath>
  82 #include <typeinfo>
  83 #include <iostream>
  84 #include <exception>
  85
  86 std::string show_versions(void);
  87
  88 //! Build a message of information about build, architecture, address model, platform, ...
  89 std::string show_versions(void)
  90 {
  91   // Some of this information can also be obtained from running with a Custom Post-build step
  92   // adding the option --build_info=yes
  93     // "$(TargetDir)$(TargetName).exe" --build_info=yes
  94
  95   std::ostringstream message;
  96
  97   message << "Program: " << __FILE__ << "\n";
  98 #ifdef __TIMESTAMP__
  99   message << __TIMESTAMP__;
 100 #endif
 101   message << "\nBuildInfo:\n" "  Platform " << BOOST_PLATFORM;
 102   // http://stackoverflow.com/questions/1505582/determining-32-vs-64-bit-in-c
 103 #if defined(__LP64__) || defined(_WIN64) || (defined(__x86_64__) && !defined(__ILP32__) ) || defined(_M_X64) || defined(__ia64) || defined (_M_IA64) || defined(__aarch64__) || defined(__powerpc64__)
 104   message << ", 64-bit.";
 105 #else
 106   message << ", 32-bit.";
 107 #endif
 108
 109   message << "\n  Compiler " BOOST_COMPILER;
 110 #ifdef BOOST_MSC_VER
 111 #ifdef _MSC_FULL_VER
 112   message << "\n  MSVC version " << BOOST_STRINGIZE(_MSC_FULL_VER) << ".";
 113 #endif
 114 #ifdef __WIN64
 115   mess age << "\n WIN64" << std::endl;
 116 #endif // __WIN64
 117 #ifdef _WIN32
 118   message << "\n WIN32" << std::endl;
 119 #endif  // __WIN32
 120 #endif
 121 #ifdef __GNUC__
 122   //PRINT_MACRO(__GNUC__);
 123   //PRINT_MACRO(__GNUC_MINOR__);
 124   //PRINT_MACRO(__GNUC_PATCH__);
 125   std::cout << "GCC " << __VERSION__ << std::endl;
 126   //PRINT_MACRO(LONG_MAX);
 127 #endif // __GNUC__
 128
 129 #ifdef __MINGW64__
 130 std::cout << "MINGW64 " << __MINGW32_MAJOR_VERSION << __MINGW32_MINOR_VERSION << std::endl;
 131 //
 132 //  << __MINGW64_MAJOR_VERSION << __MINGW64_MINOR_VERSION << std::endl; not declared in this scope???
 133 #endif // __MINGW64__
 134
 135 #ifdef __MINGW32__
 136 std::cout << "MINGW64 " << __MINGW32_MAJOR_VERSION << __MINGW32_MINOR_VERSION << std::endl;
 137 #endif // __MINGW32__
 138
 139   message << "\n  STL " << BOOST_STDLIB;
 140   message << "\n  Boost version " << BOOST_VERSION / 100000 << "." << BOOST_VERSION / 100 % 1000 << "." << BOOST_VERSION % 100;
 141
 142 #ifdef BOOST_MATH_TEST_MULTIPRECISION
 143   message << "\nBOOST_MATH_TEST_MULTIPRECISION defined for multiprecision tests. " << std::endl;
 144 #else
 145  message << "\nBOOST_MATH_TEST_MULTIPRECISION not defined so NO multiprecision tests. " << std::endl;
 146 #endif // BOOST_MATH_TEST_MULTIPRECISION
 147
 148 #ifdef BOOST_HAS_FLOAT128
 149   message << "BOOST_HAS_FLOAT128 is defined." << std::endl;
 150 #endif // ifdef BOOST_HAS_FLOAT128
 151
 152   message << std::endl;
 153   return message.str();
 154 } // std::string show_versions()
 155
 156
 157 template <class T>
 158 void wolfram_test_moderate_values()
 159 {
 160    //
 161    // Spots of moderate value http://www.wolframalpha.com/input/?i=TABLE%5B%5BN%5B-1%2Fe%2Bi,+50%5D,+N%5BLambertW%5B-1%2Fe%2Bi%5D,+50%5D%5D,+%7Bi,+1%2F8,+6,+1%2F8%7D%5D
 162    //
 163    static const boost::array<boost::array<typename table_type<T>::type, 2>, 96/2> wolfram_test_small_neg =
 164    {{
 165       {{ SC_(-0.24287944117144232159552377016146086744581113103177), SC_(-0.34187241316000572901412382650748493957063539755395) }},{{ SC_(-0.11787944117144232159552377016146086744581113103177), SC_(-0.13490446826612135454875992607636577833255418182633) }},{{ SC_(0.0071205588285576784044762298385391325541888689682322), SC_(0.0070703912528860797819274709355398032954165697080076) }},{{ SC_(0.13212055882855767840447622983853913255418886896823), SC_(0.11747650174894814471295063763686399700941650918302) }},{{ SC_(0.25712055882855767840447622983853913255418886896823), SC_(0.20869089404810562424547046857454995304964242368484) }},{{ SC_(0.38212055882855767840447622983853913255418886896823), SC_(0.28683366713002653952708635029764106993377156175310) }},{{ SC_(0.50712055882855767840447622983853913255418886896823), SC_(0.35542749308004931507852679571061486656821523044053) }},{{ SC_(0.63212055882855767840447622983853913255418886896823), SC_(0.41670399881776590750659327292575356285757792776250) }},{{ SC_(0.75712055882855767840447622983853913255418886896823), SC_(0.47217430075943420437939326812963066971059146681283) }},{{ SC_(0.88212055882855767840447622983853913255418886896823), SC_(0.52291321715862065064992942239384690347359852107504) }},{{ SC_(1.0071205588285576784044762298385391325541888689682), SC_(0.56971477154593975582335630229323210831843899740884) }},{{ SC_(1.1321205588285576784044762298385391325541888689682), SC_(0.61318350578224462394572352964726524514921241969798) }},{{ SC_(1.2571205588285576784044762298385391325541888689682), SC_(0.65379115237566259933564436658873734121781110980034) }},{{ SC_(1.3821205588285576784044762298385391325541888689682), SC_(0.69191341320406026236753559968630177636780741203666) }},{{ SC_(1.5071205588285576784044762298385391325541888689682), SC_(0.72785472286747598788295903283683432537852776142064) }},{{ SC_(1.6321205588285576784044762298385391325541888689682), SC_(0.76186544538805130363636977458614856100481979440639) }},{{ SC_(1.7571205588285576784044762298385391325541888689682), SC_(0.79415413501531119849043049331889268136479923750037) }},{{ SC_(1.8821205588285576784044762298385391325541888689682), SC_(0.82489647878345700122288701550494847447982817483512) }},{{ SC_(2.0071205588285576784044762298385391325541888689682), SC_(0.85424194939386899439722948096520865643710851410970) }},{{ SC_(2.1321205588285576784044762298385391325541888689682), SC_(0.88231884173371311472940735780441644004275449741412) }},{{ SC_(2.2571205588285576784044762298385391325541888689682), SC_(0.90923814516532488963517314558961057510689871415824) }},{{ SC_(2.3821205588285576784044762298385391325541888689682), SC_(0.93509656212104191797135657485515114635876341802516) }},{{ SC_(2.5071205588285576784044762298385391325541888689682), SC_(0.95997889061117906067636869169049106690165665554172) }},{{ SC_(2.6321205588285576784044762298385391325541888689682), SC_(0.98395992590529701946948066548039809917492328184099) }},{{ SC_(2.7571205588285576784044762298385391325541888689682), SC_(1.0071059939771381126732041109492705496242899774655) }},{{ SC_(2.8821205588285576784044762298385391325541888689682), SC_(1.0294761995723706229651673877352399077168142413723) }},{{ SC_(3.0071205588285576784044762298385391325541888689682), SC_(1.0511234507020167125769191146012321442040919222298) }},{{ SC_(3.1321205588285576784044762298385391325541888689682), SC_(1.0720953062286332723365148290552887215464891915069) }},{{ SC_(3.2571205588285576784044762298385391325541888689682), SC_(1.0924346821831089228990349517861599064007594751702) }},{{ SC_(3.3821205588285576784044762298385391325541888689682), SC_(1.1121804443118533629930276674418322662764569673766) }},{{ SC_(3.5071205588285576784044762298385391325541888689682), SC_(1.1313679082795201044696522785560810652358663683706) }},{{ SC_(3.6321205588285576784044762298385391325541888689682), SC_(1.1500292643692387775614691790201052907317404963905) }},{{ SC_(3.7571205588285576784044762298385391325541888689682), SC_(1.1681939400299161555212785901786587344721733034978) }},{{ SC_(3.8821205588285576784044762298385391325541888689682), SC_(1.1858889109341735194685896928615740804115521714257) }},{{ SC_(4.0071205588285576784044762298385391325541888689682), SC_(1.2031389691267953962289622785796365085402661808452) }},{{ SC_(4.1321205588285576784044762298385391325541888689682), SC_(1.2199669552139996161903252772502362264684476580522) }},{{ SC_(4.2571205588285576784044762298385391325541888689682), SC_(1.2363939602597347325278067608637615539794532870296) }},{{ SC_(4.3821205588285576784044762298385391325541888689682), SC_(1.2524395020361026107226019920575290018966524482736) }},{{ SC_(4.5071205588285576784044762298385391325541888689682), SC_(1.2681216794607666389159742215265331040507889789444) }},{{ SC_(4.6321205588285576784044762298385391325541888689682), SC_(1.2834573083995295018572263393035905604511320189369) }},{{ SC_(4.7571205588285576784044762298385391325541888689682), SC_(1.2984620414827281167361144981111712803667945033184) }},{{ SC_(4.8821205588285576784044762298385391325541888689682), SC_(1.3131504741533499076663954559108617687274731330916) }},{{ SC_(5.0071205588285576784044762298385391325541888689682), SC_(1.3275362388125116267199919229657120782894307415376) }},{{ SC_(5.1321205588285576784044762298385391325541888689682), SC_(1.3416320886383928057123774168081846145768561516693) }},{{ SC_(5.2571205588285576784044762298385391325541888689682), SC_(1.3554499724155634924134183248962114419200302481356) }},{{ SC_(5.3821205588285576784044762298385391325541888689682), SC_(1.3690011015132087699425938733927188719869603184010) }},{{ SC_(5.5071205588285576784044762298385391325541888689682), SC_(1.3822960099853765706075495327819109601506356054327) }},{{ SC_(5.6321205588285576784044762298385391325541888689682), SC_(1.3953446086279755263512146907828727538440007615239) }}
 166    }};
 167    T tolerance = boost::math::tools::epsilon<T>() * 3;
 168    if (std::numeric_limits<T>::digits10 > 40)
 169       tolerance *= 4;  // arbitrary precision types have lower accuracy on exp(z).
 170    for (unsigned i = 0; i < wolfram_test_small_neg.size(); ++i)
 171    {
 172       BOOST_CHECK_CLOSE_FRACTION(boost::math::lambert_w0(T(wolfram_test_small_neg[i][0])), T(wolfram_test_small_neg[i][1]), tolerance);
 173    }
 174 }
 175
 176 template <class T>
 177 void wolfram_test_small_pos()
 178 {
 179    //
 180    // Spots near zero and positive http://www.wolframalpha.com/input/?i=TABLE%5B%5BN%5BPi+*+10%5Ei,+50%5D,+N%5BLambertW%5BPi+*+10%5Ei%5D,+50%5D%5D,+%7Bi,+-25,+-1%7D%5D
 181    //
 182    static const boost::array<boost::array<typename table_type<T>::type, 2>, 25> wolfram_test_small_neg =
 183    {{
 184       {{ SC_(3.1415926535897932384626433832795028841971693993751e-25), SC_(3.1415926535897932384626423963190627752613075159265e-25) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-24), SC_(3.1415926535897932384626335136751017948385505649306e-24) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-23), SC_(3.1415926535897932384625446872354919906109810591160e-23) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-22), SC_(3.1415926535897932384616564228393939483352864153693e-22) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-21), SC_(3.1415926535897932384527737788784135255783814177903e-21) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-20), SC_(3.1415926535897932383639473392686092980134754308784e-20) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-19), SC_(3.1415926535897932374756829431705670227788144495920e-19) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-18), SC_(3.1415926535897932285930389821901443118720934199487e-18) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-17), SC_(3.1415926535897931397665993723859213467937614455864e-17) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-16), SC_(3.1415926535897922515022032743441060948982739088029e-16) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-15), SC_(3.1415926535897833688582422939673934647266189937296e-15) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-14), SC_(3.1415926535896945424186324943442560413318839066091e-14) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-13), SC_(3.1415926535888062780225349125117696393347268403158e-13) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-12), SC_(3.1415926535799236340616005340756885831699803736331e-12) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-11), SC_(3.1415926534910971944564007385929431896486546006413e-11) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-10), SC_(3.1415926526028327988188016713407935109104110982749e-10) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-9), SC_(3.1415926437201888838826995251371676507148394412103e-9) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-8), SC_(3.1415925548937538785102994823474670579278874210259e-8) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-7), SC_(3.1415916666298182234172285804275105377159084331529e-7) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e-6), SC_(3.1415827840319013043684920305205420694740106954961e-6) }},{{ SC_(0.000031415926535897932384626433832795028841971693993751), SC_(0.000031414939621964641052828244109272729597989570861172) }},{{ SC_(0.00031415926535897932384626433832795028841971693993751), SC_(0.00031406061579842362125003023838529350597159230209458) }},{{ SC_(0.0031415926535897932384626433832795028841971693993751), SC_(0.0031317693004296877733926356188004473035977501714541) }},{{ SC_(0.031415926535897932384626433832795028841971693993751), SC_(0.030473027596269883517196555192955092247613270959259) }},{{ SC_(0.31415926535897932384626433832795028841971693993751), SC_(0.24571751376320572448656753973370462139374436325987) }}
 185    }};
 186    T tolerance = boost::math::tools::epsilon<T>() * 3;
 187    if (std::numeric_limits<T>::digits10 > 40)
 188       tolerance *= 3;  // arbitrary precision types have lower accuracy on exp(z).
 189    for (unsigned i = 0; i < wolfram_test_small_neg.size(); ++i)
 190    {
 191       BOOST_CHECK_CLOSE_FRACTION(boost::math::lambert_w0(T(wolfram_test_small_neg[i][0])), T(wolfram_test_small_neg[i][1]), tolerance);
 192    }
 193 }
 194
 195 template <class T>
 196 void wolfram_test_small_neg()
 197 {
 198    //
 199    // Spots near zero and negative http://www.wolframalpha.com/input/?i=TABLE%5B%5BN%5B-Pi+*+10%5Ei,+50%5D,+N%5BLambertW%5B-Pi+*+10%5Ei%5D,+50%5D%5D,+%7Bi,+-25,+-1%7D%5D
 200    //
 201    static const boost::array<boost::array<typename table_type<T>::type, 2>, 70/2> wolfram_test_small_neg =
 202    {{
 203       {{ SC_(-3.1415926535897932384626433832795028841971693993751e-25), SC_(-3.1415926535897932384626443702399429931330312828247e-25) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-24), SC_(-3.1415926535897932384626532528839039735557882339126e-24) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-23), SC_(-3.1415926535897932384627420793235137777833577489360e-23) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-22), SC_(-3.1415926535897932384636303437196118200590533135692e-22) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-21), SC_(-3.1415926535897932384725129876805922428160503997900e-21) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-20), SC_(-3.1415926535897932385613394272903964703901652508759e-20) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-19), SC_(-3.1415926535897932394496038233884387465457126495672e-19) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-18), SC_(-3.1415926535897932483322477843688615495410754197010e-18) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-17), SC_(-3.1415926535897933371586873941730937234835814431099e-17) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-16), SC_(-3.1415926535897942254230834922158298617964738845526e-16) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-15), SC_(-3.1415926535898031080670444726846311337086192655470e-15) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-14), SC_(-3.1415926535898919345066542815166327311524009447840e-14) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-13), SC_(-3.1415926535907801989027527842355365380542172227242e-13) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-12), SC_(-3.1415926535996628428637792513133580846848848572500e-12) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-11), SC_(-3.1415926536884892824781879109701525247983589696795e-11) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-10), SC_(-3.1415926545767536790366733956272068630669876574730e-10) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-9), SC_(-3.1415926634593976860614172823213018318134944055260e-9) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-8), SC_(-3.1415927522858419002979913741894684038594384671969e-8) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-7), SC_(-3.1415936405506984418084674995072645049396296346958e-7) }},{{ SC_(-3.1415926535897932384626433832795028841971693993751e-6), SC_(-3.1416025232407040026008819016148803316716797067967e-6) }},{{ SC_(-0.000031415926535897932384626433832795028841971693993751), SC_(-0.000031416913542850054076094590477471913042739704497976) }},{{ SC_(-0.00031415926535897932384626433832795028841971693993751), SC_(-0.00031425800793839694440655801311183879569843264709852) }},{{ SC_(-0.0031415926535897932384626433832795028841971693993751), SC_(-0.0031515090287677856656576839914749012339811781712486) }},{{ SC_(-0.031415926535897932384626433832795028841971693993751), SC_(-0.032452164493239992272463616095775075564894751832128) }},{{ SC_(-0.31415926535897932384626433832795028841971693993751), SC_(-0.53804834513759287053587977755877044660611017981968) }},
 204       {{ SC_(-0.090099009900990099009900990099009900990099009900990), SC_(-0.099527797075226962190621767732039397602197803169897)}},{{ SC_(-0.080198019801980198019801980198019801980198019801980), SC_(-0.087534530933383521242151071722737877728489741787814) }},{{ SC_(-0.070297029702970297029702970297029702970297029702970), SC_(-0.075835379000403488962496062196568904002201151736290) }},{{ SC_(-0.060396039603960396039603960396039603960396039603960), SC_(-0.064414449758822413858363348099340678962612835311800) }},{{ SC_(-0.050495049504950495049504950495049504950495049504950), SC_(-0.053257171600878093079366736202964706966166164696873) }},{{ SC_(-0.040594059405940594059405940594059405940594059405941), SC_(-0.042350146588050412657332988380168720859403591863698) }},{{ SC_(-0.030693069306930693069306930693069306930693069306931), SC_(-0.031681024260949098136757222042165581145138786336298) }},{{ SC_(-0.020792079207920792079207920792079207920792079207921), SC_(-0.021238392251213645736199359110665662967213312773617) }},{{ SC_(-0.010891089108910891089108910891089108910891089108911), SC_(-0.011011681049909946810068329378571761407667575030714) }},{{ SC_(-0.00099009900990099009900990099009900990099009900990099), SC_(-0.00099108076440319890968631186785975507712384928918616) }}
 205    }};
 206    T tolerance = boost::math::tools::epsilon<T>() * 3;
 207    if (std::numeric_limits<T>::digits10 > 40)
 208       tolerance *= 3;  // arbitrary precision types have lower accuracy on exp(z).
 209    for (unsigned i = 0; i < wolfram_test_small_neg.size(); ++i)
 210    {
 211       BOOST_CHECK_CLOSE_FRACTION(boost::math::lambert_w0(T(wolfram_test_small_neg[i][0])), T(wolfram_test_small_neg[i][1]), tolerance);
 212    }
 213 }
 214
 215 template <class T>
 216 void wolfram_test_large(const boost::mpl::true_&)
 217 {
 218    //
 219    // Spots near the singularity from http://www.wolframalpha.com/input/?i=TABLE%5B%5BN%5B-1%2Fe%2B2%5E-i,+50%5D,+N%5BLambertW%5B-1%2Fe+%2B+2%5E-i%5D,+50%5D%5D,+%7Bi,+2,+40%7D%5D
 220    //
 221    static const boost::array<boost::array<typename table_type<T>::type, 2>, 28/2> wolfram_test_large_data =
 222    { {
 223       {{ SC_(3.1415926535897932384626433832795028841971693993751e350), SC_(800.36444525326526998205084284403447902093784176640) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e400), SC_(915.35945025352715923124904626896745356022974283730) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e450), SC_(1030.3703481552571717312484086444052442055003737018) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e500), SC_(1145.3937726197879355969554296951287620979399652268) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e550), SC_(1260.4273249433458391941776841900870933799293511610) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e600), SC_(1375.4692354682341092954911299903937009237749971748) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e650), SC_(1490.5181612342761763990969379122584268166707632003) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e700), SC_(1605.5730589637597079362569020729894833435943718597) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e750), SC_(1720.6331020467166402802313799793443913873949058922) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e800), SC_(1835.6976244160526737141293452999638879204852786698) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e850), SC_(1950.7660814940759743605616247252782614446819652848) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e900), SC_(2065.8380223354646200773160641407055989098916114637) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e950), SC_(2180.9130693229593212006354812037286740424563145700) }},{{ SC_(3.1415926535897932384626433832795028841971693993751e1000), SC_(2295.9909030845346718801238821248991904602625884450) }}
 224    } };
 225    T tolerance = boost::math::tools::epsilon<T>() * 3;
 226    if (std::numeric_limits<T>::digits10 > 40)
 227       tolerance *= 3;  // arbitrary precision types have lower accuracy on exp(z).
 228    for (unsigned i = 0; i < wolfram_test_large_data.size(); ++i)
 229    {
 230       BOOST_CHECK_CLOSE_FRACTION(boost::math::lambert_w0(T(wolfram_test_large_data[i][0])), T(wolfram_test_large_data[i][1]), tolerance);
 231    }
 232 }
 233 template <class T>
 234 void wolfram_test_large(const boost::mpl::false_&){}
 235
 236 template <class T>
 237 void wolfram_test_large()
 238 {
 239    wolfram_test_large<T>(boost::mpl::bool_<(std::numeric_limits<T>::max_exponent10 > 1000)>());
 240 }
 241
 242
 243 template <class T>
 244 void wolfram_test_near_singularity()
 245 {
 246    //
 247    // Spots near the singularity from http://www.wolframalpha.com/input/?i=TABLE%5B%5BN%5B-1%2Fe%2B2%5E-i,+50%5D,+N%5BLambertW%5B-1%2Fe+%2B+2%5E-i%5D,+50%5D%5D,+%7Bi,+2,+40%7D%5D
 248    //
 249    static const boost::array<boost::array<typename table_type<T>::type, 2>, 39> wolfram_test_near_singularity_data =
 250    {{
 251       { { SC_(-0.11787944117144233402427744294982403516769409179688), SC_(-0.13490446826612137099065142885543349308605449591189) } },{ { SC_(-0.24287944117144233402427744294982403516769409179688), SC_(-0.34187241316000575559631565516533717918703951393828) } },{ { SC_(-0.30537944117144233402427744294982403516769409179688), SC_(-0.50704532478540670242736394530166187052909039079642) } },{ { SC_(-0.33662944117144233402427744294982403516769409179688), SC_(-0.63562321628494791544895212508757067989859372121549) } },{ { SC_(-0.35225444117144233402427744294982403516769409179688), SC_(-0.73357201771558852140844624841371893543359405991894) } },{ { SC_(-0.36006694117144233402427744294982403516769409179688), SC_(-0.80685912552602238275976720505076149562188136941981) } },{ { SC_(-0.36397319117144233402427744294982403516769409179688), SC_(-0.86091151614390373770305184939107560322835214525382) } },{ { SC_(-0.36592631617144233402427744294982403516769409179688), SC_(-0.90033567669608907987528169545609510444951296636737) } },{ { SC_(-0.36690287867144233402427744294982403516769409179688), SC_(-0.92884889586304130900291705545970353898661233095513) } },{ { SC_(-0.36739115992144233402427744294982403516769409179688), SC_(-0.94934196763921122756108351994184213101752011076782) } },{ { SC_(-0.36763530054644233402427744294982403516769409179688), SC_(-0.96400324129495105632485735566132352543383271582526) } },{ { SC_(-0.36775737085894233402427744294982403516769409179688), SC_(-0.97445736712728703357755243595334553847237474201138) } },{ { SC_(-0.36781840601519233402427744294982403516769409179688), SC_(-0.98189372378619472154195350108189165241865132390473) } },{ { SC_(-0.36784892359331733402427744294982403516769409179688), SC_(-0.98717434434269671591894280580432721487757138768109) } },{ { SC_(-0.36786418238237983402427744294982403516769409179688), SC_(-0.99091955260257317141206161906086819616043312707614) } },{ { SC_(-0.36787181177691108402427744294982403516769409179688), SC_(-0.99357346775773151586057357459040504547191256911173) } },{ { SC_(-0.36787562647417670902427744294982403516769409179688), SC_(-0.99545290640175819861266174073519228782773422561472) } },{ { SC_(-0.36787753382280952152427744294982403516769409179688), SC_(-0.99678329264937600678258333756796350065436689760936) } },{ { SC_(-0.36787848749712592777427744294982403516769409179688), SC_(-0.99772473035978895659981485126201758865515569761514) } },{ { SC_(-0.36787896433428413089927744294982403516769409179688), SC_(-0.99839078411548014765525278348680286544429555739338) } },{ { SC_(-0.36787920275286323246177744294982403516769409179688), SC_(-0.99886193379608135520603487963907992157933985302350) } },{ { SC_(-0.36787932196215278324302744294982403516769409179688), SC_(-0.99919517626703684624524893082905669989578841060892) } },{ { SC_(-0.36787938156679755863365244294982403516769409179688), SC_(-0.99943085896775657378245957087668418410735469441835) } },{ { SC_(-0.36787941136911994632896494294982403516769409179688), SC_(-0.99959753415605033951327478977234592072050509074480) } },{ { SC_(-0.36787942627028114017662119294982403516769409179688), SC_(-0.99971540249082798050505534900918173321899800190957) } },{ { SC_(-0.36787943372086173710044931794982403516769409179688), SC_(-0.99979875358003464529770521637722571161846456343102) } },{ { SC_(-0.36787943744615203556236338044982403516769409179688), SC_(-0.99985769449598686744630754715710430111838645655608) } },{ { SC_(-0.36787943930879718479332041169982403516769409179688), SC_(-0.99989937341527312969776294577792175610005161268265) } },{ { SC_(-0.36787944024011975940879892732482403516769409179688), SC_(-0.99992884556078314715423832743355922518662235135757) } },{ { SC_(-0.36787944070578104671653818513732403516769409179688), SC_(-0.99994968586433278794146581248117772412549843583586) } },{ { SC_(-0.36787944093861169037040781404357403516769409179688), SC_(-0.99996442235919152892644019456912452486892832990114) } },{ { SC_(-0.36787944105502701219734262849669903516769409179688), SC_(-0.99997484272221444495021480907850566954322542216868) } },{ { SC_(-0.36787944111323467311081003572326153516769409179688), SC_(-0.99998221107553951227244139186618591264285119372063) } },{ { SC_(-0.36787944114233850356754373933654278516769409179688), SC_(-0.99998742131038091608107093454795869661238860012568) } },{ { SC_(-0.36787944115689041879591059114318341016769409179688), SC_(-0.99999110551424805741455916942650424910940130482916) } },{ { SC_(-0.36787944116416637641009401704650372266769409179688), SC_(-0.99999371064603396347995131962984747427523504609782) } },{ { SC_(-0.36787944116780435521718572999816387891769409179688), SC_(-0.99999555275622895023796382943893319302015254415029) } },{ { SC_(-0.36787944116962334462073158647399395704269409179688), SC_(-0.99999685532777825691586263781552103878671869687024) } },{ { SC_(-0.36787944117053283932250451471190899610519409179688), SC_(-0.99999777638786151731498560321162974199505119200634) } }
 252    }};
 253    T tolerance = boost::math::tools::epsilon<T>() * 3;
 254    if (boost::math::tools::epsilon<T>() <= boost::math::tools::epsilon<long double>())
 255       tolerance *= 5e5;
 256    T endpoint = -boost::math::constants::exp_minus_one<T>();
 257    for (unsigned i = 0; i < wolfram_test_near_singularity_data.size(); ++i)
 258    {
 259       if (wolfram_test_near_singularity_data[i][0] <= endpoint)
 260          break;
 261       else
 262          BOOST_CHECK_CLOSE_FRACTION(boost::math::lambert_w0(T(wolfram_test_near_singularity_data[i][0])), T(wolfram_test_near_singularity_data[i][1]), tolerance);
 263    }
 264 }
 265
 266 template <>
 267 void wolfram_test_near_singularity<float>()
 268 {
 269    //
 270    // Spot values near the singularity with inputs truncated to float precision,
 271    // from http://www.wolframalpha.com/input/?i=TABLE%5B%5BN%5BROUND%5B-1%2Fe%2B2%5E-i,+2%5E-23%5D,+50%5D,+N%5BLambertW%5BROUND%5B-1%2Fe+%2B+2%5E-i,+2%5E-23%5D%5D,+50%5D%5D,+%7Bi,+2,+40%7D%5D
 272    //
 273    static const boost::array<boost::array<float, 2>, 39> wolfram_test_near_singularity_data =
 274    {{
 275       {{ -0.11787939071655273437500000000000000000000000000000f, -0.13490440151978599948261696847702203722148729212591f }},{{ -0.24287939071655273437500000000000000000000000000000f, -0.34187230524883404685074938529655332889057132590877f }},{{ -0.30537939071655273437500000000000000000000000000000f, -0.50704515484245965628066570100405225451296978841169f }},{{ -0.33662939071655273437500000000000000000000000000000f, -0.63562295482810970976475066480034941107064440641758f }},{{ -0.35225439071655273437500000000000000000000000000000f, -0.73357162334066102207977288738307124189083069773180f }},{{ -0.36006689071655273437500000000000000000000000000000f, -0.80685854013946199386910756662972252220827924037205f }},{{ -0.36397314071655273437500000000000000000000000000000f, -0.86091065811941702413570870801021404654934249886505f }},{{ -0.36592626571655273437500000000000000000000000000000f, -0.90033443111682454984393817004965279949925483847744f }},{{ -0.36690282821655273437500000000000000000000000000000f, -0.92884710067602836873486989954484681592392882968841f }},{{ -0.36739110946655273437500000000000000000000000000000f, -0.94933939406123900376318336910404763737960907662666f }},{{ -0.36763525009155273437500000000000000000000000000000f, -0.96399956611859464483214118051190513364901860207328f }},{{ -0.36775732040405273437500000000000000000000000000000f, -0.97445213361280651797731195324654593603807971082292f }},{{ -0.36781835556030273437500000000000000000000000000000f, -0.98188628650256330812037232517657284107351472091741f }},{{ -0.36784887313842773437500000000000000000000000000000f, -0.98716379155663346207408852364078406478772014890806f }},{{ -0.36786413192749023437500000000000000000000000000000f, -0.99090459761086986284393759319956676727684106186028f }},{{ -0.36787176132202148437500000000000000000000000000000f, -0.99355229825129408828026714426677096743753950457546f }},{{ -0.36787557601928710937500000000000000000000000000000f, -0.99542297991285328482403963994064328331346049089419f }},{{ -0.36787748336791992187500000000000000000000000000000f, -0.99674107062291256263133271694520294422529881114769f }},{{ -0.36787843704223632812500000000000000000000000000000f, -0.99766536478294767461296564658785293377699068226332f }},{{ -0.36787891387939453125000000000000000000000000000000f, -0.99830783438342654552199009076049244789994050996944f }},{{ -0.36787915229797363281250000000000000000000000000000f, -0.99874733565614076859582844941545958416543067187493f }},{{ -0.36787927150726318359375000000000000000000000000000f, -0.99903989590053869025356285499889881633845057984872f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }},{{ -0.36787939071655273437500000000000000000000000000000f, -0.99947635367299698033494423493356278945921228277354f }}
 276    }};
 277    float tolerance = boost::math::tools::epsilon<float>() * 16;
 278    float endpoint = -boost::math::constants::exp_minus_one<float>();
 279    for (unsigned i = 0; i < wolfram_test_near_singularity_data.size(); ++i)
 280    {
 281       if (wolfram_test_near_singularity_data[i][0] <= endpoint)
 282          break;
 283       else
 284          BOOST_CHECK_CLOSE_FRACTION(boost::math::lambert_w0(wolfram_test_near_singularity_data[i][0]), wolfram_test_near_singularity_data[i][1], tolerance);
 285    }
 286 }
 287
 288 template <>
 289 void wolfram_test_near_singularity<double>()
 290 {
 291    //
 292    // Spot values near the singularity with inputs truncated to double precision,
 293    // from http://www.wolframalpha.com/input/?i=TABLE%5B%5BN%5BROUND%5B-1%2Fe%2B2%5E-i,+2%5E-23%5D,+50%5D,+N%5BLambertW%5BROUND%5B-1%2Fe+%2B+2%5E-i,+2%5E-23%5D%5D,+50%5D%5D,+%7Bi,+2,+40%7D%5D
 294    //
 295    static const boost::array<boost::array<double, 2>, 39> wolfram_test_near_singularity_data =
 296    {{
 297       {{ -0.11787944117144233402427744294982403516769409179688, -0.13490446826612137099065142885543349308605449591189 }},{{ -0.24287944117144233402427744294982403516769409179688, -0.34187241316000575559631565516533717918703951393828 }},{{ -0.30537944117144233402427744294982403516769409179688, -0.50704532478540670242736394530166187052909039079642 }},{{ -0.33662944117144233402427744294982403516769409179688, -0.63562321628494791544895212508757067989859372121549 }},{{ -0.35225444117144233402427744294982403516769409179688, -0.73357201771558852140844624841371893543359405991894 }},{{ -0.36006694117144233402427744294982403516769409179688, -0.80685912552602238275976720505076149562188136941981 }},{{ -0.36397319117144233402427744294982403516769409179688, -0.86091151614390373770305184939107560322835214525382 }},{{ -0.36592631617144233402427744294982403516769409179688, -0.90033567669608907987528169545609510444951296636737 }},{{ -0.36690287867144233402427744294982403516769409179688, -0.92884889586304130900291705545970353898661233095513 }},{{ -0.36739115992144233402427744294982403516769409179688, -0.94934196763921122756108351994184213101752011076782 }},{{ -0.36763530054644233402427744294982403516769409179688, -0.96400324129495105632485735566132352543383271582526 }},{{ -0.36775737085894233402427744294982403516769409179688, -0.97445736712728703357755243595334553847237474201138 }},{{ -0.36781840601519233402427744294982403516769409179688, -0.98189372378619472154195350108189165241865132390473 }},{{ -0.36784892359331733402427744294982403516769409179688, -0.98717434434269671591894280580432721487757138768109 }},{{ -0.36786418238237983402427744294982403516769409179688, -0.99091955260257317141206161906086819616043312707614 }},{{ -0.36787181177691108402427744294982403516769409179688, -0.99357346775773151586057357459040504547191256911173 }},{{ -0.36787562647417670902427744294982403516769409179688, -0.99545290640175819861266174073519228782773422561472 }},{{ -0.36787753382280952152427744294982403516769409179688, -0.99678329264937600678258333756796350065436689760936 }},{{ -0.36787848749712592777427744294982403516769409179688, -0.99772473035978895659981485126201758865515569761514 }},{{ -0.36787896433428413089927744294982403516769409179688, -0.99839078411548014765525278348680286544429555739338 }},{{ -0.36787920275286323246177744294982403516769409179688, -0.99886193379608135520603487963907992157933985302350 }},{{ -0.36787932196215278324302744294982403516769409179688, -0.99919517626703684624524893082905669989578841060892 }},{{ -0.36787938156679755863365244294982403516769409179688, -0.99943085896775657378245957087668418410735469441835 }},{{ -0.36787941136911994632896494294982403516769409179688, -0.99959753415605033951327478977234592072050509074480 }},{{ -0.36787942627028114017662119294982403516769409179688, -0.99971540249082798050505534900918173321899800190957 }},{{ -0.36787943372086173710044931794982403516769409179688, -0.99979875358003464529770521637722571161846456343102 }},{{ -0.36787943744615203556236338044982403516769409179688, -0.99985769449598686744630754715710430111838645655608 }},{{ -0.36787943930879718479332041169982403516769409179688, -0.99989937341527312969776294577792175610005161268265 }},{{ -0.36787944024011975940879892732482403516769409179688, -0.99992884556078314715423832743355922518662235135757 }},{{ -0.36787944070578104671653818513732403516769409179688, -0.99994968586433278794146581248117772412549843583586 }},{{ -0.36787944093861169037040781404357403516769409179688, -0.99996442235919152892644019456912452486892832990114 }},{{ -0.36787944105502701219734262849669903516769409179688, -0.99997484272221444495021480907850566954322542216868 }},{{ -0.36787944111323467311081003572326153516769409179688, -0.99998221107553951227244139186618591264285119372063 }},{{ -0.36787944114233850356754373933654278516769409179688, -0.99998742131038091608107093454795869661238860012568 }},{{ -0.36787944115689041879591059114318341016769409179688, -0.99999110551424805741455916942650424910940130482916 }},{{ -0.36787944116416637641009401704650372266769409179688, -0.99999371064603396347995131962984747427523504609782 }},{{ -0.36787944116780435521718572999816387891769409179688, -0.99999555275622895023796382943893319302015254415029 }},{{ -0.36787944116962334462073158647399395704269409179688, -0.99999685532777825691586263781552103878671869687024 }},{{ -0.36787944117053283932250451471190899610519409179688, -0.99999777638786151731498560321162974199505119200634 }}
 298    }};
 299    double tolerance = boost::math::tools::epsilon<double>() * 5;
 300    if (std::numeric_limits<double>::digits >= std::numeric_limits<long double>::digits)
 301       tolerance *= 1e5;
 302    else if (std::numeric_limits<double>::digits * 2 >= std::numeric_limits<long double>::digits)
 303       tolerance *= 5e4;
 304    double endpoint = -boost::math::constants::exp_minus_one<double>();
 305    for (unsigned i = 0; i < wolfram_test_near_singularity_data.size(); ++i)
 306    {
 307       if (wolfram_test_near_singularity_data[i][0] <= endpoint)
 308          break;
 309       else
 310          BOOST_CHECK_CLOSE_FRACTION(boost::math::lambert_w0(wolfram_test_near_singularity_data[i][0]), wolfram_test_near_singularity_data[i][1], tolerance);
 311    }
 312 }
 313
 314 template <class RealType>
 315 void test_spots(RealType)
 316 {
 317   // (Unused Parameter value, arbitrarily zero, only communicates the floating point type).
 318   // test_spots(0.F); test_spots(0.); test_spots(0.L);
 319
 320   using boost::math::lambert_w0;
 321   using boost::math::lambert_wm1;
 322   using boost::math::constants::exp_minus_one;
 323   using boost::math::constants::e;
 324   using boost::math::policies::policy;
 325
 326   /*  Example of an exception-free 'ignore_all' policy (possibly ill-advised?).
 327   */
 328   typedef policy <
 329     boost::math::policies::domain_error<boost::math::policies::ignore_error>,
 330     boost::math::policies::overflow_error<boost::math::policies::ignore_error>,
 331     boost::math::policies::underflow_error<boost::math::policies::ignore_error>,
 332     boost::math::policies::denorm_error<boost::math::policies::ignore_error>,
 333     boost::math::policies::pole_error<boost::math::policies::ignore_error>,
 334     boost::math::policies::evaluation_error<boost::math::policies::ignore_error>
 335   > ignore_all_policy;
 336
 337 //  Test some bad parameters to the function, with default policy and also with ignore_all policy.
 338 #ifndef BOOST_NO_EXCEPTIONS
 339   BOOST_CHECK_THROW(lambert_w0<RealType>(-1.), std::domain_error);
 340   BOOST_CHECK_THROW(lambert_wm1<RealType>(-1.), std::domain_error);
 341   if (std::numeric_limits<RealType>::has_quiet_NaN)
 342   {
 343      BOOST_CHECK_THROW(lambert_w0<RealType>(std::numeric_limits<RealType>::quiet_NaN()), std::domain_error); // Would be NaN.
 344      //BOOST_CHECK_EQUAL(lambert_w0<RealType>(std::numeric_limits<RealType>::quiet_NaN(), ignore_all_policy()), std::numeric_limits<RealType>::quiet_NaN()); // Should be NaN.
 345      // Fails as NaN != NaN by definition.
 346      BOOST_CHECK(boost::math::isnan(lambert_w0<RealType>(std::numeric_limits<RealType>::quiet_NaN(), ignore_all_policy())));
 347      //BOOST_MATH_CHECK_EQUAL(boost::math::lambert_w0<RealType>(std::numeric_limits<RealType>::infinity(), ignore_all_policy()), std::numeric_limits<RealType::infinity()); // infinity.
 348   }
 349
 350   // BOOST_CHECK_THROW(lambert_w0<RealType>(std::numeric_limits<RealType>::infinity()), std::domain_error); // Was if infinity should throw, now infinity.
 351   BOOST_CHECK_THROW(lambert_w0<RealType>(-static_cast<RealType>(0.4)), std::domain_error); // Would be complex.
 352
 353 #else // No exceptions, so set policy to ignore and check result is NaN.
 354   BOOST_MATH_CHECK_EQUAL(boost::math::lambert_w0<RealType>(std::numeric_limits<RealType>::quiet_NaN(), ignore_all_policy()), std::numeric_limits<RealType::quiet_NaN()); // NaN.
 355   BOOST_MATH_CHECK_EQUAL(boost::math::lambert_w0<RealType>(std::numeric_limits<RealType>::infinity(), ignore_all_policy()), std::numeric_limits<RealType::infinity()); // infinity.
 356   BOOST_MATH_CHECK_EQUAL(boost::math::lambert_w0<RealType>(std::numeric_limits<RealType>::infinity(), ignore_all_policy()), std::numeric_limits<RealType::infinity()); // infinity.
 357 #endif
 358
 359   std::cout << "\nTesting type " << typeid(RealType).name() << std::endl;
 360   int epsilons = 2;
 361   if (std::numeric_limits<RealType>::digits > 53)
 362   { // Multiprecision types.
 363     epsilons *= 8; // (Perhaps needed because need slightly longer (55) reference values?).
 364   }
 365   RealType tolerance = boost::math::tools::epsilon<RealType>() * epsilons; // 2 eps as a fraction.
 366   std::cout << "Tolerance " << epsilons << " * epsilon == " << tolerance << std::endl;
 367 #ifndef BOOST_NO_CXX11_NUMERIC_LIMITS
 368   std::cout << "Precision " << std::numeric_limits<RealType>::digits10 << " decimal digits, max_digits10 = " << std::numeric_limits <RealType>::max_digits10<< std::endl;
 369   // std::cout.precision(std::numeric_limits<RealType>::digits10);
 370   std::cout.precision(std::numeric_limits <RealType>::max_digits10);
 371 #endif
 372   std::cout.setf(std::ios_base::showpoint);  // show trailing significant zeros.
 373   std::cout << "-exp(-1) = " << -exp_minus_one<RealType>() << std::endl;
 374
 375   wolfram_test_near_singularity<RealType>();
 376   wolfram_test_large<RealType>();
 377   wolfram_test_small_neg<RealType>();
 378   wolfram_test_small_pos<RealType>();
 379   wolfram_test_moderate_values<RealType>();
 380
 381   // Test at singularity.
 382   // RealType test_value = BOOST_MATH_TEST_VALUE(RealType, -0.36787944117144232159552377016146086744581113103176783450783680169746149574489980335714727434591964374662732527);
 383   RealType singular_value = -exp_minus_one<RealType>();
 384   // -exp(-1) = -0.36787944117144232159552377016146086744581113103176783450783680169746149574489980335714727434591964374662732527
 385   // lambert_w0[-0.367879441171442321595523770161460867445811131031767834] == -1
 386   //           -0.36787945032119751
 387   RealType minus_one_value = BOOST_MATH_TEST_VALUE(RealType, -1.);
 388   //std::cout << "singular_value " << singular_value << ", expected Lambert W = " << minus_one_value << std::endl;
 389
 390   BOOST_CHECK_CLOSE_FRACTION( // Check -exp(-1) = -0.367879450 = -1max
 391     lambert_w0(singular_value),
 392     minus_one_value,
 393     tolerance);  // OK
 394
 395   BOOST_CHECK_CLOSE_FRACTION(  // Check -exp(-1) ~= -0.367879450 == -1
 396     lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -0.36787944117144232159552377016146086744581113103176783450783680169746149574489980335714727434591964374662732527)),
 397     BOOST_MATH_TEST_VALUE(RealType, -1.),
 398     tolerance);
 399
 400   BOOST_CHECK_CLOSE_FRACTION(  // Check -exp(-1) ~= -0.367879450 == -1
 401     lambert_w0<RealType>(-exp_minus_one<RealType>()),
 402     BOOST_MATH_TEST_VALUE(RealType, -1.),
 403     tolerance);
 404
 405   // Tests with some spot values computed using
 406   // https://www.wolframalpha.com/input
 407   // For example: N[lambert_w[1], 50] outputs:
 408   // 0.56714329040978387299996866221035554975381578718651
 409
 410   // At branch junction singularity.
 411   BOOST_CHECK_CLOSE_FRACTION(  // Check -exp(-1) ~= -0.367879450 == -1
 412     lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -0.36787944117144232159552377016146086744581113103176783450783680169746149574489980335714727434591964374662732527)),
 413     BOOST_MATH_TEST_VALUE(RealType, -1.),
 414     tolerance);
 415
 416     BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 0.1)),
 417     BOOST_MATH_TEST_VALUE(RealType, 0.091276527160862264299895721423179568653119224051472),
 418     // Output from https://www.wolframalpha.com/input/?i=lambert_w0(0.2)
 419     tolerance);
 420
 421   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 0.2)),
 422     BOOST_MATH_TEST_VALUE(RealType, 0.16891597349910956511647490370581839872844691351073),
 423     // Output from https://www.wolframalpha.com/input/?i=lambert_w0(0.2)
 424     tolerance);
 425
 426   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 0.5)),
 427     BOOST_MATH_TEST_VALUE(RealType, 0.351733711249195826024909300929951065171464215517111804046),
 428     // Output from https://www.wolframalpha.com/input/?i=lambert_w0(0.5)
 429     tolerance);
 430
 431   BOOST_CHECK_CLOSE_FRACTION(
 432     lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 1.)),
 433     BOOST_MATH_TEST_VALUE(RealType, 0.56714329040978387299996866221035554975381578718651),
 434    // Output from https://www.wolframalpha.com/input/?i=lambert_w0(1)
 435    tolerance);
 436
 437   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 2.)),
 438     BOOST_MATH_TEST_VALUE(RealType, 0.852605502013725491346472414695317466898453300151403508772),
 439     // Output from https://www.wolframalpha.com/input/?i=lambert_w0(2.)
 440     tolerance);
 441
 442   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 3.)),
 443     BOOST_MATH_TEST_VALUE(RealType, 1.049908894964039959988697070552897904589466943706341452932),
 444     // Output from https://www.wolframalpha.com/input/?i=lambert_w0(3.)
 445     tolerance);
 446
 447   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 5.)),
 448     BOOST_MATH_TEST_VALUE(RealType, 1.326724665242200223635099297758079660128793554638047479789),
 449     // Output from https://www.wolframalpha.com/input/?i=lambert_w0(0.5)
 450     tolerance);
 451
 452   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 6.)),
 453     BOOST_MATH_TEST_VALUE(RealType, 1.432404775898300311234078007212058694786434608804302025655),
 454     // Output from https://www.wolframalpha.com/input/?i=lambert_w0(6)
 455     tolerance);
 456
 457   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 100.)),
 458     BOOST_MATH_TEST_VALUE(RealType, 3.3856301402900501848882443645297268674916941701578),
 459     // Output from https://www.wolframalpha.com/input/?i=lambert_w0(100)
 460     tolerance);
 461
 462   if (std::numeric_limits<RealType>::has_infinity)
 463   {
 464     BOOST_CHECK_THROW(lambert_w0(std::numeric_limits<RealType>::infinity()), std::overflow_error); // If should throw exception for infinity.
 465     //BOOST_CHECK_EQUAL(lambert_w0(std::numeric_limits<RealType>::infinity()), +std::numeric_limits<RealType>::infinity()); // message is:
 466     // Error in "test_types": class boost::exception_detail::clone_impl<struct boost::exception_detail::error_info_injector<class std::overflow_error> > :
 467     // Error in function boost::math::lambert_w0<RealType>(<RealType>) : Argument z is infinite!
 468     //BOOST_CHECK_EQUAL(lambert_w0(std::numeric_limits<RealType>::infinity()), +std::numeric_limits<RealType>::infinity()); // If infinity allowed.
 469     BOOST_CHECK_THROW(lambert_wm1(std::numeric_limits<RealType>::infinity()), std::domain_error); // Infinity NOT allowed at all (not an edge case).
 470   }
 471   if (std::numeric_limits<RealType>::has_quiet_NaN)
 472   { // Argument Z == NaN is always an throwable error for both branches.
 473     // BOOST_CHECK_EQUAL(lambert_w0(std::numeric_limits<RealType>::quiet_NaN()), +std::numeric_limits<RealType>::infinity()); // message is:
 474     // Error in function boost::math::lambert_w0<RealType>(<RealType>): Argument z is NaN!
 475     BOOST_CHECK_THROW(lambert_w0(std::numeric_limits<RealType>::quiet_NaN()), std::domain_error);
 476     BOOST_CHECK_THROW(lambert_wm1(std::numeric_limits<RealType>::quiet_NaN()), std::domain_error);
 477   }
 478
 479   // denorm - but might be == min or zero?
 480   if (std::numeric_limits<RealType>::has_denorm == true)
 481   { // Might also return infinity like z == 0?
 482     BOOST_CHECK_THROW(lambert_wm1(std::numeric_limits<RealType>::denorm_min()), std::overflow_error);
 483   }
 484
 485     // Tests of Lambert W-1 branch.
 486     BOOST_CHECK_CLOSE_FRACTION(  // Check -exp(-1) ~= -0.367879450 == -1 at the singularity branch point.
 487     lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -0.36787944117144232159552377016146086744581113103176783450783680169746149574489980335714727434591964374662732527)),
 488     BOOST_MATH_TEST_VALUE(RealType, -1.),
 489     tolerance);
 490
 491     // Near singularity and using series approximation.
 492     // N[productlog(-1, -0.36), 50] = -1.2227701339785059531429380734238623131735264411311
 493     BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -0.36)),
 494       BOOST_MATH_TEST_VALUE(RealType, -1.2227701339785059531429380734238623131735264411311),
 495     10 *   tolerance); // tolerance OK for quad
 496     // -1.2227701339785059531429380734238623131735264411311
 497     // -1.222770133978505953142938073423862313173526441131033
 498
 499     // Just using series approximation (switch at -0.35).
 500     // N[productlog(-0.351), 50] = -0.72398644140937651483634596143951001600417138085814
 501     BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -0.351)),
 502       BOOST_MATH_TEST_VALUE(RealType, -0.72398644140937651483634596143951001600417138085814),
 503       // 2 * tolerance); // Note 2 * tolerance for PB fukushima
 504     // got -0.723986441409376931150560229265736446 without Halley
 505     // exp -0.72398644140937651483634596143951001
 506     // got -0.72398644140937651483634596143951029 with Halley
 507      10 * tolerance); // expect -0.72398644140937651 float -0.723987103 needs 10 * tolerance
 508      // 2 * tolerance is fine for double and up.
 509     // Float is OK
 510
 511     // Same for W-1 branch
 512     BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -0.351)),
 513       BOOST_MATH_TEST_VALUE(RealType, -1.3385736984773431852492145715526995809854973408320),
 514       10 * tolerance); // 2 tolerance OK for quad
 515
 516     // Near singularity and NOT using series approximation (switch at -0.35)
 517     // N[productlog(-1, -0.34), 50]
 518     BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -0.34)),
 519       BOOST_MATH_TEST_VALUE(RealType, -1.4512014851325470735077533710339268100722032730024),
 520      10 * tolerance); // tolerance OK for quad
 521     //
 522
 523     // Decreasing z until near zero (small z) .
 524     //N[productlog(-1, -0.3), 50] = -1.7813370234216276119741702815127452608215583564545
 525     BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -0.3)),
 526      BOOST_MATH_TEST_VALUE(RealType, -1.7813370234216276119741702815127452608215583564545),
 527     2 * tolerance);
 528     //                               -1.78133702342162761197417028151274526082155835645446
 529
 530     //N[productlog(-1, -0.2), 50] = -2.5426413577735264242938061566618482901614749075294
 531   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -0.2)),
 532     BOOST_MATH_TEST_VALUE(RealType, -2.5426413577735264242938061566618482901614749075294),
 533    2 * tolerance);
 534
 535   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -0.1)),
 536     BOOST_MATH_TEST_VALUE(RealType, -3.577152063957297218409391963511994880401796257793),
 537     tolerance);
 538
 539      //N[productlog(-1, -0.01), 50] = -6.4727751243940046947410578927244880371043455902257
 540   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -0.01)),
 541     BOOST_MATH_TEST_VALUE(RealType, -6.4727751243940046947410578927244880371043455902257),
 542     tolerance);
 543
 544   //  N[productlog(-1, -0.001), 50] = -9.1180064704027401212583371820468142742704349737639
 545   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -0.001)),
 546     BOOST_MATH_TEST_VALUE(RealType, -9.1180064704027401212583371820468142742704349737639),
 547     tolerance);
 548
 549   //  N[productlog(-1, -0.000001), 50] = -16.626508901372473387706432163984684996461726803805
 550   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -0.000001)),
 551     BOOST_MATH_TEST_VALUE(RealType, -16.626508901372473387706432163984684996461726803805),
 552     tolerance);
 553
 554   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -1e-12)),
 555     BOOST_MATH_TEST_VALUE(RealType, -31.067172842017230842039496250208586707880448763222),
 556     tolerance);
 557
 558   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -1e-25)),
 559     BOOST_MATH_TEST_VALUE(RealType, -61.686695602074505366866968627049381352503620377944),
 560     tolerance);
 561
 562   // z nearly too small.
 563   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -2e-26)),
 564     BOOST_MATH_TEST_VALUE(RealType, -63.322302839923597803393585145387854867226970485197),
 565     tolerance* 2);
 566
 567   // z very nearly too small.  G(k=64) g[63] = -1.0264389699511303e-26 to using 1.027e-26
 568   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -1.027e-26)),
 569     BOOST_MATH_TEST_VALUE(RealType, -63.999444896732265186957073549916026532499356695343),
 570     tolerance);
 571   // So -64 is the most negative value that can be determined using lookup.
 572   // N[productlog(-1, -1.0264389699511303 * 10^-26 ), 50] -63.999999999999997947255011093606206983577811736472 == -64
 573   // G[k=64] = g[63] = -1.0264389699511303e-26
 574
 575   // z too small for G(k=64) g[63] = -1.0264389699511303e-26 to using 1.027e-26
 576   //  N[productlog(-1, -10 ^ -26), 50]    = -31.067172842017230842039496250208586707880448763222
 577   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -1e-26)),
 578     BOOST_MATH_TEST_VALUE(RealType, -64.026509628385889681156090340691637712441162092868),
 579     tolerance); //                  -64.0265121
 580
 581   if (std::numeric_limits<RealType>::has_infinity)
 582   {
 583     BOOST_CHECK_EQUAL(lambert_wm1(0), -std::numeric_limits<RealType>::infinity());
 584   }
 585   if (std::numeric_limits<RealType>::has_quiet_NaN)
 586   {
 587     // BOOST_CHECK_EQUAL(lambert_w0(std::numeric_limits<RealType>::quiet_NaN()), +std::numeric_limits<RealType>::infinity()); // message is:
 588     // Error in function boost::math::lambert_w0<RealType>(<RealType>): Argument z is NaN!
 589     BOOST_CHECK_THROW(lambert_wm1(std::numeric_limits<RealType>::quiet_NaN()), std::domain_error);
 590   }
 591
 592    // W0 Tests for too big and too small to use lookup table.
 593    // Exactly W = 64, not enough to be OK for lookup.
 594     BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 3.9904954117194348050619127737142206366920907815909119e+29)),
 595     BOOST_MATH_TEST_VALUE(RealType, 64.0),
 596     tolerance);
 597
 598     // Just below z for F[64]
 599     BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 3.99045411719434e+29)),
 600      BOOST_MATH_TEST_VALUE(RealType, 63.999989810930513468726486827408823607175844852495), tolerance);
 601     // Fails for quad_float -1.22277013397850595265
 602     //                      -1.22277013397850595319
 603
 604   // Just too big, so using log approx and Halley refinement.
 605   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 4e+29)),
 606     BOOST_MATH_TEST_VALUE(RealType, 64.002342375637950350970694519073803643686041499677),
 607     tolerance);
 608
 609   // Check at reduced precision.
 610   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 4e+29), policy<digits2<11> >()),
 611     BOOST_MATH_TEST_VALUE(RealType, 64.002342375637950350970694519073803643686041499677),
 612     0.00002);  // 0.00001 fails.
 613
 614   // Tests to ensure that all JM rational polynomials are being checked.
 615
 616   // 1st polynomal if (z < 0.5)   // 0.05 < z < 0.5
 617   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 0.49)),
 618     BOOST_MATH_TEST_VALUE(RealType, 0.3465058086974944293540338951489158955895910665452626949),
 619     tolerance);
 620   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 0.051)),
 621     BOOST_MATH_TEST_VALUE(RealType, 0.04858156174600359264950777241723801201748517590507517888),
 622     tolerance);
 623
 624   // 2st polynomal if 0.5 < z < 2
 625   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 0.51)),
 626     BOOST_MATH_TEST_VALUE(RealType, 0.3569144916935871518694242462560450385494399307379277704),
 627     tolerance);
 628
 629   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 1.9)),
 630     BOOST_MATH_TEST_VALUE(RealType, 0.8291763302658400337004358009672187071638421282477162293),
 631     tolerance);
 632
 633   // 3rd polynomials 2 < z < 6
 634   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 2.1)),
 635     BOOST_MATH_TEST_VALUE(RealType, 0.8752187586805470099843211502166029752154384079916131962),
 636     tolerance);
 637
 638   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 5.9)),
 639     BOOST_MATH_TEST_VALUE(RealType, 1.422521411785098213935338853943459424120416844150520831),
 640     tolerance);
 641
 642   // 4th polynomials 6 < z < 18
 643   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 6.1)),
 644     BOOST_MATH_TEST_VALUE(RealType, 1.442152194116056579987235881273412088690824214100254315),
 645     tolerance);
 646
 647   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 17.9)),
 648     BOOST_MATH_TEST_VALUE(RealType, 2.129100923757568114366514708174691237123820852409339147),
 649     tolerance);
 650
 651   // 5th polynomials if (z < 9897.12905874)  // 2.8 < log(z) < 9.2
 652   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 18.1)),
 653     BOOST_MATH_TEST_VALUE(RealType, 2.136665501382339778305178680563584563343639180897328666),
 654     tolerance);
 655
 656   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 9897.)),
 657     BOOST_MATH_TEST_VALUE(RealType, 7.222751047988674263127929506116648714752441161828893633),
 658     tolerance);
 659
 660   // 6th polynomials if (z < 7.896296e+13)  // 9.2 < log(z) <= 32
 661   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 9999.)),
 662     BOOST_MATH_TEST_VALUE(RealType, 7.231758181708737258902175236106030961433080976032516996),
 663     tolerance);
 664
 665   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 7.7e+13)),
 666     BOOST_MATH_TEST_VALUE(RealType, 28.62069643025822480911439831021393125282095606713326376),
 667     tolerance);
 668
 669   // 7th polynomial // 32 < log(z) < 100
 670   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 8.0e+18)),
 671     BOOST_MATH_TEST_VALUE(RealType, 39.84107480517853176296156400093560722439428484537515586),
 672     tolerance);
 673
 674   // Largest 32-bit float. (Larger values for other types tested using max())
 675   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 1.e38)),
 676     BOOST_MATH_TEST_VALUE(RealType, 83.07844821316409592720410446942538465411465113447713574),
 677     tolerance);
 678
 679   // Using z small series function if z < 0.05  if (z < -0.051)  -0.27 < z < -0.051
 680
 681   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -0.28)),
 682     BOOST_MATH_TEST_VALUE(RealType, -0.4307588745271127579165306568413721388196459822705155385),
 683     tolerance);
 684
 685   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -0.25)),
 686     BOOST_MATH_TEST_VALUE(RealType, -0.3574029561813889030688111040559047533165905550760120436),
 687     tolerance);
 688
 689   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, +0.25)),
 690     BOOST_MATH_TEST_VALUE(RealType, 0.2038883547022401644431818313271398701493524772101596350),
 691     tolerance);
 692
 693   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -0.051)), // just above 0.05 cutoff.
 694     BOOST_MATH_TEST_VALUE(RealType, -0.05382002772543396036830469500362485089791914689728115249),
 695     tolerance * 4);
 696
 697   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -0.05)), // at cutoff.
 698     BOOST_MATH_TEST_VALUE(RealType, -0.05270598355154634795995650617915721289427674396592395160),
 699     tolerance * 8);
 700
 701   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 0.049)), // Just below cutoff.
 702     BOOST_MATH_TEST_VALUE(RealType, 0.04676143671340832342497289393737051868103596756298863555),
 703     tolerance * 4);
 704
 705   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 0.01)),
 706     BOOST_MATH_TEST_VALUE(RealType, 0.009901473843595011885336326816570107953627746494917415483),
 707     tolerance);
 708
 709   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -0.01)),
 710     BOOST_MATH_TEST_VALUE(RealType, -0.01010152719853875327292018767138623973670903993475235877),
 711     tolerance);
 712
 713   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -0.049)),
 714     BOOST_MATH_TEST_VALUE(RealType, -0.05159448479219405354564920228913331280713177046648170658),
 715     tolerance * 8);
 716
 717   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 1e-6)),
 718     BOOST_MATH_TEST_VALUE(RealType, 9.999990000014999973333385416558666900096702096424344715e-7),
 719     tolerance);
 720
 721   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -1e-6)),
 722     BOOST_MATH_TEST_VALUE(RealType, -1.000001000001500002666671875010800023343107568372593753e-6),
 723     tolerance);
 724
 725   // Near Smallest float.
 726   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 1e-38)),
 727     BOOST_MATH_TEST_VALUE(RealType, 9.99999999999999999999999999999999999990000000000000000e-39),
 728     tolerance);
 729
 730   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -1e-38)),
 731     BOOST_MATH_TEST_VALUE(RealType, -1.000000000000000000000000000000000000010000000000000000e-38),
 732     tolerance);
 733
 734   // Similar 'too near zero' tests for W-1 branch.
 735   // lambert_wm1(-1.0264389699511283e-26) = -64.000000000000000
 736   // Exactly z for W=-64
 737   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -1.026438969951128225904695701851094643838952857740385870e-26)),
 738     BOOST_MATH_TEST_VALUE(RealType, -64.000000000000000000000000000000000000),
 739    2 * tolerance);
 740
 741   // Just more negative than G[64 max] = wm1zs[63] so can't use lookup table.
 742   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -1.5e-27)),
 743     BOOST_MATH_TEST_VALUE(RealType, -65.953279000145077719128800110134854577850889171784),
 744     tolerance); //                  -65.9532776
 745
 746   // Just less negative than G[64 max] = wm1zs[63] so can use lookup table.
 747   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -1.1e-26)),
 748     BOOST_MATH_TEST_VALUE(RealType, -63.929686062157630858625440758283127600360210072859),
 749     tolerance);
 750
 751    // N[productlog(-1, -10 ^ -26), 50]    = -31.067172842017230842039496250208586707880448763222
 752   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -1e-26)),
 753     BOOST_MATH_TEST_VALUE(RealType, -64.026509628385889681156090340691637712441162092868),
 754     tolerance);
 755
 756   // 1e-28 is too small
 757   //  N[productlog(-1, -10 ^ -28), 50]    = -31.067172842017230842039496250208586707880448763222
 758   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -1e-28)),
 759     BOOST_MATH_TEST_VALUE(RealType, -68.702163291525429160769761667024460023336801014578),
 760     tolerance);
 761
 762   // Check for overflow when using a double (including when using for approximate value for refinement for higher precision).
 763
 764   // N[productlog(-1, -10 ^ -30), 50]    = -73.373110313822976797067478758120874529181611813766
 765   //BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -1e30)),
 766   //  BOOST_MATH_TEST_VALUE(RealType, -73.373110313822976797067478758120874529181611813766),
 767   //  tolerance);
 768   //unknown location : fatal error : in "test_types" :
 769   //class boost::exception_detail::clone_impl<struct boost::exception_detail::error_info_injector<class std::domain_error> >
 770   //  : Error in function boost::math::lambert_wm1<RealType>(<RealType>) :
 771   //  Argument z = -1.00000002e+30 out of range(z < -exp(-1) = -3.6787944) for Lambert W - 1 branch!
 772
 773   BOOST_CHECK_THROW(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -1e30)), std::domain_error);
 774
 775   // Too negative
 776   BOOST_CHECK_THROW(lambert_wm1(RealType(-0.5)), std::domain_error);
 777
 778   // This fails for fixed_point type used for other tests because out of range?
 779     //BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 1.0e6)),
 780     //BOOST_MATH_TEST_VALUE(RealType, 11.383358086140052622000156781585004289033774706019),
 781     //// Output from https://www.wolframalpha.com/input/?i=lambert_w0(1e6)
 782     //// tolerance * 1000); // fails for fixed_point type exceeds 0.00015258789063
 783     //  // 15.258789063
 784     //  // 11.383346558
 785     // tolerance * 100000);
 786   // So need to use some spot tests for specific types, or use a bigger fixed_point type.
 787
 788   // Check zero.
 789   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 0.0)),
 790     BOOST_MATH_TEST_VALUE(RealType, 0.0),
 791     tolerance);
 792   // these fail for cpp_dec_float_50
 793   // 'boost::multiprecision::detail::expression<boost::multiprecision::detail::negate,boost::multiprecision::number<boost::multiprecision::backends::cpp_dec_float<50,int32_t,void>,boost::multiprecision::et_on>,void,void,void>'
 794   // : no appropriate default constructor available
 795   // TODO ???????????
 796
 797  } // template <class RealType>void test_spots(RealType)
 798
 799 BOOST_AUTO_TEST_CASE( test_types )
 800 {
 801   BOOST_MATH_CONTROL_FP;
 802   // BOOST_TEST_MESSAGE output only appears if command line has --log_level="message"
 803   // or call set_threshold_level function:
 804   boost::unit_test_framework::unit_test_log.set_threshold_level(boost::unit_test_framework::log_messages);
 805   BOOST_TEST_MESSAGE("\nTest Lambert W function for several types.");
 806   BOOST_TEST_MESSAGE(show_versions());  // Full version of Boost, STL and compiler info.
 807 #ifndef BOOST_MATH_TEST_MULTIPRECISION
 808   // Fundamental built-in types:
 809   test_spots(0.0F); // float
 810   test_spots(0.0); // double
 811 #ifndef BOOST_MATH_NO_LONG_DOUBLE_MATH_FUNCTIONS
 812   if (sizeof(long double) > sizeof(double))
 813   { // Avoid pointless re-testing if double and long double are identical (for example, MSVC).
 814     test_spots(0.0L); // long double
 815   }
 816   test_spots(boost::math::concepts::real_concept(0));
 817 #endif
 818
 819   #else // BOOST_MATH_TEST_MULTIPRECISION
 820   // Multiprecision types:
 821 #if BOOST_MATH_TEST_MULTIPRECISION == 1
 822   test_spots(static_cast<boost::multiprecision::cpp_bin_float_double_extended>(0));
 823 #endif
 824 #if BOOST_MATH_TEST_MULTIPRECISION == 2
 825   test_spots(static_cast<boost::multiprecision::cpp_bin_float_quad>(0));
 826 #endif
 827 #if BOOST_MATH_TEST_MULTIPRECISION == 3
 828   test_spots(static_cast<boost::multiprecision::cpp_bin_float_50>(0));
 829 #endif
 830 #endif // ifdef BOOST_MATH_TEST_MULTIPRECISION
 831
 832   #ifdef BOOST_MATH_TEST_FLOAT128
 833    std::cout << "\nBOOST_MATH_TEST_FLOAT128 defined for float128 tests." << std::endl;
 834
 835 #ifdef BOOST_HAS_FLOAT128
 836   //  GCC and Intel only.
 837   // Requires link to libquadmath library, see
 838   // http://www.boost.org/doc/libs/release/libs/multiprecision/doc/html/boost_multiprecision/tut/floats/float128.html
 839   // for example:
 840   // C:\Program Files\mingw-w64\x86_64-7.2.0-win32-seh-rt_v5-rev1\mingw64\lib\gcc\x86_64-w64-mingw32\7.2.0\libquadmath.a
 841
 842   using boost::multiprecision::float128;
 843   std::cout << "BOOST_HAS_FLOAT128" << std::endl;
 844
 845   std::cout.precision(std::numeric_limits<float128>::max_digits10);
 846
 847   test_spots(static_cast<float128>(0));
 848 #endif // BOOST_HAS_FLOAT128
 849 #else
 850   std::cout << "\nBOOST_MATH_TEST_FLOAT128 NOT defined so NO float128 tests." << std::endl;
 851 #endif // #ifdef BOOST_MATH_TEST_FLOAT128
 852
 853 } // BOOST_AUTO_TEST_CASE( test_types )
 854
 855
 856 BOOST_AUTO_TEST_CASE( test_range_of_double_values )
 857 {
 858   using boost::math::constants::exp_minus_one;
 859   using boost::math::lambert_w0;
 860
 861   BOOST_TEST_MESSAGE("\nTest Lambert W function type double for range of values.");
 862
 863   // Want to test almost largest value.
 864   // test_value = (std::numeric_limits<RealType>::max)() / 4;
 865   // std::cout << std::setprecision(std::numeric_limits<RealType>::max_digits10) << "Max value = " << test_value << std::endl;
 866   // Can't use a test like this for all types because max_value depends on RealType
 867   // and thus the expected result of lambert_w0 does too.
 868   //BOOST_CHECK_CLOSE_FRACTION(lambert_w0<RealType>(test_value),
 869   //  BOOST_MATH_TEST_VALUE(RealType, ???),
 870   //  tolerance);
 871   // So this section just tests a single type, say IEEE 64-bit double, for a range of spot values.
 872
 873   typedef double RealType; // Some tests assume type is double.
 874
 875   int epsilons = 1;
 876   RealType tolerance = boost::math::tools::epsilon<RealType>() * epsilons; // 2 eps as a fraction.
 877   std::cout << "Tolerance " << epsilons << " * epsilon == " << tolerance << std::endl;
 878
 879 #ifndef BOOST_MATH_TEST_MULTIPRECISION
 880   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 1.0e-6)),
 881     BOOST_MATH_TEST_VALUE(RealType, 9.9999900000149999733333854165586669000967020964243e-7),
 882     // Output from https://www.wolframalpha.com/input/ N[lambert_w[1e-6],50])
 883     tolerance);
 884   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 0.0001)),
 885     BOOST_MATH_TEST_VALUE(RealType, 0.000099990001499733385405869000452213835767629477903460),
 886     // Output from https://www.wolframalpha.com/input/ N[lambert_w[0.001],50])
 887     tolerance);
 888   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 0.001)),
 889     BOOST_MATH_TEST_VALUE(RealType, 0.00099900149733853088995782787410778559957065467928884),
 890     // Output from https://www.wolframalpha.com/input/ N[lambert_w[0.001],50])
 891     tolerance);
 892   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 0.01)),
 893     BOOST_MATH_TEST_VALUE(RealType, 0.0099014738435950118853363268165701079536277464949174),
 894     // Output from https://www.wolframalpha.com/input/ N[lambert_w[0.01],50])
 895     tolerance * 25);  // <<< Needs a much bigger tolerance???
 896   // 0.0099014738435951096 this test max_digits10
 897   // 0.00990147384359511  digits10
 898   // 0.0099014738435950118  wolfram
 899   // 0.00990147384359501  wolfram  digits10
 900   // 0.0099014738435950119 N[lambert_w[0.01],17]
 901   // 0.00990147384359501   N[lambert_w[0.01],15] which really is more different than expected.
 902   // 0.00990728209160670  approx
 903   // 0.00990147384359511  previous
 904
 905   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 0.05)),
 906     BOOST_MATH_TEST_VALUE(RealType, 0.047672308600129374726388900514160870747062965933891),
 907     // Output from https://www.wolframalpha.com/input/ N[lambert_w[0.01],50])
 908     tolerance);
 909
 910   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 0.1)),
 911     BOOST_MATH_TEST_VALUE(RealType, 0.091276527160862264299895721423179568653119224051472),
 912     // Output from https://www.wolframalpha.com/input/ N[lambert_w[1],50])
 913     tolerance);
 914
 915   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 1.)),
 916     BOOST_MATH_TEST_VALUE(RealType, 0.56714329040978387299996866221035554975381578718651),
 917     // Output from https://www.wolframalpha.com/input/ N[lambert_w[1],50])
 918     tolerance);
 919
 920   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 2.)),
 921     BOOST_MATH_TEST_VALUE(RealType, 0.852605502013725491346472414695317466898453300151403508772),
 922     // Output from https://www.wolframalpha.com/input/?i=lambert_w0(2.)
 923     tolerance);
 924
 925   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 3.)),
 926     BOOST_MATH_TEST_VALUE(RealType, 1.049908894964039959988697070552897904589466943706341452932),
 927     // Output from https://www.wolframalpha.com/input/?i=lambert_w0(3.)
 928     tolerance);
 929
 930   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 5.)),
 931     BOOST_MATH_TEST_VALUE(RealType, 1.326724665242200223635099297758079660128793554638047479789),
 932     // Output from https://www.wolframalpha.com/input/?i=lambert_w0(0.5)
 933     tolerance);
 934
 935   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 6.)),
 936     BOOST_MATH_TEST_VALUE(RealType, 1.432404775898300311234078007212058694786434608804302025655),
 937     // Output from https://www.wolframalpha.com/input/?i=lambert_w0(6)
 938     tolerance);
 939
 940   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 10.)),
 941     BOOST_MATH_TEST_VALUE(RealType, 1.7455280027406993830743012648753899115352881290809),
 942     // Output from https://www.wolframalpha.com/input/ N[lambert_w[10],50])
 943     tolerance);
 944
 945     BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 100.)),
 946     BOOST_MATH_TEST_VALUE(RealType, 3.3856301402900501848882443645297268674916941701578),
 947     // Output from https://www.wolframalpha.com/input/?i=lambert_w0(100)
 948     tolerance);
 949
 950   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 1000.)),
 951     BOOST_MATH_TEST_VALUE(RealType, 5.2496028524015962271260563196973062825214723860596),
 952     // Output from https://www.wolframalpha.com/input/?i=lambert_w0(1000)
 953     tolerance);
 954
 955   // This fails for fixed_point type used for other tests because out of range of the type?
 956   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, 1.0e6)),
 957     BOOST_MATH_TEST_VALUE(RealType, 11.383358086140052622000156781585004289033774706019),
 958     // Output from https://www.wolframalpha.com/input/?i=lambert_w0(1e6)
 959     tolerance); //
 960
 961   // Tests for double only near the max and the singularity where Lambert_w estimates are less precise.
 962   if (std::numeric_limits<RealType>::is_specialized)
 963   { // is_specialized means that can use numeric_limits for tests.
 964     // Check near std::numeric_limits<>::max() for type.
 965     //std::cout << std::setprecision(std::numeric_limits<RealType>::max_digits10)
 966     //  << (std::numeric_limits<double>::max)()          // == 1.7976931348623157e+308
 967     //  << " " << (std::numeric_limits<double>::max)()/4 // == 4.4942328371557893e+307
 968     //  << std::endl;
 969
 970     // All these result in faulty error message
 971     // unknown location : fatal error : in "test_range_of_values": class boost::exception_detail::clone_impl<struct boost::exception_detail::error_info_injector<class std::domain_error> >: Error in function boost::math::lambert_w0<RealType>(<RealType>): Argument z = %1 too large.
 972     // I:\modular - boost\libs\math\test\test_lambert_w.cpp(456) : last checkpoint
 973
 974     BOOST_CHECK_CLOSE_FRACTION(lambert_w0(1.7976931348623157e+308 ), // max_value for IEEE 64-bit double.
 975       static_cast<double>(703.2270331047701868711791887193075929608934699575820028L),
 976       // N[productlog[0, 1.7976931348623157*10^308 /2],50] == 702.53487067487671916110655783739076368512998658347
 977       tolerance);
 978
 979     BOOST_CHECK_CLOSE_FRACTION(lambert_w0(1.7976931348623157e+308 / 2), // max_value/2 for IEEE 64-bit double.
 980       static_cast<double>(702.53487067487671916110655783739076368512998658347L),
 981       // N[productlog[0, 1.7976931348623157*10^308 /2],50] == 702.53487067487671916110655783739076368512998658347
 982       tolerance);
 983
 984     BOOST_CHECK_CLOSE_FRACTION(lambert_w0(1.7976931348623157e+308 /4), // near max_value/4 for IEEE 64-bit double.
 985       static_cast<double>(701.8427092142920014223182853764045476L),
 986       // N[productlog(0, 1.7976931348623157* 10^308 /4 ), 37] =701.8427092142920014223182853764045476
 987       // N[productlog(0, 0.25 * 1.7976931348623157*10^307), 37]
 988       tolerance);
 989
 990     BOOST_CHECK_CLOSE_FRACTION(lambert_w0(4.4942328371557893e+307), // max_value/4 for IEEE 64-bit double.
 991       static_cast<double>(701.84270921429200143342782556643059L),
 992       // N[lambert_w[4.4942328371557893e+307], 35]  == 701.8427092142920014334278255664305887
 993       // as a double == 701.83341468208209
 994       // Lambert computed 702.02379914670587
 995       0.000003); // OK Much less precise at the max edge???
 996
 997     BOOST_CHECK_CLOSE_FRACTION(lambert_w0((std::numeric_limits<double>::max)()), // max_value for IEEE 64-bit double.
 998       static_cast<double>(703.2270331047701868711791887193075930),
 999       // N[productlog(0, 1.7976931348623157* 10^308), 37] = 703.2270331047701868711791887193075930
1000       //                                                    703.22700325995515 lambert W
1001       //                                                    703.22703310477016  Wolfram
1002       tolerance * 2e8); // OK but much less accurate near max.
1003
1004   // Compare precisions very close to the singularity.
1005     // This test value is one epsilon close to the singularity at -exp(-1) * z
1006     // (below which the result has a non-zero imaginary part).
1007     RealType test_value = -exp_minus_one<RealType>();
1008     test_value += (std::numeric_limits<RealType>::epsilon() * 1);
1009     BOOST_CHECK_CLOSE_FRACTION(lambert_w0(test_value),
1010       BOOST_MATH_TEST_VALUE(RealType, -0.99999996349975895),
1011       tolerance * 1000000000);
1012     // -0.99999996788201051
1013     // -0.99999996349975895
1014     // Would not expect to get a result closer than sqrt(epsilon)?
1015   } //  if (std::numeric_limits<RealType>::is_specialized)
1016
1017   // Can only compare float_next for specific type T = double.
1018   // Comparison with Wolfram N[productlog(0,-0.36787944117144228 ), 17]
1019   // Note big loss of precision and big tolerance needed to pass.
1020   BOOST_CHECK_CLOSE_FRACTION(  // Check float_next(-exp(-1) )
1021     lambert_w0(BOOST_MATH_TEST_VALUE(double, -0.36787944117144228)),
1022     BOOST_MATH_TEST_VALUE(RealType, -0.99999998496215738),
1023     1e8 * tolerance); // diff 6.03558e-09 v 2.2204460492503131e-16
1024
1025    BOOST_CHECK_CLOSE_FRACTION(  // Check  float_next(float_next(-exp(-1) ))
1026     lambert_w0(BOOST_MATH_TEST_VALUE(double, -0.36787944117144222)),
1027     BOOST_MATH_TEST_VALUE(RealType, -0.99999997649828679),
1028     5e7 * tolerance);// diff 2.30785e-09 v 2.2204460492503131e-16
1029
1030   // Compare with previous PB/FK computations at double precision.
1031   BOOST_CHECK_CLOSE_FRACTION(  // Check float_next(-exp(-1) )
1032     lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -0.36787944117144228)),
1033     BOOST_MATH_TEST_VALUE(RealType, -0.99999997892657588),
1034     tolerance); // diff 6.03558e-09 v 2.2204460492503131e-16
1035
1036   BOOST_CHECK_CLOSE_FRACTION(  // Check  float_next(float_next(-exp(-1) ))
1037     lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -0.36787944117144222)),
1038     BOOST_MATH_TEST_VALUE(RealType, -0.99999997419043196),
1039     tolerance);// diff 2.30785e-09 v 2.2204460492503131e-16
1040
1041                // z increasingly close to singularity.
1042   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -0.36)),
1043      BOOST_MATH_TEST_VALUE(RealType, -0.8060843159708177782855213616209920019974599683466713016),
1044      2 * tolerance); // -0.806084335
1045
1046   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -0.365)),
1047      BOOST_MATH_TEST_VALUE(RealType, -0.8798200914159538111724840007674053239388642469453350954),
1048      5 * tolerance); // Note 5 * tolerance
1049
1050   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -0.3678)),
1051      BOOST_MATH_TEST_VALUE(RealType, -0.9793607149578284774761844434886481686055949229547379368),
1052      15 * tolerance); // Note 15 * tolerance when this close to singularity.
1053
1054                       // Just using series approximation (Fukushima switch at -0.35, but JM at 0.01 of singularity < -0.3679).
1055                       // N[productlog(-0.351), 50] = -0.72398644140937651483634596143951001600417138085814
1056                       // N[productlog(-0.351), 55] = -0.7239864414093765148363459614395100160041713808581379727
1057   BOOST_CHECK_CLOSE_FRACTION(lambert_w0(BOOST_MATH_TEST_VALUE(RealType, -0.351)),
1058      BOOST_MATH_TEST_VALUE(RealType, -0.72398644140937651483634596143951001600417138085814),
1059      10 * tolerance); // Note was 2 * tolerance
1060
1061                       // Check value just not using near_singularity series approximation (and using rational polynomial instead).
1062   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -0.3)),
1063      BOOST_MATH_TEST_VALUE(RealType, -1.7813370234216276119741702815127452608215583564545),
1064      // Output from https://www.wolframalpha.com/input/
1065      //N[productlog(-1, -0.3), 50] = -1.7813370234216276119741702815127452608215583564545
1066      tolerance);
1067
1068   // Using table lookup and schroeder with decreasing z to zero.
1069   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -0.2)),
1070      BOOST_MATH_TEST_VALUE(RealType, -2.5426413577735264242938061566618482901614749075294),
1071      // N[productlog[-1, -0.2],50] -2.5426413577735264242938061566618482901614749075294
1072      tolerance);
1073
1074   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -0.1)),
1075      BOOST_MATH_TEST_VALUE(RealType, -3.5771520639572972184093919635119948804017962577931),
1076      //N[productlog(-1, -0.1), 50] = -3.5771520639572972184093919635119948804017962577931
1077      tolerance);
1078
1079   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -0.001)),
1080      BOOST_MATH_TEST_VALUE(RealType, -9.1180064704027401212583371820468142742704349737639),
1081      //  N[productlog(-1, -0.001), 50] = -9.1180064704027401212583371820468142742704349737639
1082      tolerance);
1083
1084   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -0.000001)),
1085      BOOST_MATH_TEST_VALUE(RealType, -16.626508901372473387706432163984684996461726803805),
1086      //  N[productlog(-1, -0.000001), 50] = -16.626508901372473387706432163984684996461726803805
1087      tolerance);
1088
1089   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -1e-6)),
1090      BOOST_MATH_TEST_VALUE(RealType, -16.626508901372473387706432163984684996461726803805),
1091      //  N[productlog(-1, -10 ^ -6), 50] = -16.626508901372473387706432163984684996461726803805
1092      tolerance);
1093
1094   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -1.0e-26)),
1095      BOOST_MATH_TEST_VALUE(RealType, -64.026509628385889681156090340691637712441162092868),
1096      // Output from https://www.wolframalpha.com/input/
1097      // N[productlog(-1, -1 * 10^-26 ), 50] = -64.026509628385889681156090340691637712441162092868
1098      tolerance);
1099
1100   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -2e-26)),
1101      BOOST_MATH_TEST_VALUE(RealType, -63.322302839923597803393585145387854867226970485197),
1102      // N[productlog[-1, -2*10^-26],50] = -63.322302839923597803393585145387854867226970485197
1103      tolerance * 2);
1104
1105   // Smaller than lookup table, so must use approx and Halley refinements.
1106   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -1e-30)),
1107      BOOST_MATH_TEST_VALUE(RealType, -73.373110313822976797067478758120874529181611813766),
1108      //  N[productlog(-1, -10 ^ -30), 50]    = -73.373110313822976797067478758120874529181611813766
1109      tolerance);
1110
1111   // std::numeric_limits<RealType>::min
1112 #ifndef BOOST_NO_CXX11_NUMERIC_LIMITS
1113   std::cout.precision(std::numeric_limits<RealType>::max_digits10);
1114 #endif
1115   std::cout << "(std::numeric_limits<RealType>::min)() " << (std::numeric_limits<RealType>::min)() << std::endl;
1116
1117   BOOST_CHECK_CLOSE_FRACTION(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, -2.2250738585072014e-308)),
1118      BOOST_MATH_TEST_VALUE(RealType, -714.96865723796647086868547560654825435542227693935),
1119      // N[productlog[-1, -2.2250738585072014e-308],50] =  -714.96865723796647086868547560654825435542227693935
1120      tolerance);
1121
1122   // For z = 0, W = -infinity
1123   if (std::numeric_limits<RealType>::has_infinity)
1124   {
1125      BOOST_CHECK_EQUAL(lambert_wm1(BOOST_MATH_TEST_VALUE(RealType, 0.)),
1126         -std::numeric_limits<RealType>::infinity());
1127   }
1128
1129 #elif BOOST_MATH_TEST_MULTIPRECISION == 2
1130
1131   // Comparison with Wolfram N[productlog(0,-0.36787944117144228 ), 17]
1132   // Using conversion from double to higher precision cpp_bin_float_quad
1133   using boost::multiprecision::cpp_bin_float_quad;
1134   BOOST_CHECK_CLOSE_FRACTION(  // Check float_next(-exp(-1) )
1135     lambert_w0(BOOST_MATH_TEST_VALUE(cpp_bin_float_quad, -0.36787944117144228)),
1136     BOOST_MATH_TEST_VALUE(cpp_bin_float_quad, -0.99999998496215738),
1137     tolerance); // OK
1138
1139   BOOST_CHECK_CLOSE_FRACTION(  // Check  float_next(float_next(-exp(-1) ))
1140     lambert_w0(BOOST_MATH_TEST_VALUE(cpp_bin_float_quad, -0.36787944117144222)),
1141     BOOST_MATH_TEST_VALUE(cpp_bin_float_quad, -0.99999997649828679),
1142     tolerance);// OK
1143 #endif
1144 } // BOOST_AUTO_TEST_CASE(test_range_of_double_values)
1145