js/src/jit-test/tests/sunspider/check-string-fasta.js

changeset 0
6474c204b198
     1.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
     1.2 +++ b/js/src/jit-test/tests/sunspider/check-string-fasta.js	Wed Dec 31 06:09:35 2014 +0100
     1.3 @@ -0,0 +1,90 @@
     1.4 +// The Great Computer Language Shootout
     1.5 +//  http://shootout.alioth.debian.org
     1.6 +//
     1.7 +//  Contributed by Ian Osgood
     1.8 +
     1.9 +var last = 42, A = 3877, C = 29573, M = 139968;
    1.10 +
    1.11 +function rand(max) {
    1.12 +  last = (last * A + C) % M;
    1.13 +  return max * last / M;
    1.14 +}
    1.15 +
    1.16 +var ALU =
    1.17 +  "GGCCGGGCGCGGTGGCTCACGCCTGTAATCCCAGCACTTTGG" +
    1.18 +  "GAGGCCGAGGCGGGCGGATCACCTGAGGTCAGGAGTTCGAGA" +
    1.19 +  "CCAGCCTGGCCAACATGGTGAAACCCCGTCTCTACTAAAAAT" +
    1.20 +  "ACAAAAATTAGCCGGGCGTGGTGGCGCGCGCCTGTAATCCCA" +
    1.21 +  "GCTACTCGGGAGGCTGAGGCAGGAGAATCGCTTGAACCCGGG" +
    1.22 +  "AGGCGGAGGTTGCAGTGAGCCGAGATCGCGCCACTGCACTCC" +
    1.23 +  "AGCCTGGGCGACAGAGCGAGACTCCGTCTCAAAAA";
    1.24 +
    1.25 +var IUB = {
    1.26 +  a:0.27, c:0.12, g:0.12, t:0.27,
    1.27 +  B:0.02, D:0.02, H:0.02, K:0.02,
    1.28 +  M:0.02, N:0.02, R:0.02, S:0.02,
    1.29 +  V:0.02, W:0.02, Y:0.02
    1.30 +}
    1.31 +
    1.32 +var HomoSap = {
    1.33 +  a: 0.3029549426680,
    1.34 +  c: 0.1979883004921,
    1.35 +  g: 0.1975473066391,
    1.36 +  t: 0.3015094502008
    1.37 +}
    1.38 +
    1.39 +function makeCumulative(table) {
    1.40 +  var last = null;
    1.41 +  for (var c in table) {
    1.42 +    if (last) table[c] += table[last];
    1.43 +    last = c;
    1.44 +  }
    1.45 +}
    1.46 +
    1.47 +function fastaRepeat(n, seq) {
    1.48 +  var seqi = 0, lenOut = 60;
    1.49 +  while (n>0) {
    1.50 +    if (n<lenOut) lenOut = n;
    1.51 +    if (seqi + lenOut < seq.length) {
    1.52 +      ret = seq.substring(seqi, seqi+lenOut);
    1.53 +      seqi += lenOut;
    1.54 +    } else {
    1.55 +      var s = seq.substring(seqi);
    1.56 +      seqi = lenOut - s.length;
    1.57 +      ret = s + seq.substring(0, seqi);
    1.58 +    }
    1.59 +    n -= lenOut;
    1.60 +  }
    1.61 +  return ret;
    1.62 +}
    1.63 +
    1.64 +function fastaRandom(n, table) {
    1.65 +  var line = new Array(60);
    1.66 +  makeCumulative(table);
    1.67 +  while (n>0) {
    1.68 +    if (n<line.length) line = new Array(n);
    1.69 +    for (var i=0; i<line.length; i++) {
    1.70 +      var r = rand(1);
    1.71 +      for (var c in table) {
    1.72 +        if (r < table[c]) {
    1.73 +          line[i] = c;
    1.74 +          break;
    1.75 +        }
    1.76 +      }
    1.77 +    }
    1.78 +    ret = line.join('');
    1.79 +    n -= line.length;
    1.80 +  }
    1.81 +  return ret;
    1.82 +}
    1.83 +
    1.84 +var ret;
    1.85 +
    1.86 +var count = 7;
    1.87 +var actual1 = fastaRepeat(2*count*100000, ALU);
    1.88 +var actual2 = fastaRandom(3*count*1000, IUB);
    1.89 +var actual3 = fastaRandom(5*count*1000, HomoSap);
    1.90 +
    1.91 +assertEq(actual1, "CAAAAAGGCCGGGCGCGGTG");
    1.92 +assertEq(actual2, "VtttaDtKgcaaWaaaaatSccMcVatgtKgtaKgcgatatgtagtSaaaDttatacaaa");
    1.93 +assertEq(actual3, "ttggctatatttatgttgga");

mercurial