Browse Source

* src/float/transcendental/cl_F_lnx.cc: Make actuallen of type uintC.

* src/float/transcendental/cl_F_expx.cc: Likewise.
        * src/float/transcendental/cl_F_sinhx.cc: Likewise.
        * src/float/transcendental/cl_F_sinx.cc: Likewise.
master
Richard Kreckel 17 years ago
parent
commit
3c77def744
  1. 7
      ChangeLog
  2. 2
      src/float/transcendental/cl_F_expx.cc
  3. 2
      src/float/transcendental/cl_F_lnx.cc
  4. 2
      src/float/transcendental/cl_F_sinhx.cc
  5. 2
      src/float/transcendental/cl_F_sinx.cc

7
ChangeLog

@ -1,3 +1,10 @@
2007-12-18 Richard B. Kreckel <kreckel@ginac.de>
* src/float/transcendental/cl_F_lnx.cc: Make actuallen of type uintC.
* src/float/transcendental/cl_F_expx.cc: Likewise.
* src/float/transcendental/cl_F_sinhx.cc: Likewise.
* src/float/transcendental/cl_F_sinx.cc: Likewise.
2007-12-17 Richard B. Kreckel <kreckel@ginac.de>
Silly workaround for silly bug in gmp.h:

2
src/float/transcendental/cl_F_expx.cc

@ -46,7 +46,7 @@ const cl_LF expx_naive (const cl_LF& x)
// wie oben, mit adaptiver Genauigkeit während der Potenzreihen-Summation.
if (zerop(x))
return cl_float(1,x);
var uintL actuallen = TheLfloat(x)->len;
var uintC actuallen = TheLfloat(x)->len;
var uintC d = float_digits(x);
var sintE e = float_exponent(x);
if (e < -(sintC)d) // e < -d ?

2
src/float/transcendental/cl_F_lnx.cc

@ -47,7 +47,7 @@ const cl_LF lnx_naive (const cl_LF& x)
var cl_LF y = x-cl_float(1,x);
if (zerop(y)) // y=0.0 -> y als Ergebnis
return y;
var uintL actuallen = TheLfloat(x)->len;
var uintC actuallen = TheLfloat(x)->len;
var uintC d = float_digits(x);
var sintE e = float_exponent(y);
if (e <= -(sintC)d) // e <= -d ?

2
src/float/transcendental/cl_F_sinhx.cc

@ -114,7 +114,7 @@ const cl_LF sinhx_naive (const cl_LF& x)
if (zerop(x))
return x;
var uintL actuallen = TheLfloat(x)->len;
var uintC actuallen = TheLfloat(x)->len;
var uintC d = float_digits(x);
var sintE e = float_exponent(x);
if (e <= (1-(sintC)d)>>1) // e <= (1-d)/2 <==> e <= -ceiling((d-1)/2) ?

2
src/float/transcendental/cl_F_sinx.cc

@ -123,7 +123,7 @@ const cl_LF sinx_naive (const cl_LF& x)
if (zerop(x))
return x;
var uintL actuallen = TheLfloat(x)->len;
var uintC actuallen = TheLfloat(x)->len;
var uintC d = float_digits(x);
var sintE e = float_exponent(x);
if (e <= (-(sintC)d)>>1) // e <= (-d)/2 <==> e <= -ceiling(d/2) ?

Loading…
Cancel
Save