| 1 |
3a515b92
|
cagy
|
'use strict';
|
| 2 |
|
|
|
| 3 |
|
|
function genCombMulTo (alen, blen) {
|
| 4 |
|
|
var len = alen + blen - 1;
|
| 5 |
|
|
var src = [
|
| 6 |
|
|
'var a = self.words;',
|
| 7 |
|
|
'var b = num.words;',
|
| 8 |
|
|
'var o = out.words;',
|
| 9 |
|
|
'var c = 0;',
|
| 10 |
|
|
'var lo;',
|
| 11 |
|
|
'var mid;',
|
| 12 |
|
|
'var hi;'
|
| 13 |
|
|
];
|
| 14 |
|
|
for (var i = 0; i < alen; i++) {
|
| 15 |
|
|
src.push('var a' + i + ' = a[' + i + '] | 0;');
|
| 16 |
|
|
src.push('var al' + i + ' = a' + i + ' & 0x1fff;');
|
| 17 |
|
|
src.push('var ah' + i + ' = a' + i + ' >>> 13;');
|
| 18 |
|
|
}
|
| 19 |
|
|
for (i = 0; i < blen; i++) {
|
| 20 |
|
|
src.push('var b' + i + ' = b[' + i + '] | 0;');
|
| 21 |
|
|
src.push('var bl' + i + ' = b' + i + ' & 0x1fff;');
|
| 22 |
|
|
src.push('var bh' + i + ' = b' + i + ' >>> 13;');
|
| 23 |
|
|
}
|
| 24 |
|
|
src.push('');
|
| 25 |
|
|
src.push('out.negative = self.negative ^ num.negative;');
|
| 26 |
|
|
src.push('out.length = ' + len + ';');
|
| 27 |
|
|
|
| 28 |
|
|
for (var k = 0; k < len; k++) {
|
| 29 |
|
|
var minJ = Math.max(0, k - alen + 1);
|
| 30 |
|
|
var maxJ = Math.min(k, blen - 1);
|
| 31 |
|
|
|
| 32 |
|
|
src.push('\/* k = ' + k + ' *\/');
|
| 33 |
|
|
src.push('lo = Math.imul(al' + (k - minJ) + ', bl' + minJ + ');');
|
| 34 |
|
|
src.push('mid = Math.imul(al' + (k - minJ) + ', bh' + minJ + ');');
|
| 35 |
|
|
src.push(
|
| 36 |
|
|
'mid = (mid + Math.imul(ah' + (k - minJ) + ', bl' + minJ + ')) | 0;');
|
| 37 |
|
|
src.push('hi = Math.imul(ah' + (k - minJ) + ', bh' + minJ + ');');
|
| 38 |
|
|
|
| 39 |
|
|
for (var j = minJ + 1; j <= maxJ; j++) {
|
| 40 |
|
|
i = k - j;
|
| 41 |
|
|
|
| 42 |
|
|
src.push('lo = (lo + Math.imul(al' + i + ', bl' + j + ')) | 0;');
|
| 43 |
|
|
src.push('mid = (mid + Math.imul(al' + i + ', bh' + j + ')) | 0;');
|
| 44 |
|
|
src.push('mid = (mid + Math.imul(ah' + i + ', bl' + j + ')) | 0;');
|
| 45 |
|
|
src.push('hi = (hi + Math.imul(ah' + i + ', bh' + j + ')) | 0;');
|
| 46 |
|
|
}
|
| 47 |
|
|
|
| 48 |
|
|
src.push('var w' + k + ' = (((c + lo) | 0) + ((mid & 0x1fff) << 13)) | 0;');
|
| 49 |
|
|
src.push('c = (((hi + (mid >>> 13)) | 0) + (w' + k + ' >>> 26)) | 0;');
|
| 50 |
|
|
src.push('w' + k + ' &= 0x3ffffff;');
|
| 51 |
|
|
}
|
| 52 |
|
|
// Store in separate step for better memory access
|
| 53 |
|
|
for (k = 0; k < len; k++) {
|
| 54 |
|
|
src.push('o[' + k + '] = w' + k + ';');
|
| 55 |
|
|
}
|
| 56 |
|
|
src.push('if (c !== 0) {',
|
| 57 |
|
|
' o[' + k + '] = c;',
|
| 58 |
|
|
' out.length++;',
|
| 59 |
|
|
'}',
|
| 60 |
|
|
'return out;');
|
| 61 |
|
|
|
| 62 |
|
|
return src.join('\n');
|
| 63 |
|
|
}
|
| 64 |
|
|
|
| 65 |
|
|
console.log(genCombMulTo(10, 10));
|