http://sources.redhat.com/ml/newlib/2002/msg00230.html Stephen L Moshier wrote: > > pow(x,y) returns 0 when x is very close to -1.0 and y is very large. > The following test program prints > > pow(1.0000000000000002e+00 4.5035996273704970e+15) = 2.7182818284590455e+00 > pow(-1.0000000000000002e+00 4.5035996273704970e+15) =0.0000000000000000e+00 > pow(9.9999999999999978e-01 4.5035996273704970e+15) = 3.6787944117144222e-01 > pow(-9.9999999999999978e-01 4.5035996273704970e+15) = 0.0000000000000000e+00 > > which is incorrect for the negative arguments raised to an odd integer > power. > > ----- > double pow (double, double); > > int > main () > { > double x, y, z; > > x = 1.0 + pow (2.0, -52.0); > y = 1.0 + pow (2.0, 52.0); > z = pow (x, y); > printf ("pow(%.16e %.16e) = %.16e\n", x, y, z); > x = -x; > z = pow (x, y); > printf ("pow(%.16e %.16e) = %.16e\n", x, y, z); > x = 1.0 - pow (2.0, -52.0); > z = pow (x, y); > printf ("pow(%.16e %.16e) = %.16e\n", x, y, z); > x = -x; > z = pow (x, y); > printf ("pow(%.16e %.16e) = %.16e\n", x, y, z); > } > ----- > > Here is a patch for newlib/libm/math/epow.c: Patch checked in and duplicated for ef_pow.c. Thanks. -- Jeff J.

<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE plist PUBLIC "-//Apple Computer//DTD PLIST 1.0//EN" "http://www.apple.com/DTDs/PropertyList-1.0.dtd">
<plist version="1.0">
<dict>
	<key>CFBundleDevelopmentRegion</key>
	<string>English</string>
	<key>CFBundleExecutable</key>
	<string>SDL</string>
	<key>CFBundleGetInfoString</key>
	<string>http://www.libsdl.org</string>
	<key>CFBundleIconFile</key>
	<string></string>
	<key>CFBundleIdentifier</key>
	<string>SDL</string>
	<key>CFBundleInfoDictionaryVersion</key>
	<string>6.0</string>
	<key>CFBundleName</key>
	<string>Simple DirectMedia Layer</string>
	<key>CFBundlePackageType</key>
	<string>FMWK</string>
	<key>CFBundleShortVersionString</key>
	<string>1.2.9</string>
	<key>CFBundleSignature</key>
	<string>SDLX</string>
	<key>CFBundleVersion</key>
	<string>1.2.10</string>
</dict>
</plist>