Commit d4783c9a authored by Morten Welinder's avatar Morten Welinder

House keeping.

parent d5fb139f
2013-12-19 Morten Welinder <terra@gnome.org>
* src/sf-gamma.c (qfactf): Squeeze a few extra bits out of this,
especially when |x|<1.
especially when |x|<<1. Ditto for x<-1.
2013-12-17 Morten Welinder <terra@gnome.org>
......
......@@ -119,6 +119,7 @@ typedef long double gnm_float;
#define gnm_quad_mul12 go_quad_mul12l
#define gnm_quad_one go_quad_onel
#define gnm_quad_pi go_quad_pil
#define gnm_quad_2pi go_quad_2pil
#define gnm_quad_pow go_quad_powl
#define gnm_quad_sin go_quad_sinl
#define gnm_quad_sinpi go_quad_sinpil
......@@ -227,6 +228,7 @@ typedef double gnm_float;
#define gnm_quad_mul12 go_quad_mul12
#define gnm_quad_one go_quad_one
#define gnm_quad_pi go_quad_pi
#define gnm_quad_2pi go_quad_2pi
#define gnm_quad_pow go_quad_pow
#define gnm_quad_sin go_quad_sin
#define gnm_quad_sinpi go_quad_sinpi
......
......@@ -920,8 +920,7 @@ qfactf (gnm_float x, GnmQuad *mant, int *exp2)
*exp2 = 0;
/* sqrt(2Pi) */
gnm_quad_add (&f1, &gnm_quad_pi, &gnm_quad_pi);
gnm_quad_sqrt (&f1, &f1);
gnm_quad_sqrt (&f1, &gnm_quad_2pi);
if (debug) g_printerr ("f1=%.20g\n", gnm_quad_value (&f1));
/* (y/e)^y */
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment