X-Authentication-Warning: delorie.com: mail set sender to djgpp-workers-bounces using -f From: Kbwms AT aol DOT com Message-ID: <59.17546f3.2d109998@aol.com> Date: Tue, 16 Dec 2003 12:23:36 EST Subject: e_pow.c To: djgpp-workers AT delorie DOT com MIME-Version: 1.0 Content-Type: multipart/alternative; boundary="part1_59.17546f3.2d109998_boundary" X-Mailer: 8.0 for Windows sub 6021 Reply-To: djgpp-workers AT delorie DOT com --part1_59.17546f3.2d109998_boundary Content-Type: text/plain; charset="US-ASCII" Content-Transfer-Encoding: 7bit Kudos to the person who put a note in a patch for e_pow.c concerning cases where, for z=pow(x,y), when x is very near 1 and y is very large, z would vanish. It prompted me to create test cases for all three floating point precisions. The same problem existed for the float version and was easily fixed. For testing the double precision case, the parameters are: x = 1-DBL_EPSILON y = log(DBL_MIN)/log(1-DBL_EPSILON) KB Williams --part1_59.17546f3.2d109998_boundary Content-Type: text/html; charset="US-ASCII" Content-Transfer-Encoding: quoted-printable Kudos to the person who put a note in= a patch for e_pow.c concerning cases where, for z=3Dpow(x,y), when x is ver= y near 1 and y is very large, z would vanish.  It prompted me to create= test cases for all three floating point precisions.  The same problem=20= existed for the float version and was easily fixed.

For testing the double precision case, the parameters are:

x =3D 1-DBL_EPSILON
y =3D log(DBL_MIN)/log(1-DBL_EPSILON)


KB Williams --part1_59.17546f3.2d109998_boundary--