From: Lowell Johnson <ldj00@sio.midco.net>
To: picca@synchrotron-soleil.fr
Cc: gsl-discuss@sources.redhat.com
Subject: Re: cvs version: make check FAIL
Date: Thu, 25 May 2006 20:11:00 -0000 [thread overview]
Message-ID: <200605241728.02573.ldj00@sio.midco.net> (raw)
In-Reply-To: <20060524084212.78554aa6.picca@synchrotron-soleil.fr>
[-- Attachment #1: Type: text/plain, Size: 829 bytes --]
On Wednesday 24 May 2006 01:42 am, picca@synchrotron-soleil.fr wrote:
> Hello
>
> Instead of computing 2 times fabs(fn), why not storing it in a
> variable ?
>
> result->val = fn;
> - result->err = GSL_DBL_EPSILON*fabs(fn);
> + result->err = 2.0*GSL_DBL_EPSILON;
> + double factor = fabs(fn);
> + if (factor > 1.0)
> + result->err *= factor;
>
> Have a nice day.
Sure, that sounds fine to me. The cost of the extra fabs() is rather
trivial, but avoiding it will save a few clock cycles.
I've attached a patch to be used in place of the earlier patch.
Thanks.
Lowell
--
Lowell D. Johnson
Linux: Bringing stability, security, and freedom to home and business
computing since 1991. www.linux.org
Free and Open Source Software: Of the people, by the people, for the
people.
[-- Attachment #2: mathieu.patch --]
[-- Type: text/x-diff, Size: 3318 bytes --]
Index: mathieu_angfunc.c
===================================================================
RCS file: /cvs/gsl/gsl/specfunc/mathieu_angfunc.c,v
retrieving revision 1.1
diff -u -r1.1 mathieu_angfunc.c
--- mathieu_angfunc.c 18 Apr 2006 17:59:46 -0000 1.1
+++ mathieu_angfunc.c 24 May 2006 22:17:49 -0000
@@ -31,7 +31,7 @@
int gsl_sf_mathieu_c(int order, double qq, double zz, gsl_sf_result *result)
{
int even_odd, ii, status;
- double coeff[NUM_MATHIEU_COEFF], aa, norm, fn;
+ double coeff[NUM_MATHIEU_COEFF], aa, norm, fn, factor;
norm = 0.0;
@@ -54,7 +54,10 @@
fn = cos(order*zz)/norm;
result->val = fn;
- result->err = GSL_DBL_EPSILON*fabs(fn);
+ result->err = 2.0*GSL_DBL_EPSILON;
+ factor = fabs(fn);
+ if (factor > 1.0)
+ result->err *= factor;
return GSL_SUCCESS;
}
@@ -97,7 +100,10 @@
fn /= norm;
result->val = fn;
- result->err = GSL_DBL_EPSILON*fabs(fn);
+ result->err = 2.0*GSL_DBL_EPSILON;
+ factor = fabs(fn);
+ if (factor > 1.0)
+ result->err *= factor;
return GSL_SUCCESS;
}
@@ -106,7 +112,7 @@
int gsl_sf_mathieu_s(int order, double qq, double zz,gsl_sf_result *result)
{
int even_odd, ii, status;
- double coeff[NUM_MATHIEU_COEFF], aa, norm, fn;
+ double coeff[NUM_MATHIEU_COEFF], aa, norm, fn, factor;
norm = 0.0;
@@ -126,7 +132,10 @@
fn = sin(order*zz);
result->val = fn;
- result->err = GSL_DBL_EPSILON*fabs(fn);
+ result->err = 2.0*GSL_DBL_EPSILON;
+ factor = fabs(fn);
+ if (factor > 1.0)
+ result->err *= factor;
return GSL_SUCCESS;
}
@@ -167,7 +176,10 @@
fn /= norm;
result->val = fn;
- result->err = GSL_DBL_EPSILON*fabs(fn);
+ result->err = 2.0*GSL_DBL_EPSILON;
+ factor = fabs(fn);
+ if (factor > 1.0)
+ result->err *= factor;
return GSL_SUCCESS;
}
Index: mathieu_radfunc.c
===================================================================
RCS file: /cvs/gsl/gsl/specfunc/mathieu_radfunc.c,v
retrieving revision 1.1
diff -u -r1.1 mathieu_radfunc.c
--- mathieu_radfunc.c 18 Apr 2006 17:59:46 -0000 1.1
+++ mathieu_radfunc.c 24 May 2006 22:17:49 -0000
@@ -34,7 +34,7 @@
gsl_sf_result *result)
{
int even_odd, kk, mm, status;
- double maxerr = 1e-14, amax, pi = acos(-1.0), fn;
+ double maxerr = 1e-14, amax, pi = acos(-1.0), fn, factor;
double coeff[NUM_MATHIEU_COEFF], aa, fc, fj, fjp;
double j1c, j2c, j1pc, j2pc;
double u1, u2;
@@ -141,7 +141,10 @@
}
result->val = fn;
- result->err = GSL_DBL_EPSILON*fabs(fn);
+ result->err = 2.0*GSL_DBL_EPSILON;
+ factor = fabs(fn);
+ if (factor > 1.0)
+ result->err *= factor;
return GSL_SUCCESS;
}
@@ -151,7 +154,7 @@
gsl_sf_result *result)
{
int even_odd, kk, mm, status;
- double maxerr = 1e-14, amax, pi = acos(-1.0), fn;
+ double maxerr = 1e-14, amax, pi = acos(-1.0), fn, factor;
double coeff[NUM_MATHIEU_COEFF], aa, fc, fj, fjp, fjm;
double j1c, j2c, j1mc, j2mc, j1pc, j2pc;
double u1, u2;
@@ -268,7 +271,10 @@
}
result->val = fn;
- result->err = GSL_DBL_EPSILON*fabs(fn);
+ result->err = 2.0*GSL_DBL_EPSILON;
+ factor = fabs(fn);
+ if (factor > 1.0)
+ result->err *= factor;
return GSL_SUCCESS;
}
next prev parent reply other threads:[~2006-05-24 22:28 UTC|newest]
Thread overview: 9+ messages / expand[flat|nested] mbox.gz Atom feed top
2006-05-23 22:03 Giulio Bottazzi
2006-05-24 2:09 ` Lowell Johnson
2006-05-24 6:42 ` Lowell Johnson
2006-05-24 14:28 ` picca
2006-05-25 20:11 ` Lowell Johnson [this message]
2006-05-25 20:12 ` Gerard Jungman
2006-05-25 20:18 ` Gerard Jungman
2006-05-29 16:10 ` Giulio Bottazzi
2006-06-01 12:33 ` Brian Gough
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=200605241728.02573.ldj00@sio.midco.net \
--to=ldj00@sio.midco.net \
--cc=gsl-discuss@sources.redhat.com \
--cc=picca@synchrotron-soleil.fr \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).