Marc/server
1
0
Fork 0
forked from github/server
Code Issues Pull requests Projects Releases Wiki Activity

Happy Pi day!

Browse source 
Let's use M_PI, it's available in all known math.h implementations.
This commit is contained in:
Enno Rehling 2021-03-14 17:39:53 +01:00
parent 07c1dee907
commit d1aff44d7b
1 changed files with 5 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

12
src/util/rand.c
View file

@ -1,3 +1,6 @@
#ifdef _MSC_VER
#define _USE_MATH_DEFINES
#endif
#include <platform.h> #include <platform.h>
#include "rand.h" #include "rand.h"
#include "mtrand.h" #include "mtrand.h"
@ -10,11 +13,6 @@
#include <float.h> #include <float.h>
#include <ctype.h> #include <ctype.h>
/* do not use M_PI, use one of these instead: */
#define PI_F 3.1415926535897932384626433832795F
#define PI_D 3.1415926535897932384626433832795
#define PI_L 3.1415926535897932384626433832795L
int lovar(double xpct_x2) int lovar(double xpct_x2)
{ {
int n = (int)(xpct_x2 * 500) + 1; int n = (int)(xpct_x2 * 500) + 1;
@ -36,10 +34,10 @@ double normalvariate(double mu, double sigma)
if (phase == 0) { if (phase == 0) {
U = (rng_int() + 1.) / (RNG_RAND_MAX + 2.); U = (rng_int() + 1.) / (RNG_RAND_MAX + 2.);
V = rng_int() / (RNG_RAND_MAX + 1.); V = rng_int() / (RNG_RAND_MAX + 1.);
Z = sqrt(-2 * log(U)) * sin(2 * PI_D * V); Z = sqrt(-2 * log(U)) * sin(2 * M_PI * V);
} }
else { else {
Z = sqrt(-2 * log(U)) * cos(2 * PI_D * V); Z = sqrt(-2 * log(U)) * cos(2 * M_PI * V);
} }
phase = 1 - phase; phase = 1 - phase;