]> www.ginac.de Git - cln.git/commitdiff
Get rid CL_REQUIRE/CL_PROVIDE(cl_F_epspos), it is not really necessary.
authorAlexei Sheplyakov <varg@theor.jinr.ru>
Thu, 21 Aug 2008 12:06:02 +0000 (16:06 +0400)
committerAlexei Sheplyakov <varg@theor.jinr.ru>
Wed, 27 Aug 2008 04:41:07 +0000 (08:41 +0400)
Move static variables [SFDL]F_epsilon into LF_epsilon() function (which
is the only user of those variables) in order to avoid possible static
order initialization problems. While at it, make LF_epsilon() function
static.

include/cln/float.h
src/float/misc/cl_F_epspos.cc

index 5a06ff8f2f8e1eb6c32f909040cbe6a42ae7d3b0..fe8a9c8fd749ad9c1c354bf4a0a1fa3473fd75fe 100644 (file)
@@ -570,7 +570,6 @@ extern const cl_F least_negative_float (float_format_t f);
 
 // Returns the smallest floating point number e > 0 such that 1+e != 1.
 extern const cl_F float_epsilon (float_format_t f);
-//CL_REQUIRE(cl_F_epspos)
 
 // Returns the smallest floating point number e > 0 such that 1-e != 1.
 extern const cl_F float_negative_epsilon (float_format_t f);
index f3ade60722ab5950d6902beb78786f2e4ded1ef0..8342defd087c583ddb6a951159411dc3aef0f692 100644 (file)
@@ -3,8 +3,6 @@
 // General includes.
 #include "cl_sysdep.h"
 
-CL_PROVIDE(cl_F_epspos)
-
 // Specification.
 #include "cln/float.h"
 
@@ -19,23 +17,7 @@ CL_PROVIDE(cl_F_epspos)
 
 namespace cln {
 
-// Bei Floats mit d Bits (incl. Hiddem Bit, also d = ?F_mant_len+1)
-// ist ?F_epsilon = 2^-d*(1+2^(1-d)), d.h. Mantisse 10...01, Vorzeichen +.
-
-static const cl_SF SF_epsilon =
-       make_SF(0,SF_exp_mid-SF_mant_len,bit(SF_mant_len)+1);
-
-static const cl_FF FF_epsilon =
-       encode_FF(0,-FF_mant_len,bit(FF_mant_len)+1);
-
-static const cl_DF DF_epsilon =
-       #if (cl_word_size==64)
-         encode_DF(0,-DF_mant_len,bit(DF_mant_len)+1);
-       #else
-         encode_DF(0,-DF_mant_len,bit(DF_mant_len-32),1);
-       #endif
-
-inline const cl_LF LF_epsilon (uintC len)
+static inline const cl_LF LF_epsilon (uintC len)
 {
        var Lfloat erg = allocate_lfloat(len,LF_exp_mid+1-intDsize*len,0);
        var uintD* ptr = &TheLfloat(erg)->data[0];
@@ -53,6 +35,22 @@ inline const cl_LF LF_epsilon (uintC len)
 
 const cl_F float_epsilon (float_format_t f)
 {
+       // Bei Floats mit d Bits (incl. Hiddem Bit, also d = ?F_mant_len+1)
+       // ist ?F_epsilon = 2^-d*(1+2^(1-d)), d.h. Mantisse 10...01, Vorzeichen +.
+
+       static const cl_SF SF_epsilon =
+               make_SF(0,SF_exp_mid-SF_mant_len,bit(SF_mant_len)+1);
+
+       static const cl_FF FF_epsilon =
+               encode_FF(0,-FF_mant_len,bit(FF_mant_len)+1);
+
+       static const cl_DF DF_epsilon =
+       #if (cl_word_size==64)
+               encode_DF(0,-DF_mant_len,bit(DF_mant_len)+1);
+       #else
+               encode_DF(0,-DF_mant_len,bit(DF_mant_len-32),1);
+       #endif
+
        floatformatcase((uintC)f
        ,       return SF_epsilon;
        ,       return FF_epsilon;
@@ -63,4 +61,3 @@ const cl_F float_epsilon (float_format_t f)
 
 }  // namespace cln
 
-CL_PROVIDE_END(cl_F_epspos)