Исходный код языка R - что такое `NOMORE_FOR_THREADS` в коде C
Я шел через исходный код R для функции gamma()
который указывает на этот исходный файл C.
Фрагмент кода выглядит следующим образом:
#ifdef NOMORE_FOR_THREADS
static int ngam = 0;
static double xmin = 0, xmax = 0., xsml = 0., dxrel = 0.;
/* Initialize machine dependent constants, the first time gamma() is called.
FIXME for threads ! */
if (ngam == 0) {
ngam = chebyshev_init(gamcs, 42, DBL_EPSILON/20);/*was .1*d1mach(3)*/
gammalims(&xmin, &xmax);/*-> ./gammalims.c */
xsml = exp(fmax2(log(DBL_MIN), -log(DBL_MAX)) + 0.01);
/* = exp(.01)*DBL_MIN = 2.247e-308 for IEEE */
dxrel = sqrt(DBL_EPSILON);/*was sqrt(d1mach(4)) */
}
#else
/* For IEEE double precision DBL_EPSILON = 2^-52 = 2.220446049250313e-16 :
* (xmin, xmax) are non-trivial, see ./gammalims.c
* xsml = exp(.01)*DBL_MIN
* dxrel = sqrt(DBL_EPSILON) = 2 ^ -26
*/
# define ngam 22
# define xmin -170.5674972726612
# define xmax 171.61447887182298
# define xsml 2.2474362225598545e-308
# define dxrel 1.490116119384765696e-8
#endif
Теперь я понимаю, что этот кусок кода инициализирует значения переменных, таких как ngam
, xsml
, dxrel
и т.д. Но я не понимаю, что такое NOMORE_FOR_THREADS
,
Я не мог найти, где это определено, поэтому я не понимаю, когда он войдет в #ifdef
пункт или же просто пойти в #else
предложение и определение переменных в соответствии со стандартом IEEE.
Если кто-то может пролить свет на то, что NOMORE_FOR_THREADS
или где это определено, это было бы действительно полезно.
Спасибо за помощь.
1 ответ
Выглядит так #ifdef
был добавлен в ревизию 13433:
> svn log -vr13433
------------------------------------------------------------------------
r13433 | maechler | 2001-04-10 10:05:54 -0500 (Tue, 10 Apr 2001) | 2 lines
Changed paths:
M /trunk/src/nmath/beta.c
M /trunk/src/nmath/gamma.c
M /trunk/src/nmath/lgamma.c
M /trunk/src/nmath/lgammacor.c
globals now #ifdef NOMORE_FOR_THREADS
------------------------------------------------------------------------
А вот разница с предыдущей ревизией:
> svn diff -r13432:13433 src/nmath/gamma.c
Index: src/nmath/gamma.c
===================================================================
--- src/nmath/gamma.c (revision 13432)
+++ src/nmath/gamma.c (revision 13433)
@@ -1,7 +1,7 @@
/*
* Mathlib : A C Library of Special Functions
* Copyright (C) 1998 Ross Ihaka
- * Copyright (C) 2000 The R Development Core Team
+ * Copyright (C) 2000-2001 The R Development Core Team
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
@@ -40,7 +40,7 @@
double gammafn(double x)
{
- static /* const */ double gamcs[42] = {
+ const double gamcs[42] = {
+.8571195590989331421920062399942e-2,
+.4415381324841006757191315771652e-2,
+.5685043681599363378632664588789e-1,
@@ -89,6 +89,7 @@
double y;
double sinpiy, value;
+#ifdef NOMORE_FOR_THREADS
static int ngam = 0;
static double xmin = 0, xmax = 0., xsml = 0., dxrel = 0.;
@@ -101,6 +102,18 @@
/* = exp(.01)*DBL_MIN = 2.247e-308 for IEEE */
dxrel = sqrt(1/DBL_EPSILON);/*was (1/d1mach(4)) */
}
+#else
+/* For IEEE double precision DBL_EPSILON = 2^-52 = 2.220446049250313e-16 :
+ * (xmin, xmax) are non-trivial, see ./gammalims.c
+ * xsml = exp(.01)*DBL_MIN
+ * dxrel = sqrt(1/DBL_EPSILON) = 2 ^ 26
+*/
+# define ngam 22
+# define xmin -170.5674972726612
+# define xmax 171.61447887182298
+# define xsml 2.2474362225598545e-308
+# define dxrel 67108864.
+#endif
if(ISNAN(x)) return x;
Я не видела NOMORE_FOR_THREADS
определяется в любом месте исходного кода в этой ревизии, следующих двух сотен ревизий после его добавления или в текущих ревизиях. Возможно (хотя кажется маловероятным), что это может быть определено в некотором системном заголовочном файле. Кажется более вероятным, что это будет указано через аргумент командной строки.