Subject: [Boost-commit] svn:boost r77758 - in sandbox/big_number: boost/multiprecision/detail/functions libs/multiprecision/test
From: john_at_[hidden]
Date: 2012-04-04 13:04:01
Author: johnmaddock
Date: 2012-04-04 13:04:00 EDT (Wed, 04 Apr 2012)
New Revision: 77758
URL: http://svn.boost.org/trac/boost/changeset/77758
Log:
Add string versions of constants for faster (and more accurate) constant initialization.
Text files modified:
sandbox/big_number/boost/multiprecision/detail/functions/constants.hpp | 70 ++++++++++++++++++++++++++++++++++++++++
sandbox/big_number/libs/multiprecision/test/Jamfile.v2 | 12 +++++-
2 files changed, 80 insertions(+), 2 deletions(-)
Modified: sandbox/big_number/boost/multiprecision/detail/functions/constants.hpp
==============================================================================
--- sandbox/big_number/boost/multiprecision/detail/functions/constants.hpp (original)
+++ sandbox/big_number/boost/multiprecision/detail/functions/constants.hpp 2012-04-04 13:04:00 EDT (Wed, 04 Apr 2012)
@@ -17,6 +17,29 @@
typedef typename mpl::deref<iter_type>::type ui_type;
typedef typename mpl::front<typename T::signed_types>::type si_type;
+ //
+ // String value with 1100 digits:
+ //
+ static const char* string_val = "0."
+ "6931471805599453094172321214581765680755001343602552541206800094933936219696947156058633269964186875"
+ "4200148102057068573368552023575813055703267075163507596193072757082837143519030703862389167347112335"
+ "0115364497955239120475172681574932065155524734139525882950453007095326366642654104239157814952043740"
+ "4303855008019441706416715186447128399681717845469570262716310645461502572074024816377733896385506952"
+ "6066834113727387372292895649354702576265209885969320196505855476470330679365443254763274495125040606"
+ "9438147104689946506220167720424524529612687946546193165174681392672504103802546259656869144192871608"
+ "2938031727143677826548775664850856740776484514644399404614226031930967354025744460703080960850474866"
+ "3852313818167675143866747664789088143714198549423151997354880375165861275352916610007105355824987941"
+ "4729509293113897155998205654392871700072180857610252368892132449713893203784393530887748259701715591"
+ "0708823683627589842589185353024363421436706118923678919237231467232172053401649256872747782344535347"
+ "6481149418642386776774406069562657379600867076257199184734022651462837904883062033061144630073719489";
+ //
+ // Check if we can just construct from string:
+ //
+ if(digits * 1000L / 301 < 1100)
+ {
+ num = string_val;
+ return;
+ }
num = static_cast<ui_type>(1180509120u);
T denom, next_term, temp;
denom = static_cast<ui_type>(1277337600u);
@@ -51,6 +74,30 @@
{
typedef typename mpl::front<typename T::unsigned_types>::type ui_type;
typedef typename mpl::front<typename T::float_types>::type real_type;
+ //
+ // 1100 digits in string form:
+ //
+ const char* string_val = "2."
+ "7182818284590452353602874713526624977572470936999595749669676277240766303535475945713821785251664274"
+ "2746639193200305992181741359662904357290033429526059563073813232862794349076323382988075319525101901"
+ "1573834187930702154089149934884167509244761460668082264800168477411853742345442437107539077744992069"
+ "5517027618386062613313845830007520449338265602976067371132007093287091274437470472306969772093101416"
+ "9283681902551510865746377211125238978442505695369677078544996996794686445490598793163688923009879312"
+ "7736178215424999229576351482208269895193668033182528869398496465105820939239829488793320362509443117"
+ "3012381970684161403970198376793206832823764648042953118023287825098194558153017567173613320698112509"
+ "9618188159304169035159888851934580727386673858942287922849989208680582574927961048419844436346324496"
+ "8487560233624827041978623209002160990235304369941849146314093431738143640546253152096183690888707016"
+ "7683964243781405927145635490613031072085103837505101157477041718986106873969655212671546889570350354"
+ "0212340784981933432106817012100562788023519303322474501585390473041995777709350366041699732972508869";
+ //
+ // Check if we can just construct from string:
+ //
+ if(digits * 1000L / 301 < 1100)
+ {
+ result = string_val;
+ return;
+ }
+
T lim;
lim = ui_type(1);
eval_ldexp(lim, lim, digits);
@@ -73,6 +120,29 @@
{
typedef typename mpl::front<typename T::unsigned_types>::type ui_type;
typedef typename mpl::front<typename T::float_types>::type real_type;
+ //
+ // 1100 digits in string form:
+ //
+ const char* string_val = "3."
+ "1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679"
+ "8214808651328230664709384460955058223172535940812848111745028410270193852110555964462294895493038196"
+ "4428810975665933446128475648233786783165271201909145648566923460348610454326648213393607260249141273"
+ "7245870066063155881748815209209628292540917153643678925903600113305305488204665213841469519415116094"
+ "3305727036575959195309218611738193261179310511854807446237996274956735188575272489122793818301194912"
+ "9833673362440656643086021394946395224737190702179860943702770539217176293176752384674818467669405132"
+ "0005681271452635608277857713427577896091736371787214684409012249534301465495853710507922796892589235"
+ "4201995611212902196086403441815981362977477130996051870721134999999837297804995105973173281609631859"
+ "5024459455346908302642522308253344685035261931188171010003137838752886587533208381420617177669147303"
+ "5982534904287554687311595628638823537875937519577818577805321712268066130019278766111959092164201989"
+ "3809525720106548586327886593615338182796823030195203530185296899577362259941389124972177528347913152";
+ //
+ // Check if we can just construct from string:
+ //
+ if(digits * 1000L / 301 < 1100)
+ {
+ result = string_val;
+ return;
+ }
T a;
a = ui_type(1);
Modified: sandbox/big_number/libs/multiprecision/test/Jamfile.v2
==============================================================================
--- sandbox/big_number/libs/multiprecision/test/Jamfile.v2 (original)
+++ sandbox/big_number/libs/multiprecision/test/Jamfile.v2 2012-04-04 13:04:00 EDT (Wed, 04 Apr 2012)
@@ -553,6 +553,7 @@
: # input files
: # requirements
<define>TEST_CPP_DEC_FLOAT
+ release # Otherwise runtime is slow
: test_float_io_cpp_dec_float ;
run test_float_io.cpp gmp
@@ -560,6 +561,7 @@
: # input files
: # requirements
<define>TEST_MPF_50
+ release # Otherwise runtime is slow
[ check-target-builds ../config//has_gmp : : <build>no ]
: test_float_io_mpf ;
@@ -568,6 +570,7 @@
: # input files
: # requirements
<define>TEST_MPFR_50
+ release # Otherwise runtime is slow
[ check-target-builds ../config//has_mpfr : : <build>no ]
: test_float_io_mpfr ;
@@ -576,6 +579,7 @@
: # input files
: # requirements
<define>TEST_TOMMATH
+ release # Otherwise runtime is slow
[ check-target-builds ../config//has_tommath : : <build>no ]
: test_int_io_tommath ;
@@ -584,6 +588,7 @@
: # input files
: # requirements
<define>TEST_MPZ
+ release # Otherwise runtime is slow
[ check-target-builds ../config//has_gmp : : <build>no ]
: test_int_io_mpz ;
@@ -592,6 +597,7 @@
: # input files
: # requirements
<define>TEST_CPP_INT
+ release # Otherwise runtime is slow
: test_int_io_cpp_int ;
run test_cpp_int.cpp gmp
@@ -616,6 +622,7 @@
: # requirements
<define>TEST_TOMMATH
[ check-target-builds ../config//has_tommath : : <build>no ]
+ release # Otherwise runtime is slow
: test_rational_io_tommath ;
run test_rational_io.cpp gmp
@@ -624,6 +631,7 @@
: # requirements
<define>TEST_MPQ
[ check-target-builds ../config//has_gmp : : <build>no ]
+ release # Otherwise runtime is slow
: test_rational_io_mpz ;
run test_generic_conv.cpp
@@ -633,7 +641,7 @@
[ check-target-builds ../config//has_gmp : <define>TEST_GMP <source>gmp : ]
[ check-target-builds ../config//has_tommath : <define>TEST_TOMMATH <source>$(TOMMATH) : ]
[ check-target-builds ../config//has_mpfr : <define>TEST_MPFR <source>mpfr : ]
- release
+ release # Otherwise runtime is slow
;
@@ -657,7 +665,7 @@
: # input files
: # requirements
[ check-target-builds ../config//has_tommath : : <build>no ] ;
-
+run ../example/floating_point_examples.cpp ;
if $(enable-specfun)
{