Press n or j to go to the next uncovered block, b, p or k for the previous block.
| 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 | 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 376x 376x 376x 376x 376x 376x 376x 376x 376x 376x 376x 376x 376x 376x 376x 376x 376x 376x 376x 376x 376x 10x 10x 10x 10x 10x 376x 376x 376x 1x 1x 1x 1x 1x | /**
* @license Apache-2.0
*
* Copyright (c) 2025 The Stdlib Authors.
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
* ## Notice
*
* The following copyright and license were part of the original implementation available as part of [FreeBSD]{@link https://svnweb.freebsd.org/base/release/9.3.0/lib/msun/src/k_rem_pio2.c}. The implementation follows the original, but has been modified for JavaScript.
*
* ```text
* Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
*
* Developed at SunPro, a Sun Microsystems, Inc. business.
* Permission to use, copy, modify, and distribute this
* software is freely granted, provided that this notice
* is preserved.
* ```
*/
'use strict';
// MODULES //
var round = require( '@stdlib/math/base/special/round' );
var getWord = require( '@stdlib/number/float32/base/to-word' ); // Single-precision word
// VARIABLES //
// 24 bits of 2/π (single precision):
var INVPIO2 = 0.636619772; // 0x3f22f983
// First 24 bits of π/2:
var PIO2_1 = 1.570796326; // 0x3fc90fda
// PIO2_1T = π/2 - PIO2_1:
var PIO2_1T = 4.960468e-8; // 0x33522168
// Exponent mask (255 => 0xff for single precision):
var EXPONENT_MASK = 0xff|0; // asm type annotation
// MAIN //
/**
* Computes `x - nπ/2 = r` for medium-sized inputs using single precision.
*
* @private
* @param {number} x - input value
* @param {uint32} ix - word of `x` (single precision)
* @param {(Array|TypedArray|Object)} y - remainder element (single value)
* @returns {integer} factor of `π/2`
*/
function rempio2fMedium( x, ix, y ) {
var word;
var n;
var t;
var r;
var w;
var i;
var j;
n = round( x * INVPIO2 );
r = x - ( n * PIO2_1 );
w = n * PIO2_1T;
// First rounding (good to 24 bits for single precision)...
j = (ix >> 23)|0; // Exponent from 32-bit word
y[ 0 ] = r - w;
word = getWord( y[ 0 ] );
i = j - ( (word >> 23) & EXPONENT_MASK );
// Check if a second iteration is needed...
if ( i > 16 ) {
t = r;
w = n * PIO2_1T;
r = t - w;
y[ 0 ] = r;
}
// No y[1] for single precision
return n;
}
// EXPORTS //
module.exports = rempio2fMedium;
|