Mercurial > sdl-ios-xcode
view src/timer/mint/SDL_systimer.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 | e857f02e709e |
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" #ifdef SDL_TIMER_MINT /* * TOS/MiNT timer driver * based on vbl vector * * Patrice Mandin */ #include <stdio.h> #include <sys/time.h> #include <signal.h> #include <unistd.h> #include <string.h> #include <errno.h> #include <mint/cookie.h> #include <mint/sysvars.h> #include <mint/osbind.h> #include <mint/mintbind.h> #include "SDL_timer.h" #include "../SDL_timer_c.h" #include "SDL_thread.h" #include "SDL_vbltimer_s.h" /* from audio/mint */ void SDL_MintAudio_CheckFpu(void); /* The first ticks value of the application */ static Uint32 start; static SDL_bool read_hz200_from_vbl = SDL_FALSE; static int mint_present; /* can we use Syield() ? */ void SDL_StartTicks(void) { void *oldpile; unsigned long dummy; /* Set first ticks value */ oldpile = (void *) Super(0); start = *((volatile long *) _hz_200); Super(oldpile); start *= 5; /* One _hz_200 tic is 5ms */ mint_present = (Getcookie(C_MiNT, &dummy) == C_FOUND); } Uint32 SDL_GetTicks(void) { Uint32 now = start; if (read_hz200_from_vbl) { now = SDL_Atari_hz200; } else { void *old_stack = (void *) Super(0); now = *((volatile long *) _hz_200); Super(old_stack); } return ((now * 5) - start); } void SDL_Delay(Uint32 ms) { Uint32 now; now = SDL_GetTicks(); while ((SDL_GetTicks() - now) < ms) { if (mint_present) { Syield(); } } } /* Data to handle a single periodic alarm */ static SDL_bool timer_installed = SDL_FALSE; /* This is only called if the event thread is not running */ int SDL_SYS_TimerInit(void) { void *old_stack; SDL_MintAudio_CheckFpu(); /* Install RunTimer in vbl vector */ old_stack = (void *) Super(0); timer_installed = !SDL_AtariVblInstall(SDL_ThreadedTimerCheck); Super(old_stack); if (!timer_installed) { return (-1); } read_hz200_from_vbl = SDL_TRUE; return (SDL_SetTimerThreaded(0)); } void SDL_SYS_TimerQuit(void) { if (timer_installed) { /* Uninstall RunTimer vbl vector */ void *old_stack = (void *) Super(0); SDL_AtariVblUninstall(SDL_ThreadedTimerCheck); Super(old_stack); timer_installed = SDL_FALSE; } read_hz200_from_vbl = SDL_FALSE; } int SDL_SYS_StartTimer(void) { SDL_SetError("Internal logic error: MiNT uses vbl timer"); return (-1); } void SDL_SYS_StopTimer(void) { return; } #endif /* SDL_TIMER_MINT */ /* vi: set ts=4 sw=4 expandtab: */