mirror of
https://gitlab.winehq.org/wine/wine-gecko.git
synced 2024-09-13 09:24:08 -07:00
86 lines
1.8 KiB
JavaScript
86 lines
1.8 KiB
JavaScript
// The Great Computer Language Shootout
|
|
// http://shootout.alioth.debian.org
|
|
//
|
|
// Contributed by Ian Osgood
|
|
|
|
var last = 42, A = 3877, C = 29573, M = 139968;
|
|
|
|
function rand(max) {
|
|
last = (last * A + C) % M;
|
|
return max * last / M;
|
|
}
|
|
|
|
var ALU =
|
|
"GGCCGGGCGCGGTGGCTCACGCCTGTAATCCCAGCACTTTGG" +
|
|
"GAGGCCGAGGCGGGCGGATCACCTGAGGTCAGGAGTTCGAGA" +
|
|
"CCAGCCTGGCCAACATGGTGAAACCCCGTCTCTACTAAAAAT" +
|
|
"ACAAAAATTAGCCGGGCGTGGTGGCGCGCGCCTGTAATCCCA" +
|
|
"GCTACTCGGGAGGCTGAGGCAGGAGAATCGCTTGAACCCGGG" +
|
|
"AGGCGGAGGTTGCAGTGAGCCGAGATCGCGCCACTGCACTCC" +
|
|
"AGCCTGGGCGACAGAGCGAGACTCCGTCTCAAAAA";
|
|
|
|
var IUB = {
|
|
a:0.27, c:0.12, g:0.12, t:0.27,
|
|
B:0.02, D:0.02, H:0.02, K:0.02,
|
|
M:0.02, N:0.02, R:0.02, S:0.02,
|
|
V:0.02, W:0.02, Y:0.02
|
|
}
|
|
|
|
var HomoSap = {
|
|
a: 0.3029549426680,
|
|
c: 0.1979883004921,
|
|
g: 0.1975473066391,
|
|
t: 0.3015094502008
|
|
}
|
|
|
|
function makeCumulative(table) {
|
|
var last = null;
|
|
for (var c in table) {
|
|
if (last) table[c] += table[last];
|
|
last = c;
|
|
}
|
|
}
|
|
|
|
function fastaRepeat(n, seq) {
|
|
var seqi = 0, lenOut = 60;
|
|
while (n>0) {
|
|
if (n<lenOut) lenOut = n;
|
|
if (seqi + lenOut < seq.length) {
|
|
ret = seq.substring(seqi, seqi+lenOut);
|
|
seqi += lenOut;
|
|
} else {
|
|
var s = seq.substring(seqi);
|
|
seqi = lenOut - s.length;
|
|
ret = s + seq.substring(0, seqi);
|
|
}
|
|
n -= lenOut;
|
|
}
|
|
}
|
|
|
|
function fastaRandom(n, table) {
|
|
var line = new Array(60);
|
|
makeCumulative(table);
|
|
while (n>0) {
|
|
if (n<line.length) line = new Array(n);
|
|
for (var i=0; i<line.length; i++) {
|
|
var r = rand(1);
|
|
for (var c in table) {
|
|
if (r < table[c]) {
|
|
line[i] = c;
|
|
break;
|
|
}
|
|
}
|
|
}
|
|
ret = line.join('');
|
|
n -= line.length;
|
|
}
|
|
}
|
|
|
|
var ret;
|
|
|
|
var count = 7;
|
|
ret = fastaRepeat(2*count*100000, ALU);
|
|
ret = fastaRandom(3*count*1000, IUB);
|
|
ret = fastaRandom(5*count*1000, HomoSap);
|
|
|