js/src/jit-test/tests/sunspider/check-math-spectral-norm.js

Wed, 31 Dec 2014 06:09:35 +0100

author
Michael Schloh von Bennewitz <michael@schloh.com>
date
Wed, 31 Dec 2014 06:09:35 +0100
changeset 0
6474c204b198
permissions
-rw-r--r--

Cloned upstream origin tor-browser at tor-browser-31.3.0esr-4.5-1-build1
revision ID fc1c9ff7c1b2defdbc039f12214767608f46423f for hacking purpose.

     1 // The Great Computer Language Shootout
     2 // http://shootout.alioth.debian.org/
     3 //
     4 // contributed by Ian Osgood
     6 function A(i,j) {
     7   return 1/((i+j)*(i+j+1)/2+i+1);
     8 }
    10 function Au(u,v) {
    11   for (var i=0; i<u.length; ++i) {
    12     var t = 0;
    13     for (var j=0; j<u.length; ++j)
    14       t += A(i,j) * u[j];
    15     v[i] = t;
    16   }
    17 }
    19 function Atu(u,v) {
    20   for (var i=0; i<u.length; ++i) {
    21     var t = 0;
    22     for (var j=0; j<u.length; ++j)
    23       t += A(j,i) * u[j];
    24     v[i] = t;
    25   }
    26 }
    28 function AtAu(u,v,w) {
    29   Au(u,w);
    30   Atu(w,v);
    31 }
    33 function spectralnorm(n) {
    34   var i, u=[], v=[], w=[], vv=0, vBv=0;
    35   for (i=0; i<n; ++i) {
    36     u[i] = 1; v[i] = w[i] = 0;
    37   }
    38   for (i=0; i<10; ++i) {
    39     AtAu(u,v,w);
    40     AtAu(v,u,w);
    41   }
    42   for (i=0; i<n; ++i) {
    43     vBv += u[i]*v[i];
    44     vv  += v[i]*v[i];
    45   }
    46   return Math.sqrt(vBv/vv);
    47 }
    49 var actual = '';
    50 for (var i = 6; i <= 48; i *= 2) {
    51   actual += spectralnorm(i) + ',';
    52 }
    53 assertEq(actual, "1.2657786149754053,1.2727355112619148,1.273989979775574,1.274190125290389,");

mercurial