--- /dev/null
+/*\r
+ * 11/19/04 1.0 moved to LGPL.\r
+ * \r
+ * 04/01/00 Fixes for running under build 23xx Microsoft JVM. mdm.\r
+ * \r
+ * 19/12/99 Performance improvements to compute_pcm_samples(). \r
+ * Mat McGowan. mdm@techie.com. \r
+ *\r
+ * 16/02/99 Java Conversion by E.B , javalayer@javazoom.net\r
+ *\r
+ * @(#) synthesis_filter.h 1.8, last edit: 6/15/94 16:52:00\r
+ * @(#) Copyright (C) 1993, 1994 Tobias Bading (bading@cs.tu-berlin.de)\r
+ * @(#) Berlin University of Technology\r
+ *\r
+ *-----------------------------------------------------------------------\r
+ * This program is free software; you can redistribute it and/or modify\r
+ * it under the terms of the GNU Library General Public License as published\r
+ * by the Free Software Foundation; either version 2 of the License, or\r
+ * (at your option) any later version.\r
+ *\r
+ * This program is distributed in the hope that it will be useful,\r
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of\r
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the\r
+ * GNU Library General Public License for more details.\r
+ *\r
+ * You should have received a copy of the GNU Library General Public\r
+ * License along with this program; if not, write to the Free Software\r
+ * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.\r
+ *----------------------------------------------------------------------\r
+ */\r
+\r
+/**\r
+ * A class for the synthesis filter bank. This class does a fast downsampling\r
+ * from 32, 44.1 or 48 kHz to 8 kHz, if ULAW is defined. Frequencies above 4 kHz\r
+ * are removed by ignoring higher subbands.\r
+ */\r
+\r
+\r
+final class SynthesisFilter {\r
+\r
+ \r
+ public int vidx;\r
+ \r
+ public float[] v1;\r
+ \r
+ public float[] v2;\r
+\r
+ \r
+ public float[] prev1;\r
+ \r
+ public float[] prev2;\r
+\r
+ // private float[] actual_v; // v1 or v2\r
+ \r
+ public int actual_write_pos; // 0-15\r
+ \r
+ private float[] samples; // 32 new subband samples\r
+ \r
+ public final int channel;\r
+ \r
+ public final float scalefactor;\r
+ \r
+ private float[] eq;\r
+\r
+ /**\r
+ * Quality value for controlling CPU usage/quality tradeoff.\r
+ */\r
+ /*\r
+ * private int quality;\r
+ * \r
+ * private int v_inc;\r
+ * \r
+ * \r
+ * \r
+ * public static final int HIGH_QUALITY = 1; public static final int\r
+ * MEDIUM_QUALITY = 2; public static final int LOW_QUALITY = 4;\r
+ */\r
+\r
+ /**\r
+ * Contructor. The scalefactor scales the calculated float pcm samples to\r
+ * short values (raw pcm samples are in [-1.0, 1.0], if no violations occur).\r
+ */\r
+ public SynthesisFilter(int channelnumber, float factor, float[] eq0) {\r
+\r
+ vidx = 1;\r
+ d16 = splitArray(d, 16);\r
+\r
+ v1 = new float[512];\r
+ v2 = new float[512];\r
+ prev1 = new float[512];\r
+ prev2 = new float[512];\r
+ samples = new float[32];\r
+ channel = channelnumber;\r
+ scalefactor = factor;\r
+ // setEQ(eq);\r
+ // setQuality(HIGH_QUALITY);\r
+\r
+ if (eq == null) {\r
+ eq = new float[32];\r
+ for (int i = 0; i < 32; i++)\r
+ eq[i] = 1.0f;\r
+ }\r
+ if (eq.length < 32) {\r
+ throw new IllegalArgumentException("eq0");\r
+ }\r
+\r
+ // reset();\r
+\r
+ for (int p = 0; p < 512; p++)\r
+ v1[p] = v2[p] = 0.0f;\r
+\r
+ for (int p2 = 0; p2 < 32; p2++)\r
+ samples[p2] = 0.0f;\r
+\r
+ // actual_v = v1;\r
+ actual_write_pos = 15;\r
+\r
+ }\r
+\r
+ /*\r
+ * private void setQuality(int quality0) { switch (quality0) { case\r
+ * HIGH_QUALITY: case MEDIUM_QUALITY: case LOW_QUALITY: v_inc = 16 * quality0;\r
+ * quality = quality0; break; default : throw new\r
+ * IllegalArgumentException("Unknown quality value"); } }\r
+ * \r
+ * public int getQuality() { return quality; }\r
+ */\r
+\r
+ /**\r
+ * Inject Sample.\r
+ */\r
+ public void input_sample( float sample, int subbandnumber) {\r
+ samples[subbandnumber] = eq[subbandnumber] * sample;\r
+ }\r
+\r
+ public void input_samples( float[] s) {\r
+ TERMINATE: for ( int i = 31; i >= 0; i--) {\r
+ samples[i] = s[i] * eq[i];\r
+ }\r
+ }\r
+\r
+ private void compute_new_v2_v1() {\r
+\r
+ float new_v0 = 0.0f;\r
+ float new_v1 = 0.0f;\r
+ float new_v2 = 0.0f;\r
+ float new_v3 = 0.0f;\r
+ float new_v4 = 0.0f;\r
+ float new_v5 = 0.0f;\r
+ float new_v6 = 0.0f;\r
+ float new_v7 = 0.0f;\r
+ float new_v8 = 0.0f;\r
+ float new_v9 = 0.0f;\r
+ float new_v10 = 0.0f;\r
+ float new_v11 = 0.0f;\r
+ float new_v12 = 0.0f;\r
+ float new_v13 = 0.0f;\r
+ float new_v14 = 0.0f;\r
+ float new_v15 = 0.0f;\r
+ float new_v16 = 0.0f;\r
+ float new_v17 = 0.0f;\r
+ float new_v18 = 0.0f;\r
+ float new_v19 = 0.0f;\r
+ float new_v20 = 0.0f;\r
+ float new_v21 = 0.0f;\r
+ float new_v22 = 0.0f;\r
+ float new_v23 = 0.0f;\r
+ float new_v24 = 0.0f;\r
+ float new_v25 = 0.0f;\r
+ float new_v26 = 0.0f;\r
+ float new_v27 = 0.0f;\r
+ float new_v28 = 0.0f;\r
+ float new_v29 = 0.0f;\r
+ float new_v30 = 0.0f;\r
+ float new_v31 = 0.0f;\r
+\r
+ // float[] new_v = new float[32]; // new V[0-15] and V[33-48] of Figure\r
+ // 3-A.2 in ISO DIS 11172-3\r
+ // float[] p = new float[16];\r
+ // float[] pp = new float[16];\r
+\r
+ // float[] s = samples;\r
+\r
+ float s0 = samples[0];\r
+ float s1 = samples[1];\r
+ float s2 = samples[2];\r
+ float s3 = samples[3];\r
+ float s4 = samples[4];\r
+ float s5 = samples[5];\r
+ float s6 = samples[6];\r
+ float s7 = samples[7];\r
+ float s8 = samples[8];\r
+ float s9 = samples[9];\r
+ float s10 = samples[10];\r
+ float s11 = samples[11];\r
+ float s12 = samples[12];\r
+ float s13 = samples[13];\r
+ float s14 = samples[14];\r
+ float s15 = samples[15];\r
+ float s16 = samples[16];\r
+ float s17 = samples[17];\r
+ float s18 = samples[18];\r
+ float s19 = samples[19];\r
+ float s20 = samples[20];\r
+ float s21 = samples[21];\r
+ float s22 = samples[22];\r
+ float s23 = samples[23];\r
+ float s24 = samples[24];\r
+ float s25 = samples[25];\r
+ float s26 = samples[26];\r
+ float s27 = samples[27];\r
+ float s28 = samples[28];\r
+ float s29 = samples[29];\r
+ float s30 = samples[30];\r
+ float s31 = samples[31];\r
+\r
+ float p0 = s0 + s31;\r
+ float p1 = s1 + s30;\r
+ float p2 = s2 + s29;\r
+ float p3 = s3 + s28;\r
+ float p4 = s4 + s27;\r
+ float p5 = s5 + s26;\r
+ float p6 = s6 + s25;\r
+ float p7 = s7 + s24;\r
+ float p8 = s8 + s23;\r
+ float p9 = s9 + s22;\r
+ float p10 = s10 + s21;\r
+ float p11 = s11 + s20;\r
+ float p12 = s12 + s19;\r
+ float p13 = s13 + s18;\r
+ float p14 = s14 + s17;\r
+ float p15 = s15 + s16;\r
+\r
+ float pp0 = p0 + p15;\r
+ float pp1 = p1 + p14;\r
+ float pp2 = p2 + p13;\r
+ float pp3 = p3 + p12;\r
+ float pp4 = p4 + p11;\r
+ float pp5 = p5 + p10;\r
+ float pp6 = p6 + p9;\r
+ float pp7 = p7 + p8;\r
+ float pp8 = (p0 - p15) * cos1_32;\r
+ float pp9 = (p1 - p14) * cos3_32;\r
+ float pp10 = (p2 - p13) * cos5_32;\r
+ float pp11 = (p3 - p12) * cos7_32;\r
+ float pp12 = (p4 - p11) * cos9_32;\r
+ float pp13 = (p5 - p10) * cos11_32;\r
+ float pp14 = (p6 - p9) * cos13_32;\r
+ float pp15 = (p7 - p8) * cos15_32;\r
+\r
+ p0 = pp0 + pp7;\r
+ p1 = pp1 + pp6;\r
+ p2 = pp2 + pp5;\r
+ p3 = pp3 + pp4;\r
+ p4 = (pp0 - pp7) * cos1_16;\r
+ p5 = (pp1 - pp6) * cos3_16;\r
+ p6 = (pp2 - pp5) * cos5_16;\r
+ p7 = (pp3 - pp4) * cos7_16;\r
+ p8 = pp8 + pp15;\r
+ p9 = pp9 + pp14;\r
+ p10 = pp10 + pp13;\r
+ p11 = pp11 + pp12;\r
+ p12 = (pp8 - pp15) * cos1_16;\r
+ p13 = (pp9 - pp14) * cos3_16;\r
+ p14 = (pp10 - pp13) * cos5_16;\r
+ p15 = (pp11 - pp12) * cos7_16;\r
+\r
+ pp0 = p0 + p3;\r
+ pp1 = p1 + p2;\r
+ pp2 = (p0 - p3) * cos1_8;\r
+ pp3 = (p1 - p2) * cos3_8;\r
+ pp4 = p4 + p7;\r
+ pp5 = p5 + p6;\r
+ pp6 = (p4 - p7) * cos1_8;\r
+ pp7 = (p5 - p6) * cos3_8;\r
+ pp8 = p8 + p11;\r
+ pp9 = p9 + p10;\r
+ pp10 = (p8 - p11) * cos1_8;\r
+ pp11 = (p9 - p10) * cos3_8;\r
+ pp12 = p12 + p15;\r
+ pp13 = p13 + p14;\r
+ pp14 = (p12 - p15) * cos1_8;\r
+ pp15 = (p13 - p14) * cos3_8;\r
+\r
+ p0 = pp0 + pp1;\r
+ p1 = (pp0 - pp1) * cos1_4;\r
+ p2 = pp2 + pp3;\r
+ p3 = (pp2 - pp3) * cos1_4;\r
+ p4 = pp4 + pp5;\r
+ p5 = (pp4 - pp5) * cos1_4;\r
+ p6 = pp6 + pp7;\r
+ p7 = (pp6 - pp7) * cos1_4;\r
+ p8 = pp8 + pp9;\r
+ p9 = (pp8 - pp9) * cos1_4;\r
+ p10 = pp10 + pp11;\r
+ p11 = (pp10 - pp11) * cos1_4;\r
+ p12 = pp12 + pp13;\r
+ p13 = (pp12 - pp13) * cos1_4;\r
+ p14 = pp14 + pp15;\r
+ p15 = (pp14 - pp15) * cos1_4;\r
+\r
+ // this is pretty insane coding\r
+ float tmp1;\r
+ new_v19/* 36-17 */= -(new_v4 = (new_v12 = p7) + p5) - p6;\r
+ new_v27/* 44-17 */= -p6 - p7 - p4;\r
+ new_v6 = (new_v10 = (new_v14 = p15) + p11) + p13;\r
+ new_v17/* 34-17 */= -(new_v2 = p15 + p13 + p9) - p14;\r
+ new_v21/* 38-17 */= (tmp1 = -p14 - p15 - p10 - p11) - p13;\r
+ new_v29/* 46-17 */= -p14 - p15 - p12 - p8;\r
+ new_v25/* 42-17 */= tmp1 - p12;\r
+ new_v31/* 48-17 */= -p0;\r
+ new_v0 = p1;\r
+ new_v23/* 40-17 */= -(new_v8 = p3) - p2;\r
+\r
+ p0 = (s0 - s31) * cos1_64;\r
+ p1 = (s1 - s30) * cos3_64;\r
+ p2 = (s2 - s29) * cos5_64;\r
+ p3 = (s3 - s28) * cos7_64;\r
+ p4 = (s4 - s27) * cos9_64;\r
+ p5 = (s5 - s26) * cos11_64;\r
+ p6 = (s6 - s25) * cos13_64;\r
+ p7 = (s7 - s24) * cos15_64;\r
+ p8 = (s8 - s23) * cos17_64;\r
+ p9 = (s9 - s22) * cos19_64;\r
+ p10 = (s10 - s21) * cos21_64;\r
+ p11 = (s11 - s20) * cos23_64;\r
+ p12 = (s12 - s19) * cos25_64;\r
+ p13 = (s13 - s18) * cos27_64;\r
+ p14 = (s14 - s17) * cos29_64;\r
+ p15 = (s15 - s16) * cos31_64;\r
+\r
+ pp0 = p0 + p15;\r
+ pp1 = p1 + p14;\r
+ pp2 = p2 + p13;\r
+ pp3 = p3 + p12;\r
+ pp4 = p4 + p11;\r
+ pp5 = p5 + p10;\r
+ pp6 = p6 + p9;\r
+ pp7 = p7 + p8;\r
+ pp8 = (p0 - p15) * cos1_32;\r
+ pp9 = (p1 - p14) * cos3_32;\r
+ pp10 = (p2 - p13) * cos5_32;\r
+ pp11 = (p3 - p12) * cos7_32;\r
+ pp12 = (p4 - p11) * cos9_32;\r
+ pp13 = (p5 - p10) * cos11_32;\r
+ pp14 = (p6 - p9) * cos13_32;\r
+ pp15 = (p7 - p8) * cos15_32;\r
+\r
+ p0 = pp0 + pp7;\r
+ p1 = pp1 + pp6;\r
+ p2 = pp2 + pp5;\r
+ p3 = pp3 + pp4;\r
+ p4 = (pp0 - pp7) * cos1_16;\r
+ p5 = (pp1 - pp6) * cos3_16;\r
+ p6 = (pp2 - pp5) * cos5_16;\r
+ p7 = (pp3 - pp4) * cos7_16;\r
+ p8 = pp8 + pp15;\r
+ p9 = pp9 + pp14;\r
+ p10 = pp10 + pp13;\r
+ p11 = pp11 + pp12;\r
+ p12 = (pp8 - pp15) * cos1_16;\r
+ p13 = (pp9 - pp14) * cos3_16;\r
+ p14 = (pp10 - pp13) * cos5_16;\r
+ p15 = (pp11 - pp12) * cos7_16;\r
+\r
+ pp0 = p0 + p3;\r
+ pp1 = p1 + p2;\r
+ pp2 = (p0 - p3) * cos1_8;\r
+ pp3 = (p1 - p2) * cos3_8;\r
+ pp4 = p4 + p7;\r
+ pp5 = p5 + p6;\r
+ pp6 = (p4 - p7) * cos1_8;\r
+ pp7 = (p5 - p6) * cos3_8;\r
+ pp8 = p8 + p11;\r
+ pp9 = p9 + p10;\r
+ pp10 = (p8 - p11) * cos1_8;\r
+ pp11 = (p9 - p10) * cos3_8;\r
+ pp12 = p12 + p15;\r
+ pp13 = p13 + p14;\r
+ pp14 = (p12 - p15) * cos1_8;\r
+ pp15 = (p13 - p14) * cos3_8;\r
+\r
+ p0 = pp0 + pp1;\r
+ p1 = (pp0 - pp1) * cos1_4;\r
+ p2 = pp2 + pp3;\r
+ p3 = (pp2 - pp3) * cos1_4;\r
+ p4 = pp4 + pp5;\r
+ p5 = (pp4 - pp5) * cos1_4;\r
+ p6 = pp6 + pp7;\r
+ p7 = (pp6 - pp7) * cos1_4;\r
+ p8 = pp8 + pp9;\r
+ p9 = (pp8 - pp9) * cos1_4;\r
+ p10 = pp10 + pp11;\r
+ p11 = (pp10 - pp11) * cos1_4;\r
+ p12 = pp12 + pp13;\r
+ p13 = (pp12 - pp13) * cos1_4;\r
+ p14 = pp14 + pp15;\r
+ p15 = (pp14 - pp15) * cos1_4;\r
+\r
+ // manually doing something that a compiler should handle sucks\r
+ // coding like this is hard to read\r
+ float tmp2;\r
+ new_v5 = (new_v11 = (new_v13 = (new_v15 = p15) + p7) + p11) + p5 + p13;\r
+ new_v7 = (new_v9 = p15 + p11 + p3) + p13;\r
+ new_v16/* 33-17 */= -(new_v1 = (tmp1 = p13 + p15 + p9) + p1) - p14;\r
+ new_v18/* 35-17 */= -(new_v3 = tmp1 + p5 + p7) - p6 - p14;\r
+\r
+ new_v22/* 39-17 */= (tmp1 = -p10 - p11 - p14 - p15) - p13 - p2 - p3;\r
+ new_v20/* 37-17 */= tmp1 - p13 - p5 - p6 - p7;\r
+ new_v24/* 41-17 */= tmp1 - p12 - p2 - p3;\r
+ new_v26/* 43-17 */= tmp1 - p12 - (tmp2 = p4 + p6 + p7);\r
+ new_v30/* 47-17 */= (tmp1 = -p8 - p12 - p14 - p15) - p0;\r
+ new_v28/* 45-17 */= tmp1 - tmp2;\r
+\r
+ // insert V[0-15] (== new_v[0-15]) into actual v:\r
+ // float[] x2 = actual_v + actual_write_pos;\r
+ // float dest[] = actual_v; v2\r
+\r
+ int pos = actual_write_pos;\r
+\r
+ v2[0 + pos] = new_v0;\r
+ v2[16 + pos] = new_v1;\r
+ v2[32 + pos] = new_v2;\r
+ v2[48 + pos] = new_v3;\r
+ v2[64 + pos] = new_v4;\r
+ v2[80 + pos] = new_v5;\r
+ v2[96 + pos] = new_v6;\r
+ v2[112 + pos] = new_v7;\r
+ v2[128 + pos] = new_v8;\r
+ v2[144 + pos] = new_v9;\r
+ v2[160 + pos] = new_v10;\r
+ v2[176 + pos] = new_v11;\r
+ v2[192 + pos] = new_v12;\r
+ v2[208 + pos] = new_v13;\r
+ v2[224 + pos] = new_v14;\r
+ v2[240 + pos] = new_v15;\r
+\r
+ // V[16] is always 0.0:\r
+ v2[256 + pos] = 0.0f;\r
+\r
+ // insert V[17-31] (== -new_v[15-1]) into actual v:\r
+ v2[272 + pos] = -new_v15;\r
+ v2[288 + pos] = -new_v14;\r
+ v2[304 + pos] = -new_v13;\r
+ v2[320 + pos] = -new_v12;\r
+ v2[336 + pos] = -new_v11;\r
+ v2[352 + pos] = -new_v10;\r
+ v2[368 + pos] = -new_v9;\r
+ v2[384 + pos] = -new_v8;\r
+ v2[400 + pos] = -new_v7;\r
+ v2[416 + pos] = -new_v6;\r
+ v2[432 + pos] = -new_v5;\r
+ v2[448 + pos] = -new_v4;\r
+ v2[464 + pos] = -new_v3;\r
+ v2[480 + pos] = -new_v2;\r
+ v2[496 + pos] = -new_v1;\r
+\r
+ // insert V[32] (== -new_v[0]) into other v:\r
+ // dest = (actual_v == v1) ? v2 : v1;\r
+\r
+ v1[0 + pos] = -new_v0;\r
+ // insert V[33-48] (== new_v[16-31]) into other v:\r
+ v1[16 + pos] = new_v16;\r
+ v1[32 + pos] = new_v17;\r
+ v1[48 + pos] = new_v18;\r
+ v1[64 + pos] = new_v19;\r
+ v1[80 + pos] = new_v20;\r
+ v1[96 + pos] = new_v21;\r
+ v1[112 + pos] = new_v22;\r
+ v1[128 + pos] = new_v23;\r
+ v1[144 + pos] = new_v24;\r
+ v1[160 + pos] = new_v25;\r
+ v1[176 + pos] = new_v26;\r
+ v1[192 + pos] = new_v27;\r
+ v1[208 + pos] = new_v28;\r
+ v1[224 + pos] = new_v29;\r
+ v1[240 + pos] = new_v30;\r
+ v1[256 + pos] = new_v31;\r
+\r
+ // insert V[49-63] (== new_v[30-16]) into other v:\r
+ v1[272 + pos] = new_v30;\r
+ v1[288 + pos] = new_v29;\r
+ v1[304 + pos] = new_v28;\r
+ v1[320 + pos] = new_v27;\r
+ v1[336 + pos] = new_v26;\r
+ v1[352 + pos] = new_v25;\r
+ v1[368 + pos] = new_v24;\r
+ v1[384 + pos] = new_v23;\r
+ v1[400 + pos] = new_v22;\r
+ v1[416 + pos] = new_v21;\r
+ v1[432 + pos] = new_v20;\r
+ v1[448 + pos] = new_v19;\r
+ v1[464 + pos] = new_v18;\r
+ v1[480 + pos] = new_v17;\r
+ v1[496 + pos] = new_v16;\r
+\r
+ /*\r
+ * setup PREV\r
+ */\r
+\r
+ prev2[0 + pos] = new_v0;\r
+ prev2[16 + pos] = new_v1;\r
+ prev2[32 + pos] = new_v2;\r
+ prev2[48 + pos] = new_v3;\r
+ prev2[64 + pos] = new_v4;\r
+ prev2[80 + pos] = new_v5;\r
+ prev2[96 + pos] = new_v6;\r
+ prev2[112 + pos] = new_v7;\r
+ prev2[128 + pos] = new_v8;\r
+ prev2[144 + pos] = new_v9;\r
+ prev2[160 + pos] = new_v10;\r
+ prev2[176 + pos] = new_v11;\r
+ prev2[192 + pos] = new_v12;\r
+ prev2[208 + pos] = new_v13;\r
+ prev2[224 + pos] = new_v14;\r
+ prev2[240 + pos] = new_v15;\r
+\r
+ // V[16] is always 0.0:\r
+ prev2[256 + pos] = 0.0f;\r
+\r
+ // insert V[17-31] (== -new_v[15-1]) into actual v:\r
+ prev2[272 + pos] = -new_v15;\r
+ prev2[288 + pos] = -new_v14;\r
+ prev2[304 + pos] = -new_v13;\r
+ prev2[320 + pos] = -new_v12;\r
+ prev2[336 + pos] = -new_v11;\r
+ prev2[352 + pos] = -new_v10;\r
+ prev2[368 + pos] = -new_v9;\r
+ prev2[384 + pos] = -new_v8;\r
+ prev2[400 + pos] = -new_v7;\r
+ prev2[416 + pos] = -new_v6;\r
+ prev2[432 + pos] = -new_v5;\r
+ prev2[448 + pos] = -new_v4;\r
+ prev2[464 + pos] = -new_v3;\r
+ prev2[480 + pos] = -new_v2;\r
+ prev2[496 + pos] = -new_v1;\r
+\r
+ // insert V[32] (== -new_v[0]) into other v:\r
+ // dest = (actual_v == v1) ? v2 : v1;\r
+\r
+ prev1[0 + pos] = -new_v0;\r
+ // insert V[33-48] (== new_v[16-31]) into other v:\r
+ prev1[16 + pos] = new_v16;\r
+ prev1[32 + pos] = new_v17;\r
+ prev1[48 + pos] = new_v18;\r
+ prev1[64 + pos] = new_v19;\r
+ prev1[80 + pos] = new_v20;\r
+ prev1[96 + pos] = new_v21;\r
+ prev1[112 + pos] = new_v22;\r
+ prev1[128 + pos] = new_v23;\r
+ prev1[144 + pos] = new_v24;\r
+ prev1[160 + pos] = new_v25;\r
+ prev1[176 + pos] = new_v26;\r
+ prev1[192 + pos] = new_v27;\r
+ prev1[208 + pos] = new_v28;\r
+ prev1[224 + pos] = new_v29;\r
+ prev1[240 + pos] = new_v30;\r
+ prev1[256 + pos] = new_v31;\r
+\r
+ // insert V[49-63] (== new_v[30-16]) into other v:\r
+ prev1[272 + pos] = new_v30;\r
+ prev1[288 + pos] = new_v29;\r
+ prev1[304 + pos] = new_v28;\r
+ prev1[320 + pos] = new_v27;\r
+ prev1[336 + pos] = new_v26;\r
+ prev1[352 + pos] = new_v25;\r
+ prev1[368 + pos] = new_v24;\r
+ prev1[384 + pos] = new_v23;\r
+ prev1[400 + pos] = new_v22;\r
+ prev1[416 + pos] = new_v21;\r
+ prev1[432 + pos] = new_v20;\r
+ prev1[448 + pos] = new_v19;\r
+ prev1[464 + pos] = new_v18;\r
+ prev1[480 + pos] = new_v17;\r
+ prev1[496 + pos] = new_v16;\r
+ }\r
+\r
+ private void compute_new_v1_v2() {\r
+\r
+ float new_v0 = 0.0f;\r
+ float new_v1 = 0.0f;\r
+ float new_v2 = 0.0f;\r
+ float new_v3 = 0.0f;\r
+ float new_v4 = 0.0f;\r
+ float new_v5 = 0.0f;\r
+ float new_v6 = 0.0f;\r
+ float new_v7 = 0.0f;\r
+ float new_v8 = 0.0f;\r
+ float new_v9 = 0.0f;\r
+ float new_v10 = 0.0f;\r
+ float new_v11 = 0.0f;\r
+ float new_v12 = 0.0f;\r
+ float new_v13 = 0.0f;\r
+ float new_v14 = 0.0f;\r
+ float new_v15 = 0.0f;\r
+ float new_v16 = 0.0f;\r
+ float new_v17 = 0.0f;\r
+ float new_v18 = 0.0f;\r
+ float new_v19 = 0.0f;\r
+ float new_v20 = 0.0f;\r
+ float new_v21 = 0.0f;\r
+ float new_v22 = 0.0f;\r
+ float new_v23 = 0.0f;\r
+ float new_v24 = 0.0f;\r
+ float new_v25 = 0.0f;\r
+ float new_v26 = 0.0f;\r
+ float new_v27 = 0.0f;\r
+ float new_v28 = 0.0f;\r
+ float new_v29 = 0.0f;\r
+ float new_v30 = 0.0f;\r
+ float new_v31 = 0.0f;\r
+\r
+ // float[] new_v = new float[32]; // new V[0-15] and V[33-48] of Figure\r
+ // 3-A.2 in ISO DIS 11172-3\r
+ // float[] p = new float[16];\r
+ // float[] pp = new float[16];\r
+\r
+ // float[] s = samples;\r
+\r
+ float s0 = samples[0];\r
+ float s1 = samples[1];\r
+ float s2 = samples[2];\r
+ float s3 = samples[3];\r
+ float s4 = samples[4];\r
+ float s5 = samples[5];\r
+ float s6 = samples[6];\r
+ float s7 = samples[7];\r
+ float s8 = samples[8];\r
+ float s9 = samples[9];\r
+ float s10 = samples[10];\r
+ float s11 = samples[11];\r
+ float s12 = samples[12];\r
+ float s13 = samples[13];\r
+ float s14 = samples[14];\r
+ float s15 = samples[15];\r
+ float s16 = samples[16];\r
+ float s17 = samples[17];\r
+ float s18 = samples[18];\r
+ float s19 = samples[19];\r
+ float s20 = samples[20];\r
+ float s21 = samples[21];\r
+ float s22 = samples[22];\r
+ float s23 = samples[23];\r
+ float s24 = samples[24];\r
+ float s25 = samples[25];\r
+ float s26 = samples[26];\r
+ float s27 = samples[27];\r
+ float s28 = samples[28];\r
+ float s29 = samples[29];\r
+ float s30 = samples[30];\r
+ float s31 = samples[31];\r
+\r
+ float p0 = s0 + s31;\r
+ float p1 = s1 + s30;\r
+ float p2 = s2 + s29;\r
+ float p3 = s3 + s28;\r
+ float p4 = s4 + s27;\r
+ float p5 = s5 + s26;\r
+ float p6 = s6 + s25;\r
+ float p7 = s7 + s24;\r
+ float p8 = s8 + s23;\r
+ float p9 = s9 + s22;\r
+ float p10 = s10 + s21;\r
+ float p11 = s11 + s20;\r
+ float p12 = s12 + s19;\r
+ float p13 = s13 + s18;\r
+ float p14 = s14 + s17;\r
+ float p15 = s15 + s16;\r
+\r
+ float pp0 = p0 + p15;\r
+ float pp1 = p1 + p14;\r
+ float pp2 = p2 + p13;\r
+ float pp3 = p3 + p12;\r
+ float pp4 = p4 + p11;\r
+ float pp5 = p5 + p10;\r
+ float pp6 = p6 + p9;\r
+ float pp7 = p7 + p8;\r
+ float pp8 = (p0 - p15) * cos1_32;\r
+ float pp9 = (p1 - p14) * cos3_32;\r
+ float pp10 = (p2 - p13) * cos5_32;\r
+ float pp11 = (p3 - p12) * cos7_32;\r
+ float pp12 = (p4 - p11) * cos9_32;\r
+ float pp13 = (p5 - p10) * cos11_32;\r
+ float pp14 = (p6 - p9) * cos13_32;\r
+ float pp15 = (p7 - p8) * cos15_32;\r
+\r
+ p0 = pp0 + pp7;\r
+ p1 = pp1 + pp6;\r
+ p2 = pp2 + pp5;\r
+ p3 = pp3 + pp4;\r
+ p4 = (pp0 - pp7) * cos1_16;\r
+ p5 = (pp1 - pp6) * cos3_16;\r
+ p6 = (pp2 - pp5) * cos5_16;\r
+ p7 = (pp3 - pp4) * cos7_16;\r
+ p8 = pp8 + pp15;\r
+ p9 = pp9 + pp14;\r
+ p10 = pp10 + pp13;\r
+ p11 = pp11 + pp12;\r
+ p12 = (pp8 - pp15) * cos1_16;\r
+ p13 = (pp9 - pp14) * cos3_16;\r
+ p14 = (pp10 - pp13) * cos5_16;\r
+ p15 = (pp11 - pp12) * cos7_16;\r
+\r
+ pp0 = p0 + p3;\r
+ pp1 = p1 + p2;\r
+ pp2 = (p0 - p3) * cos1_8;\r
+ pp3 = (p1 - p2) * cos3_8;\r
+ pp4 = p4 + p7;\r
+ pp5 = p5 + p6;\r
+ pp6 = (p4 - p7) * cos1_8;\r
+ pp7 = (p5 - p6) * cos3_8;\r
+ pp8 = p8 + p11;\r
+ pp9 = p9 + p10;\r
+ pp10 = (p8 - p11) * cos1_8;\r
+ pp11 = (p9 - p10) * cos3_8;\r
+ pp12 = p12 + p15;\r
+ pp13 = p13 + p14;\r
+ pp14 = (p12 - p15) * cos1_8;\r
+ pp15 = (p13 - p14) * cos3_8;\r
+\r
+ p0 = pp0 + pp1;\r
+ p1 = (pp0 - pp1) * cos1_4;\r
+ p2 = pp2 + pp3;\r
+ p3 = (pp2 - pp3) * cos1_4;\r
+ p4 = pp4 + pp5;\r
+ p5 = (pp4 - pp5) * cos1_4;\r
+ p6 = pp6 + pp7;\r
+ p7 = (pp6 - pp7) * cos1_4;\r
+ p8 = pp8 + pp9;\r
+ p9 = (pp8 - pp9) * cos1_4;\r
+ p10 = pp10 + pp11;\r
+ p11 = (pp10 - pp11) * cos1_4;\r
+ p12 = pp12 + pp13;\r
+ p13 = (pp12 - pp13) * cos1_4;\r
+ p14 = pp14 + pp15;\r
+ p15 = (pp14 - pp15) * cos1_4;\r
+\r
+ // this is pretty insane coding\r
+ float tmp1;\r
+ new_v19/* 36-17 */= -(new_v4 = (new_v12 = p7) + p5) - p6;\r
+ new_v27/* 44-17 */= -p6 - p7 - p4;\r
+ new_v6 = (new_v10 = (new_v14 = p15) + p11) + p13;\r
+ new_v17/* 34-17 */= -(new_v2 = p15 + p13 + p9) - p14;\r
+ new_v21/* 38-17 */= (tmp1 = -p14 - p15 - p10 - p11) - p13;\r
+ new_v29/* 46-17 */= -p14 - p15 - p12 - p8;\r
+ new_v25/* 42-17 */= tmp1 - p12;\r
+ new_v31/* 48-17 */= -p0;\r
+ new_v0 = p1;\r
+ new_v23/* 40-17 */= -(new_v8 = p3) - p2;\r
+\r
+ p0 = (s0 - s31) * cos1_64;\r
+ p1 = (s1 - s30) * cos3_64;\r
+ p2 = (s2 - s29) * cos5_64;\r
+ p3 = (s3 - s28) * cos7_64;\r
+ p4 = (s4 - s27) * cos9_64;\r
+ p5 = (s5 - s26) * cos11_64;\r
+ p6 = (s6 - s25) * cos13_64;\r
+ p7 = (s7 - s24) * cos15_64;\r
+ p8 = (s8 - s23) * cos17_64;\r
+ p9 = (s9 - s22) * cos19_64;\r
+ p10 = (s10 - s21) * cos21_64;\r
+ p11 = (s11 - s20) * cos23_64;\r
+ p12 = (s12 - s19) * cos25_64;\r
+ p13 = (s13 - s18) * cos27_64;\r
+ p14 = (s14 - s17) * cos29_64;\r
+ p15 = (s15 - s16) * cos31_64;\r
+\r
+ pp0 = p0 + p15;\r
+ pp1 = p1 + p14;\r
+ pp2 = p2 + p13;\r
+ pp3 = p3 + p12;\r
+ pp4 = p4 + p11;\r
+ pp5 = p5 + p10;\r
+ pp6 = p6 + p9;\r
+ pp7 = p7 + p8;\r
+ pp8 = (p0 - p15) * cos1_32;\r
+ pp9 = (p1 - p14) * cos3_32;\r
+ pp10 = (p2 - p13) * cos5_32;\r
+ pp11 = (p3 - p12) * cos7_32;\r
+ pp12 = (p4 - p11) * cos9_32;\r
+ pp13 = (p5 - p10) * cos11_32;\r
+ pp14 = (p6 - p9) * cos13_32;\r
+ pp15 = (p7 - p8) * cos15_32;\r
+\r
+ p0 = pp0 + pp7;\r
+ p1 = pp1 + pp6;\r
+ p2 = pp2 + pp5;\r
+ p3 = pp3 + pp4;\r
+ p4 = (pp0 - pp7) * cos1_16;\r
+ p5 = (pp1 - pp6) * cos3_16;\r
+ p6 = (pp2 - pp5) * cos5_16;\r
+ p7 = (pp3 - pp4) * cos7_16;\r
+ p8 = pp8 + pp15;\r
+ p9 = pp9 + pp14;\r
+ p10 = pp10 + pp13;\r
+ p11 = pp11 + pp12;\r
+ p12 = (pp8 - pp15) * cos1_16;\r
+ p13 = (pp9 - pp14) * cos3_16;\r
+ p14 = (pp10 - pp13) * cos5_16;\r
+ p15 = (pp11 - pp12) * cos7_16;\r
+\r
+ pp0 = p0 + p3;\r
+ pp1 = p1 + p2;\r
+ pp2 = (p0 - p3) * cos1_8;\r
+ pp3 = (p1 - p2) * cos3_8;\r
+ pp4 = p4 + p7;\r
+ pp5 = p5 + p6;\r
+ pp6 = (p4 - p7) * cos1_8;\r
+ pp7 = (p5 - p6) * cos3_8;\r
+ pp8 = p8 + p11;\r
+ pp9 = p9 + p10;\r
+ pp10 = (p8 - p11) * cos1_8;\r
+ pp11 = (p9 - p10) * cos3_8;\r
+ pp12 = p12 + p15;\r
+ pp13 = p13 + p14;\r
+ pp14 = (p12 - p15) * cos1_8;\r
+ pp15 = (p13 - p14) * cos3_8;\r
+\r
+ p0 = pp0 + pp1;\r
+ p1 = (pp0 - pp1) * cos1_4;\r
+ p2 = pp2 + pp3;\r
+ p3 = (pp2 - pp3) * cos1_4;\r
+ p4 = pp4 + pp5;\r
+ p5 = (pp4 - pp5) * cos1_4;\r
+ p6 = pp6 + pp7;\r
+ p7 = (pp6 - pp7) * cos1_4;\r
+ p8 = pp8 + pp9;\r
+ p9 = (pp8 - pp9) * cos1_4;\r
+ p10 = pp10 + pp11;\r
+ p11 = (pp10 - pp11) * cos1_4;\r
+ p12 = pp12 + pp13;\r
+ p13 = (pp12 - pp13) * cos1_4;\r
+ p14 = pp14 + pp15;\r
+ p15 = (pp14 - pp15) * cos1_4;\r
+\r
+ // manually doing something that a compiler should handle sucks\r
+ // coding like this is hard to read\r
+ float tmp2;\r
+ new_v5 = (new_v11 = (new_v13 = (new_v15 = p15) + p7) + p11) + p5 + p13;\r
+ new_v7 = (new_v9 = p15 + p11 + p3) + p13;\r
+ new_v16/* 33-17 */= -(new_v1 = (tmp1 = p13 + p15 + p9) + p1) - p14;\r
+ new_v18/* 35-17 */= -(new_v3 = tmp1 + p5 + p7) - p6 - p14;\r
+\r
+ new_v22/* 39-17 */= (tmp1 = -p10 - p11 - p14 - p15) - p13 - p2 - p3;\r
+ new_v20/* 37-17 */= tmp1 - p13 - p5 - p6 - p7;\r
+ new_v24/* 41-17 */= tmp1 - p12 - p2 - p3;\r
+ new_v26/* 43-17 */= tmp1 - p12 - (tmp2 = p4 + p6 + p7);\r
+ new_v30/* 47-17 */= (tmp1 = -p8 - p12 - p14 - p15) - p0;\r
+ new_v28/* 45-17 */= tmp1 - tmp2;\r
+\r
+ // insert V[0-15] (== new_v[0-15]) into actual v:\r
+ // float[] x2 = actual_v + actual_write_pos;\r
+ // float dest[] = actual_v; actual_v=v1;\r
+\r
+ int pos = actual_write_pos;\r
+\r
+ v1[0 + pos] = new_v0;\r
+ v1[16 + pos] = new_v1;\r
+ v1[32 + pos] = new_v2;\r
+ v1[48 + pos] = new_v3;\r
+ v1[64 + pos] = new_v4;\r
+ v1[80 + pos] = new_v5;\r
+ v1[96 + pos] = new_v6;\r
+ v1[112 + pos] = new_v7;\r
+ v1[128 + pos] = new_v8;\r
+ v1[144 + pos] = new_v9;\r
+ v1[160 + pos] = new_v10;\r
+ v1[176 + pos] = new_v11;\r
+ v1[192 + pos] = new_v12;\r
+ v1[208 + pos] = new_v13;\r
+ v1[224 + pos] = new_v14;\r
+ v1[240 + pos] = new_v15;\r
+\r
+ // V[16] is always 0.0:\r
+ v1[256 + pos] = 0.0f;\r
+\r
+ // insert V[17-31] (== -new_v[15-1]) into actual v:\r
+ v1[272 + pos] = -new_v15;\r
+ v1[288 + pos] = -new_v14;\r
+ v1[304 + pos] = -new_v13;\r
+ v1[320 + pos] = -new_v12;\r
+ v1[336 + pos] = -new_v11;\r
+ v1[352 + pos] = -new_v10;\r
+ v1[368 + pos] = -new_v9;\r
+ v1[384 + pos] = -new_v8;\r
+ v1[400 + pos] = -new_v7;\r
+ v1[416 + pos] = -new_v6;\r
+ v1[432 + pos] = -new_v5;\r
+ v1[448 + pos] = -new_v4;\r
+ v1[464 + pos] = -new_v3;\r
+ v1[480 + pos] = -new_v2;\r
+ v1[496 + pos] = -new_v1;\r
+\r
+ // insert V[32] (== -new_v[0]) into other v:\r
+ // dest = (actual_v == v1) ? v2 : v1;\r
+\r
+ v2[0 + pos] = -new_v0;\r
+ // insert V[33-48] (== new_v[16-31]) into other v:\r
+ v2[16 + pos] = new_v16;\r
+ v2[32 + pos] = new_v17;\r
+ v2[48 + pos] = new_v18;\r
+ v2[64 + pos] = new_v19;\r
+ v2[80 + pos] = new_v20;\r
+ v2[96 + pos] = new_v21;\r
+ v2[112 + pos] = new_v22;\r
+ v2[128 + pos] = new_v23;\r
+ v2[144 + pos] = new_v24;\r
+ v2[160 + pos] = new_v25;\r
+ v2[176 + pos] = new_v26;\r
+ v2[192 + pos] = new_v27;\r
+ v2[208 + pos] = new_v28;\r
+ v2[224 + pos] = new_v29;\r
+ v2[240 + pos] = new_v30;\r
+ v2[256 + pos] = new_v31;\r
+\r
+ // insert V[49-63] (== new_v[30-16]) into other v:\r
+ v2[272 + pos] = new_v30;\r
+ v2[288 + pos] = new_v29;\r
+ v2[304 + pos] = new_v28;\r
+ v2[320 + pos] = new_v27;\r
+ v2[336 + pos] = new_v26;\r
+ v2[352 + pos] = new_v25;\r
+ v2[368 + pos] = new_v24;\r
+ v2[384 + pos] = new_v23;\r
+ v2[400 + pos] = new_v22;\r
+ v2[416 + pos] = new_v21;\r
+ v2[432 + pos] = new_v20;\r
+ v2[448 + pos] = new_v19;\r
+ v2[464 + pos] = new_v18;\r
+ v2[480 + pos] = new_v17;\r
+ v2[496 + pos] = new_v16;\r
+\r
+ /*\r
+ * setup PREV\r
+ */\r
+\r
+ prev1[0 + pos] = new_v0;\r
+ prev1[16 + pos] = new_v1;\r
+ prev1[32 + pos] = new_v2;\r
+ prev1[48 + pos] = new_v3;\r
+ prev1[64 + pos] = new_v4;\r
+ prev1[80 + pos] = new_v5;\r
+ prev1[96 + pos] = new_v6;\r
+ prev1[112 + pos] = new_v7;\r
+ prev1[128 + pos] = new_v8;\r
+ prev1[144 + pos] = new_v9;\r
+ prev1[160 + pos] = new_v10;\r
+ prev1[176 + pos] = new_v11;\r
+ prev1[192 + pos] = new_v12;\r
+ prev1[208 + pos] = new_v13;\r
+ prev1[224 + pos] = new_v14;\r
+ prev1[240 + pos] = new_v15;\r
+\r
+ // V[16] is always 0.0:\r
+ prev1[256 + pos] = 0.0f;\r
+\r
+ // insert V[17-31] (== -new_v[15-1]) into actual v:\r
+ prev1[272 + pos] = -new_v15;\r
+ prev1[288 + pos] = -new_v14;\r
+ prev1[304 + pos] = -new_v13;\r
+ prev1[320 + pos] = -new_v12;\r
+ prev1[336 + pos] = -new_v11;\r
+ prev1[352 + pos] = -new_v10;\r
+ prev1[368 + pos] = -new_v9;\r
+ prev1[384 + pos] = -new_v8;\r
+ prev1[400 + pos] = -new_v7;\r
+ prev1[416 + pos] = -new_v6;\r
+ prev1[432 + pos] = -new_v5;\r
+ prev1[448 + pos] = -new_v4;\r
+ prev1[464 + pos] = -new_v3;\r
+ prev1[480 + pos] = -new_v2;\r
+ prev1[496 + pos] = -new_v1;\r
+\r
+ // insert V[32] (== -new_v[0]) into other v:\r
+ // dest = (actual_v == v1) ? v2 : v1;\r
+\r
+ prev2[0 + pos] = -new_v0;\r
+ // insert V[33-48] (== new_v[16-31]) into other v:\r
+ prev2[16 + pos] = new_v16;\r
+ prev2[32 + pos] = new_v17;\r
+ prev2[48 + pos] = new_v18;\r
+ prev2[64 + pos] = new_v19;\r
+ prev2[80 + pos] = new_v20;\r
+ prev2[96 + pos] = new_v21;\r
+ prev2[112 + pos] = new_v22;\r
+ prev2[128 + pos] = new_v23;\r
+ prev2[144 + pos] = new_v24;\r
+ prev2[160 + pos] = new_v25;\r
+ prev2[176 + pos] = new_v26;\r
+ prev2[192 + pos] = new_v27;\r
+ prev2[208 + pos] = new_v28;\r
+ prev2[224 + pos] = new_v29;\r
+ prev2[240 + pos] = new_v30;\r
+ prev2[256 + pos] = new_v31;\r
+\r
+ // insert V[49-63] (== new_v[30-16]) into other v:\r
+ prev2[272 + pos] = new_v30;\r
+ prev2[288 + pos] = new_v29;\r
+ prev2[304 + pos] = new_v28;\r
+ prev2[320 + pos] = new_v27;\r
+ prev2[336 + pos] = new_v26;\r
+ prev2[352 + pos] = new_v25;\r
+ prev2[368 + pos] = new_v24;\r
+ prev2[384 + pos] = new_v23;\r
+ prev2[400 + pos] = new_v22;\r
+ prev2[416 + pos] = new_v21;\r
+ prev2[432 + pos] = new_v20;\r
+ prev2[448 + pos] = new_v19;\r
+ prev2[464 + pos] = new_v18;\r
+ prev2[480 + pos] = new_v17;\r
+ prev2[496 + pos] = new_v16;\r
+ }\r
+\r
+ /**\r
+ * Compute PCM Samples.\r
+ */\r
+\r
+ \r
+ private float[] _tmpOut = new float[32];\r
+\r
+ private void compute_pcm_samples0() {\r
+\r
+ if (vidx == 1) {\r
+ // final float[] vp = actual_v;\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ float pcm_sample;\r
+ // final float[] dp = d16[i];\r
+ pcm_sample =\r
+ (float) (((v1[0 + dvp] * d16[i][0]) + (v1[15 + dvp] * d16[i][1])\r
+ + (v1[14 + dvp] * d16[i][2]) + (v1[13 + dvp] * d16[i][3])\r
+ + (v1[12 + dvp] * d16[i][4]) + (v1[11 + dvp] * d16[i][5])\r
+ + (v1[10 + dvp] * d16[i][6]) + (v1[9 + dvp] * d16[i][7])\r
+ + (v1[8 + dvp] * d16[i][8]) + (v1[7 + dvp] * d16[i][9])\r
+ + (v1[6 + dvp] * d16[i][10]) + (v1[5 + dvp] * d16[i][11])\r
+ + (v1[4 + dvp] * d16[i][12]) + (v1[3 + dvp] * d16[i][13])\r
+ + (v1[2 + dvp] * d16[i][14]) + (v1[1 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ } else {\r
+ // final float[] vp = actual_v;\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ float pcm_sample;\r
+ // final float[] dp = d16[i];\r
+ pcm_sample =\r
+ (float) (((v2[0 + dvp] * d16[i][0]) + (v2[15 + dvp] * d16[i][1])\r
+ + (v2[14 + dvp] * d16[i][2]) + (v2[13 + dvp] * d16[i][3])\r
+ + (v2[12 + dvp] * d16[i][4]) + (v2[11 + dvp] * d16[i][5])\r
+ + (v2[10 + dvp] * d16[i][6]) + (v2[9 + dvp] * d16[i][7])\r
+ + (v2[8 + dvp] * d16[i][8]) + (v2[7 + dvp] * d16[i][9])\r
+ + (v2[6 + dvp] * d16[i][10]) + (v2[5 + dvp] * d16[i][11])\r
+ + (v2[4 + dvp] * d16[i][12]) + (v2[3 + dvp] * d16[i][13])\r
+ + (v2[2 + dvp] * d16[i][14]) + (v2[1 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ }\r
+\r
+ }\r
+\r
+ private void compute_pcm_samples1() {\r
+\r
+ if (vidx == 1) {\r
+ // final float[] vp = actual_v;\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v1[1 + dvp] * d16[i][0]) + (v1[0 + dvp] * d16[i][1])\r
+ + (v1[15 + dvp] * d16[i][2]) + (v1[14 + dvp] * d16[i][3])\r
+ + (v1[13 + dvp] * d16[i][4]) + (v1[12 + dvp] * d16[i][5])\r
+ + (v1[11 + dvp] * d16[i][6]) + (v1[10 + dvp] * d16[i][7])\r
+ + (v1[9 + dvp] * d16[i][8]) + (v1[8 + dvp] * d16[i][9])\r
+ + (v1[7 + dvp] * d16[i][10]) + (v1[6 + dvp] * d16[i][11])\r
+ + (v1[5 + dvp] * d16[i][12]) + (v1[4 + dvp] * d16[i][13])\r
+ + (v1[3 + dvp] * d16[i][14]) + (v1[2 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ } else {\r
+ // final float[] vp = actual_v;\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v2[1 + dvp] * d16[i][0]) + (v2[0 + dvp] * d16[i][1])\r
+ + (v2[15 + dvp] * d16[i][2]) + (v2[14 + dvp] * d16[i][3])\r
+ + (v2[13 + dvp] * d16[i][4]) + (v2[12 + dvp] * d16[i][5])\r
+ + (v2[11 + dvp] * d16[i][6]) + (v2[10 + dvp] * d16[i][7])\r
+ + (v2[9 + dvp] * d16[i][8]) + (v2[8 + dvp] * d16[i][9])\r
+ + (v2[7 + dvp] * d16[i][10]) + (v2[6 + dvp] * d16[i][11])\r
+ + (v2[5 + dvp] * d16[i][12]) + (v2[4 + dvp] * d16[i][13])\r
+ + (v2[3 + dvp] * d16[i][14]) + (v2[2 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ }\r
+\r
+ }\r
+\r
+ private void compute_pcm_samples2() {\r
+\r
+ if (vidx == 1) {\r
+ // final float[] vp = actual_v;\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v1[2 + dvp] * d16[i][0]) + (v1[1 + dvp] * d16[i][1])\r
+ + (v1[0 + dvp] * d16[i][2]) + (v1[15 + dvp] * d16[i][3])\r
+ + (v1[14 + dvp] * d16[i][4]) + (v1[13 + dvp] * d16[i][5])\r
+ + (v1[12 + dvp] * d16[i][6]) + (v1[11 + dvp] * d16[i][7])\r
+ + (v1[10 + dvp] * d16[i][8]) + (v1[9 + dvp] * d16[i][9])\r
+ + (v1[8 + dvp] * d16[i][10]) + (v1[7 + dvp] * d16[i][11])\r
+ + (v1[6 + dvp] * d16[i][12]) + (v1[5 + dvp] * d16[i][13])\r
+ + (v1[4 + dvp] * d16[i][14]) + (v1[3 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ } else {\r
+ // final float[] vp = actual_v;\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v2[2 + dvp] * d16[i][0]) + (v2[1 + dvp] * d16[i][1])\r
+ + (v2[0 + dvp] * d16[i][2]) + (v2[15 + dvp] * d16[i][3])\r
+ + (v2[14 + dvp] * d16[i][4]) + (v2[13 + dvp] * d16[i][5])\r
+ + (v2[12 + dvp] * d16[i][6]) + (v2[11 + dvp] * d16[i][7])\r
+ + (v2[10 + dvp] * d16[i][8]) + (v2[9 + dvp] * d16[i][9])\r
+ + (v2[8 + dvp] * d16[i][10]) + (v2[7 + dvp] * d16[i][11])\r
+ + (v2[6 + dvp] * d16[i][12]) + (v2[5 + dvp] * d16[i][13])\r
+ + (v2[4 + dvp] * d16[i][14]) + (v2[3 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ }\r
+\r
+ }\r
+\r
+ private void compute_pcm_samples3() {\r
+\r
+ if (vidx == 1) {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v1[3 + dvp] * d16[i][0]) + (v1[2 + dvp] * d16[i][1])\r
+ + (v1[1 + dvp] * d16[i][2]) + (v1[0 + dvp] * d16[i][3])\r
+ + (v1[15 + dvp] * d16[i][4]) + (v1[14 + dvp] * d16[i][5])\r
+ + (v1[13 + dvp] * d16[i][6]) + (v1[12 + dvp] * d16[i][7])\r
+ + (v1[11 + dvp] * d16[i][8]) + (v1[10 + dvp] * d16[i][9])\r
+ + (v1[9 + dvp] * d16[i][10]) + (v1[8 + dvp] * d16[i][11])\r
+ + (v1[7 + dvp] * d16[i][12]) + (v1[6 + dvp] * d16[i][13])\r
+ + (v1[5 + dvp] * d16[i][14]) + (v1[4 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ } else {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v2[3 + dvp] * d16[i][0]) + (v2[2 + dvp] * d16[i][1])\r
+ + (v2[1 + dvp] * d16[i][2]) + (v2[0 + dvp] * d16[i][3])\r
+ + (v2[15 + dvp] * d16[i][4]) + (v2[14 + dvp] * d16[i][5])\r
+ + (v2[13 + dvp] * d16[i][6]) + (v2[12 + dvp] * d16[i][7])\r
+ + (v2[11 + dvp] * d16[i][8]) + (v2[10 + dvp] * d16[i][9])\r
+ + (v2[9 + dvp] * d16[i][10]) + (v2[8 + dvp] * d16[i][11])\r
+ + (v2[7 + dvp] * d16[i][12]) + (v2[6 + dvp] * d16[i][13])\r
+ + (v2[5 + dvp] * d16[i][14]) + (v2[4 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ }\r
+\r
+ }\r
+\r
+ private void compute_pcm_samples4() {\r
+\r
+ if (vidx == 1) {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v1[4 + dvp] * d16[i][0]) + (v1[3 + dvp] * d16[i][1])\r
+ + (v1[2 + dvp] * d16[i][2]) + (v1[1 + dvp] * d16[i][3]) + (v1[0 + dvp] * d16[i][4])\r
+ + (v1[15 + dvp] * d16[i][5]) + (v1[14 + dvp] * d16[i][6])\r
+ + (v1[13 + dvp] * d16[i][7]) + (v1[12 + dvp] * d16[i][8])\r
+ + (v1[11 + dvp] * d16[i][9]) + (v1[10 + dvp] * d16[i][10])\r
+ + (v1[9 + dvp] * d16[i][11]) + (v1[8 + dvp] * d16[i][12])\r
+ + (v1[7 + dvp] * d16[i][13]) + (v1[6 + dvp] * d16[i][14]) + (v1[5 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ } else {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v2[4 + dvp] * d16[i][0]) + (v2[3 + dvp] * d16[i][1])\r
+ + (v2[2 + dvp] * d16[i][2]) + (v2[1 + dvp] * d16[i][3]) + (v2[0 + dvp] * d16[i][4])\r
+ + (v2[15 + dvp] * d16[i][5]) + (v2[14 + dvp] * d16[i][6])\r
+ + (v2[13 + dvp] * d16[i][7]) + (v2[12 + dvp] * d16[i][8])\r
+ + (v2[11 + dvp] * d16[i][9]) + (v2[10 + dvp] * d16[i][10])\r
+ + (v2[9 + dvp] * d16[i][11]) + (v2[8 + dvp] * d16[i][12])\r
+ + (v2[7 + dvp] * d16[i][13]) + (v2[6 + dvp] * d16[i][14]) + (v2[5 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ }\r
+\r
+ }\r
+\r
+ private void compute_pcm_samples5() {\r
+\r
+ if (vidx == 1) {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v1[5 + dvp] * d16[i][0]) + (v1[4 + dvp] * d16[i][1])\r
+ + (v1[3 + dvp] * d16[i][2]) + (v1[2 + dvp] * d16[i][3]) + (v1[1 + dvp] * d16[i][4])\r
+ + (v1[0 + dvp] * d16[i][5]) + (v1[15 + dvp] * d16[i][6])\r
+ + (v1[14 + dvp] * d16[i][7]) + (v1[13 + dvp] * d16[i][8])\r
+ + (v1[12 + dvp] * d16[i][9]) + (v1[11 + dvp] * d16[i][10])\r
+ + (v1[10 + dvp] * d16[i][11]) + (v1[9 + dvp] * d16[i][12])\r
+ + (v1[8 + dvp] * d16[i][13]) + (v1[7 + dvp] * d16[i][14]) + (v1[6 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ } else {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v2[5 + dvp] * d16[i][0]) + (v2[4 + dvp] * d16[i][1])\r
+ + (v2[3 + dvp] * d16[i][2]) + (v2[2 + dvp] * d16[i][3]) + (v2[1 + dvp] * d16[i][4])\r
+ + (v2[0 + dvp] * d16[i][5]) + (v2[15 + dvp] * d16[i][6])\r
+ + (v2[14 + dvp] * d16[i][7]) + (v2[13 + dvp] * d16[i][8])\r
+ + (v2[12 + dvp] * d16[i][9]) + (v2[11 + dvp] * d16[i][10])\r
+ + (v2[10 + dvp] * d16[i][11]) + (v2[9 + dvp] * d16[i][12])\r
+ + (v2[8 + dvp] * d16[i][13]) + (v2[7 + dvp] * d16[i][14]) + (v2[6 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ }\r
+\r
+ }\r
+\r
+ private void compute_pcm_samples6() {\r
+\r
+ if (vidx == 1) {\r
+\r
+ // final float[] vp = actual_v;\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v1[6 + dvp] * d16[i][0]) + (v1[5 + dvp] * d16[i][1])\r
+ + (v1[4 + dvp] * d16[i][2]) + (v1[3 + dvp] * d16[i][3]) + (v1[2 + dvp] * d16[i][4])\r
+ + (v1[1 + dvp] * d16[i][5]) + (v1[0 + dvp] * d16[i][6])\r
+ + (v1[15 + dvp] * d16[i][7]) + (v1[14 + dvp] * d16[i][8])\r
+ + (v1[13 + dvp] * d16[i][9]) + (v1[12 + dvp] * d16[i][10])\r
+ + (v1[11 + dvp] * d16[i][11]) + (v1[10 + dvp] * d16[i][12])\r
+ + (v1[9 + dvp] * d16[i][13]) + (v1[8 + dvp] * d16[i][14]) + (v1[7 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ } else {\r
+\r
+ // final float[] vp = actual_v;\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v2[6 + dvp] * d16[i][0]) + (v2[5 + dvp] * d16[i][1])\r
+ + (v2[4 + dvp] * d16[i][2]) + (v2[3 + dvp] * d16[i][3]) + (v2[2 + dvp] * d16[i][4])\r
+ + (v2[1 + dvp] * d16[i][5]) + (v2[0 + dvp] * d16[i][6])\r
+ + (v2[15 + dvp] * d16[i][7]) + (v2[14 + dvp] * d16[i][8])\r
+ + (v2[13 + dvp] * d16[i][9]) + (v2[12 + dvp] * d16[i][10])\r
+ + (v2[11 + dvp] * d16[i][11]) + (v2[10 + dvp] * d16[i][12])\r
+ + (v2[9 + dvp] * d16[i][13]) + (v2[8 + dvp] * d16[i][14]) + (v2[7 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ }\r
+\r
+ }\r
+\r
+ private void compute_pcm_samples7() {\r
+\r
+ if (vidx == 1) {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v1[7 + dvp] * d16[i][0]) + (v1[6 + dvp] * d16[i][1])\r
+ + (v1[5 + dvp] * d16[i][2]) + (v1[4 + dvp] * d16[i][3]) + (v1[3 + dvp] * d16[i][4])\r
+ + (v1[2 + dvp] * d16[i][5]) + (v1[1 + dvp] * d16[i][6]) + (v1[0 + dvp] * d16[i][7])\r
+ + (v1[15 + dvp] * d16[i][8]) + (v1[14 + dvp] * d16[i][9])\r
+ + (v1[13 + dvp] * d16[i][10]) + (v1[12 + dvp] * d16[i][11])\r
+ + (v1[11 + dvp] * d16[i][12]) + (v1[10 + dvp] * d16[i][13])\r
+ + (v1[9 + dvp] * d16[i][14]) + (v1[8 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ } else {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v2[7 + dvp] * d16[i][0]) + (v2[6 + dvp] * d16[i][1])\r
+ + (v2[5 + dvp] * d16[i][2]) + (v2[4 + dvp] * d16[i][3]) + (v2[3 + dvp] * d16[i][4])\r
+ + (v2[2 + dvp] * d16[i][5]) + (v2[1 + dvp] * d16[i][6]) + (v2[0 + dvp] * d16[i][7])\r
+ + (v2[15 + dvp] * d16[i][8]) + (v2[14 + dvp] * d16[i][9])\r
+ + (v2[13 + dvp] * d16[i][10]) + (v2[12 + dvp] * d16[i][11])\r
+ + (v2[11 + dvp] * d16[i][12]) + (v2[10 + dvp] * d16[i][13])\r
+ + (v2[9 + dvp] * d16[i][14]) + (v2[8 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ }\r
+\r
+ }\r
+\r
+ private void compute_pcm_samples8() {\r
+\r
+ if (vidx == 1) {\r
+\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v1[8 + dvp] * d16[i][0]) + (v1[7 + dvp] * d16[i][1])\r
+ + (v1[6 + dvp] * d16[i][2]) + (v1[5 + dvp] * d16[i][3]) + (v1[4 + dvp] * d16[i][4])\r
+ + (v1[3 + dvp] * d16[i][5]) + (v1[2 + dvp] * d16[i][6]) + (v1[1 + dvp] * d16[i][7])\r
+ + (v1[0 + dvp] * d16[i][8]) + (v1[15 + dvp] * d16[i][9])\r
+ + (v1[14 + dvp] * d16[i][10]) + (v1[13 + dvp] * d16[i][11])\r
+ + (v1[12 + dvp] * d16[i][12]) + (v1[11 + dvp] * d16[i][13])\r
+ + (v1[10 + dvp] * d16[i][14]) + (v1[9 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ } else {\r
+\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v2[8 + dvp] * d16[i][0]) + (v2[7 + dvp] * d16[i][1])\r
+ + (v2[6 + dvp] * d16[i][2]) + (v2[5 + dvp] * d16[i][3]) + (v2[4 + dvp] * d16[i][4])\r
+ + (v2[3 + dvp] * d16[i][5]) + (v2[2 + dvp] * d16[i][6]) + (v2[1 + dvp] * d16[i][7])\r
+ + (v2[0 + dvp] * d16[i][8]) + (v2[15 + dvp] * d16[i][9])\r
+ + (v2[14 + dvp] * d16[i][10]) + (v2[13 + dvp] * d16[i][11])\r
+ + (v2[12 + dvp] * d16[i][12]) + (v2[11 + dvp] * d16[i][13])\r
+ + (v2[10 + dvp] * d16[i][14]) + (v2[9 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ }\r
+\r
+ }\r
+\r
+ private void compute_pcm_samples9() {\r
+\r
+ if (vidx == 1) {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v1[9 + dvp] * d16[i][0]) + (v1[8 + dvp] * d16[i][1])\r
+ + (v1[7 + dvp] * d16[i][2]) + (v1[6 + dvp] * d16[i][3]) + (v1[5 + dvp] * d16[i][4])\r
+ + (v1[4 + dvp] * d16[i][5]) + (v1[3 + dvp] * d16[i][6]) + (v1[2 + dvp] * d16[i][7])\r
+ + (v1[1 + dvp] * d16[i][8]) + (v1[0 + dvp] * d16[i][9])\r
+ + (v1[15 + dvp] * d16[i][10]) + (v1[14 + dvp] * d16[i][11])\r
+ + (v1[13 + dvp] * d16[i][12]) + (v1[12 + dvp] * d16[i][13])\r
+ + (v1[11 + dvp] * d16[i][14]) + (v1[10 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ } else {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v2[9 + dvp] * d16[i][0]) + (v2[8 + dvp] * d16[i][1])\r
+ + (v2[7 + dvp] * d16[i][2]) + (v2[6 + dvp] * d16[i][3]) + (v2[5 + dvp] * d16[i][4])\r
+ + (v2[4 + dvp] * d16[i][5]) + (v2[3 + dvp] * d16[i][6]) + (v2[2 + dvp] * d16[i][7])\r
+ + (v2[1 + dvp] * d16[i][8]) + (v2[0 + dvp] * d16[i][9])\r
+ + (v2[15 + dvp] * d16[i][10]) + (v2[14 + dvp] * d16[i][11])\r
+ + (v2[13 + dvp] * d16[i][12]) + (v2[12 + dvp] * d16[i][13])\r
+ + (v2[11 + dvp] * d16[i][14]) + (v2[10 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ }\r
+\r
+ }\r
+\r
+ private void compute_pcm_samples10() {\r
+ if (vidx == 1) {\r
+ // final float[] vp = actual_v;\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v1[10 + dvp] * d16[i][0]) + (v1[9 + dvp] * d16[i][1])\r
+ + (v1[8 + dvp] * d16[i][2]) + (v1[7 + dvp] * d16[i][3]) + (v1[6 + dvp] * d16[i][4])\r
+ + (v1[5 + dvp] * d16[i][5]) + (v1[4 + dvp] * d16[i][6]) + (v1[3 + dvp] * d16[i][7])\r
+ + (v1[2 + dvp] * d16[i][8]) + (v1[1 + dvp] * d16[i][9])\r
+ + (v1[0 + dvp] * d16[i][10]) + (v1[15 + dvp] * d16[i][11])\r
+ + (v1[14 + dvp] * d16[i][12]) + (v1[13 + dvp] * d16[i][13])\r
+ + (v1[12 + dvp] * d16[i][14]) + (v1[11 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ } else {\r
+ // final float[] vp = actual_v;\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v2[10 + dvp] * d16[i][0]) + (v2[9 + dvp] * d16[i][1])\r
+ + (v2[8 + dvp] * d16[i][2]) + (v2[7 + dvp] * d16[i][3]) + (v2[6 + dvp] * d16[i][4])\r
+ + (v2[5 + dvp] * d16[i][5]) + (v2[4 + dvp] * d16[i][6]) + (v2[3 + dvp] * d16[i][7])\r
+ + (v2[2 + dvp] * d16[i][8]) + (v2[1 + dvp] * d16[i][9])\r
+ + (v2[0 + dvp] * d16[i][10]) + (v2[15 + dvp] * d16[i][11])\r
+ + (v2[14 + dvp] * d16[i][12]) + (v2[13 + dvp] * d16[i][13])\r
+ + (v2[12 + dvp] * d16[i][14]) + (v2[11 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ }\r
+\r
+ }\r
+\r
+ private void compute_pcm_samples11() {\r
+\r
+ if (vidx == 1) {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v1[11 + dvp] * d16[i][0]) + (v1[10 + dvp] * d16[i][1])\r
+ + (v1[9 + dvp] * d16[i][2]) + (v1[8 + dvp] * d16[i][3]) + (v1[7 + dvp] * d16[i][4])\r
+ + (v1[6 + dvp] * d16[i][5]) + (v1[5 + dvp] * d16[i][6]) + (v1[4 + dvp] * d16[i][7])\r
+ + (v1[3 + dvp] * d16[i][8]) + (v1[2 + dvp] * d16[i][9])\r
+ + (v1[1 + dvp] * d16[i][10]) + (v1[0 + dvp] * d16[i][11])\r
+ + (v1[15 + dvp] * d16[i][12]) + (v1[14 + dvp] * d16[i][13])\r
+ + (v1[13 + dvp] * d16[i][14]) + (v1[12 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ } else {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v2[11 + dvp] * d16[i][0]) + (v2[10 + dvp] * d16[i][1])\r
+ + (v2[9 + dvp] * d16[i][2]) + (v2[8 + dvp] * d16[i][3]) + (v2[7 + dvp] * d16[i][4])\r
+ + (v2[6 + dvp] * d16[i][5]) + (v2[5 + dvp] * d16[i][6]) + (v2[4 + dvp] * d16[i][7])\r
+ + (v2[3 + dvp] * d16[i][8]) + (v2[2 + dvp] * d16[i][9])\r
+ + (v2[1 + dvp] * d16[i][10]) + (v2[0 + dvp] * d16[i][11])\r
+ + (v2[15 + dvp] * d16[i][12]) + (v2[14 + dvp] * d16[i][13])\r
+ + (v2[13 + dvp] * d16[i][14]) + (v2[12 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ }\r
+\r
+ }\r
+\r
+ private void compute_pcm_samples12() {\r
+\r
+ if (vidx == 1) {\r
+ // final float[] vp = actual_v;\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v1[12 + dvp] * d16[i][0]) + (v1[11 + dvp] * d16[i][1])\r
+ + (v1[10 + dvp] * d16[i][2]) + (v1[9 + dvp] * d16[i][3])\r
+ + (v1[8 + dvp] * d16[i][4]) + (v1[7 + dvp] * d16[i][5]) + (v1[6 + dvp] * d16[i][6])\r
+ + (v1[5 + dvp] * d16[i][7]) + (v1[4 + dvp] * d16[i][8]) + (v1[3 + dvp] * d16[i][9])\r
+ + (v1[2 + dvp] * d16[i][10]) + (v1[1 + dvp] * d16[i][11])\r
+ + (v1[0 + dvp] * d16[i][12]) + (v1[15 + dvp] * d16[i][13])\r
+ + (v1[14 + dvp] * d16[i][14]) + (v1[13 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ } else {\r
+ // final float[] vp = actual_v;\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v2[12 + dvp] * d16[i][0]) + (v2[11 + dvp] * d16[i][1])\r
+ + (v2[10 + dvp] * d16[i][2]) + (v2[9 + dvp] * d16[i][3])\r
+ + (v2[8 + dvp] * d16[i][4]) + (v2[7 + dvp] * d16[i][5]) + (v2[6 + dvp] * d16[i][6])\r
+ + (v2[5 + dvp] * d16[i][7]) + (v2[4 + dvp] * d16[i][8]) + (v2[3 + dvp] * d16[i][9])\r
+ + (v2[2 + dvp] * d16[i][10]) + (v2[1 + dvp] * d16[i][11])\r
+ + (v2[0 + dvp] * d16[i][12]) + (v2[15 + dvp] * d16[i][13])\r
+ + (v2[14 + dvp] * d16[i][14]) + (v2[13 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ }\r
+\r
+ }\r
+\r
+ private void compute_pcm_samples13() {\r
+\r
+ if (vidx == 1) {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v1[13 + dvp] * d16[i][0]) + (v1[12 + dvp] * d16[i][1])\r
+ + (v1[11 + dvp] * d16[i][2]) + (v1[10 + dvp] * d16[i][3])\r
+ + (v1[9 + dvp] * d16[i][4]) + (v1[8 + dvp] * d16[i][5]) + (v1[7 + dvp] * d16[i][6])\r
+ + (v1[6 + dvp] * d16[i][7]) + (v1[5 + dvp] * d16[i][8]) + (v1[4 + dvp] * d16[i][9])\r
+ + (v1[3 + dvp] * d16[i][10]) + (v1[2 + dvp] * d16[i][11])\r
+ + (v1[1 + dvp] * d16[i][12]) + (v1[0 + dvp] * d16[i][13])\r
+ + (v1[15 + dvp] * d16[i][14]) + (v1[14 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ } else {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v2[13 + dvp] * d16[i][0]) + (v2[12 + dvp] * d16[i][1])\r
+ + (v2[11 + dvp] * d16[i][2]) + (v2[10 + dvp] * d16[i][3])\r
+ + (v2[9 + dvp] * d16[i][4]) + (v2[8 + dvp] * d16[i][5]) + (v2[7 + dvp] * d16[i][6])\r
+ + (v2[6 + dvp] * d16[i][7]) + (v2[5 + dvp] * d16[i][8]) + (v2[4 + dvp] * d16[i][9])\r
+ + (v2[3 + dvp] * d16[i][10]) + (v2[2 + dvp] * d16[i][11])\r
+ + (v2[1 + dvp] * d16[i][12]) + (v2[0 + dvp] * d16[i][13])\r
+ + (v2[15 + dvp] * d16[i][14]) + (v2[14 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ }\r
+\r
+ }\r
+\r
+ private void compute_pcm_samples14() {\r
+\r
+ if (vidx == 1) {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v1[14 + dvp] * d16[i][0]) + (v1[13 + dvp] * d16[i][1])\r
+ + (v1[12 + dvp] * d16[i][2]) + (v1[11 + dvp] * d16[i][3])\r
+ + (v1[10 + dvp] * d16[i][4]) + (v1[9 + dvp] * d16[i][5])\r
+ + (v1[8 + dvp] * d16[i][6]) + (v1[7 + dvp] * d16[i][7]) + (v1[6 + dvp] * d16[i][8])\r
+ + (v1[5 + dvp] * d16[i][9]) + (v1[4 + dvp] * d16[i][10])\r
+ + (v1[3 + dvp] * d16[i][11]) + (v1[2 + dvp] * d16[i][12])\r
+ + (v1[1 + dvp] * d16[i][13]) + (v1[0 + dvp] * d16[i][14]) + (v1[15 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ } else {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ // final float[] dp = d16[i];\r
+ float pcm_sample;\r
+\r
+ pcm_sample =\r
+ (float) (((v2[14 + dvp] * d16[i][0]) + (v2[13 + dvp] * d16[i][1])\r
+ + (v2[12 + dvp] * d16[i][2]) + (v2[11 + dvp] * d16[i][3])\r
+ + (v2[10 + dvp] * d16[i][4]) + (v2[9 + dvp] * d16[i][5])\r
+ + (v2[8 + dvp] * d16[i][6]) + (v2[7 + dvp] * d16[i][7]) + (v2[6 + dvp] * d16[i][8])\r
+ + (v2[5 + dvp] * d16[i][9]) + (v2[4 + dvp] * d16[i][10])\r
+ + (v2[3 + dvp] * d16[i][11]) + (v2[2 + dvp] * d16[i][12])\r
+ + (v2[1 + dvp] * d16[i][13]) + (v2[0 + dvp] * d16[i][14]) + (v2[15 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+\r
+ dvp += 16;\r
+ } // for\r
+ }\r
+\r
+ }\r
+\r
+ private void compute_pcm_samples15() {\r
+ if (vidx == 1) {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ float pcm_sample;\r
+ // final float d16[i][] = d16[i];\r
+ pcm_sample =\r
+ (float) (((v1[15 + dvp] * d16[i][0]) + (v1[14 + dvp] * d16[i][1])\r
+ + (v1[13 + dvp] * d16[i][2]) + (v1[12 + dvp] * d16[i][3])\r
+ + (v1[11 + dvp] * d16[i][4]) + (v1[10 + dvp] * d16[i][5])\r
+ + (v1[9 + dvp] * d16[i][6]) + (v1[8 + dvp] * d16[i][7]) + (v1[7 + dvp] * d16[i][8])\r
+ + (v1[6 + dvp] * d16[i][9]) + (v1[5 + dvp] * d16[i][10])\r
+ + (v1[4 + dvp] * d16[i][11]) + (v1[3 + dvp] * d16[i][12])\r
+ + (v1[2 + dvp] * d16[i][13]) + (v1[1 + dvp] * d16[i][14]) + (v1[0 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+ dvp += 16;\r
+ } // for\r
+ } else {\r
+ // final float[] vp = actual_v;\r
+\r
+ // int inc = v_inc;\r
+ // final float[] tmpOut = _tmpOut;\r
+ int dvp = 0;\r
+\r
+ // fat chance of having this loop unroll\r
+ for ( int i = 0; i < 32; i++) {\r
+ float pcm_sample;\r
+ // final float d16[i][] = d16[i];\r
+ pcm_sample =\r
+ (float) (((v2[15 + dvp] * d16[i][0]) + (v2[14 + dvp] * d16[i][1])\r
+ + (v2[13 + dvp] * d16[i][2]) + (v2[12 + dvp] * d16[i][3])\r
+ + (v2[11 + dvp] * d16[i][4]) + (v2[10 + dvp] * d16[i][5])\r
+ + (v2[9 + dvp] * d16[i][6]) + (v2[8 + dvp] * d16[i][7]) + (v2[7 + dvp] * d16[i][8])\r
+ + (v2[6 + dvp] * d16[i][9]) + (v2[5 + dvp] * d16[i][10])\r
+ + (v2[4 + dvp] * d16[i][11]) + (v2[3 + dvp] * d16[i][12])\r
+ + (v2[2 + dvp] * d16[i][13]) + (v2[1 + dvp] * d16[i][14]) + (v2[0 + dvp] * d16[i][15])) * scalefactor);\r
+\r
+ _tmpOut[i] = pcm_sample;\r
+ dvp += 16;\r
+ } // for\r
+ }\r
+\r
+ }\r
+\r
+ private void compute_pcm_samples() {\r
+\r
+ switch (actual_write_pos) {\r
+ case 0:\r
+ compute_pcm_samples0();\r
+ break;\r
+ case 1:\r
+ compute_pcm_samples1();\r
+ break;\r
+ case 2:\r
+ compute_pcm_samples2();\r
+ break;\r
+ case 3:\r
+ compute_pcm_samples3();\r
+ break;\r
+ case 4:\r
+ compute_pcm_samples4();\r
+ break;\r
+ case 5:\r
+ compute_pcm_samples5();\r
+ break;\r
+ case 6:\r
+ compute_pcm_samples6();\r
+ break;\r
+ case 7:\r
+ compute_pcm_samples7();\r
+ break;\r
+ case 8:\r
+ compute_pcm_samples8();\r
+ break;\r
+ case 9:\r
+ compute_pcm_samples9();\r
+ break;\r
+ case 10:\r
+ compute_pcm_samples10();\r
+ break;\r
+ case 11:\r
+ compute_pcm_samples11();\r
+ break;\r
+ case 12:\r
+ compute_pcm_samples12();\r
+ break;\r
+ case 13:\r
+ compute_pcm_samples13();\r
+ break;\r
+ case 14:\r
+ compute_pcm_samples14();\r
+ break;\r
+ case 15:\r
+ compute_pcm_samples15();\r
+ break;\r
+ }\r
+\r
+ // if (buffer != null) {\r
+ // buffer.appendSamples(channel, _tmpOut);\r
+ // }\r
+ SampleBufferWrapper.appendSamples(channel, _tmpOut);\r
+\r
+ /*\r
+ * // MDM: I was considering putting in quality control for // low-spec\r
+ * CPUs, but the performance gain (about 10-15%) // did not justify the\r
+ * considerable drop in audio quality. switch (inc) { case 16:\r
+ * buffer.appendSamples(channel, tmpOut); break; case 32: for (int i=0;\r
+ * i<16; i++) { buffer.append(channel, (short)tmpOut[i]);\r
+ * buffer.append(channel, (short)tmpOut[i]); } break; case 64: for (int i=0;\r
+ * i<8; i++) { buffer.append(channel, (short)tmpOut[i]);\r
+ * buffer.append(channel, (short)tmpOut[i]); buffer.append(channel,\r
+ * (short)tmpOut[i]); buffer.append(channel, (short)tmpOut[i]); } break;\r
+ * \r
+ * }\r
+ */\r
+ }\r
+\r
+ \r
+ public void clear() {\r
+ // clear out v1,v2\r
+ SSJAVA.arrayinit(v1, 0);\r
+ SSJAVA.arrayinit(v2, 0);\r
+ // copy previous v1,v2\r
+\r
+ for ( int i = 0; i < prev1.length; i++) {\r
+ v1[i] = prev1[i];\r
+ }\r
+ for ( int i = 0; i < prev2.length; i++) {\r
+ v2[i] = prev2[i];\r
+ }\r
+ // clear out previous buffer\r
+ SSJAVA.arrayinit(prev1, 0);\r
+ SSJAVA.arrayinit(prev2, 0);\r
+ SSJAVA.arrayinit(samples, 0);\r
+ SSJAVA.arrayinit(_tmpOut, 0);\r
+ }\r
+\r
+ /**\r
+ * Calculate 32 PCM samples and put the into the Obuffer-object.\r
+ */\r
+\r
+ public void calculate_pcm_samples() {\r
+\r
+// System.out.println("#calculate_pcm_samples::actual_write_pos=" + actual_write_pos);\r
+\r
+ if (vidx == 1) {\r
+ compute_new_v1_v2();\r
+ } else {\r
+ compute_new_v2_v1();\r
+ }\r
+\r
+ // System.out.println("1.actual_v=" + (actual_v == v1) + " vidx=" + vidx);\r
+ // compute_new_v();\r
+ // System.out.println("2.actual_v=" + (actual_v == v1) + " vidx=" + vidx);\r
+ compute_pcm_samples();\r
+ // System.out.println("3.actual_v=" + (actual_v == v1) + " vidx=" + vidx);\r
+\r
+ actual_write_pos = (actual_write_pos + 1) & 0xf;\r
+// System.out.println("actual_write_pos="+actual_write_pos);\r
+ // actual_v = (actual_v == v1) ? v2 : v1;\r
+\r
+ if (vidx == 1) {\r
+ vidx = 2;\r
+ } else {\r
+ vidx = 1;\r
+ }\r
+\r
+ // initialize samples[]:\r
+ // for (register float *floatp = samples + 32; floatp > samples; )\r
+ // *--floatp = 0.0f;\r
+\r
+ // MDM: this may not be necessary. The Layer III decoder always\r
+ // outputs 32 subband samples, but I haven't checked layer I & II.\r
+ // for ( int p = 0; p < 32; p++){\r
+ // samples[p] = 0.0f;\r
+ // }\r
+ SSJAVA.arrayinit(samples, 0);\r
+ }\r
+\r
+ private static final double MY_PI = 3.14159265358979323846;\r
+ private static final float cos1_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI / 64.0)));\r
+ private static final float cos3_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 3.0 / 64.0)));\r
+ private static final float cos5_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 5.0 / 64.0)));\r
+ private static final float cos7_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 7.0 / 64.0)));\r
+ private static final float cos9_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 9.0 / 64.0)));\r
+ private static final float cos11_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 11.0 / 64.0)));\r
+ private static final float cos13_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 13.0 / 64.0)));\r
+ private static final float cos15_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 15.0 / 64.0)));\r
+ private static final float cos17_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 17.0 / 64.0)));\r
+ private static final float cos19_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 19.0 / 64.0)));\r
+ private static final float cos21_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 21.0 / 64.0)));\r
+ private static final float cos23_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 23.0 / 64.0)));\r
+ private static final float cos25_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 25.0 / 64.0)));\r
+ private static final float cos27_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 27.0 / 64.0)));\r
+ private static final float cos29_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 29.0 / 64.0)));\r
+ private static final float cos31_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 31.0 / 64.0)));\r
+ private static final float cos1_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI / 32.0)));\r
+ private static final float cos3_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 3.0 / 32.0)));\r
+ private static final float cos5_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 5.0 / 32.0)));\r
+ private static final float cos7_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 7.0 / 32.0)));\r
+ private static final float cos9_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 9.0 / 32.0)));\r
+ private static final float cos11_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 11.0 / 32.0)));\r
+ private static final float cos13_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 13.0 / 32.0)));\r
+ private static final float cos15_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 15.0 / 32.0)));\r
+ private static final float cos1_16 = (float) (1.0 / (2.0 * Math.cos(MY_PI / 16.0)));\r
+ private static final float cos3_16 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 3.0 / 16.0)));\r
+ private static final float cos5_16 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 5.0 / 16.0)));\r
+ private static final float cos7_16 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 7.0 / 16.0)));\r
+ private static final float cos1_8 = (float) (1.0 / (2.0 * Math.cos(MY_PI / 8.0)));\r
+ private static final float cos3_8 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 3.0 / 8.0)));\r
+ private static final float cos1_4 = (float) (1.0 / (2.0 * Math.cos(MY_PI / 4.0)));\r
+\r
+ // Note: These values are not in the same order\r
+ // as in Annex 3-B.3 of the ISO/IEC DIS 11172-3\r
+ // private float d[] = {0.000000000, -4.000442505};\r
+\r
+ /**\r
+ * d[] split into subarrays of length 16. This provides for more faster access\r
+ * by allowing a block of 16 to be addressed with constant offset.\r
+ **/\r
+ // TODO CONST\r
+ \r
+ private final static float d16[][] = null;\r
+\r
+ /**\r
+ * Converts a 1D array into a number of smaller arrays. This is used to\r
+ * achieve offset + constant indexing into an array. Each sub-array represents\r
+ * a block of values of the original array.\r
+ * \r
+ * @param array\r
+ * The array to split up into blocks.\r
+ * @param blockSize\r
+ * The size of the blocks to split the array into. This must be an\r
+ * exact divisor of the length of the array, or some data will be\r
+ * lost from the main array.\r
+ * \r
+ * @return An array of arrays in which each element in the returned array will\r
+ * be of length <code>blockSize</code>.\r
+ */\r
+ static private float[][] splitArray(final float[] array, final int blockSize) {\r
+ int size = array.length / blockSize;\r
+ float[][] split = new float[size][];\r
+ for (int i = 0; i < size; i++) {\r
+ split[i] = subArray(array, i * blockSize, blockSize);\r
+ }\r
+ return split;\r
+ }\r
+\r
+ /**\r
+ * Returns a subarray of an existing array.\r
+ * \r
+ * @param array\r
+ * The array to retrieve a subarra from.\r
+ * @param offs\r
+ * The offset in the array that corresponds to the first index of the\r
+ * subarray.\r
+ * @param len\r
+ * The number of indeces in the subarray.\r
+ * @return The subarray, which may be of length 0.\r
+ */\r
+ static private float[] subArray(final float[] array, final int offs, int len) {\r
+ if (offs + len > array.length) {\r
+ len = array.length - offs;\r
+ }\r
+\r
+ if (len < 0)\r
+ len = 0;\r
+\r
+ float[] subarray = new float[len];\r
+ for (int i = 0; i < len; i++) {\r
+ subarray[i] = array[offs + i];\r
+ }\r
+\r
+ return subarray;\r
+ }\r
+\r
+ // The original data for d[]. This data is loaded from a file\r
+ // to reduce the overall package size and to improve performance.\r
+\r
+ static final float d[] = { 0.000000000f, -0.000442505f, 0.003250122f, -0.007003784f,\r
+ 0.031082153f, -0.078628540f, 0.100311279f, -0.572036743f, 1.144989014f, 0.572036743f,\r
+ 0.100311279f, 0.078628540f, 0.031082153f, 0.007003784f, 0.003250122f, 0.000442505f,\r
+ -0.000015259f, -0.000473022f, 0.003326416f, -0.007919312f, 0.030517578f, -0.084182739f,\r
+ 0.090927124f, -0.600219727f, 1.144287109f, 0.543823242f, 0.108856201f, 0.073059082f,\r
+ 0.031478882f, 0.006118774f, 0.003173828f, 0.000396729f, -0.000015259f, -0.000534058f,\r
+ 0.003387451f, -0.008865356f, 0.029785156f, -0.089706421f, 0.080688477f, -0.628295898f,\r
+ 1.142211914f, 0.515609741f, 0.116577148f, 0.067520142f, 0.031738281f, 0.005294800f,\r
+ 0.003082275f, 0.000366211f, -0.000015259f, -0.000579834f, 0.003433228f, -0.009841919f,\r
+ 0.028884888f, -0.095169067f, 0.069595337f, -0.656219482f, 1.138763428f, 0.487472534f,\r
+ 0.123474121f, 0.061996460f, 0.031845093f, 0.004486084f, 0.002990723f, 0.000320435f,\r
+ -0.000015259f, -0.000625610f, 0.003463745f, -0.010848999f, 0.027801514f, -0.100540161f,\r
+ 0.057617188f, -0.683914185f, 1.133926392f, 0.459472656f, 0.129577637f, 0.056533813f,\r
+ 0.031814575f, 0.003723145f, 0.002899170f, 0.000289917f, -0.000015259f, -0.000686646f,\r
+ 0.003479004f, -0.011886597f, 0.026535034f, -0.105819702f, 0.044784546f, -0.711318970f,\r
+ 1.127746582f, 0.431655884f, 0.134887695f, 0.051132202f, 0.031661987f, 0.003005981f,\r
+ 0.002792358f, 0.000259399f, -0.000015259f, -0.000747681f, 0.003479004f, -0.012939453f,\r
+ 0.025085449f, -0.110946655f, 0.031082153f, -0.738372803f, 1.120223999f, 0.404083252f,\r
+ 0.139450073f, 0.045837402f, 0.031387329f, 0.002334595f, 0.002685547f, 0.000244141f,\r
+ -0.000030518f, -0.000808716f, 0.003463745f, -0.014022827f, 0.023422241f, -0.115921021f,\r
+ 0.016510010f, -0.765029907f, 1.111373901f, 0.376800537f, 0.143264771f, 0.040634155f,\r
+ 0.031005859f, 0.001693726f, 0.002578735f, 0.000213623f, -0.000030518f, -0.000885010f,\r
+ 0.003417969f, -0.015121460f, 0.021575928f, -0.120697021f, 0.001068115f, -0.791213989f,\r
+ 1.101211548f, 0.349868774f, 0.146362305f, 0.035552979f, 0.030532837f, 0.001098633f,\r
+ 0.002456665f, 0.000198364f, -0.000030518f, -0.000961304f, 0.003372192f, -0.016235352f,\r
+ 0.019531250f, -0.125259399f, -0.015228271f, -0.816864014f, 1.089782715f, 0.323318481f,\r
+ 0.148773193f, 0.030609131f, 0.029937744f, 0.000549316f, 0.002349854f, 0.000167847f,\r
+ -0.000030518f, -0.001037598f, 0.003280640f, -0.017349243f, 0.017257690f, -0.129562378f,\r
+ -0.032379150f, -0.841949463f, 1.077117920f, 0.297210693f, 0.150497437f, 0.025817871f,\r
+ 0.029281616f, 0.000030518f, 0.002243042f, 0.000152588f, -0.000045776f, -0.001113892f,\r
+ 0.003173828f, -0.018463135f, 0.014801025f, -0.133590698f, -0.050354004f, -0.866363525f,\r
+ 1.063217163f, 0.271591187f, 0.151596069f, 0.021179199f, 0.028533936f, -0.000442505f,\r
+ 0.002120972f, 0.000137329f, -0.000045776f, -0.001205444f, 0.003051758f, -0.019577026f,\r
+ 0.012115479f, -0.137298584f, -0.069168091f, -0.890090942f, 1.048156738f, 0.246505737f,\r
+ 0.152069092f, 0.016708374f, 0.027725220f, -0.000869751f, 0.002014160f, 0.000122070f,\r
+ -0.000061035f, -0.001296997f, 0.002883911f, -0.020690918f, 0.009231567f, -0.140670776f,\r
+ -0.088775635f, -0.913055420f, 1.031936646f, 0.221984863f, 0.151962280f, 0.012420654f,\r
+ 0.026840210f, -0.001266479f, 0.001907349f, 0.000106812f, -0.000061035f, -0.001388550f,\r
+ 0.002700806f, -0.021789551f, 0.006134033f, -0.143676758f, -0.109161377f, -0.935195923f,\r
+ 1.014617920f, 0.198059082f, 0.151306152f, 0.008316040f, 0.025909424f, -0.001617432f,\r
+ 0.001785278f, 0.000106812f, -0.000076294f, -0.001480103f, 0.002487183f, -0.022857666f,\r
+ 0.002822876f, -0.146255493f, -0.130310059f, -0.956481934f, 0.996246338f, 0.174789429f,\r
+ 0.150115967f, 0.004394531f, 0.024932861f, -0.001937866f, 0.001693726f, 0.000091553f,\r
+ -0.000076294f, -0.001586914f, 0.002227783f, -0.023910522f, -0.000686646f, -0.148422241f,\r
+ -0.152206421f, -0.976852417f, 0.976852417f, 0.152206421f, 0.148422241f, 0.000686646f,\r
+ 0.023910522f, -0.002227783f, 0.001586914f, 0.000076294f, -0.000091553f, -0.001693726f,\r
+ 0.001937866f, -0.024932861f, -0.004394531f, -0.150115967f, -0.174789429f, -0.996246338f,\r
+ 0.956481934f, 0.130310059f, 0.146255493f, -0.002822876f, 0.022857666f, -0.002487183f,\r
+ 0.001480103f, 0.000076294f, -0.000106812f, -0.001785278f, 0.001617432f, -0.025909424f,\r
+ -0.008316040f, -0.151306152f, -0.198059082f, -1.014617920f, 0.935195923f, 0.109161377f,\r
+ 0.143676758f, -0.006134033f, 0.021789551f, -0.002700806f, 0.001388550f, 0.000061035f,\r
+ -0.000106812f, -0.001907349f, 0.001266479f, -0.026840210f, -0.012420654f, -0.151962280f,\r
+ -0.221984863f, -1.031936646f, 0.913055420f, 0.088775635f, 0.140670776f, -0.009231567f,\r
+ 0.020690918f, -0.002883911f, 0.001296997f, 0.000061035f, -0.000122070f, -0.002014160f,\r
+ 0.000869751f, -0.027725220f, -0.016708374f, -0.152069092f, -0.246505737f, -1.048156738f,\r
+ 0.890090942f, 0.069168091f, 0.137298584f, -0.012115479f, 0.019577026f, -0.003051758f,\r
+ 0.001205444f, 0.000045776f, -0.000137329f, -0.002120972f, 0.000442505f, -0.028533936f,\r
+ -0.021179199f, -0.151596069f, -0.271591187f, -1.063217163f, 0.866363525f, 0.050354004f,\r
+ 0.133590698f, -0.014801025f, 0.018463135f, -0.003173828f, 0.001113892f, 0.000045776f,\r
+ -0.000152588f, -0.002243042f, -0.000030518f, -0.029281616f, -0.025817871f, -0.150497437f,\r
+ -0.297210693f, -1.077117920f, 0.841949463f, 0.032379150f, 0.129562378f, -0.017257690f,\r
+ 0.017349243f, -0.003280640f, 0.001037598f, 0.000030518f, -0.000167847f, -0.002349854f,\r
+ -0.000549316f, -0.029937744f, -0.030609131f, -0.148773193f, -0.323318481f, -1.089782715f,\r
+ 0.816864014f, 0.015228271f, 0.125259399f, -0.019531250f, 0.016235352f, -0.003372192f,\r
+ 0.000961304f, 0.000030518f, -0.000198364f, -0.002456665f, -0.001098633f, -0.030532837f,\r
+ -0.035552979f, -0.146362305f, -0.349868774f, -1.101211548f, 0.791213989f, -0.001068115f,\r
+ 0.120697021f, -0.021575928f, 0.015121460f, -0.003417969f, 0.000885010f, 0.000030518f,\r
+ -0.000213623f, -0.002578735f, -0.001693726f, -0.031005859f, -0.040634155f, -0.143264771f,\r
+ -0.376800537f, -1.111373901f, 0.765029907f, -0.016510010f, 0.115921021f, -0.023422241f,\r
+ 0.014022827f, -0.003463745f, 0.000808716f, 0.000030518f, -0.000244141f, -0.002685547f,\r
+ -0.002334595f, -0.031387329f, -0.045837402f, -0.139450073f, -0.404083252f, -1.120223999f,\r
+ 0.738372803f, -0.031082153f, 0.110946655f, -0.025085449f, 0.012939453f, -0.003479004f,\r
+ 0.000747681f, 0.000015259f, -0.000259399f, -0.002792358f, -0.003005981f, -0.031661987f,\r
+ -0.051132202f, -0.134887695f, -0.431655884f, -1.127746582f, 0.711318970f, -0.044784546f,\r
+ 0.105819702f, -0.026535034f, 0.011886597f, -0.003479004f, 0.000686646f, 0.000015259f,\r
+ -0.000289917f, -0.002899170f, -0.003723145f, -0.031814575f, -0.056533813f, -0.129577637f,\r
+ -0.459472656f, -1.133926392f, 0.683914185f, -0.057617188f, 0.100540161f, -0.027801514f,\r
+ 0.010848999f, -0.003463745f, 0.000625610f, 0.000015259f, -0.000320435f, -0.002990723f,\r
+ -0.004486084f, -0.031845093f, -0.061996460f, -0.123474121f, -0.487472534f, -1.138763428f,\r
+ 0.656219482f, -0.069595337f, 0.095169067f, -0.028884888f, 0.009841919f, -0.003433228f,\r
+ 0.000579834f, 0.000015259f, -0.000366211f, -0.003082275f, -0.005294800f, -0.031738281f,\r
+ -0.067520142f, -0.116577148f, -0.515609741f, -1.142211914f, 0.628295898f, -0.080688477f,\r
+ 0.089706421f, -0.029785156f, 0.008865356f, -0.003387451f, 0.000534058f, 0.000015259f,\r
+ -0.000396729f, -0.003173828f, -0.006118774f, -0.031478882f, -0.073059082f, -0.108856201f,\r
+ -0.543823242f, -1.144287109f, 0.600219727f, -0.090927124f, 0.084182739f, -0.030517578f,\r
+ 0.007919312f, -0.003326416f, 0.000473022f, 0.000015259f };\r
+\r
+}\r