1
+
2
+ <!doctype html>
3
+ < html lang ="en ">
4
+
5
+ < head >
6
+ < title > Code coverage report for math/base/special/gcdf/lib/binary_gcd.js</ title >
7
+ < meta charset ="utf-8 " />
8
+ < link rel ="stylesheet " href ="../../../../../prettify.css " />
9
+ < link rel ="stylesheet " href ="../../../../../base.css " />
10
+ < link rel ="shortcut icon " type ="image/x-icon " href ="../../../../../favicon.png " />
11
+ < meta name ="viewport " content ="width=device-width, initial-scale=1 " />
12
+ < style type ='text/css '>
13
+ .coverage-summary .sorter {
14
+ background-image : url(../ ../ ../ ../ ../ sor t-arrow-sprite.png);
15
+ }
16
+ </ style >
17
+ </ head >
18
+
19
+ < body >
20
+ < div class ='wrapper '>
21
+ < div class ='pad1 '>
22
+ < h1 > < a href ="../../../../../index.html "> All files</ a > / < a href ="index.html "> math/base/special/gcdf/lib</ a > binary_gcd.js</ h1 >
23
+ < div class ='clearfix '>
24
+
25
+ < div class ='fl pad1y space-right2 '>
26
+ < span class ="strong "> 100% </ span >
27
+ < span class ="quiet "> Statements</ span >
28
+ < span class ='fraction '> 81/81</ span >
29
+ </ div >
30
+
31
+
32
+ < div class ='fl pad1y space-right2 '>
33
+ < span class ="strong "> 100% </ span >
34
+ < span class ="quiet "> Branches</ span >
35
+ < span class ='fraction '> 15/15</ span >
36
+ </ div >
37
+
38
+
39
+ < div class ='fl pad1y space-right2 '>
40
+ < span class ="strong "> 100% </ span >
41
+ < span class ="quiet "> Functions</ span >
42
+ < span class ='fraction '> 1/1</ span >
43
+ </ div >
44
+
45
+
46
+ < div class ='fl pad1y space-right2 '>
47
+ < span class ="strong "> 100% </ span >
48
+ < span class ="quiet "> Lines</ span >
49
+ < span class ='fraction '> 81/81</ span >
50
+ </ div >
51
+
52
+
53
+ </ div >
54
+ < p class ="quiet ">
55
+ Press < em > n</ em > or < em > j</ em > to go to the next uncovered block, < em > b</ em > , < em > p</ em > or < em > k</ em > for the previous block.
56
+ </ p >
57
+ < template id ="filterTemplate ">
58
+ < div class ="quiet ">
59
+ Filter:
60
+ < input type ="search " id ="fileSearch ">
61
+ </ div >
62
+ </ template >
63
+ </ div >
64
+ < div class ='status-line high '> </ div >
65
+ < pre > < table class ="coverage ">
66
+ < tr > < td class ="line-count quiet "> < a name ='L1 '> </ a > < a href ='#L1 '> 1</ a >
67
+ < a name ='L2 '> </ a > < a href ='#L2 '> 2</ a >
68
+ < a name ='L3 '> </ a > < a href ='#L3 '> 3</ a >
69
+ < a name ='L4 '> </ a > < a href ='#L4 '> 4</ a >
70
+ < a name ='L5 '> </ a > < a href ='#L5 '> 5</ a >
71
+ < a name ='L6 '> </ a > < a href ='#L6 '> 6</ a >
72
+ < a name ='L7 '> </ a > < a href ='#L7 '> 7</ a >
73
+ < a name ='L8 '> </ a > < a href ='#L8 '> 8</ a >
74
+ < a name ='L9 '> </ a > < a href ='#L9 '> 9</ a >
75
+ < a name ='L10 '> </ a > < a href ='#L10 '> 10</ a >
76
+ < a name ='L11 '> </ a > < a href ='#L11 '> 11</ a >
77
+ < a name ='L12 '> </ a > < a href ='#L12 '> 12</ a >
78
+ < a name ='L13 '> </ a > < a href ='#L13 '> 13</ a >
79
+ < a name ='L14 '> </ a > < a href ='#L14 '> 14</ a >
80
+ < a name ='L15 '> </ a > < a href ='#L15 '> 15</ a >
81
+ < a name ='L16 '> </ a > < a href ='#L16 '> 16</ a >
82
+ < a name ='L17 '> </ a > < a href ='#L17 '> 17</ a >
83
+ < a name ='L18 '> </ a > < a href ='#L18 '> 18</ a >
84
+ < a name ='L19 '> </ a > < a href ='#L19 '> 19</ a >
85
+ < a name ='L20 '> </ a > < a href ='#L20 '> 20</ a >
86
+ < a name ='L21 '> </ a > < a href ='#L21 '> 21</ a >
87
+ < a name ='L22 '> </ a > < a href ='#L22 '> 22</ a >
88
+ < a name ='L23 '> </ a > < a href ='#L23 '> 23</ a >
89
+ < a name ='L24 '> </ a > < a href ='#L24 '> 24</ a >
90
+ < a name ='L25 '> </ a > < a href ='#L25 '> 25</ a >
91
+ < a name ='L26 '> </ a > < a href ='#L26 '> 26</ a >
92
+ < a name ='L27 '> </ a > < a href ='#L27 '> 27</ a >
93
+ < a name ='L28 '> </ a > < a href ='#L28 '> 28</ a >
94
+ < a name ='L29 '> </ a > < a href ='#L29 '> 29</ a >
95
+ < a name ='L30 '> </ a > < a href ='#L30 '> 30</ a >
96
+ < a name ='L31 '> </ a > < a href ='#L31 '> 31</ a >
97
+ < a name ='L32 '> </ a > < a href ='#L32 '> 32</ a >
98
+ < a name ='L33 '> </ a > < a href ='#L33 '> 33</ a >
99
+ < a name ='L34 '> </ a > < a href ='#L34 '> 34</ a >
100
+ < a name ='L35 '> </ a > < a href ='#L35 '> 35</ a >
101
+ < a name ='L36 '> </ a > < a href ='#L36 '> 36</ a >
102
+ < a name ='L37 '> </ a > < a href ='#L37 '> 37</ a >
103
+ < a name ='L38 '> </ a > < a href ='#L38 '> 38</ a >
104
+ < a name ='L39 '> </ a > < a href ='#L39 '> 39</ a >
105
+ < a name ='L40 '> </ a > < a href ='#L40 '> 40</ a >
106
+ < a name ='L41 '> </ a > < a href ='#L41 '> 41</ a >
107
+ < a name ='L42 '> </ a > < a href ='#L42 '> 42</ a >
108
+ < a name ='L43 '> </ a > < a href ='#L43 '> 43</ a >
109
+ < a name ='L44 '> </ a > < a href ='#L44 '> 44</ a >
110
+ < a name ='L45 '> </ a > < a href ='#L45 '> 45</ a >
111
+ < a name ='L46 '> </ a > < a href ='#L46 '> 46</ a >
112
+ < a name ='L47 '> </ a > < a href ='#L47 '> 47</ a >
113
+ < a name ='L48 '> </ a > < a href ='#L48 '> 48</ a >
114
+ < a name ='L49 '> </ a > < a href ='#L49 '> 49</ a >
115
+ < a name ='L50 '> </ a > < a href ='#L50 '> 50</ a >
116
+ < a name ='L51 '> </ a > < a href ='#L51 '> 51</ a >
117
+ < a name ='L52 '> </ a > < a href ='#L52 '> 52</ a >
118
+ < a name ='L53 '> </ a > < a href ='#L53 '> 53</ a >
119
+ < a name ='L54 '> </ a > < a href ='#L54 '> 54</ a >
120
+ < a name ='L55 '> </ a > < a href ='#L55 '> 55</ a >
121
+ < a name ='L56 '> </ a > < a href ='#L56 '> 56</ a >
122
+ < a name ='L57 '> </ a > < a href ='#L57 '> 57</ a >
123
+ < a name ='L58 '> </ a > < a href ='#L58 '> 58</ a >
124
+ < a name ='L59 '> </ a > < a href ='#L59 '> 59</ a >
125
+ < a name ='L60 '> </ a > < a href ='#L60 '> 60</ a >
126
+ < a name ='L61 '> </ a > < a href ='#L61 '> 61</ a >
127
+ < a name ='L62 '> </ a > < a href ='#L62 '> 62</ a >
128
+ < a name ='L63 '> </ a > < a href ='#L63 '> 63</ a >
129
+ < a name ='L64 '> </ a > < a href ='#L64 '> 64</ a >
130
+ < a name ='L65 '> </ a > < a href ='#L65 '> 65</ a >
131
+ < a name ='L66 '> </ a > < a href ='#L66 '> 66</ a >
132
+ < a name ='L67 '> </ a > < a href ='#L67 '> 67</ a >
133
+ < a name ='L68 '> </ a > < a href ='#L68 '> 68</ a >
134
+ < a name ='L69 '> </ a > < a href ='#L69 '> 69</ a >
135
+ < a name ='L70 '> </ a > < a href ='#L70 '> 70</ a >
136
+ < a name ='L71 '> </ a > < a href ='#L71 '> 71</ a >
137
+ < a name ='L72 '> </ a > < a href ='#L72 '> 72</ a >
138
+ < a name ='L73 '> </ a > < a href ='#L73 '> 73</ a >
139
+ < a name ='L74 '> </ a > < a href ='#L74 '> 74</ a >
140
+ < a name ='L75 '> </ a > < a href ='#L75 '> 75</ a >
141
+ < a name ='L76 '> </ a > < a href ='#L76 '> 76</ a >
142
+ < a name ='L77 '> </ a > < a href ='#L77 '> 77</ a >
143
+ < a name ='L78 '> </ a > < a href ='#L78 '> 78</ a >
144
+ < a name ='L79 '> </ a > < a href ='#L79 '> 79</ a >
145
+ < a name ='L80 '> </ a > < a href ='#L80 '> 80</ a >
146
+ < a name ='L81 '> </ a > < a href ='#L81 '> 81</ a >
147
+ < a name ='L82 '> </ a > < a href ='#L82 '> 82</ a > </ td > < td class ="line-coverage quiet "> < span class ="cline-any cline-yes "> 1x</ span >
148
+ < span class ="cline-any cline-yes "> 1x</ span >
149
+ < span class ="cline-any cline-yes "> 1x</ span >
150
+ < span class ="cline-any cline-yes "> 1x</ span >
151
+ < span class ="cline-any cline-yes "> 1x</ span >
152
+ < span class ="cline-any cline-yes "> 1x</ span >
153
+ < span class ="cline-any cline-yes "> 1x</ span >
154
+ < span class ="cline-any cline-yes "> 1x</ span >
155
+ < span class ="cline-any cline-yes "> 1x</ span >
156
+ < span class ="cline-any cline-yes "> 1x</ span >
157
+ < span class ="cline-any cline-yes "> 1x</ span >
158
+ < span class ="cline-any cline-yes "> 1x</ span >
159
+ < span class ="cline-any cline-yes "> 1x</ span >
160
+ < span class ="cline-any cline-yes "> 1x</ span >
161
+ < span class ="cline-any cline-yes "> 1x</ span >
162
+ < span class ="cline-any cline-yes "> 1x</ span >
163
+ < span class ="cline-any cline-yes "> 1x</ span >
164
+ < span class ="cline-any cline-yes "> 1x</ span >
165
+ < span class ="cline-any cline-yes "> 1x</ span >
166
+ < span class ="cline-any cline-yes "> 1x</ span >
167
+ < span class ="cline-any cline-yes "> 1x</ span >
168
+ < span class ="cline-any cline-yes "> 1x</ span >
169
+ < span class ="cline-any cline-yes "> 1x</ span >
170
+ < span class ="cline-any cline-yes "> 1x</ span >
171
+ < span class ="cline-any cline-yes "> 1x</ span >
172
+ < span class ="cline-any cline-yes "> 1x</ span >
173
+ < span class ="cline-any cline-yes "> 1x</ span >
174
+ < span class ="cline-any cline-yes "> 1x</ span >
175
+ < span class ="cline-any cline-yes "> 1x</ span >
176
+ < span class ="cline-any cline-yes "> 1x</ span >
177
+ < span class ="cline-any cline-yes "> 1x</ span >
178
+ < span class ="cline-any cline-yes "> 1x</ span >
179
+ < span class ="cline-any cline-yes "> 1x</ span >
180
+ < span class ="cline-any cline-yes "> 1x</ span >
181
+ < span class ="cline-any cline-yes "> 1x</ span >
182
+ < span class ="cline-any cline-yes "> 1x</ span >
183
+ < span class ="cline-any cline-yes "> 1x</ span >
184
+ < span class ="cline-any cline-yes "> 1x</ span >
185
+ < span class ="cline-any cline-yes "> 6x</ span >
186
+ < span class ="cline-any cline-yes "> 6x</ span >
187
+ < span class ="cline-any cline-yes "> 6x</ span >
188
+ < span class ="cline-any cline-yes "> 6x</ span >
189
+ < span class ="cline-any cline-yes "> 6x</ span >
190
+ < span class ="cline-any cline-yes "> 6x</ span >
191
+ < span class ="cline-any cline-yes "> 1x</ span >
192
+ < span class ="cline-any cline-yes "> 1x</ span >
193
+ < span class ="cline-any cline-yes "> 6x</ span >
194
+ < span class ="cline-any cline-yes "> 1x</ span >
195
+ < span class ="cline-any cline-yes "> 1x</ span >
196
+ < span class ="cline-any cline-yes "> 4x</ span >
197
+ < span class ="cline-any cline-yes "> 6x</ span >
198
+ < span class ="cline-any cline-yes "> 59x</ span >
199
+ < span class ="cline-any cline-yes "> 59x</ span >
200
+ < span class ="cline-any cline-yes "> 59x</ span >
201
+ < span class ="cline-any cline-yes "> 59x</ span >
202
+ < span class ="cline-any cline-yes "> 4x</ span >
203
+ < span class ="cline-any cline-yes "> 6x</ span >
204
+ < span class ="cline-any cline-yes "> 244x</ span >
205
+ < span class ="cline-any cline-yes "> 244x</ span >
206
+ < span class ="cline-any cline-yes "> 4x</ span >
207
+ < span class ="cline-any cline-yes "> 6x</ span >
208
+ < span class ="cline-any cline-yes "> 38x</ span >
209
+ < span class ="cline-any cline-yes "> 38x</ span >
210
+ < span class ="cline-any cline-yes "> 149x</ span >
211
+ < span class ="cline-any cline-yes "> 149x</ span >
212
+ < span class ="cline-any cline-yes "> 38x</ span >
213
+ < span class ="cline-any cline-yes "> 38x</ span >
214
+ < span class ="cline-any cline-yes "> 1x</ span >
215
+ < span class ="cline-any cline-yes "> 1x</ span >
216
+ < span class ="cline-any cline-yes "> 1x</ span >
217
+ < span class ="cline-any cline-yes "> 1x</ span >
218
+ < span class ="cline-any cline-yes "> 38x</ span >
219
+ < span class ="cline-any cline-yes "> 38x</ span >
220
+ < span class ="cline-any cline-yes "> 4x</ span >
221
+ < span class ="cline-any cline-yes "> 4x</ span >
222
+ < span class ="cline-any cline-yes "> 6x</ span >
223
+ < span class ="cline-any cline-yes "> 1x</ span >
224
+ < span class ="cline-any cline-yes "> 1x</ span >
225
+ < span class ="cline-any cline-yes "> 1x</ span >
226
+ < span class ="cline-any cline-yes "> 1x</ span >
227
+ < span class ="cline-any cline-yes "> 1x</ span >
228
+ < span class ="cline-any cline-neutral "> </ span > </ td > < td class ="text "> < pre class ="prettyprint lang-js "> /**
229
+ * @license Apache-2.0
230
+ *
231
+ * Copyright (c) 2024 The Stdlib Authors.
232
+ *
233
+ * Licensed under the Apache License, Version 2.0 (the "License");
234
+ * you may not use this file except in compliance with the License.
235
+ * You may obtain a copy of the License at
236
+ *
237
+ * http://www.apache.org/licenses/LICENSE-2.0
238
+ *
239
+ * Unless required by applicable law or agreed to in writing, software
240
+ * distributed under the License is distributed on an "AS IS" BASIS,
241
+ * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
242
+ * See the License for the specific language governing permissions and
243
+ * limitations under the License.
244
+ */
245
+
246
+ 'use strict';
247
+
248
+ /**
249
+ * Computes the greatest common divisor (gcd) of two single-precision floating point numbers using the binary GCD algorithm.
250
+ *
251
+ * ## References
252
+ *
253
+ * - Stein, Josef. 1967. "Computational problems associated with Racah algebra." _Journal of Computational Physics_ 1 (3): 397–405. doi:[10.1016/0021-9991(67)90047-2][@stein:1967].
254
+ *
255
+ * [@stein:1967]: https://doi.org/10.1016/0021-9991(67)90047-2
256
+ *
257
+ * @private
258
+ * @param {integer} a - first number
259
+ * @param {integer} b - second number
260
+ * @returns {integer} greatest common divisor
261
+ *
262
+ * @example
263
+ * var v = gcdf( 16777216.0, 65536.0 );
264
+ * // returns 65536
265
+ */
266
+ function gcdf( a, b ) {
267
+ var k = 1;
268
+ var t;
269
+
270
+ // Simple cases:
271
+ if ( a === 0 ) {
272
+ return b;
273
+ }
274
+ if ( b === 0 ) {
275
+ return a;
276
+ }
277
+ // Reduce `a` and/or `b` to odd numbers and keep track of the greatest power of 2 dividing both `a` and `b`...
278
+ while ( a%2 === 0 && b%2 === 0 ) {
279
+ a /= 2; // right shift
280
+ b /= 2; // right shift
281
+ k *= 2; // left shift
282
+ }
283
+ // Reduce `a` to an odd number...
284
+ while ( a%2 === 0 ) {
285
+ a /= 2; // right shift
286
+ }
287
+ // Henceforth, `a` is always odd...
288
+ while ( b ) {
289
+ // Remove all factors of 2 in `b`, as they are not common...
290
+ while ( b%2 === 0 ) {
291
+ b /= 2; // right shift
292
+ }
293
+ // `a` and `b` are both odd. Swap values such that `b` is the larger of the two values, and then set `b` to the difference (which is even)...
294
+ if ( a > b ) {
295
+ t = b;
296
+ b = a;
297
+ a = t;
298
+ }
299
+ b -= a; // b=0 iff b=a
300
+ }
301
+ // Restore common factors of 2...
302
+ return k * a;
303
+ }
304
+
305
+
306
+ // EXPORTS //
307
+
308
+ module.exports = gcdf;
309
+ </ pre > </ td > </ tr > </ table > </ pre >
310
+
311
+ < div class ='push '> </ div > <!-- for sticky footer -->
312
+ </ div > <!-- /wrapper -->
313
+ < div class ='footer quiet pad2 space-top1 center small '>
314
+ Code coverage generated by
315
+ < a href ="https://istanbul.js.org/ " target ="_blank " rel ="noopener noreferrer "> istanbul</ a >
316
+ at 2024-11-11T15:20:25.784Z
317
+ </ div >
318
+ < script src ="../../../../../prettify.js "> </ script >
319
+ < script >
320
+ window . onload = function ( ) {
321
+ prettyPrint ( ) ;
322
+ } ;
323
+ </ script >
324
+ < script src ="../../../../../sorter.js "> </ script >
325
+ < script src ="../../../../../block-navigation.js "> </ script >
326
+ </ body >
327
+ </ html >
328
+
0 commit comments