829N/A/*
2362N/A * Copyright (c) 2007, Oracle and/or its affiliates. All rights reserved.
829N/A * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS FILE HEADER.
829N/A *
829N/A * This code is free software; you can redistribute it and/or modify it
829N/A * under the terms of the GNU General Public License version 2 only, as
2362N/A * published by the Free Software Foundation. Oracle designates this
829N/A * particular file as subject to the "Classpath" exception as provided
2362N/A * by Oracle in the LICENSE file that accompanied this code.
829N/A *
829N/A * This code is distributed in the hope that it will be useful, but WITHOUT
829N/A * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
829N/A * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
829N/A * version 2 for more details (a copy is included in the LICENSE file that
829N/A * accompanied this code).
829N/A *
829N/A * You should have received a copy of the GNU General Public License version
829N/A * 2 along with this work; if not, write to the Free Software Foundation,
829N/A * Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA.
829N/A *
2362N/A * Please contact Oracle, 500 Oracle Parkway, Redwood Shores, CA 94065 USA
2362N/A * or visit www.oracle.com if you need additional information or have any
2362N/A * questions.
829N/A */
829N/A
829N/A/* @test
829N/A @summary Test SoftCubicResampler interpolate method */
829N/A
829N/Aimport java.io.File;
829N/Aimport java.io.FileOutputStream;
829N/Aimport java.io.IOException;
829N/A
829N/Aimport javax.sound.sampled.*;
829N/A
829N/Aimport com.sun.media.sound.*;
829N/A
829N/Apublic class Interpolate {
829N/A
829N/A private static float getResamplerTestValue(double i)
829N/A {
829N/A return (float)Math.sin(i / 10.0);
829N/A }
829N/A
829N/A private static void perfectInterpolation(float[] in_offset, float in_end,
829N/A float[] startpitch, float pitchstep, float[] out, int[] out_offset,
829N/A int out_end) {
829N/A
829N/A float pitch = startpitch[0];
829N/A float ix = in_offset[0];
829N/A int ox = out_offset[0];
829N/A float ix_end = in_end;
829N/A int ox_end = out_end;
829N/A if (pitchstep == 0f) {
829N/A while (ix < ix_end && ox < ox_end) {
829N/A out[ox++] = getResamplerTestValue(ix);
829N/A ix += pitch;
829N/A }
829N/A } else {
829N/A while (ix < ix_end && ox < ox_end) {
829N/A out[ox++] = getResamplerTestValue(ix);
829N/A ix += pitch;
829N/A pitch += pitchstep;
829N/A }
829N/A }
829N/A in_offset[0] = ix;
829N/A out_offset[0] = ox;
829N/A startpitch[0] = pitch;
829N/A
829N/A }
829N/A
829N/A private static float testResampler(SoftAbstractResampler resampler, float p_pitch, float p_pitch2)
829N/A {
829N/A float[] testbuffer = new float[4096];
829N/A float[] testbuffer2 = new float[1024];
829N/A float[] testbuffer3 = new float[1024];
829N/A for (int i = 0; i < testbuffer.length; i++)
829N/A testbuffer[i] = getResamplerTestValue(i);
829N/A int pads = resampler.getPadding();
829N/A float pitchstep = (p_pitch2 - p_pitch)/1024f;
829N/A int[] out_offset2 = {0};
829N/A int[] out_offset3 = {0};
829N/A resampler.interpolate(testbuffer, new float[] {pads}, testbuffer.length - pads, new float[] {p_pitch}, pitchstep, testbuffer2, out_offset2, testbuffer2.length);
829N/A perfectInterpolation(new float[] {pads}, testbuffer.length - pads, new float[] {p_pitch}, pitchstep, testbuffer3, out_offset3, testbuffer3.length);
829N/A int out_off = out_offset2[0];
829N/A if(out_offset3[0] < out_off)
829N/A out_off = out_offset3[0];
829N/A float ac_error = 0;
829N/A int counter = 0;
829N/A for (int i = pads; i < out_off; i++) {
829N/A ac_error += Math.abs(testbuffer2[i] - testbuffer3[i]);
829N/A counter++;
829N/A }
829N/A return ac_error / ((float)counter);
829N/A }
829N/A
829N/A private static void fail(String error) throws Exception
829N/A {
829N/A throw new RuntimeException(error);
829N/A }
829N/A
829N/A public static void main(String[] args) throws Exception {
829N/A SoftCubicResampler resampler = new SoftCubicResampler();
829N/A float max = testResampler(resampler, 0.3f, 0.3f);
829N/A if(max > 0.005)
829N/A fail("Interpolation failed, error="+max);
829N/A max = testResampler(resampler, 0.3f, 0.01f);
829N/A if(max > 0.005)
829N/A fail("Interpolation failed, error="+max);
829N/A max = testResampler(resampler, 1.0f, 0.00f);
829N/A if(max > 0.005)
829N/A fail("Interpolation failed, error="+max);
829N/A }
829N/A}