michael@0: /*
michael@0:  * Copyright (C) 2007 Apple Inc.  All rights reserved.
michael@0:  *
michael@0:  * Redistribution and use in source and binary forms, with or without
michael@0:  * modification, are permitted provided that the following conditions
michael@0:  * are met:
michael@0:  * 1. Redistributions of source code must retain the above copyright
michael@0:  *    notice, this list of conditions and the following disclaimer.
michael@0:  * 2. Redistributions in binary form must reproduce the above copyright
michael@0:  *    notice, this list of conditions and the following disclaimer in the
michael@0:  *    documentation and/or other materials provided with the distribution.
michael@0:  *
michael@0:  * THIS SOFTWARE IS PROVIDED BY APPLE COMPUTER, INC. ``AS IS'' AND ANY
michael@0:  * EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
michael@0:  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
michael@0:  * PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL APPLE COMPUTER, INC. OR
michael@0:  * CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
michael@0:  * EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
michael@0:  * PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
michael@0:  * PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY
michael@0:  * OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
michael@0:  * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
michael@0:  * OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 
michael@0:  */
michael@0: 
michael@0: var loops = 15
michael@0: var nx = 120
michael@0: var nz = 120
michael@0: 
michael@0: function morph(a, f) {
michael@0:     var PI2nx = Math.PI * 8/nx
michael@0:     var sin = Math.sin
michael@0:     var f30 = -(50 * sin(f*Math.PI*2))
michael@0:     
michael@0:     for (var i = 0; i < nz; ++i) {
michael@0:         for (var j = 0; j < nx; ++j) {
michael@0:             a[3*(i*nx+j)+1]    = sin((j-1) * PI2nx ) * -f30
michael@0:         }
michael@0:     }
michael@0: }
michael@0: 
michael@0:     
michael@0: var a = Array()
michael@0: for (var i=0; i < nx*nz*3; ++i) 
michael@0:     a[i] = 0
michael@0: 
michael@0: for (var i = 0; i < loops; ++i) {
michael@0:     morph(a, i/loops)
michael@0: }
michael@0: 
michael@0: testOutput = 0;
michael@0: for (var i = 0; i < nx; i++)
michael@0:     testOutput += a[3*(i*nx+i)+1];
michael@0: a = null;
michael@0: 
michael@0: /* not based on any mathematical error calculation.*/
michael@0: acceptableDelta = 4e-15
michael@0: 
michael@0: assertEq((testOutput - 6.394884621840902e-14) < acceptableDelta, true);