testing/performance/talos/page_load_test/jss/real-spectral-norm-2.html

   1 <html><head>
   2 <!-- MOZ_INSERT_CONTENT_HOOK -->
   3 <script src = runner.js></script>
   4 <script>
   5 var onlyName = 'Spectral Norm', onlyNum = 16;
   6 function thisTest() {
   7
   8 // The Great Computer Language Shootout
   9 // http://shootout.alioth.debian.org/
  10 //
  11 // contributed by Ian Osgood
  12
  13 function A(i,j) {
  14   return 1/((i+j)*(i+j+1)/2+i+1);
  15 }
  16
  17 function Au(u,v) {
  18   for (var i=0; i<u.length; ++i) {
  19     var t = 0;
  20     for (var j=0; j<u.length; ++j)
  21       t += A(i,j) * u[j];
  22     v[i] = t;
  23   }
  24 }
  25
  26 function Atu(u,v) {
  27   for (var i=0; i<u.length; ++i) {
  28     var t = 0;
  29     for (var j=0; j<u.length; ++j)
  30       t += A(j,i) * u[j];
  31     v[i] = t;
  32   }
  33 }
  34
  35 function AtAu(u,v,w) {
  36   Au(u,w);
  37   Atu(w,v);
  38 }
  39
  40 function spectralnorm(n) {
  41   var i, u=[], v=[], w=[], vv=0, vBv=0;
  42   for (i=0; i<n; ++i) {
  43     u[i] = 1; v[i] = w[i] = 0;
  44   }
  45   for (i=0; i<10; ++i) {
  46     AtAu(u,v,w);
  47     AtAu(v,u,w);
  48   }
  49   for (i=0; i<n; ++i) {
  50     vBv += u[i]*v[i];
  51     vv  += v[i]*v[i];
  52   }
  53   return Math.sqrt(vBv/vv);
  54 }
  55
  56 startTest("real-spectral-norm");
  57
  58 for ( var i = 8; i <= 64; i *= 2 ) (function(i){
  59   test( "Spectral Norm", i, function(){
  60     spectralnorm(i);
  61   });
  62 })(i);
  63
  64 endTest();
  65 }
  66 </script>
  67 <body onload="thisTest()"></body></html>