| <!DOCTYPE html> |
| <head> |
| <!-- |
| Copyright (C) 2007 Apple Inc. All rights reserved. |
| |
| Redistribution and use in source and binary forms, with or without |
| modification, are permitted provided that the following conditions |
| are met: |
| 1. Redistributions of source code must retain the above copyright |
| notice, this list of conditions and the following disclaimer. |
| 2. Redistributions in binary form must reproduce the above copyright |
| notice, this list of conditions and the following disclaimer in the |
| documentation and/or other materials provided with the distribution. |
| |
| THIS SOFTWARE IS PROVIDED BY APPLE COMPUTER, INC. ``AS IS'' AND ANY |
| EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE |
| IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR |
| PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL APPLE COMPUTER, INC. OR |
| CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, |
| EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, |
| PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR |
| PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY |
| OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT |
| (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE |
| OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. |
| --> |
| |
| <title>SunSpider math-spectral-norm</title> |
| |
| </head> |
| |
| <body> |
| <h3>math-spectral-norm</h3> |
| <div id="console"> |
| </div> |
| |
| <script> |
| |
| var _sunSpiderStartDate = new Date(); |
| |
| // The Great Computer Language Shootout |
| // http://shootout.alioth.debian.org/ |
| // |
| // contributed by Ian Osgood |
| |
| function A(i,j) { |
| return 1/((i+j)*(i+j+1)/2+i+1); |
| } |
| |
| function Au(u,v) { |
| for (var i=0; i<u.length; ++i) { |
| var t = 0; |
| for (var j=0; j<u.length; ++j) |
| t += A(i,j) * u[j]; |
| v[i] = t; |
| } |
| } |
| |
| function Atu(u,v) { |
| for (var i=0; i<u.length; ++i) { |
| var t = 0; |
| for (var j=0; j<u.length; ++j) |
| t += A(j,i) * u[j]; |
| v[i] = t; |
| } |
| } |
| |
| function AtAu(u,v,w) { |
| Au(u,w); |
| Atu(w,v); |
| } |
| |
| function spectralnorm(n) { |
| var i, u=[], v=[], w=[], vv=0, vBv=0; |
| for (i=0; i<n; ++i) { |
| u[i] = 1; v[i] = w[i] = 0; |
| } |
| for (i=0; i<10; ++i) { |
| AtAu(u,v,w); |
| AtAu(v,u,w); |
| } |
| for (i=0; i<n; ++i) { |
| vBv += u[i]*v[i]; |
| vv += v[i]*v[i]; |
| } |
| return Math.sqrt(vBv/vv); |
| } |
| |
| for (var i = 6; i <= 48; i *= 2) { |
| spectralnorm(i); |
| } |
| |
| |
| var _sunSpiderInterval = new Date() - _sunSpiderStartDate; |
| |
| document.getElementById("console").innerHTML = _sunSpiderInterval; |
| </script> |
| |
| |
| </body> |
| </html> |