Skip to content

Commit 1aed4e0

Browse files
committed
chore: fix comment for jk
--- type: pre_commit_static_analysis_report description: Results of running static analysis checks when committing changes. report: - task: lint_filenames status: passed - task: lint_editorconfig status: passed - task: lint_markdown status: na - task: lint_package_json status: na - task: lint_repl_help status: na - task: lint_javascript_src status: passed - task: lint_javascript_cli status: na - task: lint_javascript_examples status: na - task: lint_javascript_tests status: na - task: lint_javascript_benchmarks status: na - task: lint_python status: na - task: lint_r status: na - task: lint_c_src status: passed - task: lint_c_examples status: na - task: lint_c_benchmarks status: na - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: na - task: lint_typescript_tests status: na - task: lint_license_headers status: passed ---
1 parent 54ab542 commit 1aed4e0

File tree

2 files changed

+8
-26
lines changed

2 files changed

+8
-26
lines changed

lib/node_modules/@stdlib/math/base/special/rempio2f/lib/kernel_rempio2f.js

Lines changed: 4 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -52,20 +52,11 @@ var zeros = require( '@stdlib/array/base/zeros' );
5252
* \operatorname{ipio2}[i] \cdot 2^{-24(i+1)}
5353
* ```
5454
*
55-
* This table must have at least `(e0-3)/24 + jk` terms. For quad precision (`e0 <= 16360`, `jk = 6`), this is `686`.
55+
* This table must have at least `(e0-3)/24 + jk` terms. For single precision (`e0 <= 127`, `jk = 3`), this is `9`.
5656
*/
5757
var IPIO2 = [
5858
0xA2F983, 0x6E4E44, 0x1529FC, 0x2757D1, 0xF534DD, 0xC0DB62,
59-
0x95993C, 0x439041, 0xFE5163, 0xABDEBB, 0xC561B7, 0x246E3A,
60-
0x424DD2, 0xE00649, 0x2EEA09, 0xD1921C, 0xFE1DEB, 0x1CB129,
61-
0xA73EE8, 0x8235F5, 0x2EBB44, 0x84E99C, 0x7026B4, 0x5F7E41,
62-
0x3991D6, 0x398353, 0x39F49C, 0x845F8B, 0xBDF928, 0x3B1FF8,
63-
0x97FFDE, 0x05980F, 0xEF2F11, 0x8B5A0A, 0x6D1F6D, 0x367ECF,
64-
0x27CB09, 0xB74F46, 0x3F669E, 0x5FEA2D, 0x7527BA, 0xC7EBE5,
65-
0xF17B3D, 0x0739F7, 0x8A5292, 0xEA6BFB, 0x5FB11F, 0x8D5D08,
66-
0x560330, 0x46FC7B, 0x6BABF0, 0xCFBC20, 0x9AF436, 0x1DA9E3,
67-
0x91615E, 0xE61B08, 0x659985, 0x5F14A0, 0x68408D, 0xFFD880,
68-
0x4D7327, 0x310606, 0x1556CA, 0x73A8C9, 0x60E27B, 0xC08C6B
59+
0x95993C, 0x439041, 0xFE5163
6960
];
7061

7162
// Double precision array, obtained by cutting `π/2` into `24` bits chunks...
@@ -364,7 +355,7 @@ function compute( x, y, jz, q, q0, jk, jv, jx, f ) {
364355
* @private
365356
* @param {PositiveNumber} x - input value
366357
* @param {(Array|TypedArray|Object)} y - remainder element
367-
* @param {PositiveInteger} e0 - the exponent of `x[0]` (must be <= 16360)
358+
* @param {PositiveInteger} e0 - the exponent of `x[0]` (must be <= 127)
368359
* @param {PositiveInteger} nx - dimension of `x[]`
369360
* @returns {number} last three binary digits of `N`
370361
*/
@@ -379,7 +370,7 @@ function kernelRempio2f( x, y, e0, nx ) {
379370
var j;
380371
var m;
381372

382-
// Initialize `jk` for double-precision floating-point numbers:
373+
// Initialize `jk` for single-precision floating-point numbers:
383374
jk = 3;
384375

385376
// Determine `jx`, `jv`, `q0` (note that `q0 < 3`):

lib/node_modules/@stdlib/math/base/special/rempio2f/src/main.c

Lines changed: 4 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -60,20 +60,11 @@ static const int32_t MEDIUM = 0x4dc90fdb;
6060
* \operatorname{ipio2}[i] \cdot 2^{-24(i+1)}
6161
* ```
6262
*
63-
* This table must have at least `(e0-3)/24 + jk` terms. For quad precision (`e0 <= 16360`, `jk = 6`), this is `686`.
63+
* This table must have at least `(e0-3)/24 + jk` terms. For single precision (`e0 <= 127`, `jk = 3`), this is `9`.
6464
*/
6565
static const int32_t IPIO2[] = {
6666
0xA2F983, 0x6E4E44, 0x1529FC, 0x2757D1, 0xF534DD, 0xC0DB62,
67-
0x95993C, 0x439041, 0xFE5163, 0xABDEBB, 0xC561B7, 0x246E3A,
68-
0x424DD2, 0xE00649, 0x2EEA09, 0xD1921C, 0xFE1DEB, 0x1CB129,
69-
0xA73EE8, 0x8235F5, 0x2EBB44, 0x84E99C, 0x7026B4, 0x5F7E41,
70-
0x3991D6, 0x398353, 0x39F49C, 0x845F8B, 0xBDF928, 0x3B1FF8,
71-
0x97FFDE, 0x05980F, 0xEF2F11, 0x8B5A0A, 0x6D1F6D, 0x367ECF,
72-
0x27CB09, 0xB74F46, 0x3F669E, 0x5FEA2D, 0x7527BA, 0xC7EBE5,
73-
0xF17B3D, 0x0739F7, 0x8A5292, 0xEA6BFB, 0x5FB11F, 0x8D5D08,
74-
0x560330, 0x46FC7B, 0x6BABF0, 0xCFBC20, 0x9AF436, 0x1DA9E3,
75-
0x91615E, 0xE61B08, 0x659985, 0x5F14A0, 0x68408D, 0xFFD880,
76-
0x4D7327, 0x310606, 0x1556CA, 0x73A8C9, 0x60E27B, 0xC08C6B
67+
0x95993C, 0x439041, 0xFE5163
7768
};
7869

7970
// Double precision array, obtained by cutting `π/2` into `24` bits chunks...
@@ -195,7 +186,7 @@ static const double PIO2[] = {
195186
*
196187
* @param x input value
197188
* @param y output array to store the remainder
198-
* @param e0 the exponent of `x[0]` (must be <= 16360)
189+
* @param e0 the exponent of `x[0]` (must be <= 127)
199190
* @param nx dimension of `x[]`
200191
* @return last 3 binary digits
201192
*/
@@ -220,7 +211,7 @@ int32_t rempio2fKernel( const double* x, double *y, int32_t e0, int32_t nx ) {
220211
double fw;
221212
double z;
222213

223-
// Initialize `jk` for double-precision floating-point numbers:
214+
// Initialize `jk` for single-precision floating-point numbers:
224215
jk = 3;
225216
jp = jk;
226217

0 commit comments

Comments
 (0)