Mercurial > sdl-ios-xcode
view src/SDL_fatal.c @ 2752:edd2839b36f7
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.
author | Sam Lantinga <slouken@libsdl.org> |
---|---|
date | Mon, 15 Sep 2008 04:31:30 +0000 |
parents | c121d94672cb |
children | 99210400e8b9 |
line wrap: on
line source
/* SDL - Simple DirectMedia Layer Copyright (C) 1997-2006 Sam Lantinga This library is free software; you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation; either version 2.1 of the License, or (at your option) any later version. This library is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more details. You should have received a copy of the GNU Lesser General Public License along with this library; if not, write to the Free Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA Sam Lantinga slouken@libsdl.org */ #include "SDL_config.h" /* General fatal signal handling code for SDL */ #ifdef HAVE_SIGNAL_H #include <signal.h> #include "SDL.h" #include "SDL_fatal.h" /* This installs some signal handlers for the more common fatal signals, so that if the programmer is lazy, the app doesn't die so horribly if the program crashes. */ static void SDL_Parachute(int sig) { signal(sig, SIG_DFL); SDL_Quit(); raise(sig); } static int SDL_fatal_signals[] = { SIGSEGV, #ifdef SIGBUS SIGBUS, #endif #ifdef SIGFPE SIGFPE, #endif #ifdef SIGQUIT SIGQUIT, #endif 0 }; void SDL_InstallParachute(void) { /* Set a handler for any fatal signal not already handled */ int i; #ifdef HAVE_SIGACTION struct sigaction action; for (i = 0; SDL_fatal_signals[i]; ++i) { sigaction(SDL_fatal_signals[i], NULL, &action); if (action.sa_handler == SIG_DFL) { action.sa_handler = SDL_Parachute; sigaction(SDL_fatal_signals[i], &action, NULL); } } #ifdef SIGALRM /* Set SIGALRM to be ignored -- necessary on Solaris */ sigaction(SIGALRM, NULL, &action); if (action.sa_handler == SIG_DFL) { action.sa_handler = SIG_IGN; sigaction(SIGALRM, &action, NULL); } #endif #else void (*ohandler) (int); for (i = 0; SDL_fatal_signals[i]; ++i) { ohandler = signal(SDL_fatal_signals[i], SDL_Parachute); if (ohandler != SIG_DFL) { signal(SDL_fatal_signals[i], ohandler); } } #endif /* HAVE_SIGACTION */ return; } void SDL_UninstallParachute(void) { /* Remove a handler for any fatal signal handled */ int i; #ifdef HAVE_SIGACTION struct sigaction action; for (i = 0; SDL_fatal_signals[i]; ++i) { sigaction(SDL_fatal_signals[i], NULL, &action); if (action.sa_handler == SDL_Parachute) { action.sa_handler = SIG_DFL; sigaction(SDL_fatal_signals[i], &action, NULL); } } #else void (*ohandler) (int); for (i = 0; SDL_fatal_signals[i]; ++i) { ohandler = signal(SDL_fatal_signals[i], SIG_DFL); if (ohandler != SDL_Parachute) { signal(SDL_fatal_signals[i], ohandler); } } #endif /* HAVE_SIGACTION */ } #else /* No signals on this platform, nothing to do.. */ void SDL_InstallParachute(void) { return; } void SDL_UninstallParachute(void) { return; } #endif /* HAVE_SIGNAL_H */ /* vi: set ts=4 sw=4 expandtab: */