|
1 <!DOCTYPE html> |
|
2 <head> |
|
3 <!-- |
|
4 Copyright (C) 2007 Apple Inc. All rights reserved. |
|
5 |
|
6 Redistribution and use in source and binary forms, with or without |
|
7 modification, are permitted provided that the following conditions |
|
8 are met: |
|
9 1. Redistributions of source code must retain the above copyright |
|
10 notice, this list of conditions and the following disclaimer. |
|
11 2. Redistributions in binary form must reproduce the above copyright |
|
12 notice, this list of conditions and the following disclaimer in the |
|
13 documentation and/or other materials provided with the distribution. |
|
14 |
|
15 THIS SOFTWARE IS PROVIDED BY APPLE COMPUTER, INC. ``AS IS'' AND ANY |
|
16 EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE |
|
17 IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR |
|
18 PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL APPLE COMPUTER, INC. OR |
|
19 CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, |
|
20 EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, |
|
21 PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR |
|
22 PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY |
|
23 OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT |
|
24 (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE |
|
25 OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. |
|
26 --> |
|
27 |
|
28 <title>SunSpider math-spectral-norm</title> |
|
29 |
|
30 </head> |
|
31 |
|
32 <body> |
|
33 <h3>math-spectral-norm</h3> |
|
34 <div id="console"> |
|
35 </div> |
|
36 |
|
37 <script> |
|
38 |
|
39 var _sunSpiderStartDate = new Date(); |
|
40 |
|
41 // The Great Computer Language Shootout |
|
42 // http://shootout.alioth.debian.org/ |
|
43 // |
|
44 // contributed by Ian Osgood |
|
45 |
|
46 function A(i,j) { |
|
47 return 1/((i+j)*(i+j+1)/2+i+1); |
|
48 } |
|
49 |
|
50 function Au(u,v) { |
|
51 for (var i=0; i<u.length; ++i) { |
|
52 var t = 0; |
|
53 for (var j=0; j<u.length; ++j) |
|
54 t += A(i,j) * u[j]; |
|
55 v[i] = t; |
|
56 } |
|
57 } |
|
58 |
|
59 function Atu(u,v) { |
|
60 for (var i=0; i<u.length; ++i) { |
|
61 var t = 0; |
|
62 for (var j=0; j<u.length; ++j) |
|
63 t += A(j,i) * u[j]; |
|
64 v[i] = t; |
|
65 } |
|
66 } |
|
67 |
|
68 function AtAu(u,v,w) { |
|
69 Au(u,w); |
|
70 Atu(w,v); |
|
71 } |
|
72 |
|
73 function spectralnorm(n) { |
|
74 var i, u=[], v=[], w=[], vv=0, vBv=0; |
|
75 for (i=0; i<n; ++i) { |
|
76 u[i] = 1; v[i] = w[i] = 0; |
|
77 } |
|
78 for (i=0; i<10; ++i) { |
|
79 AtAu(u,v,w); |
|
80 AtAu(v,u,w); |
|
81 } |
|
82 for (i=0; i<n; ++i) { |
|
83 vBv += u[i]*v[i]; |
|
84 vv += v[i]*v[i]; |
|
85 } |
|
86 return Math.sqrt(vBv/vv); |
|
87 } |
|
88 |
|
89 for (var i = 6; i <= 48; i *= 2) { |
|
90 spectralnorm(i); |
|
91 } |
|
92 |
|
93 |
|
94 var _sunSpiderInterval = new Date() - _sunSpiderStartDate; |
|
95 |
|
96 document.getElementById("console").innerHTML = _sunSpiderInterval; |
|
97 </script> |
|
98 |
|
99 |
|
100 </body> |
|
101 </html> |