js/src/devtools/jint/sunspider/string-fasta.js

branch
TOR_BUG_3246
changeset 7
129ffea94266
equal deleted inserted replaced
-1:000000000000 0:7a447934739d
1 // The Great Computer Language Shootout
2 // http://shootout.alioth.debian.org
3 //
4 // Contributed by Ian Osgood
5
6 var last = 42, A = 3877, C = 29573, M = 139968;
7
8 function rand(max) {
9 last = (last * A + C) % M;
10 return max * last / M;
11 }
12
13 var ALU =
14 "GGCCGGGCGCGGTGGCTCACGCCTGTAATCCCAGCACTTTGG" +
15 "GAGGCCGAGGCGGGCGGATCACCTGAGGTCAGGAGTTCGAGA" +
16 "CCAGCCTGGCCAACATGGTGAAACCCCGTCTCTACTAAAAAT" +
17 "ACAAAAATTAGCCGGGCGTGGTGGCGCGCGCCTGTAATCCCA" +
18 "GCTACTCGGGAGGCTGAGGCAGGAGAATCGCTTGAACCCGGG" +
19 "AGGCGGAGGTTGCAGTGAGCCGAGATCGCGCCACTGCACTCC" +
20 "AGCCTGGGCGACAGAGCGAGACTCCGTCTCAAAAA";
21
22 var IUB = {
23 a:0.27, c:0.12, g:0.12, t:0.27,
24 B:0.02, D:0.02, H:0.02, K:0.02,
25 M:0.02, N:0.02, R:0.02, S:0.02,
26 V:0.02, W:0.02, Y:0.02
27 }
28
29 var HomoSap = {
30 a: 0.3029549426680,
31 c: 0.1979883004921,
32 g: 0.1975473066391,
33 t: 0.3015094502008
34 }
35
36 function makeCumulative(table) {
37 var last = null;
38 /* BEGIN LOOP */
39 for (var c in table) {
40 if (last) table[c] += table[last];
41 last = c;
42 }
43 /* END LOOP */
44 }
45
46 function fastaRepeat(n, seq) {
47 var seqi = 0, lenOut = 60;
48 /* BEGIN LOOP */
49 while (n>0) {
50 if (n<lenOut) lenOut = n;
51 if (seqi + lenOut < seq.length) {
52 ret = seq.substring(seqi, seqi+lenOut);
53 seqi += lenOut;
54 } else {
55 var s = seq.substring(seqi);
56 seqi = lenOut - s.length;
57 ret = s + seq.substring(0, seqi);
58 }
59 n -= lenOut;
60 }
61 /* END LOOP */
62 }
63
64 function fastaRandom(n, table) {
65 var line = new Array(60);
66 makeCumulative(table);
67 /* BEGIN LOOP */
68 while (n>0) {
69 if (n<line.length) line = new Array(n);
70 /* BEGIN LOOP */
71 for (var i=0; i<line.length; i++) {
72 var r = rand(1);
73 /* BEGIN LOOP */
74 for (var c in table) {
75 if (r < table[c]) {
76 line[i] = c;
77 break;
78 }
79 }
80 /* END LOOP */
81 }
82 /* END LOOP */
83 ret = line.join('');
84 n -= line.length;
85 }
86 /* END LOOP */
87 }
88
89 var ret;
90
91 var count = 7;
92 ret = fastaRepeat(2*count*100000, ALU);
93 ret = fastaRandom(3*count*1000, IUB);
94 ret = fastaRandom(5*count*1000, HomoSap);
95

mercurial