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