OLD | NEW |
| (Empty) |
1 <!DOCTYPE html> | |
2 <head> | |
3 <!-- | |
4 Copyright (C) 2007 Apple Inc. All rights reserved. | |
5 | |
6 Redistribution and use in source and binary forms, with or without | |
7 modification, are permitted provided that the following conditions | |
8 are met: | |
9 1. Redistributions of source code must retain the above copyright | |
10 notice, this list of conditions and the following disclaimer. | |
11 2. Redistributions in binary form must reproduce the above copyright | |
12 notice, this list of conditions and the following disclaimer in the | |
13 documentation and/or other materials provided with the distribution. | |
14 | |
15 THIS SOFTWARE IS PROVIDED BY APPLE COMPUTER, INC. ``AS IS'' AND ANY | |
16 EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE | |
17 IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR | |
18 PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL APPLE COMPUTER, INC. OR | |
19 CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, | |
20 EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, | |
21 PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR | |
22 PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY | |
23 OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT | |
24 (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE | |
25 OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. | |
26 --> | |
27 | |
28 <title>SunSpider string-fasta</title> | |
29 <link rel="stylesheet" href="sunspider.css"> | |
30 </head> | |
31 | |
32 <body> | |
33 <h3>string-fasta</h3> | |
34 <div id="console"> | |
35 </div> | |
36 <script src="sunspider-record-result.js"></script> | |
37 <script> | |
38 | |
39 var _sunSpiderStartDate = new Date(); | |
40 | |
41 // The Great Computer Language Shootout | |
42 // http://shootout.alioth.debian.org | |
43 // | |
44 // Contributed by Ian Osgood | |
45 | |
46 var last = 42, A = 3877, C = 29573, M = 139968; | |
47 | |
48 function rand(max) { | |
49 last = (last * A + C) % M; | |
50 return max * last / M; | |
51 } | |
52 | |
53 var ALU = | |
54 "GGCCGGGCGCGGTGGCTCACGCCTGTAATCCCAGCACTTTGG" + | |
55 "GAGGCCGAGGCGGGCGGATCACCTGAGGTCAGGAGTTCGAGA" + | |
56 "CCAGCCTGGCCAACATGGTGAAACCCCGTCTCTACTAAAAAT" + | |
57 "ACAAAAATTAGCCGGGCGTGGTGGCGCGCGCCTGTAATCCCA" + | |
58 "GCTACTCGGGAGGCTGAGGCAGGAGAATCGCTTGAACCCGGG" + | |
59 "AGGCGGAGGTTGCAGTGAGCCGAGATCGCGCCACTGCACTCC" + | |
60 "AGCCTGGGCGACAGAGCGAGACTCCGTCTCAAAAA"; | |
61 | |
62 var IUB = { | |
63 a:0.27, c:0.12, g:0.12, t:0.27, | |
64 B:0.02, D:0.02, H:0.02, K:0.02, | |
65 M:0.02, N:0.02, R:0.02, S:0.02, | |
66 V:0.02, W:0.02, Y:0.02 | |
67 } | |
68 | |
69 var HomoSap = { | |
70 a: 0.3029549426680, | |
71 c: 0.1979883004921, | |
72 g: 0.1975473066391, | |
73 t: 0.3015094502008 | |
74 } | |
75 | |
76 function makeCumulative(table) { | |
77 var last = null; | |
78 for (var c in table) { | |
79 if (last) table[c] += table[last]; | |
80 last = c; | |
81 } | |
82 } | |
83 | |
84 function fastaRepeat(n, seq) { | |
85 var seqi = 0, lenOut = 60; | |
86 while (n>0) { | |
87 if (n<lenOut) lenOut = n; | |
88 if (seqi + lenOut < seq.length) { | |
89 ret = seq.substring(seqi, seqi+lenOut); | |
90 seqi += lenOut; | |
91 } else { | |
92 var s = seq.substring(seqi); | |
93 seqi = lenOut - s.length; | |
94 ret = s + seq.substring(0, seqi); | |
95 } | |
96 n -= lenOut; | |
97 } | |
98 } | |
99 | |
100 function fastaRandom(n, table) { | |
101 var line = new Array(60); | |
102 makeCumulative(table); | |
103 while (n>0) { | |
104 if (n<line.length) line = new Array(n); | |
105 for (var i=0; i<line.length; i++) { | |
106 var r = rand(1); | |
107 for (var c in table) { | |
108 if (r < table[c]) { | |
109 line[i] = c; | |
110 break; | |
111 } | |
112 } | |
113 } | |
114 ret = line.join(''); | |
115 n -= line.length; | |
116 } | |
117 } | |
118 | |
119 var ret; | |
120 | |
121 var count = 7; | |
122 ret = fastaRepeat(2*count*100000, ALU); | |
123 ret = fastaRandom(3*count*1000, IUB); | |
124 ret = fastaRandom(5*count*1000, HomoSap); | |
125 | |
126 | |
127 | |
128 var _sunSpiderInterval = new Date() - _sunSpiderStartDate; | |
129 | |
130 record(_sunSpiderInterval); | |
131 </script> | |
132 | |
133 | |
134 </body> | |
135 </html> | |
OLD | NEW |