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