Robust/src/Benchmarks/SSJava/MP3DecoderInfer/SynthesisFilter.java

   1 /*\r
   2  * 11/19/04 1.0 moved to LGPL.\r
   3  * \r
   4  * 04/01/00 Fixes for running under build 23xx Microsoft JVM. mdm.\r
   5  * \r
   6  * 19/12/99 Performance improvements to compute_pcm_samples().  \r
   7  *                      Mat McGowan. mdm@techie.com. \r
   8  *\r
   9  * 16/02/99 Java Conversion by E.B , javalayer@javazoom.net\r
  10  *\r
  11  *  @(#) synthesis_filter.h 1.8, last edit: 6/15/94 16:52:00\r
  12  *  @(#) Copyright (C) 1993, 1994 Tobias Bading (bading@cs.tu-berlin.de)\r
  13  *  @(#) Berlin University of Technology\r
  14  *\r
  15  *-----------------------------------------------------------------------\r
  16  *   This program is free software; you can redistribute it and/or modify\r
  17  *   it under the terms of the GNU Library General Public License as published\r
  18  *   by the Free Software Foundation; either version 2 of the License, or\r
  19  *   (at your option) any later version.\r
  20  *\r
  21  *   This program is distributed in the hope that it will be useful,\r
  22  *   but WITHOUT ANY WARRANTY; without even the implied warranty of\r
  23  *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the\r
  24  *   GNU Library General Public License for more details.\r
  25  *\r
  26  *   You should have received a copy of the GNU Library General Public\r
  27  *   License along with this program; if not, write to the Free Software\r
  28  *   Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.\r
  29  *----------------------------------------------------------------------\r
  30  */\r
  31 \r
  32 /**\r
  33  * A class for the synthesis filter bank. This class does a fast downsampling\r
  34  * from 32, 44.1 or 48 kHz to 8 kHz, if ULAW is defined. Frequencies above 4 kHz\r
  35  * are removed by ignoring higher subbands.\r
  36  */\r
  37 \r
  38 final class SynthesisFilter {\r
  39 \r
  40   public int vidx;\r
  41 \r
  42   public float[] v1;\r
  43 \r
  44   public float[] v2;\r
  45 \r
  46   public float[] prev1;\r
  47 \r
  48   public float[] prev2;\r
  49 \r
  50   // private float[] actual_v; // v1 or v2\r
  51 \r
  52   public int actual_write_pos; // 0-15\r
  53 \r
  54   private float[] samples; // 32 new subband samples\r
  55 \r
  56   public final int channel;\r
  57 \r
  58   public final float scalefactor;\r
  59 \r
  60   private float[] eq;\r
  61 \r
  62   /**\r
  63    * Quality value for controlling CPU usage/quality tradeoff.\r
  64    */\r
  65   /*\r
  66    * private int quality;\r
  67    * \r
  68    * private int v_inc;\r
  69    * \r
  70    * \r
  71    * \r
  72    * public static final int HIGH_QUALITY = 1; public static final int\r
  73    * MEDIUM_QUALITY = 2; public static final int LOW_QUALITY = 4;\r
  74    */\r
  75 \r
  76   /**\r
  77    * Contructor. The scalefactor scales the calculated float pcm samples to\r
  78    * short values (raw pcm samples are in [-1.0, 1.0], if no violations occur).\r
  79    */\r
  80   public SynthesisFilter(int channelnumber, float factor, float[] eq0) {\r
  81 \r
  82     vidx = 1;\r
  83     d16 = splitArray(d, 16);\r
  84 \r
  85     v1 = new float[512];\r
  86     v2 = new float[512];\r
  87     prev1 = new float[512];\r
  88     prev2 = new float[512];\r
  89     samples = new float[32];\r
  90     channel = channelnumber;\r
  91     scalefactor = factor;\r
  92     // setEQ(eq);\r
  93     // setQuality(HIGH_QUALITY);\r
  94 \r
  95     if (eq == null) {\r
  96       eq = new float[32];\r
  97       for (int i = 0; i < 32; i++)\r
  98         eq[i] = 1.0f;\r
  99     }\r
 100     if (eq.length < 32) {\r
 101       throw new IllegalArgumentException("eq0");\r
 102     }\r
 103 \r
 104     // reset();\r
 105 \r
 106     for (int p = 0; p < 512; p++)\r
 107       v1[p] = v2[p] = 0.0f;\r
 108 \r
 109     for (int p2 = 0; p2 < 32; p2++)\r
 110       samples[p2] = 0.0f;\r
 111 \r
 112     // actual_v = v1;\r
 113     actual_write_pos = 15;\r
 114 \r
 115   }\r
 116 \r
 117   /*\r
 118    * private void setQuality(int quality0) { switch (quality0) { case\r
 119    * HIGH_QUALITY: case MEDIUM_QUALITY: case LOW_QUALITY: v_inc = 16 * quality0;\r
 120    * quality = quality0; break; default : throw new\r
 121    * IllegalArgumentException("Unknown quality value"); } }\r
 122    * \r
 123    * public int getQuality() { return quality; }\r
 124    */\r
 125 \r
 126   /**\r
 127    * Inject Sample.\r
 128    */\r
 129   public void input_sample(float sample, int subbandnumber) {\r
 130     samples[subbandnumber] = eq[subbandnumber] * sample;\r
 131   }\r
 132 \r
 133   public void input_samples(float[] s) {\r
 134     TERMINATE: for (int i = 31; i >= 0; i--) {\r
 135       samples[i] = s[i] * eq[i];\r
 136     }\r
 137   }\r
 138 \r
 139   private void compute_new_v2_v1() {\r
 140 \r
 141     float new_v0 = 0.0f;\r
 142     float new_v1 = 0.0f;\r
 143     float new_v2 = 0.0f;\r
 144     float new_v3 = 0.0f;\r
 145     float new_v4 = 0.0f;\r
 146     float new_v5 = 0.0f;\r
 147     float new_v6 = 0.0f;\r
 148     float new_v7 = 0.0f;\r
 149     float new_v8 = 0.0f;\r
 150     float new_v9 = 0.0f;\r
 151     float new_v10 = 0.0f;\r
 152     float new_v11 = 0.0f;\r
 153     float new_v12 = 0.0f;\r
 154     float new_v13 = 0.0f;\r
 155     float new_v14 = 0.0f;\r
 156     float new_v15 = 0.0f;\r
 157     float new_v16 = 0.0f;\r
 158     float new_v17 = 0.0f;\r
 159     float new_v18 = 0.0f;\r
 160     float new_v19 = 0.0f;\r
 161     float new_v20 = 0.0f;\r
 162     float new_v21 = 0.0f;\r
 163     float new_v22 = 0.0f;\r
 164     float new_v23 = 0.0f;\r
 165     float new_v24 = 0.0f;\r
 166     float new_v25 = 0.0f;\r
 167     float new_v26 = 0.0f;\r
 168     float new_v27 = 0.0f;\r
 169     float new_v28 = 0.0f;\r
 170     float new_v29 = 0.0f;\r
 171     float new_v30 = 0.0f;\r
 172     float new_v31 = 0.0f;\r
 173 \r
 174     // float[] new_v = new float[32]; // new V[0-15] and V[33-48] of Figure\r
 175     // 3-A.2 in ISO DIS 11172-3\r
 176     // float[] p = new float[16];\r
 177     // float[] pp = new float[16];\r
 178 \r
 179     // float[] s = samples;\r
 180 \r
 181     float s0 = samples[0];\r
 182     float s1 = samples[1];\r
 183     float s2 = samples[2];\r
 184     float s3 = samples[3];\r
 185     float s4 = samples[4];\r
 186     float s5 = samples[5];\r
 187     float s6 = samples[6];\r
 188     float s7 = samples[7];\r
 189     float s8 = samples[8];\r
 190     float s9 = samples[9];\r
 191     float s10 = samples[10];\r
 192     float s11 = samples[11];\r
 193     float s12 = samples[12];\r
 194     float s13 = samples[13];\r
 195     float s14 = samples[14];\r
 196     float s15 = samples[15];\r
 197     float s16 = samples[16];\r
 198     float s17 = samples[17];\r
 199     float s18 = samples[18];\r
 200     float s19 = samples[19];\r
 201     float s20 = samples[20];\r
 202     float s21 = samples[21];\r
 203     float s22 = samples[22];\r
 204     float s23 = samples[23];\r
 205     float s24 = samples[24];\r
 206     float s25 = samples[25];\r
 207     float s26 = samples[26];\r
 208     float s27 = samples[27];\r
 209     float s28 = samples[28];\r
 210     float s29 = samples[29];\r
 211     float s30 = samples[30];\r
 212     float s31 = samples[31];\r
 213 \r
 214     float p0 = s0 + s31;\r
 215     float p1 = s1 + s30;\r
 216     float p2 = s2 + s29;\r
 217     float p3 = s3 + s28;\r
 218     float p4 = s4 + s27;\r
 219     float p5 = s5 + s26;\r
 220     float p6 = s6 + s25;\r
 221     float p7 = s7 + s24;\r
 222     float p8 = s8 + s23;\r
 223     float p9 = s9 + s22;\r
 224     float p10 = s10 + s21;\r
 225     float p11 = s11 + s20;\r
 226     float p12 = s12 + s19;\r
 227     float p13 = s13 + s18;\r
 228     float p14 = s14 + s17;\r
 229     float p15 = s15 + s16;\r
 230 \r
 231     float pp0 = p0 + p15;\r
 232     float pp1 = p1 + p14;\r
 233     float pp2 = p2 + p13;\r
 234     float pp3 = p3 + p12;\r
 235     float pp4 = p4 + p11;\r
 236     float pp5 = p5 + p10;\r
 237     float pp6 = p6 + p9;\r
 238     float pp7 = p7 + p8;\r
 239     float pp8 = (p0 - p15) * cos1_32;\r
 240     float pp9 = (p1 - p14) * cos3_32;\r
 241     float pp10 = (p2 - p13) * cos5_32;\r
 242     float pp11 = (p3 - p12) * cos7_32;\r
 243     float pp12 = (p4 - p11) * cos9_32;\r
 244     float pp13 = (p5 - p10) * cos11_32;\r
 245     float pp14 = (p6 - p9) * cos13_32;\r
 246     float pp15 = (p7 - p8) * cos15_32;\r
 247 \r
 248     p0 = pp0 + pp7;\r
 249     p1 = pp1 + pp6;\r
 250     p2 = pp2 + pp5;\r
 251     p3 = pp3 + pp4;\r
 252     p4 = (pp0 - pp7) * cos1_16;\r
 253     p5 = (pp1 - pp6) * cos3_16;\r
 254     p6 = (pp2 - pp5) * cos5_16;\r
 255     p7 = (pp3 - pp4) * cos7_16;\r
 256     p8 = pp8 + pp15;\r
 257     p9 = pp9 + pp14;\r
 258     p10 = pp10 + pp13;\r
 259     p11 = pp11 + pp12;\r
 260     p12 = (pp8 - pp15) * cos1_16;\r
 261     p13 = (pp9 - pp14) * cos3_16;\r
 262     p14 = (pp10 - pp13) * cos5_16;\r
 263     p15 = (pp11 - pp12) * cos7_16;\r
 264 \r
 265     pp0 = p0 + p3;\r
 266     pp1 = p1 + p2;\r
 267     pp2 = (p0 - p3) * cos1_8;\r
 268     pp3 = (p1 - p2) * cos3_8;\r
 269     pp4 = p4 + p7;\r
 270     pp5 = p5 + p6;\r
 271     pp6 = (p4 - p7) * cos1_8;\r
 272     pp7 = (p5 - p6) * cos3_8;\r
 273     pp8 = p8 + p11;\r
 274     pp9 = p9 + p10;\r
 275     pp10 = (p8 - p11) * cos1_8;\r
 276     pp11 = (p9 - p10) * cos3_8;\r
 277     pp12 = p12 + p15;\r
 278     pp13 = p13 + p14;\r
 279     pp14 = (p12 - p15) * cos1_8;\r
 280     pp15 = (p13 - p14) * cos3_8;\r
 281 \r
 282     p0 = pp0 + pp1;\r
 283     p1 = (pp0 - pp1) * cos1_4;\r
 284     p2 = pp2 + pp3;\r
 285     p3 = (pp2 - pp3) * cos1_4;\r
 286     p4 = pp4 + pp5;\r
 287     p5 = (pp4 - pp5) * cos1_4;\r
 288     p6 = pp6 + pp7;\r
 289     p7 = (pp6 - pp7) * cos1_4;\r
 290     p8 = pp8 + pp9;\r
 291     p9 = (pp8 - pp9) * cos1_4;\r
 292     p10 = pp10 + pp11;\r
 293     p11 = (pp10 - pp11) * cos1_4;\r
 294     p12 = pp12 + pp13;\r
 295     p13 = (pp12 - pp13) * cos1_4;\r
 296     p14 = pp14 + pp15;\r
 297     p15 = (pp14 - pp15) * cos1_4;\r
 298 \r
 299     // this is pretty insane coding\r
 300     float tmp1;\r
 301     new_v19/* 36-17 */= -(new_v4 = (new_v12 = p7) + p5) - p6;\r
 302     new_v27/* 44-17 */= -p6 - p7 - p4;\r
 303     new_v6 = (new_v10 = (new_v14 = p15) + p11) + p13;\r
 304     new_v17/* 34-17 */= -(new_v2 = p15 + p13 + p9) - p14;\r
 305     new_v21/* 38-17 */= (tmp1 = -p14 - p15 - p10 - p11) - p13;\r
 306     new_v29/* 46-17 */= -p14 - p15 - p12 - p8;\r
 307     new_v25/* 42-17 */= tmp1 - p12;\r
 308     new_v31/* 48-17 */= -p0;\r
 309     new_v0 = p1;\r
 310     new_v23/* 40-17 */= -(new_v8 = p3) - p2;\r
 311 \r
 312     p0 = (s0 - s31) * cos1_64;\r
 313     p1 = (s1 - s30) * cos3_64;\r
 314     p2 = (s2 - s29) * cos5_64;\r
 315     p3 = (s3 - s28) * cos7_64;\r
 316     p4 = (s4 - s27) * cos9_64;\r
 317     p5 = (s5 - s26) * cos11_64;\r
 318     p6 = (s6 - s25) * cos13_64;\r
 319     p7 = (s7 - s24) * cos15_64;\r
 320     p8 = (s8 - s23) * cos17_64;\r
 321     p9 = (s9 - s22) * cos19_64;\r
 322     p10 = (s10 - s21) * cos21_64;\r
 323     p11 = (s11 - s20) * cos23_64;\r
 324     p12 = (s12 - s19) * cos25_64;\r
 325     p13 = (s13 - s18) * cos27_64;\r
 326     p14 = (s14 - s17) * cos29_64;\r
 327     p15 = (s15 - s16) * cos31_64;\r
 328 \r
 329     pp0 = p0 + p15;\r
 330     pp1 = p1 + p14;\r
 331     pp2 = p2 + p13;\r
 332     pp3 = p3 + p12;\r
 333     pp4 = p4 + p11;\r
 334     pp5 = p5 + p10;\r
 335     pp6 = p6 + p9;\r
 336     pp7 = p7 + p8;\r
 337     pp8 = (p0 - p15) * cos1_32;\r
 338     pp9 = (p1 - p14) * cos3_32;\r
 339     pp10 = (p2 - p13) * cos5_32;\r
 340     pp11 = (p3 - p12) * cos7_32;\r
 341     pp12 = (p4 - p11) * cos9_32;\r
 342     pp13 = (p5 - p10) * cos11_32;\r
 343     pp14 = (p6 - p9) * cos13_32;\r
 344     pp15 = (p7 - p8) * cos15_32;\r
 345 \r
 346     p0 = pp0 + pp7;\r
 347     p1 = pp1 + pp6;\r
 348     p2 = pp2 + pp5;\r
 349     p3 = pp3 + pp4;\r
 350     p4 = (pp0 - pp7) * cos1_16;\r
 351     p5 = (pp1 - pp6) * cos3_16;\r
 352     p6 = (pp2 - pp5) * cos5_16;\r
 353     p7 = (pp3 - pp4) * cos7_16;\r
 354     p8 = pp8 + pp15;\r
 355     p9 = pp9 + pp14;\r
 356     p10 = pp10 + pp13;\r
 357     p11 = pp11 + pp12;\r
 358     p12 = (pp8 - pp15) * cos1_16;\r
 359     p13 = (pp9 - pp14) * cos3_16;\r
 360     p14 = (pp10 - pp13) * cos5_16;\r
 361     p15 = (pp11 - pp12) * cos7_16;\r
 362 \r
 363     pp0 = p0 + p3;\r
 364     pp1 = p1 + p2;\r
 365     pp2 = (p0 - p3) * cos1_8;\r
 366     pp3 = (p1 - p2) * cos3_8;\r
 367     pp4 = p4 + p7;\r
 368     pp5 = p5 + p6;\r
 369     pp6 = (p4 - p7) * cos1_8;\r
 370     pp7 = (p5 - p6) * cos3_8;\r
 371     pp8 = p8 + p11;\r
 372     pp9 = p9 + p10;\r
 373     pp10 = (p8 - p11) * cos1_8;\r
 374     pp11 = (p9 - p10) * cos3_8;\r
 375     pp12 = p12 + p15;\r
 376     pp13 = p13 + p14;\r
 377     pp14 = (p12 - p15) * cos1_8;\r
 378     pp15 = (p13 - p14) * cos3_8;\r
 379 \r
 380     p0 = pp0 + pp1;\r
 381     p1 = (pp0 - pp1) * cos1_4;\r
 382     p2 = pp2 + pp3;\r
 383     p3 = (pp2 - pp3) * cos1_4;\r
 384     p4 = pp4 + pp5;\r
 385     p5 = (pp4 - pp5) * cos1_4;\r
 386     p6 = pp6 + pp7;\r
 387     p7 = (pp6 - pp7) * cos1_4;\r
 388     p8 = pp8 + pp9;\r
 389     p9 = (pp8 - pp9) * cos1_4;\r
 390     p10 = pp10 + pp11;\r
 391     p11 = (pp10 - pp11) * cos1_4;\r
 392     p12 = pp12 + pp13;\r
 393     p13 = (pp12 - pp13) * cos1_4;\r
 394     p14 = pp14 + pp15;\r
 395     p15 = (pp14 - pp15) * cos1_4;\r
 396 \r
 397     // manually doing something that a compiler should handle sucks\r
 398     // coding like this is hard to read\r
 399     float tmp2;\r
 400     new_v5 = (new_v11 = (new_v13 = (new_v15 = p15) + p7) + p11) + p5 + p13;\r
 401     new_v7 = (new_v9 = p15 + p11 + p3) + p13;\r
 402     new_v16/* 33-17 */= -(new_v1 = (tmp1 = p13 + p15 + p9) + p1) - p14;\r
 403     new_v18/* 35-17 */= -(new_v3 = tmp1 + p5 + p7) - p6 - p14;\r
 404 \r
 405     new_v22/* 39-17 */= (tmp1 = -p10 - p11 - p14 - p15) - p13 - p2 - p3;\r
 406     new_v20/* 37-17 */= tmp1 - p13 - p5 - p6 - p7;\r
 407     new_v24/* 41-17 */= tmp1 - p12 - p2 - p3;\r
 408     new_v26/* 43-17 */= tmp1 - p12 - (tmp2 = p4 + p6 + p7);\r
 409     new_v30/* 47-17 */= (tmp1 = -p8 - p12 - p14 - p15) - p0;\r
 410     new_v28/* 45-17 */= tmp1 - tmp2;\r
 411 \r
 412     // insert V[0-15] (== new_v[0-15]) into actual v:\r
 413     // float[] x2 = actual_v + actual_write_pos;\r
 414     // float dest[] = actual_v; v2\r
 415 \r
 416     int pos = actual_write_pos;\r
 417 \r
 418     v2[0 + pos] = new_v0;\r
 419     v2[16 + pos] = new_v1;\r
 420     v2[32 + pos] = new_v2;\r
 421     v2[48 + pos] = new_v3;\r
 422     v2[64 + pos] = new_v4;\r
 423     v2[80 + pos] = new_v5;\r
 424     v2[96 + pos] = new_v6;\r
 425     v2[112 + pos] = new_v7;\r
 426     v2[128 + pos] = new_v8;\r
 427     v2[144 + pos] = new_v9;\r
 428     v2[160 + pos] = new_v10;\r
 429     v2[176 + pos] = new_v11;\r
 430     v2[192 + pos] = new_v12;\r
 431     v2[208 + pos] = new_v13;\r
 432     v2[224 + pos] = new_v14;\r
 433     v2[240 + pos] = new_v15;\r
 434 \r
 435     // V[16] is always 0.0:\r
 436     v2[256 + pos] = 0.0f;\r
 437 \r
 438     // insert V[17-31] (== -new_v[15-1]) into actual v:\r
 439     v2[272 + pos] = -new_v15;\r
 440     v2[288 + pos] = -new_v14;\r
 441     v2[304 + pos] = -new_v13;\r
 442     v2[320 + pos] = -new_v12;\r
 443     v2[336 + pos] = -new_v11;\r
 444     v2[352 + pos] = -new_v10;\r
 445     v2[368 + pos] = -new_v9;\r
 446     v2[384 + pos] = -new_v8;\r
 447     v2[400 + pos] = -new_v7;\r
 448     v2[416 + pos] = -new_v6;\r
 449     v2[432 + pos] = -new_v5;\r
 450     v2[448 + pos] = -new_v4;\r
 451     v2[464 + pos] = -new_v3;\r
 452     v2[480 + pos] = -new_v2;\r
 453     v2[496 + pos] = -new_v1;\r
 454 \r
 455     // insert V[32] (== -new_v[0]) into other v:\r
 456     // dest = (actual_v == v1) ? v2 : v1;\r
 457 \r
 458     v1[0 + pos] = -new_v0;\r
 459     // insert V[33-48] (== new_v[16-31]) into other v:\r
 460     v1[16 + pos] = new_v16;\r
 461     v1[32 + pos] = new_v17;\r
 462     v1[48 + pos] = new_v18;\r
 463     v1[64 + pos] = new_v19;\r
 464     v1[80 + pos] = new_v20;\r
 465     v1[96 + pos] = new_v21;\r
 466     v1[112 + pos] = new_v22;\r
 467     v1[128 + pos] = new_v23;\r
 468     v1[144 + pos] = new_v24;\r
 469     v1[160 + pos] = new_v25;\r
 470     v1[176 + pos] = new_v26;\r
 471     v1[192 + pos] = new_v27;\r
 472     v1[208 + pos] = new_v28;\r
 473     v1[224 + pos] = new_v29;\r
 474     v1[240 + pos] = new_v30;\r
 475     v1[256 + pos] = new_v31;\r
 476 \r
 477     // insert V[49-63] (== new_v[30-16]) into other v:\r
 478     v1[272 + pos] = new_v30;\r
 479     v1[288 + pos] = new_v29;\r
 480     v1[304 + pos] = new_v28;\r
 481     v1[320 + pos] = new_v27;\r
 482     v1[336 + pos] = new_v26;\r
 483     v1[352 + pos] = new_v25;\r
 484     v1[368 + pos] = new_v24;\r
 485     v1[384 + pos] = new_v23;\r
 486     v1[400 + pos] = new_v22;\r
 487     v1[416 + pos] = new_v21;\r
 488     v1[432 + pos] = new_v20;\r
 489     v1[448 + pos] = new_v19;\r
 490     v1[464 + pos] = new_v18;\r
 491     v1[480 + pos] = new_v17;\r
 492     v1[496 + pos] = new_v16;\r
 493 \r
 494     /*\r
 495      * setup PREV\r
 496      */\r
 497 \r
 498     prev2[0 + pos] = new_v0;\r
 499     prev2[16 + pos] = new_v1;\r
 500     prev2[32 + pos] = new_v2;\r
 501     prev2[48 + pos] = new_v3;\r
 502     prev2[64 + pos] = new_v4;\r
 503     prev2[80 + pos] = new_v5;\r
 504     prev2[96 + pos] = new_v6;\r
 505     prev2[112 + pos] = new_v7;\r
 506     prev2[128 + pos] = new_v8;\r
 507     prev2[144 + pos] = new_v9;\r
 508     prev2[160 + pos] = new_v10;\r
 509     prev2[176 + pos] = new_v11;\r
 510     prev2[192 + pos] = new_v12;\r
 511     prev2[208 + pos] = new_v13;\r
 512     prev2[224 + pos] = new_v14;\r
 513     prev2[240 + pos] = new_v15;\r
 514 \r
 515     // V[16] is always 0.0:\r
 516     prev2[256 + pos] = 0.0f;\r
 517 \r
 518     // insert V[17-31] (== -new_v[15-1]) into actual v:\r
 519     prev2[272 + pos] = -new_v15;\r
 520     prev2[288 + pos] = -new_v14;\r
 521     prev2[304 + pos] = -new_v13;\r
 522     prev2[320 + pos] = -new_v12;\r
 523     prev2[336 + pos] = -new_v11;\r
 524     prev2[352 + pos] = -new_v10;\r
 525     prev2[368 + pos] = -new_v9;\r
 526     prev2[384 + pos] = -new_v8;\r
 527     prev2[400 + pos] = -new_v7;\r
 528     prev2[416 + pos] = -new_v6;\r
 529     prev2[432 + pos] = -new_v5;\r
 530     prev2[448 + pos] = -new_v4;\r
 531     prev2[464 + pos] = -new_v3;\r
 532     prev2[480 + pos] = -new_v2;\r
 533     prev2[496 + pos] = -new_v1;\r
 534 \r
 535     // insert V[32] (== -new_v[0]) into other v:\r
 536     // dest = (actual_v == v1) ? v2 : v1;\r
 537 \r
 538     prev1[0 + pos] = -new_v0;\r
 539     // insert V[33-48] (== new_v[16-31]) into other v:\r
 540     prev1[16 + pos] = new_v16;\r
 541     prev1[32 + pos] = new_v17;\r
 542     prev1[48 + pos] = new_v18;\r
 543     prev1[64 + pos] = new_v19;\r
 544     prev1[80 + pos] = new_v20;\r
 545     prev1[96 + pos] = new_v21;\r
 546     prev1[112 + pos] = new_v22;\r
 547     prev1[128 + pos] = new_v23;\r
 548     prev1[144 + pos] = new_v24;\r
 549     prev1[160 + pos] = new_v25;\r
 550     prev1[176 + pos] = new_v26;\r
 551     prev1[192 + pos] = new_v27;\r
 552     prev1[208 + pos] = new_v28;\r
 553     prev1[224 + pos] = new_v29;\r
 554     prev1[240 + pos] = new_v30;\r
 555     prev1[256 + pos] = new_v31;\r
 556 \r
 557     // insert V[49-63] (== new_v[30-16]) into other v:\r
 558     prev1[272 + pos] = new_v30;\r
 559     prev1[288 + pos] = new_v29;\r
 560     prev1[304 + pos] = new_v28;\r
 561     prev1[320 + pos] = new_v27;\r
 562     prev1[336 + pos] = new_v26;\r
 563     prev1[352 + pos] = new_v25;\r
 564     prev1[368 + pos] = new_v24;\r
 565     prev1[384 + pos] = new_v23;\r
 566     prev1[400 + pos] = new_v22;\r
 567     prev1[416 + pos] = new_v21;\r
 568     prev1[432 + pos] = new_v20;\r
 569     prev1[448 + pos] = new_v19;\r
 570     prev1[464 + pos] = new_v18;\r
 571     prev1[480 + pos] = new_v17;\r
 572     prev1[496 + pos] = new_v16;\r
 573   }\r
 574 \r
 575   private void compute_new_v1_v2() {\r
 576 \r
 577     float new_v0 = 0.0f;\r
 578     float new_v1 = 0.0f;\r
 579     float new_v2 = 0.0f;\r
 580     float new_v3 = 0.0f;\r
 581     float new_v4 = 0.0f;\r
 582     float new_v5 = 0.0f;\r
 583     float new_v6 = 0.0f;\r
 584     float new_v7 = 0.0f;\r
 585     float new_v8 = 0.0f;\r
 586     float new_v9 = 0.0f;\r
 587     float new_v10 = 0.0f;\r
 588     float new_v11 = 0.0f;\r
 589     float new_v12 = 0.0f;\r
 590     float new_v13 = 0.0f;\r
 591     float new_v14 = 0.0f;\r
 592     float new_v15 = 0.0f;\r
 593     float new_v16 = 0.0f;\r
 594     float new_v17 = 0.0f;\r
 595     float new_v18 = 0.0f;\r
 596     float new_v19 = 0.0f;\r
 597     float new_v20 = 0.0f;\r
 598     float new_v21 = 0.0f;\r
 599     float new_v22 = 0.0f;\r
 600     float new_v23 = 0.0f;\r
 601     float new_v24 = 0.0f;\r
 602     float new_v25 = 0.0f;\r
 603     float new_v26 = 0.0f;\r
 604     float new_v27 = 0.0f;\r
 605     float new_v28 = 0.0f;\r
 606     float new_v29 = 0.0f;\r
 607     float new_v30 = 0.0f;\r
 608     float new_v31 = 0.0f;\r
 609 \r
 610     // float[] new_v = new float[32]; // new V[0-15] and V[33-48] of Figure\r
 611     // 3-A.2 in ISO DIS 11172-3\r
 612     // float[] p = new float[16];\r
 613     // float[] pp = new float[16];\r
 614 \r
 615     // float[] s = samples;\r
 616 \r
 617     float s0 = samples[0];\r
 618     float s1 = samples[1];\r
 619     float s2 = samples[2];\r
 620     float s3 = samples[3];\r
 621     float s4 = samples[4];\r
 622     float s5 = samples[5];\r
 623     float s6 = samples[6];\r
 624     float s7 = samples[7];\r
 625     float s8 = samples[8];\r
 626     float s9 = samples[9];\r
 627     float s10 = samples[10];\r
 628     float s11 = samples[11];\r
 629     float s12 = samples[12];\r
 630     float s13 = samples[13];\r
 631     float s14 = samples[14];\r
 632     float s15 = samples[15];\r
 633     float s16 = samples[16];\r
 634     float s17 = samples[17];\r
 635     float s18 = samples[18];\r
 636     float s19 = samples[19];\r
 637     float s20 = samples[20];\r
 638     float s21 = samples[21];\r
 639     float s22 = samples[22];\r
 640     float s23 = samples[23];\r
 641     float s24 = samples[24];\r
 642     float s25 = samples[25];\r
 643     float s26 = samples[26];\r
 644     float s27 = samples[27];\r
 645     float s28 = samples[28];\r
 646     float s29 = samples[29];\r
 647     float s30 = samples[30];\r
 648     float s31 = samples[31];\r
 649 \r
 650     float p0 = s0 + s31;\r
 651     float p1 = s1 + s30;\r
 652     float p2 = s2 + s29;\r
 653     float p3 = s3 + s28;\r
 654     float p4 = s4 + s27;\r
 655     float p5 = s5 + s26;\r
 656     float p6 = s6 + s25;\r
 657     float p7 = s7 + s24;\r
 658     float p8 = s8 + s23;\r
 659     float p9 = s9 + s22;\r
 660     float p10 = s10 + s21;\r
 661     float p11 = s11 + s20;\r
 662     float p12 = s12 + s19;\r
 663     float p13 = s13 + s18;\r
 664     float p14 = s14 + s17;\r
 665     float p15 = s15 + s16;\r
 666 \r
 667     float pp0 = p0 + p15;\r
 668     float pp1 = p1 + p14;\r
 669     float pp2 = p2 + p13;\r
 670     float pp3 = p3 + p12;\r
 671     float pp4 = p4 + p11;\r
 672     float pp5 = p5 + p10;\r
 673     float pp6 = p6 + p9;\r
 674     float pp7 = p7 + p8;\r
 675     float pp8 = (p0 - p15) * cos1_32;\r
 676     float pp9 = (p1 - p14) * cos3_32;\r
 677     float pp10 = (p2 - p13) * cos5_32;\r
 678     float pp11 = (p3 - p12) * cos7_32;\r
 679     float pp12 = (p4 - p11) * cos9_32;\r
 680     float pp13 = (p5 - p10) * cos11_32;\r
 681     float pp14 = (p6 - p9) * cos13_32;\r
 682     float pp15 = (p7 - p8) * cos15_32;\r
 683 \r
 684     p0 = pp0 + pp7;\r
 685     p1 = pp1 + pp6;\r
 686     p2 = pp2 + pp5;\r
 687     p3 = pp3 + pp4;\r
 688     p4 = (pp0 - pp7) * cos1_16;\r
 689     p5 = (pp1 - pp6) * cos3_16;\r
 690     p6 = (pp2 - pp5) * cos5_16;\r
 691     p7 = (pp3 - pp4) * cos7_16;\r
 692     p8 = pp8 + pp15;\r
 693     p9 = pp9 + pp14;\r
 694     p10 = pp10 + pp13;\r
 695     p11 = pp11 + pp12;\r
 696     p12 = (pp8 - pp15) * cos1_16;\r
 697     p13 = (pp9 - pp14) * cos3_16;\r
 698     p14 = (pp10 - pp13) * cos5_16;\r
 699     p15 = (pp11 - pp12) * cos7_16;\r
 700 \r
 701     pp0 = p0 + p3;\r
 702     pp1 = p1 + p2;\r
 703     pp2 = (p0 - p3) * cos1_8;\r
 704     pp3 = (p1 - p2) * cos3_8;\r
 705     pp4 = p4 + p7;\r
 706     pp5 = p5 + p6;\r
 707     pp6 = (p4 - p7) * cos1_8;\r
 708     pp7 = (p5 - p6) * cos3_8;\r
 709     pp8 = p8 + p11;\r
 710     pp9 = p9 + p10;\r
 711     pp10 = (p8 - p11) * cos1_8;\r
 712     pp11 = (p9 - p10) * cos3_8;\r
 713     pp12 = p12 + p15;\r
 714     pp13 = p13 + p14;\r
 715     pp14 = (p12 - p15) * cos1_8;\r
 716     pp15 = (p13 - p14) * cos3_8;\r
 717 \r
 718     p0 = pp0 + pp1;\r
 719     p1 = (pp0 - pp1) * cos1_4;\r
 720     p2 = pp2 + pp3;\r
 721     p3 = (pp2 - pp3) * cos1_4;\r
 722     p4 = pp4 + pp5;\r
 723     p5 = (pp4 - pp5) * cos1_4;\r
 724     p6 = pp6 + pp7;\r
 725     p7 = (pp6 - pp7) * cos1_4;\r
 726     p8 = pp8 + pp9;\r
 727     p9 = (pp8 - pp9) * cos1_4;\r
 728     p10 = pp10 + pp11;\r
 729     p11 = (pp10 - pp11) * cos1_4;\r
 730     p12 = pp12 + pp13;\r
 731     p13 = (pp12 - pp13) * cos1_4;\r
 732     p14 = pp14 + pp15;\r
 733     p15 = (pp14 - pp15) * cos1_4;\r
 734 \r
 735     // this is pretty insane coding\r
 736     float tmp1;\r
 737     new_v19/* 36-17 */= -(new_v4 = (new_v12 = p7) + p5) - p6;\r
 738     new_v27/* 44-17 */= -p6 - p7 - p4;\r
 739     new_v6 = (new_v10 = (new_v14 = p15) + p11) + p13;\r
 740     new_v17/* 34-17 */= -(new_v2 = p15 + p13 + p9) - p14;\r
 741     new_v21/* 38-17 */= (tmp1 = -p14 - p15 - p10 - p11) - p13;\r
 742     new_v29/* 46-17 */= -p14 - p15 - p12 - p8;\r
 743     new_v25/* 42-17 */= tmp1 - p12;\r
 744     new_v31/* 48-17 */= -p0;\r
 745     new_v0 = p1;\r
 746     new_v23/* 40-17 */= -(new_v8 = p3) - p2;\r
 747 \r
 748     p0 = (s0 - s31) * cos1_64;\r
 749     p1 = (s1 - s30) * cos3_64;\r
 750     p2 = (s2 - s29) * cos5_64;\r
 751     p3 = (s3 - s28) * cos7_64;\r
 752     p4 = (s4 - s27) * cos9_64;\r
 753     p5 = (s5 - s26) * cos11_64;\r
 754     p6 = (s6 - s25) * cos13_64;\r
 755     p7 = (s7 - s24) * cos15_64;\r
 756     p8 = (s8 - s23) * cos17_64;\r
 757     p9 = (s9 - s22) * cos19_64;\r
 758     p10 = (s10 - s21) * cos21_64;\r
 759     p11 = (s11 - s20) * cos23_64;\r
 760     p12 = (s12 - s19) * cos25_64;\r
 761     p13 = (s13 - s18) * cos27_64;\r
 762     p14 = (s14 - s17) * cos29_64;\r
 763     p15 = (s15 - s16) * cos31_64;\r
 764 \r
 765     pp0 = p0 + p15;\r
 766     pp1 = p1 + p14;\r
 767     pp2 = p2 + p13;\r
 768     pp3 = p3 + p12;\r
 769     pp4 = p4 + p11;\r
 770     pp5 = p5 + p10;\r
 771     pp6 = p6 + p9;\r
 772     pp7 = p7 + p8;\r
 773     pp8 = (p0 - p15) * cos1_32;\r
 774     pp9 = (p1 - p14) * cos3_32;\r
 775     pp10 = (p2 - p13) * cos5_32;\r
 776     pp11 = (p3 - p12) * cos7_32;\r
 777     pp12 = (p4 - p11) * cos9_32;\r
 778     pp13 = (p5 - p10) * cos11_32;\r
 779     pp14 = (p6 - p9) * cos13_32;\r
 780     pp15 = (p7 - p8) * cos15_32;\r
 781 \r
 782     p0 = pp0 + pp7;\r
 783     p1 = pp1 + pp6;\r
 784     p2 = pp2 + pp5;\r
 785     p3 = pp3 + pp4;\r
 786     p4 = (pp0 - pp7) * cos1_16;\r
 787     p5 = (pp1 - pp6) * cos3_16;\r
 788     p6 = (pp2 - pp5) * cos5_16;\r
 789     p7 = (pp3 - pp4) * cos7_16;\r
 790     p8 = pp8 + pp15;\r
 791     p9 = pp9 + pp14;\r
 792     p10 = pp10 + pp13;\r
 793     p11 = pp11 + pp12;\r
 794     p12 = (pp8 - pp15) * cos1_16;\r
 795     p13 = (pp9 - pp14) * cos3_16;\r
 796     p14 = (pp10 - pp13) * cos5_16;\r
 797     p15 = (pp11 - pp12) * cos7_16;\r
 798 \r
 799     pp0 = p0 + p3;\r
 800     pp1 = p1 + p2;\r
 801     pp2 = (p0 - p3) * cos1_8;\r
 802     pp3 = (p1 - p2) * cos3_8;\r
 803     pp4 = p4 + p7;\r
 804     pp5 = p5 + p6;\r
 805     pp6 = (p4 - p7) * cos1_8;\r
 806     pp7 = (p5 - p6) * cos3_8;\r
 807     pp8 = p8 + p11;\r
 808     pp9 = p9 + p10;\r
 809     pp10 = (p8 - p11) * cos1_8;\r
 810     pp11 = (p9 - p10) * cos3_8;\r
 811     pp12 = p12 + p15;\r
 812     pp13 = p13 + p14;\r
 813     pp14 = (p12 - p15) * cos1_8;\r
 814     pp15 = (p13 - p14) * cos3_8;\r
 815 \r
 816     p0 = pp0 + pp1;\r
 817     p1 = (pp0 - pp1) * cos1_4;\r
 818     p2 = pp2 + pp3;\r
 819     p3 = (pp2 - pp3) * cos1_4;\r
 820     p4 = pp4 + pp5;\r
 821     p5 = (pp4 - pp5) * cos1_4;\r
 822     p6 = pp6 + pp7;\r
 823     p7 = (pp6 - pp7) * cos1_4;\r
 824     p8 = pp8 + pp9;\r
 825     p9 = (pp8 - pp9) * cos1_4;\r
 826     p10 = pp10 + pp11;\r
 827     p11 = (pp10 - pp11) * cos1_4;\r
 828     p12 = pp12 + pp13;\r
 829     p13 = (pp12 - pp13) * cos1_4;\r
 830     p14 = pp14 + pp15;\r
 831     p15 = (pp14 - pp15) * cos1_4;\r
 832 \r
 833     // manually doing something that a compiler should handle sucks\r
 834     // coding like this is hard to read\r
 835     float tmp2;\r
 836     new_v5 = (new_v11 = (new_v13 = (new_v15 = p15) + p7) + p11) + p5 + p13;\r
 837     new_v7 = (new_v9 = p15 + p11 + p3) + p13;\r
 838     new_v16/* 33-17 */= -(new_v1 = (tmp1 = p13 + p15 + p9) + p1) - p14;\r
 839     new_v18/* 35-17 */= -(new_v3 = tmp1 + p5 + p7) - p6 - p14;\r
 840 \r
 841     new_v22/* 39-17 */= (tmp1 = -p10 - p11 - p14 - p15) - p13 - p2 - p3;\r
 842     new_v20/* 37-17 */= tmp1 - p13 - p5 - p6 - p7;\r
 843     new_v24/* 41-17 */= tmp1 - p12 - p2 - p3;\r
 844     new_v26/* 43-17 */= tmp1 - p12 - (tmp2 = p4 + p6 + p7);\r
 845     new_v30/* 47-17 */= (tmp1 = -p8 - p12 - p14 - p15) - p0;\r
 846     new_v28/* 45-17 */= tmp1 - tmp2;\r
 847 \r
 848     // insert V[0-15] (== new_v[0-15]) into actual v:\r
 849     // float[] x2 = actual_v + actual_write_pos;\r
 850     // float dest[] = actual_v; actual_v=v1;\r
 851 \r
 852     int pos = actual_write_pos;\r
 853 \r
 854     v1[0 + pos] = new_v0;\r
 855     v1[16 + pos] = new_v1;\r
 856     v1[32 + pos] = new_v2;\r
 857     v1[48 + pos] = new_v3;\r
 858     v1[64 + pos] = new_v4;\r
 859     v1[80 + pos] = new_v5;\r
 860     v1[96 + pos] = new_v6;\r
 861     v1[112 + pos] = new_v7;\r
 862     v1[128 + pos] = new_v8;\r
 863     v1[144 + pos] = new_v9;\r
 864     v1[160 + pos] = new_v10;\r
 865     v1[176 + pos] = new_v11;\r
 866     v1[192 + pos] = new_v12;\r
 867     v1[208 + pos] = new_v13;\r
 868     v1[224 + pos] = new_v14;\r
 869     v1[240 + pos] = new_v15;\r
 870 \r
 871     // V[16] is always 0.0:\r
 872     v1[256 + pos] = 0.0f;\r
 873 \r
 874     // insert V[17-31] (== -new_v[15-1]) into actual v:\r
 875     v1[272 + pos] = -new_v15;\r
 876     v1[288 + pos] = -new_v14;\r
 877     v1[304 + pos] = -new_v13;\r
 878     v1[320 + pos] = -new_v12;\r
 879     v1[336 + pos] = -new_v11;\r
 880     v1[352 + pos] = -new_v10;\r
 881     v1[368 + pos] = -new_v9;\r
 882     v1[384 + pos] = -new_v8;\r
 883     v1[400 + pos] = -new_v7;\r
 884     v1[416 + pos] = -new_v6;\r
 885     v1[432 + pos] = -new_v5;\r
 886     v1[448 + pos] = -new_v4;\r
 887     v1[464 + pos] = -new_v3;\r
 888     v1[480 + pos] = -new_v2;\r
 889     v1[496 + pos] = -new_v1;\r
 890 \r
 891     // insert V[32] (== -new_v[0]) into other v:\r
 892     // dest = (actual_v == v1) ? v2 : v1;\r
 893 \r
 894     v2[0 + pos] = -new_v0;\r
 895     // insert V[33-48] (== new_v[16-31]) into other v:\r
 896     v2[16 + pos] = new_v16;\r
 897     v2[32 + pos] = new_v17;\r
 898     v2[48 + pos] = new_v18;\r
 899     v2[64 + pos] = new_v19;\r
 900     v2[80 + pos] = new_v20;\r
 901     v2[96 + pos] = new_v21;\r
 902     v2[112 + pos] = new_v22;\r
 903     v2[128 + pos] = new_v23;\r
 904     v2[144 + pos] = new_v24;\r
 905     v2[160 + pos] = new_v25;\r
 906     v2[176 + pos] = new_v26;\r
 907     v2[192 + pos] = new_v27;\r
 908     v2[208 + pos] = new_v28;\r
 909     v2[224 + pos] = new_v29;\r
 910     v2[240 + pos] = new_v30;\r
 911     v2[256 + pos] = new_v31;\r
 912 \r
 913     // insert V[49-63] (== new_v[30-16]) into other v:\r
 914     v2[272 + pos] = new_v30;\r
 915     v2[288 + pos] = new_v29;\r
 916     v2[304 + pos] = new_v28;\r
 917     v2[320 + pos] = new_v27;\r
 918     v2[336 + pos] = new_v26;\r
 919     v2[352 + pos] = new_v25;\r
 920     v2[368 + pos] = new_v24;\r
 921     v2[384 + pos] = new_v23;\r
 922     v2[400 + pos] = new_v22;\r
 923     v2[416 + pos] = new_v21;\r
 924     v2[432 + pos] = new_v20;\r
 925     v2[448 + pos] = new_v19;\r
 926     v2[464 + pos] = new_v18;\r
 927     v2[480 + pos] = new_v17;\r
 928     v2[496 + pos] = new_v16;\r
 929 \r
 930     /*\r
 931      * setup PREV\r
 932      */\r
 933 \r
 934     prev1[0 + pos] = new_v0;\r
 935     prev1[16 + pos] = new_v1;\r
 936     prev1[32 + pos] = new_v2;\r
 937     prev1[48 + pos] = new_v3;\r
 938     prev1[64 + pos] = new_v4;\r
 939     prev1[80 + pos] = new_v5;\r
 940     prev1[96 + pos] = new_v6;\r
 941     prev1[112 + pos] = new_v7;\r
 942     prev1[128 + pos] = new_v8;\r
 943     prev1[144 + pos] = new_v9;\r
 944     prev1[160 + pos] = new_v10;\r
 945     prev1[176 + pos] = new_v11;\r
 946     prev1[192 + pos] = new_v12;\r
 947     prev1[208 + pos] = new_v13;\r
 948     prev1[224 + pos] = new_v14;\r
 949     prev1[240 + pos] = new_v15;\r
 950 \r
 951     // V[16] is always 0.0:\r
 952     prev1[256 + pos] = 0.0f;\r
 953 \r
 954     // insert V[17-31] (== -new_v[15-1]) into actual v:\r
 955     prev1[272 + pos] = -new_v15;\r
 956     prev1[288 + pos] = -new_v14;\r
 957     prev1[304 + pos] = -new_v13;\r
 958     prev1[320 + pos] = -new_v12;\r
 959     prev1[336 + pos] = -new_v11;\r
 960     prev1[352 + pos] = -new_v10;\r
 961     prev1[368 + pos] = -new_v9;\r
 962     prev1[384 + pos] = -new_v8;\r
 963     prev1[400 + pos] = -new_v7;\r
 964     prev1[416 + pos] = -new_v6;\r
 965     prev1[432 + pos] = -new_v5;\r
 966     prev1[448 + pos] = -new_v4;\r
 967     prev1[464 + pos] = -new_v3;\r
 968     prev1[480 + pos] = -new_v2;\r
 969     prev1[496 + pos] = -new_v1;\r
 970 \r
 971     // insert V[32] (== -new_v[0]) into other v:\r
 972     // dest = (actual_v == v1) ? v2 : v1;\r
 973 \r
 974     prev2[0 + pos] = -new_v0;\r
 975     // insert V[33-48] (== new_v[16-31]) into other v:\r
 976     prev2[16 + pos] = new_v16;\r
 977     prev2[32 + pos] = new_v17;\r
 978     prev2[48 + pos] = new_v18;\r
 979     prev2[64 + pos] = new_v19;\r
 980     prev2[80 + pos] = new_v20;\r
 981     prev2[96 + pos] = new_v21;\r
 982     prev2[112 + pos] = new_v22;\r
 983     prev2[128 + pos] = new_v23;\r
 984     prev2[144 + pos] = new_v24;\r
 985     prev2[160 + pos] = new_v25;\r
 986     prev2[176 + pos] = new_v26;\r
 987     prev2[192 + pos] = new_v27;\r
 988     prev2[208 + pos] = new_v28;\r
 989     prev2[224 + pos] = new_v29;\r
 990     prev2[240 + pos] = new_v30;\r
 991     prev2[256 + pos] = new_v31;\r
 992 \r
 993     // insert V[49-63] (== new_v[30-16]) into other v:\r
 994     prev2[272 + pos] = new_v30;\r
 995     prev2[288 + pos] = new_v29;\r
 996     prev2[304 + pos] = new_v28;\r
 997     prev2[320 + pos] = new_v27;\r
 998     prev2[336 + pos] = new_v26;\r
 999     prev2[352 + pos] = new_v25;\r
1000     prev2[368 + pos] = new_v24;\r
1001     prev2[384 + pos] = new_v23;\r
1002     prev2[400 + pos] = new_v22;\r
1003     prev2[416 + pos] = new_v21;\r
1004     prev2[432 + pos] = new_v20;\r
1005     prev2[448 + pos] = new_v19;\r
1006     prev2[464 + pos] = new_v18;\r
1007     prev2[480 + pos] = new_v17;\r
1008     prev2[496 + pos] = new_v16;\r
1009   }\r
1010 \r
1011   /**\r
1012    * Compute PCM Samples.\r
1013    */\r
1014 \r
1015   private float[] _tmpOut = new float[32];\r
1016 \r
1017   private void compute_pcm_samples0() {\r
1018 \r
1019     if (vidx == 1) {\r
1020       // final float[] vp = actual_v;\r
1021       // int inc = v_inc;\r
1022       // final float[] tmpOut = _tmpOut;\r
1023       int dvp = 0;\r
1024 \r
1025       // fat chance of having this loop unroll\r
1026       for (int i = 0; i < 32; i++) {\r
1027         float pcm_sample;\r
1028         // final float[] dp = d16[i];\r
1029         pcm_sample = (float) (((v1[0 + dvp] * d16[i][0]) + (v1[15 + dvp] * d16[i][1]) + (v1[14 + dvp] * d16[i][2]) + (v1[13 + dvp] * d16[i][3]) + (v1[12 + dvp] * d16[i][4]) + (v1[11 + dvp] * d16[i][5]) + (v1[10 + dvp] * d16[i][6]) + (v1[9 + dvp] * d16[i][7]) + (v1[8 + dvp] * d16[i][8]) + (v1[7 + dvp] * d16[i][9]) + (v1[6 + dvp] * d16[i][10]) + (v1[5 + dvp] * d16[i][11]) + (v1[4 + dvp] * d16[i][12]) + (v1[3 + dvp] * d16[i][13]) + (v1[2 + dvp] * d16[i][14]) + (v1[1 + dvp] * d16[i][15])) * scalefactor);\r
1030 \r
1031         _tmpOut[i] = pcm_sample;\r
1032 \r
1033         dvp += 16;\r
1034       } // for\r
1035     } else {\r
1036       // final float[] vp = actual_v;\r
1037       // int inc = v_inc;\r
1038       // final float[] tmpOut = _tmpOut;\r
1039       int dvp = 0;\r
1040 \r
1041       // fat chance of having this loop unroll\r
1042       for (int i = 0; i < 32; i++) {\r
1043         float pcm_sample;\r
1044         // final float[] dp = d16[i];\r
1045         pcm_sample = (float) (((v2[0 + dvp] * d16[i][0]) + (v2[15 + dvp] * d16[i][1]) + (v2[14 + dvp] * d16[i][2]) + (v2[13 + dvp] * d16[i][3]) + (v2[12 + dvp] * d16[i][4]) + (v2[11 + dvp] * d16[i][5]) + (v2[10 + dvp] * d16[i][6]) + (v2[9 + dvp] * d16[i][7]) + (v2[8 + dvp] * d16[i][8]) + (v2[7 + dvp] * d16[i][9]) + (v2[6 + dvp] * d16[i][10]) + (v2[5 + dvp] * d16[i][11]) + (v2[4 + dvp] * d16[i][12]) + (v2[3 + dvp] * d16[i][13]) + (v2[2 + dvp] * d16[i][14]) + (v2[1 + dvp] * d16[i][15])) * scalefactor);\r
1046 \r
1047         _tmpOut[i] = pcm_sample;\r
1048 \r
1049         dvp += 16;\r
1050       } // for\r
1051     }\r
1052 \r
1053   }\r
1054 \r
1055   private void compute_pcm_samples1() {\r
1056 \r
1057     if (vidx == 1) {\r
1058       // final float[] vp = actual_v;\r
1059       // int inc = v_inc;\r
1060       // final float[] tmpOut = _tmpOut;\r
1061       int dvp = 0;\r
1062 \r
1063       // fat chance of having this loop unroll\r
1064       for (int i = 0; i < 32; i++) {\r
1065         // final float[] dp = d16[i];\r
1066         float pcm_sample;\r
1067 \r
1068         pcm_sample = (float) (((v1[1 + dvp] * d16[i][0]) + (v1[0 + dvp] * d16[i][1]) + (v1[15 + dvp] * d16[i][2]) + (v1[14 + dvp] * d16[i][3]) + (v1[13 + dvp] * d16[i][4]) + (v1[12 + dvp] * d16[i][5]) + (v1[11 + dvp] * d16[i][6]) + (v1[10 + dvp] * d16[i][7]) + (v1[9 + dvp] * d16[i][8]) + (v1[8 + dvp] * d16[i][9]) + (v1[7 + dvp] * d16[i][10]) + (v1[6 + dvp] * d16[i][11]) + (v1[5 + dvp] * d16[i][12]) + (v1[4 + dvp] * d16[i][13]) + (v1[3 + dvp] * d16[i][14]) + (v1[2 + dvp] * d16[i][15])) * scalefactor);\r
1069 \r
1070         _tmpOut[i] = pcm_sample;\r
1071 \r
1072         dvp += 16;\r
1073       } // for\r
1074     } else {\r
1075       // final float[] vp = actual_v;\r
1076       // int inc = v_inc;\r
1077       // final float[] tmpOut = _tmpOut;\r
1078       int dvp = 0;\r
1079 \r
1080       // fat chance of having this loop unroll\r
1081       for (int i = 0; i < 32; i++) {\r
1082         // final float[] dp = d16[i];\r
1083         float pcm_sample;\r
1084 \r
1085         pcm_sample = (float) (((v2[1 + dvp] * d16[i][0]) + (v2[0 + dvp] * d16[i][1]) + (v2[15 + dvp] * d16[i][2]) + (v2[14 + dvp] * d16[i][3]) + (v2[13 + dvp] * d16[i][4]) + (v2[12 + dvp] * d16[i][5]) + (v2[11 + dvp] * d16[i][6]) + (v2[10 + dvp] * d16[i][7]) + (v2[9 + dvp] * d16[i][8]) + (v2[8 + dvp] * d16[i][9]) + (v2[7 + dvp] * d16[i][10]) + (v2[6 + dvp] * d16[i][11]) + (v2[5 + dvp] * d16[i][12]) + (v2[4 + dvp] * d16[i][13]) + (v2[3 + dvp] * d16[i][14]) + (v2[2 + dvp] * d16[i][15])) * scalefactor);\r
1086 \r
1087         _tmpOut[i] = pcm_sample;\r
1088 \r
1089         dvp += 16;\r
1090       } // for\r
1091     }\r
1092 \r
1093   }\r
1094 \r
1095   private void compute_pcm_samples2() {\r
1096 \r
1097     if (vidx == 1) {\r
1098       // final float[] vp = actual_v;\r
1099       // int inc = v_inc;\r
1100       // final float[] tmpOut = _tmpOut;\r
1101       int dvp = 0;\r
1102 \r
1103       // fat chance of having this loop unroll\r
1104       for (int i = 0; i < 32; i++) {\r
1105         // final float[] dp = d16[i];\r
1106         float pcm_sample;\r
1107 \r
1108         pcm_sample = (float) (((v1[2 + dvp] * d16[i][0]) + (v1[1 + dvp] * d16[i][1]) + (v1[0 + dvp] * d16[i][2]) + (v1[15 + dvp] * d16[i][3]) + (v1[14 + dvp] * d16[i][4]) + (v1[13 + dvp] * d16[i][5]) + (v1[12 + dvp] * d16[i][6]) + (v1[11 + dvp] * d16[i][7]) + (v1[10 + dvp] * d16[i][8]) + (v1[9 + dvp] * d16[i][9]) + (v1[8 + dvp] * d16[i][10]) + (v1[7 + dvp] * d16[i][11]) + (v1[6 + dvp] * d16[i][12]) + (v1[5 + dvp] * d16[i][13]) + (v1[4 + dvp] * d16[i][14]) + (v1[3 + dvp] * d16[i][15])) * scalefactor);\r
1109 \r
1110         _tmpOut[i] = pcm_sample;\r
1111 \r
1112         dvp += 16;\r
1113       } // for\r
1114     } else {\r
1115       // final float[] vp = actual_v;\r
1116       // int inc = v_inc;\r
1117       // final float[] tmpOut = _tmpOut;\r
1118       int dvp = 0;\r
1119 \r
1120       // fat chance of having this loop unroll\r
1121       for (int i = 0; i < 32; i++) {\r
1122         // final float[] dp = d16[i];\r
1123         float pcm_sample;\r
1124 \r
1125         pcm_sample = (float) (((v2[2 + dvp] * d16[i][0]) + (v2[1 + dvp] * d16[i][1]) + (v2[0 + dvp] * d16[i][2]) + (v2[15 + dvp] * d16[i][3]) + (v2[14 + dvp] * d16[i][4]) + (v2[13 + dvp] * d16[i][5]) + (v2[12 + dvp] * d16[i][6]) + (v2[11 + dvp] * d16[i][7]) + (v2[10 + dvp] * d16[i][8]) + (v2[9 + dvp] * d16[i][9]) + (v2[8 + dvp] * d16[i][10]) + (v2[7 + dvp] * d16[i][11]) + (v2[6 + dvp] * d16[i][12]) + (v2[5 + dvp] * d16[i][13]) + (v2[4 + dvp] * d16[i][14]) + (v2[3 + dvp] * d16[i][15])) * scalefactor);\r
1126 \r
1127         _tmpOut[i] = pcm_sample;\r
1128 \r
1129         dvp += 16;\r
1130       } // for\r
1131     }\r
1132 \r
1133   }\r
1134 \r
1135   private void compute_pcm_samples3() {\r
1136 \r
1137     if (vidx == 1) {\r
1138       // final float[] vp = actual_v;\r
1139 \r
1140       // int inc = v_inc;\r
1141       // final float[] tmpOut = _tmpOut;\r
1142       int dvp = 0;\r
1143 \r
1144       // fat chance of having this loop unroll\r
1145       for (int i = 0; i < 32; i++) {\r
1146         // final float[] dp = d16[i];\r
1147         float pcm_sample;\r
1148 \r
1149         pcm_sample = (float) (((v1[3 + dvp] * d16[i][0]) + (v1[2 + dvp] * d16[i][1]) + (v1[1 + dvp] * d16[i][2]) + (v1[0 + dvp] * d16[i][3]) + (v1[15 + dvp] * d16[i][4]) + (v1[14 + dvp] * d16[i][5]) + (v1[13 + dvp] * d16[i][6]) + (v1[12 + dvp] * d16[i][7]) + (v1[11 + dvp] * d16[i][8]) + (v1[10 + dvp] * d16[i][9]) + (v1[9 + dvp] * d16[i][10]) + (v1[8 + dvp] * d16[i][11]) + (v1[7 + dvp] * d16[i][12]) + (v1[6 + dvp] * d16[i][13]) + (v1[5 + dvp] * d16[i][14]) + (v1[4 + dvp] * d16[i][15])) * scalefactor);\r
1150 \r
1151         _tmpOut[i] = pcm_sample;\r
1152 \r
1153         dvp += 16;\r
1154       } // for\r
1155     } else {\r
1156       // final float[] vp = actual_v;\r
1157 \r
1158       // int inc = v_inc;\r
1159       // final float[] tmpOut = _tmpOut;\r
1160       int dvp = 0;\r
1161 \r
1162       // fat chance of having this loop unroll\r
1163       for (int i = 0; i < 32; i++) {\r
1164         // final float[] dp = d16[i];\r
1165         float pcm_sample;\r
1166 \r
1167         pcm_sample = (float) (((v2[3 + dvp] * d16[i][0]) + (v2[2 + dvp] * d16[i][1]) + (v2[1 + dvp] * d16[i][2]) + (v2[0 + dvp] * d16[i][3]) + (v2[15 + dvp] * d16[i][4]) + (v2[14 + dvp] * d16[i][5]) + (v2[13 + dvp] * d16[i][6]) + (v2[12 + dvp] * d16[i][7]) + (v2[11 + dvp] * d16[i][8]) + (v2[10 + dvp] * d16[i][9]) + (v2[9 + dvp] * d16[i][10]) + (v2[8 + dvp] * d16[i][11]) + (v2[7 + dvp] * d16[i][12]) + (v2[6 + dvp] * d16[i][13]) + (v2[5 + dvp] * d16[i][14]) + (v2[4 + dvp] * d16[i][15])) * scalefactor);\r
1168 \r
1169         _tmpOut[i] = pcm_sample;\r
1170 \r
1171         dvp += 16;\r
1172       } // for\r
1173     }\r
1174 \r
1175   }\r
1176 \r
1177   private void compute_pcm_samples4() {\r
1178 \r
1179     if (vidx == 1) {\r
1180       // final float[] vp = actual_v;\r
1181 \r
1182       // int inc = v_inc;\r
1183       // final float[] tmpOut = _tmpOut;\r
1184       int dvp = 0;\r
1185 \r
1186       // fat chance of having this loop unroll\r
1187       for (int i = 0; i < 32; i++) {\r
1188         // final float[] dp = d16[i];\r
1189         float pcm_sample;\r
1190 \r
1191         pcm_sample = (float) (((v1[4 + dvp] * d16[i][0]) + (v1[3 + dvp] * d16[i][1]) + (v1[2 + dvp] * d16[i][2]) + (v1[1 + dvp] * d16[i][3]) + (v1[0 + dvp] * d16[i][4]) + (v1[15 + dvp] * d16[i][5]) + (v1[14 + dvp] * d16[i][6]) + (v1[13 + dvp] * d16[i][7]) + (v1[12 + dvp] * d16[i][8]) + (v1[11 + dvp] * d16[i][9]) + (v1[10 + dvp] * d16[i][10]) + (v1[9 + dvp] * d16[i][11]) + (v1[8 + dvp] * d16[i][12]) + (v1[7 + dvp] * d16[i][13]) + (v1[6 + dvp] * d16[i][14]) + (v1[5 + dvp] * d16[i][15])) * scalefactor);\r
1192 \r
1193         _tmpOut[i] = pcm_sample;\r
1194 \r
1195         dvp += 16;\r
1196       } // for\r
1197     } else {\r
1198       // final float[] vp = actual_v;\r
1199 \r
1200       // int inc = v_inc;\r
1201       // final float[] tmpOut = _tmpOut;\r
1202       int dvp = 0;\r
1203 \r
1204       // fat chance of having this loop unroll\r
1205       for (int i = 0; i < 32; i++) {\r
1206         // final float[] dp = d16[i];\r
1207         float pcm_sample;\r
1208 \r
1209         pcm_sample = (float) (((v2[4 + dvp] * d16[i][0]) + (v2[3 + dvp] * d16[i][1]) + (v2[2 + dvp] * d16[i][2]) + (v2[1 + dvp] * d16[i][3]) + (v2[0 + dvp] * d16[i][4]) + (v2[15 + dvp] * d16[i][5]) + (v2[14 + dvp] * d16[i][6]) + (v2[13 + dvp] * d16[i][7]) + (v2[12 + dvp] * d16[i][8]) + (v2[11 + dvp] * d16[i][9]) + (v2[10 + dvp] * d16[i][10]) + (v2[9 + dvp] * d16[i][11]) + (v2[8 + dvp] * d16[i][12]) + (v2[7 + dvp] * d16[i][13]) + (v2[6 + dvp] * d16[i][14]) + (v2[5 + dvp] * d16[i][15])) * scalefactor);\r
1210 \r
1211         _tmpOut[i] = pcm_sample;\r
1212 \r
1213         dvp += 16;\r
1214       } // for\r
1215     }\r
1216 \r
1217   }\r
1218 \r
1219   private void compute_pcm_samples5() {\r
1220 \r
1221     if (vidx == 1) {\r
1222       // final float[] vp = actual_v;\r
1223 \r
1224       // int inc = v_inc;\r
1225       // final float[] tmpOut = _tmpOut;\r
1226       int dvp = 0;\r
1227 \r
1228       // fat chance of having this loop unroll\r
1229       for (int i = 0; i < 32; i++) {\r
1230         // final float[] dp = d16[i];\r
1231         float pcm_sample;\r
1232 \r
1233         pcm_sample = (float) (((v1[5 + dvp] * d16[i][0]) + (v1[4 + dvp] * d16[i][1]) + (v1[3 + dvp] * d16[i][2]) + (v1[2 + dvp] * d16[i][3]) + (v1[1 + dvp] * d16[i][4]) + (v1[0 + dvp] * d16[i][5]) + (v1[15 + dvp] * d16[i][6]) + (v1[14 + dvp] * d16[i][7]) + (v1[13 + dvp] * d16[i][8]) + (v1[12 + dvp] * d16[i][9]) + (v1[11 + dvp] * d16[i][10]) + (v1[10 + dvp] * d16[i][11]) + (v1[9 + dvp] * d16[i][12]) + (v1[8 + dvp] * d16[i][13]) + (v1[7 + dvp] * d16[i][14]) + (v1[6 + dvp] * d16[i][15])) * scalefactor);\r
1234 \r
1235         _tmpOut[i] = pcm_sample;\r
1236 \r
1237         dvp += 16;\r
1238       } // for\r
1239     } else {\r
1240       // final float[] vp = actual_v;\r
1241 \r
1242       // int inc = v_inc;\r
1243       // final float[] tmpOut = _tmpOut;\r
1244       int dvp = 0;\r
1245 \r
1246       // fat chance of having this loop unroll\r
1247       for (int i = 0; i < 32; i++) {\r
1248         // final float[] dp = d16[i];\r
1249         float pcm_sample;\r
1250 \r
1251         pcm_sample = (float) (((v2[5 + dvp] * d16[i][0]) + (v2[4 + dvp] * d16[i][1]) + (v2[3 + dvp] * d16[i][2]) + (v2[2 + dvp] * d16[i][3]) + (v2[1 + dvp] * d16[i][4]) + (v2[0 + dvp] * d16[i][5]) + (v2[15 + dvp] * d16[i][6]) + (v2[14 + dvp] * d16[i][7]) + (v2[13 + dvp] * d16[i][8]) + (v2[12 + dvp] * d16[i][9]) + (v2[11 + dvp] * d16[i][10]) + (v2[10 + dvp] * d16[i][11]) + (v2[9 + dvp] * d16[i][12]) + (v2[8 + dvp] * d16[i][13]) + (v2[7 + dvp] * d16[i][14]) + (v2[6 + dvp] * d16[i][15])) * scalefactor);\r
1252 \r
1253         _tmpOut[i] = pcm_sample;\r
1254 \r
1255         dvp += 16;\r
1256       } // for\r
1257     }\r
1258 \r
1259   }\r
1260 \r
1261   private void compute_pcm_samples6() {\r
1262 \r
1263     if (vidx == 1) {\r
1264 \r
1265       // final float[] vp = actual_v;\r
1266       // int inc = v_inc;\r
1267       // final float[] tmpOut = _tmpOut;\r
1268       int dvp = 0;\r
1269 \r
1270       // fat chance of having this loop unroll\r
1271       for (int i = 0; i < 32; i++) {\r
1272         // final float[] dp = d16[i];\r
1273         float pcm_sample;\r
1274 \r
1275         pcm_sample = (float) (((v1[6 + dvp] * d16[i][0]) + (v1[5 + dvp] * d16[i][1]) + (v1[4 + dvp] * d16[i][2]) + (v1[3 + dvp] * d16[i][3]) + (v1[2 + dvp] * d16[i][4]) + (v1[1 + dvp] * d16[i][5]) + (v1[0 + dvp] * d16[i][6]) + (v1[15 + dvp] * d16[i][7]) + (v1[14 + dvp] * d16[i][8]) + (v1[13 + dvp] * d16[i][9]) + (v1[12 + dvp] * d16[i][10]) + (v1[11 + dvp] * d16[i][11]) + (v1[10 + dvp] * d16[i][12]) + (v1[9 + dvp] * d16[i][13]) + (v1[8 + dvp] * d16[i][14]) + (v1[7 + dvp] * d16[i][15])) * scalefactor);\r
1276 \r
1277         _tmpOut[i] = pcm_sample;\r
1278 \r
1279         dvp += 16;\r
1280       } // for\r
1281     } else {\r
1282 \r
1283       // final float[] vp = actual_v;\r
1284       // int inc = v_inc;\r
1285       // final float[] tmpOut = _tmpOut;\r
1286       int dvp = 0;\r
1287 \r
1288       // fat chance of having this loop unroll\r
1289       for (int i = 0; i < 32; i++) {\r
1290         // final float[] dp = d16[i];\r
1291         float pcm_sample;\r
1292 \r
1293         pcm_sample = (float) (((v2[6 + dvp] * d16[i][0]) + (v2[5 + dvp] * d16[i][1]) + (v2[4 + dvp] * d16[i][2]) + (v2[3 + dvp] * d16[i][3]) + (v2[2 + dvp] * d16[i][4]) + (v2[1 + dvp] * d16[i][5]) + (v2[0 + dvp] * d16[i][6]) + (v2[15 + dvp] * d16[i][7]) + (v2[14 + dvp] * d16[i][8]) + (v2[13 + dvp] * d16[i][9]) + (v2[12 + dvp] * d16[i][10]) + (v2[11 + dvp] * d16[i][11]) + (v2[10 + dvp] * d16[i][12]) + (v2[9 + dvp] * d16[i][13]) + (v2[8 + dvp] * d16[i][14]) + (v2[7 + dvp] * d16[i][15])) * scalefactor);\r
1294 \r
1295         _tmpOut[i] = pcm_sample;\r
1296 \r
1297         dvp += 16;\r
1298       } // for\r
1299     }\r
1300 \r
1301   }\r
1302 \r
1303   private void compute_pcm_samples7() {\r
1304 \r
1305     if (vidx == 1) {\r
1306       // final float[] vp = actual_v;\r
1307 \r
1308       // int inc = v_inc;\r
1309       // final float[] tmpOut = _tmpOut;\r
1310       int dvp = 0;\r
1311 \r
1312       // fat chance of having this loop unroll\r
1313       for (int i = 0; i < 32; i++) {\r
1314         // final float[] dp = d16[i];\r
1315         float pcm_sample;\r
1316 \r
1317         pcm_sample = (float) (((v1[7 + dvp] * d16[i][0]) + (v1[6 + dvp] * d16[i][1]) + (v1[5 + dvp] * d16[i][2]) + (v1[4 + dvp] * d16[i][3]) + (v1[3 + dvp] * d16[i][4]) + (v1[2 + dvp] * d16[i][5]) + (v1[1 + dvp] * d16[i][6]) + (v1[0 + dvp] * d16[i][7]) + (v1[15 + dvp] * d16[i][8]) + (v1[14 + dvp] * d16[i][9]) + (v1[13 + dvp] * d16[i][10]) + (v1[12 + dvp] * d16[i][11]) + (v1[11 + dvp] * d16[i][12]) + (v1[10 + dvp] * d16[i][13]) + (v1[9 + dvp] * d16[i][14]) + (v1[8 + dvp] * d16[i][15])) * scalefactor);\r
1318 \r
1319         _tmpOut[i] = pcm_sample;\r
1320 \r
1321         dvp += 16;\r
1322       } // for\r
1323     } else {\r
1324       // final float[] vp = actual_v;\r
1325 \r
1326       // int inc = v_inc;\r
1327       // final float[] tmpOut = _tmpOut;\r
1328       int dvp = 0;\r
1329 \r
1330       // fat chance of having this loop unroll\r
1331       for (int i = 0; i < 32; i++) {\r
1332         // final float[] dp = d16[i];\r
1333         float pcm_sample;\r
1334 \r
1335         pcm_sample = (float) (((v2[7 + dvp] * d16[i][0]) + (v2[6 + dvp] * d16[i][1]) + (v2[5 + dvp] * d16[i][2]) + (v2[4 + dvp] * d16[i][3]) + (v2[3 + dvp] * d16[i][4]) + (v2[2 + dvp] * d16[i][5]) + (v2[1 + dvp] * d16[i][6]) + (v2[0 + dvp] * d16[i][7]) + (v2[15 + dvp] * d16[i][8]) + (v2[14 + dvp] * d16[i][9]) + (v2[13 + dvp] * d16[i][10]) + (v2[12 + dvp] * d16[i][11]) + (v2[11 + dvp] * d16[i][12]) + (v2[10 + dvp] * d16[i][13]) + (v2[9 + dvp] * d16[i][14]) + (v2[8 + dvp] * d16[i][15])) * scalefactor);\r
1336 \r
1337         _tmpOut[i] = pcm_sample;\r
1338 \r
1339         dvp += 16;\r
1340       } // for\r
1341     }\r
1342 \r
1343   }\r
1344 \r
1345   private void compute_pcm_samples8() {\r
1346 \r
1347     if (vidx == 1) {\r
1348 \r
1349       // final float[] vp = actual_v;\r
1350 \r
1351       // int inc = v_inc;\r
1352       // final float[] tmpOut = _tmpOut;\r
1353       int dvp = 0;\r
1354 \r
1355       // fat chance of having this loop unroll\r
1356       for (int i = 0; i < 32; i++) {\r
1357         // final float[] dp = d16[i];\r
1358         float pcm_sample;\r
1359 \r
1360         pcm_sample = (float) (((v1[8 + dvp] * d16[i][0]) + (v1[7 + dvp] * d16[i][1]) + (v1[6 + dvp] * d16[i][2]) + (v1[5 + dvp] * d16[i][3]) + (v1[4 + dvp] * d16[i][4]) + (v1[3 + dvp] * d16[i][5]) + (v1[2 + dvp] * d16[i][6]) + (v1[1 + dvp] * d16[i][7]) + (v1[0 + dvp] * d16[i][8]) + (v1[15 + dvp] * d16[i][9]) + (v1[14 + dvp] * d16[i][10]) + (v1[13 + dvp] * d16[i][11]) + (v1[12 + dvp] * d16[i][12]) + (v1[11 + dvp] * d16[i][13]) + (v1[10 + dvp] * d16[i][14]) + (v1[9 + dvp] * d16[i][15])) * scalefactor);\r
1361 \r
1362         _tmpOut[i] = pcm_sample;\r
1363 \r
1364         dvp += 16;\r
1365       } // for\r
1366     } else {\r
1367 \r
1368       // final float[] vp = actual_v;\r
1369 \r
1370       // int inc = v_inc;\r
1371       // final float[] tmpOut = _tmpOut;\r
1372       int dvp = 0;\r
1373 \r
1374       // fat chance of having this loop unroll\r
1375       for (int i = 0; i < 32; i++) {\r
1376         // final float[] dp = d16[i];\r
1377         float pcm_sample;\r
1378 \r
1379         pcm_sample = (float) (((v2[8 + dvp] * d16[i][0]) + (v2[7 + dvp] * d16[i][1]) + (v2[6 + dvp] * d16[i][2]) + (v2[5 + dvp] * d16[i][3]) + (v2[4 + dvp] * d16[i][4]) + (v2[3 + dvp] * d16[i][5]) + (v2[2 + dvp] * d16[i][6]) + (v2[1 + dvp] * d16[i][7]) + (v2[0 + dvp] * d16[i][8]) + (v2[15 + dvp] * d16[i][9]) + (v2[14 + dvp] * d16[i][10]) + (v2[13 + dvp] * d16[i][11]) + (v2[12 + dvp] * d16[i][12]) + (v2[11 + dvp] * d16[i][13]) + (v2[10 + dvp] * d16[i][14]) + (v2[9 + dvp] * d16[i][15])) * scalefactor);\r
1380 \r
1381         _tmpOut[i] = pcm_sample;\r
1382 \r
1383         dvp += 16;\r
1384       } // for\r
1385     }\r
1386 \r
1387   }\r
1388 \r
1389   private void compute_pcm_samples9() {\r
1390 \r
1391     if (vidx == 1) {\r
1392       // final float[] vp = actual_v;\r
1393 \r
1394       // int inc = v_inc;\r
1395       // final float[] tmpOut = _tmpOut;\r
1396       int dvp = 0;\r
1397 \r
1398       // fat chance of having this loop unroll\r
1399       for (int i = 0; i < 32; i++) {\r
1400         // final float[] dp = d16[i];\r
1401         float pcm_sample;\r
1402 \r
1403         pcm_sample = (float) (((v1[9 + dvp] * d16[i][0]) + (v1[8 + dvp] * d16[i][1]) + (v1[7 + dvp] * d16[i][2]) + (v1[6 + dvp] * d16[i][3]) + (v1[5 + dvp] * d16[i][4]) + (v1[4 + dvp] * d16[i][5]) + (v1[3 + dvp] * d16[i][6]) + (v1[2 + dvp] * d16[i][7]) + (v1[1 + dvp] * d16[i][8]) + (v1[0 + dvp] * d16[i][9]) + (v1[15 + dvp] * d16[i][10]) + (v1[14 + dvp] * d16[i][11]) + (v1[13 + dvp] * d16[i][12]) + (v1[12 + dvp] * d16[i][13]) + (v1[11 + dvp] * d16[i][14]) + (v1[10 + dvp] * d16[i][15])) * scalefactor);\r
1404 \r
1405         _tmpOut[i] = pcm_sample;\r
1406 \r
1407         dvp += 16;\r
1408       } // for\r
1409     } else {\r
1410       // final float[] vp = actual_v;\r
1411 \r
1412       // int inc = v_inc;\r
1413       // final float[] tmpOut = _tmpOut;\r
1414       int dvp = 0;\r
1415 \r
1416       // fat chance of having this loop unroll\r
1417       for (int i = 0; i < 32; i++) {\r
1418         // final float[] dp = d16[i];\r
1419         float pcm_sample;\r
1420 \r
1421         pcm_sample = (float) (((v2[9 + dvp] * d16[i][0]) + (v2[8 + dvp] * d16[i][1]) + (v2[7 + dvp] * d16[i][2]) + (v2[6 + dvp] * d16[i][3]) + (v2[5 + dvp] * d16[i][4]) + (v2[4 + dvp] * d16[i][5]) + (v2[3 + dvp] * d16[i][6]) + (v2[2 + dvp] * d16[i][7]) + (v2[1 + dvp] * d16[i][8]) + (v2[0 + dvp] * d16[i][9]) + (v2[15 + dvp] * d16[i][10]) + (v2[14 + dvp] * d16[i][11]) + (v2[13 + dvp] * d16[i][12]) + (v2[12 + dvp] * d16[i][13]) + (v2[11 + dvp] * d16[i][14]) + (v2[10 + dvp] * d16[i][15])) * scalefactor);\r
1422 \r
1423         _tmpOut[i] = pcm_sample;\r
1424 \r
1425         dvp += 16;\r
1426       } // for\r
1427     }\r
1428 \r
1429   }\r
1430 \r
1431   private void compute_pcm_samples10() {\r
1432     if (vidx == 1) {\r
1433       // final float[] vp = actual_v;\r
1434       // int inc = v_inc;\r
1435       // final float[] tmpOut = _tmpOut;\r
1436       int dvp = 0;\r
1437 \r
1438       // fat chance of having this loop unroll\r
1439       for (int i = 0; i < 32; i++) {\r
1440         // final float[] dp = d16[i];\r
1441         float pcm_sample;\r
1442 \r
1443         pcm_sample = (float) (((v1[10 + dvp] * d16[i][0]) + (v1[9 + dvp] * d16[i][1]) + (v1[8 + dvp] * d16[i][2]) + (v1[7 + dvp] * d16[i][3]) + (v1[6 + dvp] * d16[i][4]) + (v1[5 + dvp] * d16[i][5]) + (v1[4 + dvp] * d16[i][6]) + (v1[3 + dvp] * d16[i][7]) + (v1[2 + dvp] * d16[i][8]) + (v1[1 + dvp] * d16[i][9]) + (v1[0 + dvp] * d16[i][10]) + (v1[15 + dvp] * d16[i][11]) + (v1[14 + dvp] * d16[i][12]) + (v1[13 + dvp] * d16[i][13]) + (v1[12 + dvp] * d16[i][14]) + (v1[11 + dvp] * d16[i][15])) * scalefactor);\r
1444 \r
1445         _tmpOut[i] = pcm_sample;\r
1446 \r
1447         dvp += 16;\r
1448       } // for\r
1449     } else {\r
1450       // final float[] vp = actual_v;\r
1451       // int inc = v_inc;\r
1452       // final float[] tmpOut = _tmpOut;\r
1453       int dvp = 0;\r
1454 \r
1455       // fat chance of having this loop unroll\r
1456       for (int i = 0; i < 32; i++) {\r
1457         // final float[] dp = d16[i];\r
1458         float pcm_sample;\r
1459 \r
1460         pcm_sample = (float) (((v2[10 + dvp] * d16[i][0]) + (v2[9 + dvp] * d16[i][1]) + (v2[8 + dvp] * d16[i][2]) + (v2[7 + dvp] * d16[i][3]) + (v2[6 + dvp] * d16[i][4]) + (v2[5 + dvp] * d16[i][5]) + (v2[4 + dvp] * d16[i][6]) + (v2[3 + dvp] * d16[i][7]) + (v2[2 + dvp] * d16[i][8]) + (v2[1 + dvp] * d16[i][9]) + (v2[0 + dvp] * d16[i][10]) + (v2[15 + dvp] * d16[i][11]) + (v2[14 + dvp] * d16[i][12]) + (v2[13 + dvp] * d16[i][13]) + (v2[12 + dvp] * d16[i][14]) + (v2[11 + dvp] * d16[i][15])) * scalefactor);\r
1461 \r
1462         _tmpOut[i] = pcm_sample;\r
1463 \r
1464         dvp += 16;\r
1465       } // for\r
1466     }\r
1467 \r
1468   }\r
1469 \r
1470   private void compute_pcm_samples11() {\r
1471 \r
1472     if (vidx == 1) {\r
1473       // final float[] vp = actual_v;\r
1474 \r
1475       // int inc = v_inc;\r
1476       // final float[] tmpOut = _tmpOut;\r
1477       int dvp = 0;\r
1478 \r
1479       // fat chance of having this loop unroll\r
1480       for (int i = 0; i < 32; i++) {\r
1481         // final float[] dp = d16[i];\r
1482         float pcm_sample;\r
1483 \r
1484         pcm_sample = (float) (((v1[11 + dvp] * d16[i][0]) + (v1[10 + dvp] * d16[i][1]) + (v1[9 + dvp] * d16[i][2]) + (v1[8 + dvp] * d16[i][3]) + (v1[7 + dvp] * d16[i][4]) + (v1[6 + dvp] * d16[i][5]) + (v1[5 + dvp] * d16[i][6]) + (v1[4 + dvp] * d16[i][7]) + (v1[3 + dvp] * d16[i][8]) + (v1[2 + dvp] * d16[i][9]) + (v1[1 + dvp] * d16[i][10]) + (v1[0 + dvp] * d16[i][11]) + (v1[15 + dvp] * d16[i][12]) + (v1[14 + dvp] * d16[i][13]) + (v1[13 + dvp] * d16[i][14]) + (v1[12 + dvp] * d16[i][15])) * scalefactor);\r
1485 \r
1486         _tmpOut[i] = pcm_sample;\r
1487 \r
1488         dvp += 16;\r
1489       } // for\r
1490     } else {\r
1491       // final float[] vp = actual_v;\r
1492 \r
1493       // int inc = v_inc;\r
1494       // final float[] tmpOut = _tmpOut;\r
1495       int dvp = 0;\r
1496 \r
1497       // fat chance of having this loop unroll\r
1498       for (int i = 0; i < 32; i++) {\r
1499         // final float[] dp = d16[i];\r
1500         float pcm_sample;\r
1501 \r
1502         pcm_sample = (float) (((v2[11 + dvp] * d16[i][0]) + (v2[10 + dvp] * d16[i][1]) + (v2[9 + dvp] * d16[i][2]) + (v2[8 + dvp] * d16[i][3]) + (v2[7 + dvp] * d16[i][4]) + (v2[6 + dvp] * d16[i][5]) + (v2[5 + dvp] * d16[i][6]) + (v2[4 + dvp] * d16[i][7]) + (v2[3 + dvp] * d16[i][8]) + (v2[2 + dvp] * d16[i][9]) + (v2[1 + dvp] * d16[i][10]) + (v2[0 + dvp] * d16[i][11]) + (v2[15 + dvp] * d16[i][12]) + (v2[14 + dvp] * d16[i][13]) + (v2[13 + dvp] * d16[i][14]) + (v2[12 + dvp] * d16[i][15])) * scalefactor);\r
1503 \r
1504         _tmpOut[i] = pcm_sample;\r
1505 \r
1506         dvp += 16;\r
1507       } // for\r
1508     }\r
1509 \r
1510   }\r
1511 \r
1512   private void compute_pcm_samples12() {\r
1513 \r
1514     if (vidx == 1) {\r
1515       // final float[] vp = actual_v;\r
1516       // int inc = v_inc;\r
1517       // final float[] tmpOut = _tmpOut;\r
1518       int dvp = 0;\r
1519 \r
1520       // fat chance of having this loop unroll\r
1521       for (int i = 0; i < 32; i++) {\r
1522         // final float[] dp = d16[i];\r
1523         float pcm_sample;\r
1524 \r
1525         pcm_sample = (float) (((v1[12 + dvp] * d16[i][0]) + (v1[11 + dvp] * d16[i][1]) + (v1[10 + dvp] * d16[i][2]) + (v1[9 + dvp] * d16[i][3]) + (v1[8 + dvp] * d16[i][4]) + (v1[7 + dvp] * d16[i][5]) + (v1[6 + dvp] * d16[i][6]) + (v1[5 + dvp] * d16[i][7]) + (v1[4 + dvp] * d16[i][8]) + (v1[3 + dvp] * d16[i][9]) + (v1[2 + dvp] * d16[i][10]) + (v1[1 + dvp] * d16[i][11]) + (v1[0 + dvp] * d16[i][12]) + (v1[15 + dvp] * d16[i][13]) + (v1[14 + dvp] * d16[i][14]) + (v1[13 + dvp] * d16[i][15])) * scalefactor);\r
1526 \r
1527         _tmpOut[i] = pcm_sample;\r
1528 \r
1529         dvp += 16;\r
1530       } // for\r
1531     } else {\r
1532       // final float[] vp = actual_v;\r
1533       // int inc = v_inc;\r
1534       // final float[] tmpOut = _tmpOut;\r
1535       int dvp = 0;\r
1536 \r
1537       // fat chance of having this loop unroll\r
1538       for (int i = 0; i < 32; i++) {\r
1539         // final float[] dp = d16[i];\r
1540         float pcm_sample;\r
1541 \r
1542         pcm_sample = (float) (((v2[12 + dvp] * d16[i][0]) + (v2[11 + dvp] * d16[i][1]) + (v2[10 + dvp] * d16[i][2]) + (v2[9 + dvp] * d16[i][3]) + (v2[8 + dvp] * d16[i][4]) + (v2[7 + dvp] * d16[i][5]) + (v2[6 + dvp] * d16[i][6]) + (v2[5 + dvp] * d16[i][7]) + (v2[4 + dvp] * d16[i][8]) + (v2[3 + dvp] * d16[i][9]) + (v2[2 + dvp] * d16[i][10]) + (v2[1 + dvp] * d16[i][11]) + (v2[0 + dvp] * d16[i][12]) + (v2[15 + dvp] * d16[i][13]) + (v2[14 + dvp] * d16[i][14]) + (v2[13 + dvp] * d16[i][15])) * scalefactor);\r
1543 \r
1544         _tmpOut[i] = pcm_sample;\r
1545 \r
1546         dvp += 16;\r
1547       } // for\r
1548     }\r
1549 \r
1550   }\r
1551 \r
1552   private void compute_pcm_samples13() {\r
1553 \r
1554     if (vidx == 1) {\r
1555       // final float[] vp = actual_v;\r
1556 \r
1557       // int inc = v_inc;\r
1558       // final float[] tmpOut = _tmpOut;\r
1559       int dvp = 0;\r
1560 \r
1561       // fat chance of having this loop unroll\r
1562       for (int i = 0; i < 32; i++) {\r
1563         // final float[] dp = d16[i];\r
1564         float pcm_sample;\r
1565 \r
1566         pcm_sample = (float) (((v1[13 + dvp] * d16[i][0]) + (v1[12 + dvp] * d16[i][1]) + (v1[11 + dvp] * d16[i][2]) + (v1[10 + dvp] * d16[i][3]) + (v1[9 + dvp] * d16[i][4]) + (v1[8 + dvp] * d16[i][5]) + (v1[7 + dvp] * d16[i][6]) + (v1[6 + dvp] * d16[i][7]) + (v1[5 + dvp] * d16[i][8]) + (v1[4 + dvp] * d16[i][9]) + (v1[3 + dvp] * d16[i][10]) + (v1[2 + dvp] * d16[i][11]) + (v1[1 + dvp] * d16[i][12]) + (v1[0 + dvp] * d16[i][13]) + (v1[15 + dvp] * d16[i][14]) + (v1[14 + dvp] * d16[i][15])) * scalefactor);\r
1567 \r
1568         _tmpOut[i] = pcm_sample;\r
1569 \r
1570         dvp += 16;\r
1571       } // for\r
1572     } else {\r
1573       // final float[] vp = actual_v;\r
1574 \r
1575       // int inc = v_inc;\r
1576       // final float[] tmpOut = _tmpOut;\r
1577       int dvp = 0;\r
1578 \r
1579       // fat chance of having this loop unroll\r
1580       for (int i = 0; i < 32; i++) {\r
1581         // final float[] dp = d16[i];\r
1582         float pcm_sample;\r
1583 \r
1584         pcm_sample = (float) (((v2[13 + dvp] * d16[i][0]) + (v2[12 + dvp] * d16[i][1]) + (v2[11 + dvp] * d16[i][2]) + (v2[10 + dvp] * d16[i][3]) + (v2[9 + dvp] * d16[i][4]) + (v2[8 + dvp] * d16[i][5]) + (v2[7 + dvp] * d16[i][6]) + (v2[6 + dvp] * d16[i][7]) + (v2[5 + dvp] * d16[i][8]) + (v2[4 + dvp] * d16[i][9]) + (v2[3 + dvp] * d16[i][10]) + (v2[2 + dvp] * d16[i][11]) + (v2[1 + dvp] * d16[i][12]) + (v2[0 + dvp] * d16[i][13]) + (v2[15 + dvp] * d16[i][14]) + (v2[14 + dvp] * d16[i][15])) * scalefactor);\r
1585 \r
1586         _tmpOut[i] = pcm_sample;\r
1587 \r
1588         dvp += 16;\r
1589       } // for\r
1590     }\r
1591 \r
1592   }\r
1593 \r
1594   private void compute_pcm_samples14() {\r
1595 \r
1596     if (vidx == 1) {\r
1597       // final float[] vp = actual_v;\r
1598 \r
1599       // int inc = v_inc;\r
1600       // final float[] tmpOut = _tmpOut;\r
1601       int dvp = 0;\r
1602 \r
1603       // fat chance of having this loop unroll\r
1604       for (int i = 0; i < 32; i++) {\r
1605         // final float[] dp = d16[i];\r
1606         float pcm_sample;\r
1607 \r
1608         pcm_sample = (float) (((v1[14 + dvp] * d16[i][0]) + (v1[13 + dvp] * d16[i][1]) + (v1[12 + dvp] * d16[i][2]) + (v1[11 + dvp] * d16[i][3]) + (v1[10 + dvp] * d16[i][4]) + (v1[9 + dvp] * d16[i][5]) + (v1[8 + dvp] * d16[i][6]) + (v1[7 + dvp] * d16[i][7]) + (v1[6 + dvp] * d16[i][8]) + (v1[5 + dvp] * d16[i][9]) + (v1[4 + dvp] * d16[i][10]) + (v1[3 + dvp] * d16[i][11]) + (v1[2 + dvp] * d16[i][12]) + (v1[1 + dvp] * d16[i][13]) + (v1[0 + dvp] * d16[i][14]) + (v1[15 + dvp] * d16[i][15])) * scalefactor);\r
1609 \r
1610         _tmpOut[i] = pcm_sample;\r
1611 \r
1612         dvp += 16;\r
1613       } // for\r
1614     } else {\r
1615       // final float[] vp = actual_v;\r
1616 \r
1617       // int inc = v_inc;\r
1618       // final float[] tmpOut = _tmpOut;\r
1619       int dvp = 0;\r
1620 \r
1621       // fat chance of having this loop unroll\r
1622       for (int i = 0; i < 32; i++) {\r
1623         // final float[] dp = d16[i];\r
1624         float pcm_sample;\r
1625 \r
1626         pcm_sample = (float) (((v2[14 + dvp] * d16[i][0]) + (v2[13 + dvp] * d16[i][1]) + (v2[12 + dvp] * d16[i][2]) + (v2[11 + dvp] * d16[i][3]) + (v2[10 + dvp] * d16[i][4]) + (v2[9 + dvp] * d16[i][5]) + (v2[8 + dvp] * d16[i][6]) + (v2[7 + dvp] * d16[i][7]) + (v2[6 + dvp] * d16[i][8]) + (v2[5 + dvp] * d16[i][9]) + (v2[4 + dvp] * d16[i][10]) + (v2[3 + dvp] * d16[i][11]) + (v2[2 + dvp] * d16[i][12]) + (v2[1 + dvp] * d16[i][13]) + (v2[0 + dvp] * d16[i][14]) + (v2[15 + dvp] * d16[i][15])) * scalefactor);\r
1627 \r
1628         _tmpOut[i] = pcm_sample;\r
1629 \r
1630         dvp += 16;\r
1631       } // for\r
1632     }\r
1633 \r
1634   }\r
1635 \r
1636   private void compute_pcm_samples15() {\r
1637     if (vidx == 1) {\r
1638       // final float[] vp = actual_v;\r
1639 \r
1640       // int inc = v_inc;\r
1641       // final float[] tmpOut = _tmpOut;\r
1642       int dvp = 0;\r
1643 \r
1644       // fat chance of having this loop unroll\r
1645       for (int i = 0; i < 32; i++) {\r
1646         float pcm_sample;\r
1647         // final float d16[i][] = d16[i];\r
1648         pcm_sample = (float) (((v1[15 + dvp] * d16[i][0]) + (v1[14 + dvp] * d16[i][1]) + (v1[13 + dvp] * d16[i][2]) + (v1[12 + dvp] * d16[i][3]) + (v1[11 + dvp] * d16[i][4]) + (v1[10 + dvp] * d16[i][5]) + (v1[9 + dvp] * d16[i][6]) + (v1[8 + dvp] * d16[i][7]) + (v1[7 + dvp] * d16[i][8]) + (v1[6 + dvp] * d16[i][9]) + (v1[5 + dvp] * d16[i][10]) + (v1[4 + dvp] * d16[i][11]) + (v1[3 + dvp] * d16[i][12]) + (v1[2 + dvp] * d16[i][13]) + (v1[1 + dvp] * d16[i][14]) + (v1[0 + dvp] * d16[i][15])) * scalefactor);\r
1649 \r
1650         _tmpOut[i] = pcm_sample;\r
1651         dvp += 16;\r
1652       } // for\r
1653     } else {\r
1654       // final float[] vp = actual_v;\r
1655 \r
1656       // int inc = v_inc;\r
1657       // final float[] tmpOut = _tmpOut;\r
1658       int dvp = 0;\r
1659 \r
1660       // fat chance of having this loop unroll\r
1661       for (int i = 0; i < 32; i++) {\r
1662         float pcm_sample;\r
1663         // final float d16[i][] = d16[i];\r
1664         pcm_sample = (float) (((v2[15 + dvp] * d16[i][0]) + (v2[14 + dvp] * d16[i][1]) + (v2[13 + dvp] * d16[i][2]) + (v2[12 + dvp] * d16[i][3]) + (v2[11 + dvp] * d16[i][4]) + (v2[10 + dvp] * d16[i][5]) + (v2[9 + dvp] * d16[i][6]) + (v2[8 + dvp] * d16[i][7]) + (v2[7 + dvp] * d16[i][8]) + (v2[6 + dvp] * d16[i][9]) + (v2[5 + dvp] * d16[i][10]) + (v2[4 + dvp] * d16[i][11]) + (v2[3 + dvp] * d16[i][12]) + (v2[2 + dvp] * d16[i][13]) + (v2[1 + dvp] * d16[i][14]) + (v2[0 + dvp] * d16[i][15])) * scalefactor);\r
1665 \r
1666         _tmpOut[i] = pcm_sample;\r
1667         dvp += 16;\r
1668       } // for\r
1669     }\r
1670 \r
1671   }\r
1672 \r
1673   private void compute_pcm_samples() {\r
1674 \r
1675     switch (actual_write_pos) {\r
1676     case 0:\r
1677       compute_pcm_samples0();\r
1678       break;\r
1679     case 1:\r
1680       compute_pcm_samples1();\r
1681       break;\r
1682     case 2:\r
1683       compute_pcm_samples2();\r
1684       break;\r
1685     case 3:\r
1686       compute_pcm_samples3();\r
1687       break;\r
1688     case 4:\r
1689       compute_pcm_samples4();\r
1690       break;\r
1691     case 5:\r
1692       compute_pcm_samples5();\r
1693       break;\r
1694     case 6:\r
1695       compute_pcm_samples6();\r
1696       break;\r
1697     case 7:\r
1698       compute_pcm_samples7();\r
1699       break;\r
1700     case 8:\r
1701       compute_pcm_samples8();\r
1702       break;\r
1703     case 9:\r
1704       compute_pcm_samples9();\r
1705       break;\r
1706     case 10:\r
1707       compute_pcm_samples10();\r
1708       break;\r
1709     case 11:\r
1710       compute_pcm_samples11();\r
1711       break;\r
1712     case 12:\r
1713       compute_pcm_samples12();\r
1714       break;\r
1715     case 13:\r
1716       compute_pcm_samples13();\r
1717       break;\r
1718     case 14:\r
1719       compute_pcm_samples14();\r
1720       break;\r
1721     case 15:\r
1722       compute_pcm_samples15();\r
1723       break;\r
1724     }\r
1725 \r
1726     // if (buffer != null) {\r
1727     // buffer.appendSamples(channel, _tmpOut);\r
1728     // }\r
1729     SampleBufferWrapper.appendSamples(channel, _tmpOut);\r
1730 \r
1731     /*\r
1732      * // MDM: I was considering putting in quality control for // low-spec\r
1733      * CPUs, but the performance gain (about 10-15%) // did not justify the\r
1734      * considerable drop in audio quality. switch (inc) { case 16:\r
1735      * buffer.appendSamples(channel, tmpOut); break; case 32: for (int i=0;\r
1736      * i<16; i++) { buffer.append(channel, (short)tmpOut[i]);\r
1737      * buffer.append(channel, (short)tmpOut[i]); } break; case 64: for (int i=0;\r
1738      * i<8; i++) { buffer.append(channel, (short)tmpOut[i]);\r
1739      * buffer.append(channel, (short)tmpOut[i]); buffer.append(channel,\r
1740      * (short)tmpOut[i]); buffer.append(channel, (short)tmpOut[i]); } break;\r
1741      * \r
1742      * }\r
1743      */\r
1744   }\r
1745 \r
1746   public void clear() {\r
1747     // clear out v1,v2\r
1748     SSJAVA.arrayinit(v1, 0);\r
1749     SSJAVA.arrayinit(v2, 0);\r
1750     // copy previous v1,v2\r
1751 \r
1752     for (int i = 0; i < prev1.length; i++) {\r
1753       v1[i] = prev1[i];\r
1754     }\r
1755     for (int i = 0; i < prev2.length; i++) {\r
1756       v2[i] = prev2[i];\r
1757     }\r
1758     // clear out previous buffer\r
1759     SSJAVA.arrayinit(prev1, 0);\r
1760     SSJAVA.arrayinit(prev2, 0);\r
1761     SSJAVA.arrayinit(samples, 0);\r
1762     SSJAVA.arrayinit(_tmpOut, 0);\r
1763   }\r
1764 \r
1765   /**\r
1766    * Calculate 32 PCM samples and put the into the Obuffer-object.\r
1767    */\r
1768 \r
1769   public void calculate_pcm_samples() {\r
1770 \r
1771     // System.out.println("#calculate_pcm_samples::actual_write_pos=" +\r
1772     // actual_write_pos);\r
1773 \r
1774     if (vidx == 1) {\r
1775       compute_new_v1_v2();\r
1776     } else {\r
1777       compute_new_v2_v1();\r
1778     }\r
1779 \r
1780     // System.out.println("1.actual_v=" + (actual_v == v1) + " vidx=" + vidx);\r
1781     // compute_new_v();\r
1782     // System.out.println("2.actual_v=" + (actual_v == v1) + " vidx=" + vidx);\r
1783     compute_pcm_samples();\r
1784     // System.out.println("3.actual_v=" + (actual_v == v1) + " vidx=" + vidx);\r
1785 \r
1786     actual_write_pos = (actual_write_pos + 1) & 0xf;\r
1787     // System.out.println("actual_write_pos="+actual_write_pos);\r
1788     // actual_v = (actual_v == v1) ? v2 : v1;\r
1789 \r
1790     if (vidx == 1) {\r
1791       vidx = 2;\r
1792     } else {\r
1793       vidx = 1;\r
1794     }\r
1795 \r
1796     // initialize samples[]:\r
1797     // for (register float *floatp = samples + 32; floatp > samples; )\r
1798     // *--floatp = 0.0f;\r
1799 \r
1800     // MDM: this may not be necessary. The Layer III decoder always\r
1801     // outputs 32 subband samples, but I haven't checked layer I & II.\r
1802     // for ( int p = 0; p < 32; p++){\r
1803     // samples[p] = 0.0f;\r
1804     // }\r
1805     SSJAVA.arrayinit(samples, 0);\r
1806   }\r
1807 \r
1808   private static final double MY_PI = 3.14159265358979323846;\r
1809   private static final float cos1_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI / 64.0)));\r
1810   private static final float cos3_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 3.0 / 64.0)));\r
1811   private static final float cos5_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 5.0 / 64.0)));\r
1812   private static final float cos7_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 7.0 / 64.0)));\r
1813   private static final float cos9_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 9.0 / 64.0)));\r
1814   private static final float cos11_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 11.0 / 64.0)));\r
1815   private static final float cos13_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 13.0 / 64.0)));\r
1816   private static final float cos15_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 15.0 / 64.0)));\r
1817   private static final float cos17_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 17.0 / 64.0)));\r
1818   private static final float cos19_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 19.0 / 64.0)));\r
1819   private static final float cos21_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 21.0 / 64.0)));\r
1820   private static final float cos23_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 23.0 / 64.0)));\r
1821   private static final float cos25_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 25.0 / 64.0)));\r
1822   private static final float cos27_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 27.0 / 64.0)));\r
1823   private static final float cos29_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 29.0 / 64.0)));\r
1824   private static final float cos31_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 31.0 / 64.0)));\r
1825   private static final float cos1_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI / 32.0)));\r
1826   private static final float cos3_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 3.0 / 32.0)));\r
1827   private static final float cos5_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 5.0 / 32.0)));\r
1828   private static final float cos7_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 7.0 / 32.0)));\r
1829   private static final float cos9_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 9.0 / 32.0)));\r
1830   private static final float cos11_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 11.0 / 32.0)));\r
1831   private static final float cos13_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 13.0 / 32.0)));\r
1832   private static final float cos15_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 15.0 / 32.0)));\r
1833   private static final float cos1_16 = (float) (1.0 / (2.0 * Math.cos(MY_PI / 16.0)));\r
1834   private static final float cos3_16 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 3.0 / 16.0)));\r
1835   private static final float cos5_16 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 5.0 / 16.0)));\r
1836   private static final float cos7_16 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 7.0 / 16.0)));\r
1837   private static final float cos1_8 = (float) (1.0 / (2.0 * Math.cos(MY_PI / 8.0)));\r
1838   private static final float cos3_8 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 3.0 / 8.0)));\r
1839   private static final float cos1_4 = (float) (1.0 / (2.0 * Math.cos(MY_PI / 4.0)));\r
1840 \r
1841   // Note: These values are not in the same order\r
1842   // as in Annex 3-B.3 of the ISO/IEC DIS 11172-3\r
1843   // private float d[] = {0.000000000, -4.000442505};\r
1844 \r
1845   /**\r
1846    * d[] split into subarrays of length 16. This provides for more faster access\r
1847    * by allowing a block of 16 to be addressed with constant offset.\r
1848    **/\r
1849   // TODO CONST\r
1850 \r
1851   private final static float d16[][] = null;\r
1852 \r
1853   /**\r
1854    * Converts a 1D array into a number of smaller arrays. This is used to\r
1855    * achieve offset + constant indexing into an array. Each sub-array represents\r
1856    * a block of values of the original array.\r
1857    * \r
1858    * @param array\r
1859    *          The array to split up into blocks.\r
1860    * @param blockSize\r
1861    *          The size of the blocks to split the array into. This must be an\r
1862    *          exact divisor of the length of the array, or some data will be\r
1863    *          lost from the main array.\r
1864    * \r
1865    * @return An array of arrays in which each element in the returned array will\r
1866    *         be of length <code>blockSize</code>.\r
1867    */\r
1868   static private float[][] splitArray(final float[] array, final int blockSize) {\r
1869     int size = array.length / blockSize;\r
1870     float[][] split = new float[size][];\r
1871     for (int i = 0; i < size; i++) {\r
1872       split[i] = subArray(array, i * blockSize, blockSize);\r
1873     }\r
1874     return split;\r
1875   }\r
1876 \r
1877   /**\r
1878    * Returns a subarray of an existing array.\r
1879    * \r
1880    * @param array\r
1881    *          The array to retrieve a subarra from.\r
1882    * @param offs\r
1883    *          The offset in the array that corresponds to the first index of the\r
1884    *          subarray.\r
1885    * @param len\r
1886    *          The number of indeces in the subarray.\r
1887    * @return The subarray, which may be of length 0.\r
1888    */\r
1889   static private float[] subArray(final float[] array, final int offs, int len) {\r
1890     if (offs + len > array.length) {\r
1891       len = array.length - offs;\r
1892     }\r
1893 \r
1894     if (len < 0)\r
1895       len = 0;\r
1896 \r
1897     float[] subarray = new float[len];\r
1898     for (int i = 0; i < len; i++) {\r
1899       subarray[i] = array[offs + i];\r
1900     }\r
1901 \r
1902     return subarray;\r
1903   }\r
1904 \r
1905   // The original data for d[]. This data is loaded from a file\r
1906   // to reduce the overall package size and to improve performance.\r
1907 \r
1908   static final float d[] = { 0.000000000f, -0.000442505f, 0.003250122f, -0.007003784f, 0.031082153f, -0.078628540f, 0.100311279f, -0.572036743f, 1.144989014f, 0.572036743f, 0.100311279f, 0.078628540f, 0.031082153f, 0.007003784f, 0.003250122f, 0.000442505f, -0.000015259f, -0.000473022f, 0.003326416f, -0.007919312f, 0.030517578f, -0.084182739f, 0.090927124f, -0.600219727f, 1.144287109f, 0.543823242f, 0.108856201f, 0.073059082f, 0.031478882f, 0.006118774f, 0.003173828f, 0.000396729f, -0.000015259f, -0.000534058f, 0.003387451f, -0.008865356f, 0.029785156f, -0.089706421f, 0.080688477f, -0.628295898f, 1.142211914f, 0.515609741f, 0.116577148f, 0.067520142f, 0.031738281f, 0.005294800f, 0.003082275f, 0.000366211f, -0.000015259f, -0.000579834f, 0.003433228f, -0.009841919f, 0.028884888f, -0.095169067f, 0.069595337f, -0.656219482f, 1.138763428f, 0.487472534f, 0.123474121f, 0.061996460f, 0.031845093f, 0.004486084f, 0.002990723f, 0.000320435f, -0.000015259f, -0.000625610f, 0.003463745f, -0.010848999f, 0.027801514f, -0.100540161f, 0.057617188f, -0.683914185f, 1.133926392f, 0.459472656f, 0.129577637f, 0.056533813f, 0.031814575f, 0.003723145f, 0.002899170f, 0.000289917f, -0.000015259f, -0.000686646f, 0.003479004f, -0.011886597f, 0.026535034f, -0.105819702f, 0.044784546f, -0.711318970f, 1.127746582f, 0.431655884f, 0.134887695f, 0.051132202f, 0.031661987f, 0.003005981f, 0.002792358f, 0.000259399f, -0.000015259f, -0.000747681f, 0.003479004f, -0.012939453f, 0.025085449f, -0.110946655f, 0.031082153f, -0.738372803f, 1.120223999f, 0.404083252f, 0.139450073f, 0.045837402f, 0.031387329f, 0.002334595f, 0.002685547f, 0.000244141f, -0.000030518f, -0.000808716f, 0.003463745f, -0.014022827f, 0.023422241f, -0.115921021f, 0.016510010f, -0.765029907f, 1.111373901f, 0.376800537f, 0.143264771f, 0.040634155f, 0.031005859f, 0.001693726f, 0.002578735f, 0.000213623f, -0.000030518f, -0.000885010f, 0.003417969f, -0.015121460f, 0.021575928f, -0.120697021f, 0.001068115f, -0.791213989f, 1.101211548f, 0.349868774f, 0.146362305f, 0.035552979f, 0.030532837f, 0.001098633f, 0.002456665f, 0.000198364f, -0.000030518f, -0.000961304f, 0.003372192f, -0.016235352f, 0.019531250f, -0.125259399f, -0.015228271f, -0.816864014f, 1.089782715f, 0.323318481f, 0.148773193f, 0.030609131f, 0.029937744f, 0.000549316f, 0.002349854f, 0.000167847f, -0.000030518f, -0.001037598f, 0.003280640f, -0.017349243f, 0.017257690f, -0.129562378f, -0.032379150f, -0.841949463f, 1.077117920f, 0.297210693f, 0.150497437f, 0.025817871f, 0.029281616f, 0.000030518f, 0.002243042f, 0.000152588f, -0.000045776f, -0.001113892f, 0.003173828f, -0.018463135f, 0.014801025f, -0.133590698f, -0.050354004f, -0.866363525f, 1.063217163f, 0.271591187f, 0.151596069f, 0.021179199f, 0.028533936f, -0.000442505f, 0.002120972f, 0.000137329f, -0.000045776f, -0.001205444f, 0.003051758f, -0.019577026f, 0.012115479f, -0.137298584f, -0.069168091f, -0.890090942f, 1.048156738f, 0.246505737f, 0.152069092f, 0.016708374f, 0.027725220f, -0.000869751f, 0.002014160f, 0.000122070f, -0.000061035f, -0.001296997f, 0.002883911f, -0.020690918f, 0.009231567f, -0.140670776f, -0.088775635f, -0.913055420f, 1.031936646f, 0.221984863f, 0.151962280f, 0.012420654f, 0.026840210f, -0.001266479f, 0.001907349f, 0.000106812f, -0.000061035f, -0.001388550f, 0.002700806f, -0.021789551f, 0.006134033f, -0.143676758f, -0.109161377f, -0.935195923f, 1.014617920f, 0.198059082f, 0.151306152f, 0.008316040f, 0.025909424f, -0.001617432f, 0.001785278f, 0.000106812f, -0.000076294f, -0.001480103f, 0.002487183f, -0.022857666f, 0.002822876f, -0.146255493f, -0.130310059f, -0.956481934f, 0.996246338f, 0.174789429f, 0.150115967f, 0.004394531f, 0.024932861f, -0.001937866f, 0.001693726f, 0.000091553f, -0.000076294f, -0.001586914f, 0.002227783f, -0.023910522f, -0.000686646f, -0.148422241f, -0.152206421f, -0.976852417f, 0.976852417f, 0.152206421f, 0.148422241f, 0.000686646f, 0.023910522f, -0.002227783f, 0.001586914f, 0.000076294f, -0.000091553f, -0.001693726f, 0.001937866f, -0.024932861f, -0.004394531f, -0.150115967f, -0.174789429f, -0.996246338f, 0.956481934f, 0.130310059f, 0.146255493f, -0.002822876f, 0.022857666f, -0.002487183f, 0.001480103f, 0.000076294f, -0.000106812f, -0.001785278f, 0.001617432f, -0.025909424f, -0.008316040f, -0.151306152f, -0.198059082f, -1.014617920f, 0.935195923f, 0.109161377f, 0.143676758f, -0.006134033f, 0.021789551f, -0.002700806f, 0.001388550f, 0.000061035f, -0.000106812f, -0.001907349f, 0.001266479f, -0.026840210f, -0.012420654f, -0.151962280f, -0.221984863f, -1.031936646f, 0.913055420f, 0.088775635f, 0.140670776f, -0.009231567f, 0.020690918f, -0.002883911f, 0.001296997f, 0.000061035f, -0.000122070f, -0.002014160f, 0.000869751f, -0.027725220f, -0.016708374f, -0.152069092f, -0.246505737f, -1.048156738f, 0.890090942f, 0.069168091f, 0.137298584f, -0.012115479f, 0.019577026f, -0.003051758f, 0.001205444f, 0.000045776f, -0.000137329f, -0.002120972f, 0.000442505f, -0.028533936f, -0.021179199f, -0.151596069f, -0.271591187f, -1.063217163f, 0.866363525f, 0.050354004f, 0.133590698f, -0.014801025f, 0.018463135f, -0.003173828f, 0.001113892f, 0.000045776f, -0.000152588f, -0.002243042f, -0.000030518f, -0.029281616f, -0.025817871f, -0.150497437f, -0.297210693f, -1.077117920f, 0.841949463f, 0.032379150f, 0.129562378f, -0.017257690f, 0.017349243f, -0.003280640f, 0.001037598f, 0.000030518f, -0.000167847f, -0.002349854f, -0.000549316f, -0.029937744f, -0.030609131f, -0.148773193f, -0.323318481f, -1.089782715f, 0.816864014f, 0.015228271f, 0.125259399f, -0.019531250f, 0.016235352f, -0.003372192f, 0.000961304f, 0.000030518f, -0.000198364f, -0.002456665f, -0.001098633f, -0.030532837f, -0.035552979f, -0.146362305f, -0.349868774f, -1.101211548f, 0.791213989f, -0.001068115f, 0.120697021f, -0.021575928f, 0.015121460f, -0.003417969f, 0.000885010f, 0.000030518f, -0.000213623f, -0.002578735f, -0.001693726f, -0.031005859f, -0.040634155f, -0.143264771f, -0.376800537f, -1.111373901f, 0.765029907f, -0.016510010f, 0.115921021f, -0.023422241f, 0.014022827f, -0.003463745f, 0.000808716f, 0.000030518f, -0.000244141f, -0.002685547f, -0.002334595f, -0.031387329f, -0.045837402f, -0.139450073f, -0.404083252f, -1.120223999f, 0.738372803f, -0.031082153f, 0.110946655f, -0.025085449f, 0.012939453f, -0.003479004f, 0.000747681f, 0.000015259f, -0.000259399f, -0.002792358f, -0.003005981f, -0.031661987f, -0.051132202f, -0.134887695f, -0.431655884f, -1.127746582f, 0.711318970f, -0.044784546f, 0.105819702f, -0.026535034f, 0.011886597f, -0.003479004f, 0.000686646f, 0.000015259f, -0.000289917f, -0.002899170f, -0.003723145f, -0.031814575f, -0.056533813f, -0.129577637f, -0.459472656f, -1.133926392f, 0.683914185f, -0.057617188f, 0.100540161f, -0.027801514f, 0.010848999f, -0.003463745f, 0.000625610f, 0.000015259f, -0.000320435f, -0.002990723f, -0.004486084f, -0.031845093f, -0.061996460f, -0.123474121f, -0.487472534f, -1.138763428f, 0.656219482f, -0.069595337f, 0.095169067f, -0.028884888f, 0.009841919f, -0.003433228f, 0.000579834f, 0.000015259f, -0.000366211f, -0.003082275f, -0.005294800f, -0.031738281f, -0.067520142f, -0.116577148f, -0.515609741f, -1.142211914f, 0.628295898f, -0.080688477f, 0.089706421f, -0.029785156f, 0.008865356f, -0.003387451f, 0.000534058f, 0.000015259f, -0.000396729f, -0.003173828f, -0.006118774f, -0.031478882f, -0.073059082f, -0.108856201f, -0.543823242f, -1.144287109f, 0.600219727f, -0.090927124f, 0.084182739f, -0.030517578f, 0.007919312f, -0.003326416f, 0.000473022f, 0.000015259f };\r
1909 \r
1910 }\r