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

   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