Chris@3: 
Chris@3: /* for a phase vocoder, we probably want 2048-point real-to-complex
Chris@3:  * FFTs (if available) */
Chris@3: 
Chris@3: function inputReals(size) {
Chris@3:     var result = new Array(size);
Chris@3:     for (var i = 0; i < result.length; i++)
Chris@3: 	result[i] = (i % 20) / 10.0 - 1.0;
Chris@3:     return result;
Chris@3: }
Chris@3: 
Chris@3: function zeroReals(size) {
Chris@3:     var result = new Array(size);
Chris@3:     for (var i = 0; i < result.length; i++)
Chris@3: 	result[i] = 0.0;
Chris@3:     return result;
Chris@3: }
Chris@3: 
Chris@3: function inputReal64s(size) {
Chris@3:     var result = new Float64Array(size);
Chris@3:     for (var i = 0; i < result.length; i++)
Chris@3: 	result[i] = (i % 20) / 10.0 - 1.0;
Chris@3:     return result;
Chris@3: }
Chris@3: 
Chris@3: function inputComplexArray(size) {
Chris@3:     var result = new complex_array.ComplexArray(size);
Chris@3:     for (var i = 0; i < size; i++) {
Chris@3: 	result.real[i] = (i % 20) / 10.0 - 1.0;
Chris@3: 	result.imag[i] = 0.0;
Chris@3:     }
Chris@3:     return result;
Chris@3: }
Chris@3: 
Chris@3: var iterations = 2150;
Chris@3: var size = 2048;
Chris@3: 
Chris@3: function report(name, start, middle, end, total) {
Chris@3:     document.getElementById(name + "-result").innerHTML = total;
Chris@3:     document.getElementById(name + "-1").innerHTML =
Chris@3: 	Math.round(middle - start) + " ms";
Chris@3:     document.getElementById(name + "-2").innerHTML =
Chris@3: 	Math.round(end - middle) + " ms";
Chris@3:     document.getElementById(name + "-itr").innerHTML =
Chris@3: 	Math.round((1000.0 / ((end - middle) / iterations))) + " itr/sec";
Chris@3: }
Chris@3: 
Chris@3: function testNayuki() {
Chris@3: 
Chris@3:     var start = performance.now();
Chris@3:     var middle = start;
Chris@3:     var end = start;
Chris@3: 
Chris@3:     var total = 0.0;
Chris@3: 
Chris@3:     for (var i = 0; i < 2*iterations; ++i) {
Chris@3: 	if (i == iterations) {
Chris@3: 	    middle = performance.now();
Chris@3: 	}
Chris@3: 	var real = inputReals(size);
Chris@3: 	var imag = zeroReals(size);
Chris@3: 	transform(real, imag);
Chris@3: 	for (var j = 0; j < size; ++j) {
Chris@3: 	    total += real[j] + imag[j];
Chris@3: 	}
Chris@3:     }
Chris@3: 
Chris@3:     var end = performance.now();
Chris@3: 
Chris@3:     report("nayuki", start, middle, end, total);
Chris@3: }
Chris@3: 
Chris@3: function testNockert() {
Chris@3:     
Chris@3:     var fft = new FFT.complex(size, false);
Chris@3:     
Chris@3:     var start = performance.now();
Chris@3:     var middle = start;
Chris@3:     var end = start;
Chris@3: 
Chris@3:     total = 0.0;
Chris@3: 
Chris@3:     for (var i = 0; i < 2*iterations; ++i) {
Chris@3: 	if (i == iterations) {
Chris@3: 	    middle = performance.now();
Chris@3: 	}
Chris@3: 	var ri = inputReal64s(size);
Chris@3: 	var co = new Float64Array(2 * size);
Chris@3: 	fft.simple(co, ri, 'real');
Chris@3: 	for (var j = 0; j < 2 * size; ++j) {
Chris@3: 	    total += co[j];
Chris@3: 	}
Chris@3:     }
Chris@3: 
Chris@3:     var end = performance.now();
Chris@3: 
Chris@3:     report("nockert", start, middle, end, total);
Chris@3: }
Chris@3: 
Chris@3: function testDntj() {
Chris@3: 
Chris@3:     var start = performance.now();
Chris@3:     var middle = start;
Chris@3:     var end = start;
Chris@3: 
Chris@3:     total = 0.0;
Chris@3: 
Chris@3:     for (var i = 0; i < 2*iterations; ++i) {
Chris@3: 	if (i == iterations) {
Chris@3: 	    middle = performance.now();
Chris@3: 	}
Chris@3: 	var ci = inputComplexArray(size);
Chris@3: 	var co = ci.FFT();
Chris@3: 	for (var j = 0; j < size; ++j) {
Chris@3: 	    total += co.real[j] + co.imag[j];
Chris@3: 	}
Chris@3:     }
Chris@3: 
Chris@3:     var end = performance.now();
Chris@3: 
Chris@3:     report("dntj", start, middle, end, total * Math.sqrt(size));
Chris@3: }
Chris@3: 
Chris@3: function test() {
Chris@3: 
Chris@3:     document.getElementById("test-description").innerHTML =
Chris@3: 	"For " + iterations + " iterations per half (" + 2*iterations
Chris@3: 	+ " total) of size " + size;
Chris@3: 
Chris@3:     testNayuki();
Chris@3:     testNockert();
Chris@3:     testDntj();
Chris@3: }
Chris@3: 
Chris@3: window.onload = test;
Chris@3: