OLD | NEW |
(Empty) | |
| 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 <link rel="stylesheet" href="sunspider.css"> |
| 30 </head> |
| 31 |
| 32 <body> |
| 33 <h3>math-spectral-norm</h3> |
| 34 <div id="console"> |
| 35 </div> |
| 36 <script src="sunspider-record-result.js"></script> |
| 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 record(_sunSpiderInterval); |
| 97 </script> |
| 98 |
| 99 |
| 100 </body> |
| 101 </html> |
OLD | NEW |