Projekt

Obecné

Profil

Stáhnout (2.07 KB) Statistiky
| Větev: | Revize:
1
'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));
(2-2/2)