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 final class SynthesisFilter {\r
  38   private float[] v1;\r
  39   private float[] v2;\r
  40   private float[] actual_v; // v1 or v2\r
  41   private int actual_write_pos; // 0-15\r
  42   private float[] samples; // 32 new subband samples\r
  43   private int channel;\r
  44   private float scalefactor;\r
  45   private float[] eq;\r
  46 \r
  47   /**\r
  48    * Quality value for controlling CPU usage/quality tradeoff.\r
  49    */\r
  50   /*\r
  51    * private int quality;\r
  52    * \r
  53    * private int v_inc;\r
  54    * \r
  55    * \r
  56    * \r
  57    * public static final int HIGH_QUALITY = 1; public static final int\r
  58    * MEDIUM_QUALITY = 2; public static final int LOW_QUALITY = 4;\r
  59    */\r
  60 \r
  61   /**\r
  62    * Contructor. The scalefactor scales the calculated float pcm samples to\r
  63    * short values (raw pcm samples are in [-1.0, 1.0], if no violations occur).\r
  64    */\r
  65   public SynthesisFilter(int channelnumber, float factor, float[] eq0) {\r
  66     \r
  67     d16 = splitArray(d, 16);\r
  68 \r
  69     v1 = new float[512];\r
  70     v2 = new float[512];\r
  71     samples = new float[32];\r
  72     channel = channelnumber;\r
  73     scalefactor = factor;\r
  74     setEQ(eq);\r
  75     // setQuality(HIGH_QUALITY);\r
  76 \r
  77     reset();\r
  78   }\r
  79 \r
  80   public void setEQ(float[] eq0) {\r
  81     this.eq = eq0;\r
  82     if (eq == null) {\r
  83       eq = new float[32];\r
  84       for (int i = 0; i < 32; i++)\r
  85         eq[i] = 1.0f;\r
  86     }\r
  87     if (eq.length < 32) {\r
  88       throw new IllegalArgumentException("eq0");\r
  89     }\r
  90 \r
  91   }\r
  92 \r
  93   /*\r
  94    * private void setQuality(int quality0) { switch (quality0) { case\r
  95    * HIGH_QUALITY: case MEDIUM_QUALITY: case LOW_QUALITY: v_inc = 16 * quality0;\r
  96    * quality = quality0; break; default : throw new\r
  97    * IllegalArgumentException("Unknown quality value"); } }\r
  98    * \r
  99    * public int getQuality() { return quality; }\r
 100    */\r
 101 \r
 102   /**\r
 103    * Reset the synthesis filter.\r
 104    */\r
 105   public void reset() {\r
 106     // float[] floatp;\r
 107     // float[] floatp2;\r
 108 \r
 109     // initialize v1[] and v2[]:\r
 110     // for (floatp = v1 + 512, floatp2 = v2 + 512; floatp > v1; )\r
 111     // *--floatp = *--floatp2 = 0.0;\r
 112     for (int p = 0; p < 512; p++)\r
 113       v1[p] = v2[p] = 0.0f;\r
 114 \r
 115     // initialize samples[]:\r
 116     // for (floatp = samples + 32; floatp > samples; )\r
 117     // *--floatp = 0.0;\r
 118     for (int p2 = 0; p2 < 32; p2++)\r
 119       samples[p2] = 0.0f;\r
 120 \r
 121     actual_v = v1;\r
 122     actual_write_pos = 15;\r
 123   }\r
 124 \r
 125   /**\r
 126    * Inject Sample.\r
 127    */\r
 128   public void input_sample(float sample, int subbandnumber) {\r
 129     samples[subbandnumber] = eq[subbandnumber] * sample;\r
 130   }\r
 131 \r
 132   public void input_samples(float[] s) {\r
 133     for (int i = 31; i >= 0; i--) {\r
 134       samples[i] = s[i] * eq[i];\r
 135     }\r
 136   }\r
 137 \r
 138   /**\r
 139    * Compute new values via a fast cosine transform.\r
 140    */\r
 141   private void compute_new_v() {\r
 142     // p is fully initialized from x1\r
 143     // float[] p = _p;\r
 144     // pp is fully initialized from p\r
 145     // float[] pp = _pp;\r
 146 \r
 147     // float[] new_v = _new_v;\r
 148 \r
 149     // float[] new_v = new float[32]; // new V[0-15] and V[33-48] of Figure\r
 150     // 3-A.2 in ISO DIS 11172-3\r
 151     // float[] p = new float[16];\r
 152     // float[] pp = new float[16];\r
 153 \r
 154     /*\r
 155      * for (int i=31; i>=0; i--) { new_v[i] = 0.0f; }\r
 156      */\r
 157 \r
 158     float new_v0, new_v1, new_v2, new_v3, new_v4, new_v5, new_v6, new_v7, new_v8, new_v9;\r
 159     float new_v10, new_v11, new_v12, new_v13, new_v14, new_v15, new_v16, new_v17, new_v18, new_v19;\r
 160     float new_v20, new_v21, new_v22, new_v23, new_v24, new_v25, new_v26, new_v27, new_v28, new_v29;\r
 161     float new_v30, new_v31;\r
 162 \r
 163     new_v0 =\r
 164         new_v1 =\r
 165             new_v2 =\r
 166                 new_v3 =\r
 167                     new_v4 =\r
 168                         new_v5 =\r
 169                             new_v6 =\r
 170                                 new_v7 =\r
 171                                     new_v8 =\r
 172                                         new_v9 =\r
 173                                             new_v10 =\r
 174                                                 new_v11 =\r
 175                                                     new_v12 =\r
 176                                                         new_v13 =\r
 177                                                             new_v14 =\r
 178                                                                 new_v15 =\r
 179                                                                     new_v16 =\r
 180                                                                         new_v17 =\r
 181                                                                             new_v18 =\r
 182                                                                                 new_v19 =\r
 183                                                                                     new_v20 =\r
 184                                                                                         new_v21 =\r
 185                                                                                             new_v22 =\r
 186                                                                                                 new_v23 =\r
 187                                                                                                     new_v24 =\r
 188                                                                                                         new_v25 =\r
 189                                                                                                             new_v26 =\r
 190                                                                                                                 new_v27 =\r
 191                                                                                                                     new_v28 =\r
 192                                                                                                                         new_v29 =\r
 193                                                                                                                             new_v30 =\r
 194                                                                                                                                 new_v31 =\r
 195                                                                                                                                     0.0f;\r
 196 \r
 197     // float[] new_v = new float[32]; // new V[0-15] and V[33-48] of Figure\r
 198     // 3-A.2 in ISO DIS 11172-3\r
 199     // float[] p = new float[16];\r
 200     // float[] pp = new float[16];\r
 201 \r
 202     float[] s = samples;\r
 203 \r
 204     float s0 = s[0];\r
 205     float s1 = s[1];\r
 206     float s2 = s[2];\r
 207     float s3 = s[3];\r
 208     float s4 = s[4];\r
 209     float s5 = s[5];\r
 210     float s6 = s[6];\r
 211     float s7 = s[7];\r
 212     float s8 = s[8];\r
 213     float s9 = s[9];\r
 214     float s10 = s[10];\r
 215     float s11 = s[11];\r
 216     float s12 = s[12];\r
 217     float s13 = s[13];\r
 218     float s14 = s[14];\r
 219     float s15 = s[15];\r
 220     float s16 = s[16];\r
 221     float s17 = s[17];\r
 222     float s18 = s[18];\r
 223     float s19 = s[19];\r
 224     float s20 = s[20];\r
 225     float s21 = s[21];\r
 226     float s22 = s[22];\r
 227     float s23 = s[23];\r
 228     float s24 = s[24];\r
 229     float s25 = s[25];\r
 230     float s26 = s[26];\r
 231     float s27 = s[27];\r
 232     float s28 = s[28];\r
 233     float s29 = s[29];\r
 234     float s30 = s[30];\r
 235     float s31 = s[31];\r
 236 \r
 237     float p0 = s0 + s31;\r
 238     float p1 = s1 + s30;\r
 239     float p2 = s2 + s29;\r
 240     float p3 = s3 + s28;\r
 241     float p4 = s4 + s27;\r
 242     float p5 = s5 + s26;\r
 243     float p6 = s6 + s25;\r
 244     float p7 = s7 + s24;\r
 245     float p8 = s8 + s23;\r
 246     float p9 = s9 + s22;\r
 247     float p10 = s10 + s21;\r
 248     float p11 = s11 + s20;\r
 249     float p12 = s12 + s19;\r
 250     float p13 = s13 + s18;\r
 251     float p14 = s14 + s17;\r
 252     float p15 = s15 + s16;\r
 253 \r
 254     float pp0 = p0 + p15;\r
 255     float pp1 = p1 + p14;\r
 256     float pp2 = p2 + p13;\r
 257     float pp3 = p3 + p12;\r
 258     float pp4 = p4 + p11;\r
 259     float pp5 = p5 + p10;\r
 260     float pp6 = p6 + p9;\r
 261     float pp7 = p7 + p8;\r
 262     float pp8 = (p0 - p15) * cos1_32;\r
 263     float pp9 = (p1 - p14) * cos3_32;\r
 264     float pp10 = (p2 - p13) * cos5_32;\r
 265     float pp11 = (p3 - p12) * cos7_32;\r
 266     float pp12 = (p4 - p11) * cos9_32;\r
 267     float pp13 = (p5 - p10) * cos11_32;\r
 268     float pp14 = (p6 - p9) * cos13_32;\r
 269     float pp15 = (p7 - p8) * cos15_32;\r
 270 \r
 271     p0 = pp0 + pp7;\r
 272     p1 = pp1 + pp6;\r
 273     p2 = pp2 + pp5;\r
 274     p3 = pp3 + pp4;\r
 275     p4 = (pp0 - pp7) * cos1_16;\r
 276     p5 = (pp1 - pp6) * cos3_16;\r
 277     p6 = (pp2 - pp5) * cos5_16;\r
 278     p7 = (pp3 - pp4) * cos7_16;\r
 279     p8 = pp8 + pp15;\r
 280     p9 = pp9 + pp14;\r
 281     p10 = pp10 + pp13;\r
 282     p11 = pp11 + pp12;\r
 283     p12 = (pp8 - pp15) * cos1_16;\r
 284     p13 = (pp9 - pp14) * cos3_16;\r
 285     p14 = (pp10 - pp13) * cos5_16;\r
 286     p15 = (pp11 - pp12) * cos7_16;\r
 287 \r
 288     pp0 = p0 + p3;\r
 289     pp1 = p1 + p2;\r
 290     pp2 = (p0 - p3) * cos1_8;\r
 291     pp3 = (p1 - p2) * cos3_8;\r
 292     pp4 = p4 + p7;\r
 293     pp5 = p5 + p6;\r
 294     pp6 = (p4 - p7) * cos1_8;\r
 295     pp7 = (p5 - p6) * cos3_8;\r
 296     pp8 = p8 + p11;\r
 297     pp9 = p9 + p10;\r
 298     pp10 = (p8 - p11) * cos1_8;\r
 299     pp11 = (p9 - p10) * cos3_8;\r
 300     pp12 = p12 + p15;\r
 301     pp13 = p13 + p14;\r
 302     pp14 = (p12 - p15) * cos1_8;\r
 303     pp15 = (p13 - p14) * cos3_8;\r
 304 \r
 305     p0 = pp0 + pp1;\r
 306     p1 = (pp0 - pp1) * cos1_4;\r
 307     p2 = pp2 + pp3;\r
 308     p3 = (pp2 - pp3) * cos1_4;\r
 309     p4 = pp4 + pp5;\r
 310     p5 = (pp4 - pp5) * cos1_4;\r
 311     p6 = pp6 + pp7;\r
 312     p7 = (pp6 - pp7) * cos1_4;\r
 313     p8 = pp8 + pp9;\r
 314     p9 = (pp8 - pp9) * cos1_4;\r
 315     p10 = pp10 + pp11;\r
 316     p11 = (pp10 - pp11) * cos1_4;\r
 317     p12 = pp12 + pp13;\r
 318     p13 = (pp12 - pp13) * cos1_4;\r
 319     p14 = pp14 + pp15;\r
 320     p15 = (pp14 - pp15) * cos1_4;\r
 321 \r
 322     // this is pretty insane coding\r
 323     float tmp1;\r
 324     new_v19/* 36-17 */= -(new_v4 = (new_v12 = p7) + p5) - p6;\r
 325     new_v27/* 44-17 */= -p6 - p7 - p4;\r
 326     new_v6 = (new_v10 = (new_v14 = p15) + p11) + p13;\r
 327     new_v17/* 34-17 */= -(new_v2 = p15 + p13 + p9) - p14;\r
 328     new_v21/* 38-17 */= (tmp1 = -p14 - p15 - p10 - p11) - p13;\r
 329     new_v29/* 46-17 */= -p14 - p15 - p12 - p8;\r
 330     new_v25/* 42-17 */= tmp1 - p12;\r
 331     new_v31/* 48-17 */= -p0;\r
 332     new_v0 = p1;\r
 333     new_v23/* 40-17 */= -(new_v8 = p3) - p2;\r
 334 \r
 335     p0 = (s0 - s31) * cos1_64;\r
 336     p1 = (s1 - s30) * cos3_64;\r
 337     p2 = (s2 - s29) * cos5_64;\r
 338     p3 = (s3 - s28) * cos7_64;\r
 339     p4 = (s4 - s27) * cos9_64;\r
 340     p5 = (s5 - s26) * cos11_64;\r
 341     p6 = (s6 - s25) * cos13_64;\r
 342     p7 = (s7 - s24) * cos15_64;\r
 343     p8 = (s8 - s23) * cos17_64;\r
 344     p9 = (s9 - s22) * cos19_64;\r
 345     p10 = (s10 - s21) * cos21_64;\r
 346     p11 = (s11 - s20) * cos23_64;\r
 347     p12 = (s12 - s19) * cos25_64;\r
 348     p13 = (s13 - s18) * cos27_64;\r
 349     p14 = (s14 - s17) * cos29_64;\r
 350     p15 = (s15 - s16) * cos31_64;\r
 351 \r
 352     pp0 = p0 + p15;\r
 353     pp1 = p1 + p14;\r
 354     pp2 = p2 + p13;\r
 355     pp3 = p3 + p12;\r
 356     pp4 = p4 + p11;\r
 357     pp5 = p5 + p10;\r
 358     pp6 = p6 + p9;\r
 359     pp7 = p7 + p8;\r
 360     pp8 = (p0 - p15) * cos1_32;\r
 361     pp9 = (p1 - p14) * cos3_32;\r
 362     pp10 = (p2 - p13) * cos5_32;\r
 363     pp11 = (p3 - p12) * cos7_32;\r
 364     pp12 = (p4 - p11) * cos9_32;\r
 365     pp13 = (p5 - p10) * cos11_32;\r
 366     pp14 = (p6 - p9) * cos13_32;\r
 367     pp15 = (p7 - p8) * cos15_32;\r
 368 \r
 369     p0 = pp0 + pp7;\r
 370     p1 = pp1 + pp6;\r
 371     p2 = pp2 + pp5;\r
 372     p3 = pp3 + pp4;\r
 373     p4 = (pp0 - pp7) * cos1_16;\r
 374     p5 = (pp1 - pp6) * cos3_16;\r
 375     p6 = (pp2 - pp5) * cos5_16;\r
 376     p7 = (pp3 - pp4) * cos7_16;\r
 377     p8 = pp8 + pp15;\r
 378     p9 = pp9 + pp14;\r
 379     p10 = pp10 + pp13;\r
 380     p11 = pp11 + pp12;\r
 381     p12 = (pp8 - pp15) * cos1_16;\r
 382     p13 = (pp9 - pp14) * cos3_16;\r
 383     p14 = (pp10 - pp13) * cos5_16;\r
 384     p15 = (pp11 - pp12) * cos7_16;\r
 385 \r
 386     pp0 = p0 + p3;\r
 387     pp1 = p1 + p2;\r
 388     pp2 = (p0 - p3) * cos1_8;\r
 389     pp3 = (p1 - p2) * cos3_8;\r
 390     pp4 = p4 + p7;\r
 391     pp5 = p5 + p6;\r
 392     pp6 = (p4 - p7) * cos1_8;\r
 393     pp7 = (p5 - p6) * cos3_8;\r
 394     pp8 = p8 + p11;\r
 395     pp9 = p9 + p10;\r
 396     pp10 = (p8 - p11) * cos1_8;\r
 397     pp11 = (p9 - p10) * cos3_8;\r
 398     pp12 = p12 + p15;\r
 399     pp13 = p13 + p14;\r
 400     pp14 = (p12 - p15) * cos1_8;\r
 401     pp15 = (p13 - p14) * cos3_8;\r
 402 \r
 403     p0 = pp0 + pp1;\r
 404     p1 = (pp0 - pp1) * cos1_4;\r
 405     p2 = pp2 + pp3;\r
 406     p3 = (pp2 - pp3) * cos1_4;\r
 407     p4 = pp4 + pp5;\r
 408     p5 = (pp4 - pp5) * cos1_4;\r
 409     p6 = pp6 + pp7;\r
 410     p7 = (pp6 - pp7) * cos1_4;\r
 411     p8 = pp8 + pp9;\r
 412     p9 = (pp8 - pp9) * cos1_4;\r
 413     p10 = pp10 + pp11;\r
 414     p11 = (pp10 - pp11) * cos1_4;\r
 415     p12 = pp12 + pp13;\r
 416     p13 = (pp12 - pp13) * cos1_4;\r
 417     p14 = pp14 + pp15;\r
 418     p15 = (pp14 - pp15) * cos1_4;\r
 419 \r
 420     // manually doing something that a compiler should handle sucks\r
 421     // coding like this is hard to read\r
 422     float tmp2;\r
 423     new_v5 = (new_v11 = (new_v13 = (new_v15 = p15) + p7) + p11) + p5 + p13;\r
 424     new_v7 = (new_v9 = p15 + p11 + p3) + p13;\r
 425     new_v16/* 33-17 */= -(new_v1 = (tmp1 = p13 + p15 + p9) + p1) - p14;\r
 426     new_v18/* 35-17 */= -(new_v3 = tmp1 + p5 + p7) - p6 - p14;\r
 427 \r
 428     new_v22/* 39-17 */= (tmp1 = -p10 - p11 - p14 - p15) - p13 - p2 - p3;\r
 429     new_v20/* 37-17 */= tmp1 - p13 - p5 - p6 - p7;\r
 430     new_v24/* 41-17 */= tmp1 - p12 - p2 - p3;\r
 431     new_v26/* 43-17 */= tmp1 - p12 - (tmp2 = p4 + p6 + p7);\r
 432     new_v30/* 47-17 */= (tmp1 = -p8 - p12 - p14 - p15) - p0;\r
 433     new_v28/* 45-17 */= tmp1 - tmp2;\r
 434 \r
 435     // insert V[0-15] (== new_v[0-15]) into actual v:\r
 436     // float[] x2 = actual_v + actual_write_pos;\r
 437     float dest[] = actual_v;\r
 438 \r
 439     int pos = actual_write_pos;\r
 440 \r
 441     dest[0 + pos] = new_v0;\r
 442     dest[16 + pos] = new_v1;\r
 443     dest[32 + pos] = new_v2;\r
 444     dest[48 + pos] = new_v3;\r
 445     dest[64 + pos] = new_v4;\r
 446     dest[80 + pos] = new_v5;\r
 447     dest[96 + pos] = new_v6;\r
 448     dest[112 + pos] = new_v7;\r
 449     dest[128 + pos] = new_v8;\r
 450     dest[144 + pos] = new_v9;\r
 451     dest[160 + pos] = new_v10;\r
 452     dest[176 + pos] = new_v11;\r
 453     dest[192 + pos] = new_v12;\r
 454     dest[208 + pos] = new_v13;\r
 455     dest[224 + pos] = new_v14;\r
 456     dest[240 + pos] = new_v15;\r
 457 \r
 458     // V[16] is always 0.0:\r
 459     dest[256 + pos] = 0.0f;\r
 460 \r
 461     // insert V[17-31] (== -new_v[15-1]) into actual v:\r
 462     dest[272 + pos] = -new_v15;\r
 463     dest[288 + pos] = -new_v14;\r
 464     dest[304 + pos] = -new_v13;\r
 465     dest[320 + pos] = -new_v12;\r
 466     dest[336 + pos] = -new_v11;\r
 467     dest[352 + pos] = -new_v10;\r
 468     dest[368 + pos] = -new_v9;\r
 469     dest[384 + pos] = -new_v8;\r
 470     dest[400 + pos] = -new_v7;\r
 471     dest[416 + pos] = -new_v6;\r
 472     dest[432 + pos] = -new_v5;\r
 473     dest[448 + pos] = -new_v4;\r
 474     dest[464 + pos] = -new_v3;\r
 475     dest[480 + pos] = -new_v2;\r
 476     dest[496 + pos] = -new_v1;\r
 477 \r
 478     // insert V[32] (== -new_v[0]) into other v:\r
 479     dest = (actual_v == v1) ? v2 : v1;\r
 480 \r
 481     dest[0 + pos] = -new_v0;\r
 482     // insert V[33-48] (== new_v[16-31]) into other v:\r
 483     dest[16 + pos] = new_v16;\r
 484     dest[32 + pos] = new_v17;\r
 485     dest[48 + pos] = new_v18;\r
 486     dest[64 + pos] = new_v19;\r
 487     dest[80 + pos] = new_v20;\r
 488     dest[96 + pos] = new_v21;\r
 489     dest[112 + pos] = new_v22;\r
 490     dest[128 + pos] = new_v23;\r
 491     dest[144 + pos] = new_v24;\r
 492     dest[160 + pos] = new_v25;\r
 493     dest[176 + pos] = new_v26;\r
 494     dest[192 + pos] = new_v27;\r
 495     dest[208 + pos] = new_v28;\r
 496     dest[224 + pos] = new_v29;\r
 497     dest[240 + pos] = new_v30;\r
 498     dest[256 + pos] = new_v31;\r
 499 \r
 500     // insert V[49-63] (== new_v[30-16]) into other v:\r
 501     dest[272 + pos] = new_v30;\r
 502     dest[288 + pos] = new_v29;\r
 503     dest[304 + pos] = new_v28;\r
 504     dest[320 + pos] = new_v27;\r
 505     dest[336 + pos] = new_v26;\r
 506     dest[352 + pos] = new_v25;\r
 507     dest[368 + pos] = new_v24;\r
 508     dest[384 + pos] = new_v23;\r
 509     dest[400 + pos] = new_v22;\r
 510     dest[416 + pos] = new_v21;\r
 511     dest[432 + pos] = new_v20;\r
 512     dest[448 + pos] = new_v19;\r
 513     dest[464 + pos] = new_v18;\r
 514     dest[480 + pos] = new_v17;\r
 515     dest[496 + pos] = new_v16;\r
 516     /*\r
 517      * } else { v1[0 + actual_write_pos] = -new_v0; // insert V[33-48] (==\r
 518      * new_v[16-31]) into other v: v1[16 + actual_write_pos] = new_v16; v1[32 +\r
 519      * actual_write_pos] = new_v17; v1[48 + actual_write_pos] = new_v18; v1[64 +\r
 520      * actual_write_pos] = new_v19; v1[80 + actual_write_pos] = new_v20; v1[96 +\r
 521      * actual_write_pos] = new_v21; v1[112 + actual_write_pos] = new_v22; v1[128\r
 522      * + actual_write_pos] = new_v23; v1[144 + actual_write_pos] = new_v24;\r
 523      * v1[160 + actual_write_pos] = new_v25; v1[176 + actual_write_pos] =\r
 524      * new_v26; v1[192 + actual_write_pos] = new_v27; v1[208 + actual_write_pos]\r
 525      * = new_v28; v1[224 + actual_write_pos] = new_v29; v1[240 +\r
 526      * actual_write_pos] = new_v30; v1[256 + actual_write_pos] = new_v31;\r
 527      * \r
 528      * // insert V[49-63] (== new_v[30-16]) into other v: v1[272 +\r
 529      * actual_write_pos] = new_v30; v1[288 + actual_write_pos] = new_v29; v1[304\r
 530      * + actual_write_pos] = new_v28; v1[320 + actual_write_pos] = new_v27;\r
 531      * v1[336 + actual_write_pos] = new_v26; v1[352 + actual_write_pos] =\r
 532      * new_v25; v1[368 + actual_write_pos] = new_v24; v1[384 + actual_write_pos]\r
 533      * = new_v23; v1[400 + actual_write_pos] = new_v22; v1[416 +\r
 534      * actual_write_pos] = new_v21; v1[432 + actual_write_pos] = new_v20; v1[448\r
 535      * + actual_write_pos] = new_v19; v1[464 + actual_write_pos] = new_v18;\r
 536      * v1[480 + actual_write_pos] = new_v17; v1[496 + actual_write_pos] =\r
 537      * new_v16; }\r
 538      */\r
 539   }\r
 540 \r
 541   /**\r
 542    * Compute new values via a fast cosine transform.\r
 543    */\r
 544   private void compute_new_v_old() {\r
 545     // p is fully initialized from x1\r
 546     // float[] p = _p;\r
 547     // pp is fully initialized from p\r
 548     // float[] pp = _pp;\r
 549 \r
 550     // float[] new_v = _new_v;\r
 551 \r
 552     float[] new_v = new float[32]; // new V[0-15] and V[33-48] of Figure 3-A.2\r
 553                                    // in ISO DIS 11172-3\r
 554     float[] p = new float[16];\r
 555     float[] pp = new float[16];\r
 556 \r
 557     for (int i = 31; i >= 0; i--) {\r
 558       new_v[i] = 0.0f;\r
 559     }\r
 560 \r
 561     // float[] new_v = new float[32]; // new V[0-15] and V[33-48] of Figure\r
 562     // 3-A.2 in ISO DIS 11172-3\r
 563     // float[] p = new float[16];\r
 564     // float[] pp = new float[16];\r
 565 \r
 566     float[] x1 = samples;\r
 567 \r
 568     p[0] = x1[0] + x1[31];\r
 569     p[1] = x1[1] + x1[30];\r
 570     p[2] = x1[2] + x1[29];\r
 571     p[3] = x1[3] + x1[28];\r
 572     p[4] = x1[4] + x1[27];\r
 573     p[5] = x1[5] + x1[26];\r
 574     p[6] = x1[6] + x1[25];\r
 575     p[7] = x1[7] + x1[24];\r
 576     p[8] = x1[8] + x1[23];\r
 577     p[9] = x1[9] + x1[22];\r
 578     p[10] = x1[10] + x1[21];\r
 579     p[11] = x1[11] + x1[20];\r
 580     p[12] = x1[12] + x1[19];\r
 581     p[13] = x1[13] + x1[18];\r
 582     p[14] = x1[14] + x1[17];\r
 583     p[15] = x1[15] + x1[16];\r
 584 \r
 585     pp[0] = p[0] + p[15];\r
 586     pp[1] = p[1] + p[14];\r
 587     pp[2] = p[2] + p[13];\r
 588     pp[3] = p[3] + p[12];\r
 589     pp[4] = p[4] + p[11];\r
 590     pp[5] = p[5] + p[10];\r
 591     pp[6] = p[6] + p[9];\r
 592     pp[7] = p[7] + p[8];\r
 593     pp[8] = (p[0] - p[15]) * cos1_32;\r
 594     pp[9] = (p[1] - p[14]) * cos3_32;\r
 595     pp[10] = (p[2] - p[13]) * cos5_32;\r
 596     pp[11] = (p[3] - p[12]) * cos7_32;\r
 597     pp[12] = (p[4] - p[11]) * cos9_32;\r
 598     pp[13] = (p[5] - p[10]) * cos11_32;\r
 599     pp[14] = (p[6] - p[9]) * cos13_32;\r
 600     pp[15] = (p[7] - p[8]) * cos15_32;\r
 601 \r
 602     p[0] = pp[0] + pp[7];\r
 603     p[1] = pp[1] + pp[6];\r
 604     p[2] = pp[2] + pp[5];\r
 605     p[3] = pp[3] + pp[4];\r
 606     p[4] = (pp[0] - pp[7]) * cos1_16;\r
 607     p[5] = (pp[1] - pp[6]) * cos3_16;\r
 608     p[6] = (pp[2] - pp[5]) * cos5_16;\r
 609     p[7] = (pp[3] - pp[4]) * cos7_16;\r
 610     p[8] = pp[8] + pp[15];\r
 611     p[9] = pp[9] + pp[14];\r
 612     p[10] = pp[10] + pp[13];\r
 613     p[11] = pp[11] + pp[12];\r
 614     p[12] = (pp[8] - pp[15]) * cos1_16;\r
 615     p[13] = (pp[9] - pp[14]) * cos3_16;\r
 616     p[14] = (pp[10] - pp[13]) * cos5_16;\r
 617     p[15] = (pp[11] - pp[12]) * cos7_16;\r
 618 \r
 619     pp[0] = p[0] + p[3];\r
 620     pp[1] = p[1] + p[2];\r
 621     pp[2] = (p[0] - p[3]) * cos1_8;\r
 622     pp[3] = (p[1] - p[2]) * cos3_8;\r
 623     pp[4] = p[4] + p[7];\r
 624     pp[5] = p[5] + p[6];\r
 625     pp[6] = (p[4] - p[7]) * cos1_8;\r
 626     pp[7] = (p[5] - p[6]) * cos3_8;\r
 627     pp[8] = p[8] + p[11];\r
 628     pp[9] = p[9] + p[10];\r
 629     pp[10] = (p[8] - p[11]) * cos1_8;\r
 630     pp[11] = (p[9] - p[10]) * cos3_8;\r
 631     pp[12] = p[12] + p[15];\r
 632     pp[13] = p[13] + p[14];\r
 633     pp[14] = (p[12] - p[15]) * cos1_8;\r
 634     pp[15] = (p[13] - p[14]) * cos3_8;\r
 635 \r
 636     p[0] = pp[0] + pp[1];\r
 637     p[1] = (pp[0] - pp[1]) * cos1_4;\r
 638     p[2] = pp[2] + pp[3];\r
 639     p[3] = (pp[2] - pp[3]) * cos1_4;\r
 640     p[4] = pp[4] + pp[5];\r
 641     p[5] = (pp[4] - pp[5]) * cos1_4;\r
 642     p[6] = pp[6] + pp[7];\r
 643     p[7] = (pp[6] - pp[7]) * cos1_4;\r
 644     p[8] = pp[8] + pp[9];\r
 645     p[9] = (pp[8] - pp[9]) * cos1_4;\r
 646     p[10] = pp[10] + pp[11];\r
 647     p[11] = (pp[10] - pp[11]) * cos1_4;\r
 648     p[12] = pp[12] + pp[13];\r
 649     p[13] = (pp[12] - pp[13]) * cos1_4;\r
 650     p[14] = pp[14] + pp[15];\r
 651     p[15] = (pp[14] - pp[15]) * cos1_4;\r
 652 \r
 653     // this is pretty insane coding\r
 654     float tmp1;\r
 655     new_v[36 - 17] = -(new_v[4] = (new_v[12] = p[7]) + p[5]) - p[6];\r
 656     new_v[44 - 17] = -p[6] - p[7] - p[4];\r
 657     new_v[6] = (new_v[10] = (new_v[14] = p[15]) + p[11]) + p[13];\r
 658     new_v[34 - 17] = -(new_v[2] = p[15] + p[13] + p[9]) - p[14];\r
 659     new_v[38 - 17] = (tmp1 = -p[14] - p[15] - p[10] - p[11]) - p[13];\r
 660     new_v[46 - 17] = -p[14] - p[15] - p[12] - p[8];\r
 661     new_v[42 - 17] = tmp1 - p[12];\r
 662     new_v[48 - 17] = -p[0];\r
 663     new_v[0] = p[1];\r
 664     new_v[40 - 17] = -(new_v[8] = p[3]) - p[2];\r
 665 \r
 666     p[0] = (x1[0] - x1[31]) * cos1_64;\r
 667     p[1] = (x1[1] - x1[30]) * cos3_64;\r
 668     p[2] = (x1[2] - x1[29]) * cos5_64;\r
 669     p[3] = (x1[3] - x1[28]) * cos7_64;\r
 670     p[4] = (x1[4] - x1[27]) * cos9_64;\r
 671     p[5] = (x1[5] - x1[26]) * cos11_64;\r
 672     p[6] = (x1[6] - x1[25]) * cos13_64;\r
 673     p[7] = (x1[7] - x1[24]) * cos15_64;\r
 674     p[8] = (x1[8] - x1[23]) * cos17_64;\r
 675     p[9] = (x1[9] - x1[22]) * cos19_64;\r
 676     p[10] = (x1[10] - x1[21]) * cos21_64;\r
 677     p[11] = (x1[11] - x1[20]) * cos23_64;\r
 678     p[12] = (x1[12] - x1[19]) * cos25_64;\r
 679     p[13] = (x1[13] - x1[18]) * cos27_64;\r
 680     p[14] = (x1[14] - x1[17]) * cos29_64;\r
 681     p[15] = (x1[15] - x1[16]) * cos31_64;\r
 682 \r
 683     pp[0] = p[0] + p[15];\r
 684     pp[1] = p[1] + p[14];\r
 685     pp[2] = p[2] + p[13];\r
 686     pp[3] = p[3] + p[12];\r
 687     pp[4] = p[4] + p[11];\r
 688     pp[5] = p[5] + p[10];\r
 689     pp[6] = p[6] + p[9];\r
 690     pp[7] = p[7] + p[8];\r
 691     pp[8] = (p[0] - p[15]) * cos1_32;\r
 692     pp[9] = (p[1] - p[14]) * cos3_32;\r
 693     pp[10] = (p[2] - p[13]) * cos5_32;\r
 694     pp[11] = (p[3] - p[12]) * cos7_32;\r
 695     pp[12] = (p[4] - p[11]) * cos9_32;\r
 696     pp[13] = (p[5] - p[10]) * cos11_32;\r
 697     pp[14] = (p[6] - p[9]) * cos13_32;\r
 698     pp[15] = (p[7] - p[8]) * cos15_32;\r
 699 \r
 700     p[0] = pp[0] + pp[7];\r
 701     p[1] = pp[1] + pp[6];\r
 702     p[2] = pp[2] + pp[5];\r
 703     p[3] = pp[3] + pp[4];\r
 704     p[4] = (pp[0] - pp[7]) * cos1_16;\r
 705     p[5] = (pp[1] - pp[6]) * cos3_16;\r
 706     p[6] = (pp[2] - pp[5]) * cos5_16;\r
 707     p[7] = (pp[3] - pp[4]) * cos7_16;\r
 708     p[8] = pp[8] + pp[15];\r
 709     p[9] = pp[9] + pp[14];\r
 710     p[10] = pp[10] + pp[13];\r
 711     p[11] = pp[11] + pp[12];\r
 712     p[12] = (pp[8] - pp[15]) * cos1_16;\r
 713     p[13] = (pp[9] - pp[14]) * cos3_16;\r
 714     p[14] = (pp[10] - pp[13]) * cos5_16;\r
 715     p[15] = (pp[11] - pp[12]) * cos7_16;\r
 716 \r
 717     pp[0] = p[0] + p[3];\r
 718     pp[1] = p[1] + p[2];\r
 719     pp[2] = (p[0] - p[3]) * cos1_8;\r
 720     pp[3] = (p[1] - p[2]) * cos3_8;\r
 721     pp[4] = p[4] + p[7];\r
 722     pp[5] = p[5] + p[6];\r
 723     pp[6] = (p[4] - p[7]) * cos1_8;\r
 724     pp[7] = (p[5] - p[6]) * cos3_8;\r
 725     pp[8] = p[8] + p[11];\r
 726     pp[9] = p[9] + p[10];\r
 727     pp[10] = (p[8] - p[11]) * cos1_8;\r
 728     pp[11] = (p[9] - p[10]) * cos3_8;\r
 729     pp[12] = p[12] + p[15];\r
 730     pp[13] = p[13] + p[14];\r
 731     pp[14] = (p[12] - p[15]) * cos1_8;\r
 732     pp[15] = (p[13] - p[14]) * cos3_8;\r
 733 \r
 734     p[0] = pp[0] + pp[1];\r
 735     p[1] = (pp[0] - pp[1]) * cos1_4;\r
 736     p[2] = pp[2] + pp[3];\r
 737     p[3] = (pp[2] - pp[3]) * cos1_4;\r
 738     p[4] = pp[4] + pp[5];\r
 739     p[5] = (pp[4] - pp[5]) * cos1_4;\r
 740     p[6] = pp[6] + pp[7];\r
 741     p[7] = (pp[6] - pp[7]) * cos1_4;\r
 742     p[8] = pp[8] + pp[9];\r
 743     p[9] = (pp[8] - pp[9]) * cos1_4;\r
 744     p[10] = pp[10] + pp[11];\r
 745     p[11] = (pp[10] - pp[11]) * cos1_4;\r
 746     p[12] = pp[12] + pp[13];\r
 747     p[13] = (pp[12] - pp[13]) * cos1_4;\r
 748     p[14] = pp[14] + pp[15];\r
 749     p[15] = (pp[14] - pp[15]) * cos1_4;\r
 750 \r
 751     // manually doing something that a compiler should handle sucks\r
 752     // coding like this is hard to read\r
 753     float tmp2;\r
 754     new_v[5] = (new_v[11] = (new_v[13] = (new_v[15] = p[15]) + p[7]) + p[11]) + p[5] + p[13];\r
 755     new_v[7] = (new_v[9] = p[15] + p[11] + p[3]) + p[13];\r
 756     new_v[33 - 17] = -(new_v[1] = (tmp1 = p[13] + p[15] + p[9]) + p[1]) - p[14];\r
 757     new_v[35 - 17] = -(new_v[3] = tmp1 + p[5] + p[7]) - p[6] - p[14];\r
 758 \r
 759     new_v[39 - 17] = (tmp1 = -p[10] - p[11] - p[14] - p[15]) - p[13] - p[2] - p[3];\r
 760     new_v[37 - 17] = tmp1 - p[13] - p[5] - p[6] - p[7];\r
 761     new_v[41 - 17] = tmp1 - p[12] - p[2] - p[3];\r
 762     new_v[43 - 17] = tmp1 - p[12] - (tmp2 = p[4] + p[6] + p[7]);\r
 763     new_v[47 - 17] = (tmp1 = -p[8] - p[12] - p[14] - p[15]) - p[0];\r
 764     new_v[45 - 17] = tmp1 - tmp2;\r
 765 \r
 766     // insert V[0-15] (== new_v[0-15]) into actual v:\r
 767     x1 = new_v;\r
 768     // float[] x2 = actual_v + actual_write_pos;\r
 769     float[] dest = actual_v;\r
 770 \r
 771     dest[0 + actual_write_pos] = x1[0];\r
 772     dest[16 + actual_write_pos] = x1[1];\r
 773     dest[32 + actual_write_pos] = x1[2];\r
 774     dest[48 + actual_write_pos] = x1[3];\r
 775     dest[64 + actual_write_pos] = x1[4];\r
 776     dest[80 + actual_write_pos] = x1[5];\r
 777     dest[96 + actual_write_pos] = x1[6];\r
 778     dest[112 + actual_write_pos] = x1[7];\r
 779     dest[128 + actual_write_pos] = x1[8];\r
 780     dest[144 + actual_write_pos] = x1[9];\r
 781     dest[160 + actual_write_pos] = x1[10];\r
 782     dest[176 + actual_write_pos] = x1[11];\r
 783     dest[192 + actual_write_pos] = x1[12];\r
 784     dest[208 + actual_write_pos] = x1[13];\r
 785     dest[224 + actual_write_pos] = x1[14];\r
 786     dest[240 + actual_write_pos] = x1[15];\r
 787 \r
 788     // V[16] is always 0.0:\r
 789     dest[256 + actual_write_pos] = 0.0f;\r
 790 \r
 791     // insert V[17-31] (== -new_v[15-1]) into actual v:\r
 792     dest[272 + actual_write_pos] = -x1[15];\r
 793     dest[288 + actual_write_pos] = -x1[14];\r
 794     dest[304 + actual_write_pos] = -x1[13];\r
 795     dest[320 + actual_write_pos] = -x1[12];\r
 796     dest[336 + actual_write_pos] = -x1[11];\r
 797     dest[352 + actual_write_pos] = -x1[10];\r
 798     dest[368 + actual_write_pos] = -x1[9];\r
 799     dest[384 + actual_write_pos] = -x1[8];\r
 800     dest[400 + actual_write_pos] = -x1[7];\r
 801     dest[416 + actual_write_pos] = -x1[6];\r
 802     dest[432 + actual_write_pos] = -x1[5];\r
 803     dest[448 + actual_write_pos] = -x1[4];\r
 804     dest[464 + actual_write_pos] = -x1[3];\r
 805     dest[480 + actual_write_pos] = -x1[2];\r
 806     dest[496 + actual_write_pos] = -x1[1];\r
 807 \r
 808     // insert V[32] (== -new_v[0]) into other v:\r
 809 \r
 810   }\r
 811 \r
 812   /**\r
 813    * Compute PCM Samples.\r
 814    */\r
 815 \r
 816   private float[] _tmpOut = new float[32];\r
 817 \r
 818   private void compute_pcm_samples0(Obuffer buffer) {\r
 819     final float[] vp = actual_v;\r
 820     // int inc = v_inc;\r
 821     final float[] tmpOut = _tmpOut;\r
 822     int dvp = 0;\r
 823 \r
 824     // fat chance of having this loop unroll\r
 825     for (int i = 0; i < 32; i++) {\r
 826       float pcm_sample;\r
 827       final float[] dp = d16[i];\r
 828       pcm_sample =\r
 829           (float) (((vp[0 + dvp] * dp[0]) + (vp[15 + dvp] * dp[1]) + (vp[14 + dvp] * dp[2])\r
 830               + (vp[13 + dvp] * dp[3]) + (vp[12 + dvp] * dp[4]) + (vp[11 + dvp] * dp[5])\r
 831               + (vp[10 + dvp] * dp[6]) + (vp[9 + dvp] * dp[7]) + (vp[8 + dvp] * dp[8])\r
 832               + (vp[7 + dvp] * dp[9]) + (vp[6 + dvp] * dp[10]) + (vp[5 + dvp] * dp[11])\r
 833               + (vp[4 + dvp] * dp[12]) + (vp[3 + dvp] * dp[13]) + (vp[2 + dvp] * dp[14]) + (vp[1 + dvp] * dp[15])) * scalefactor);\r
 834 \r
 835       tmpOut[i] = pcm_sample;\r
 836 \r
 837       dvp += 16;\r
 838     } // for\r
 839   }\r
 840 \r
 841   private void compute_pcm_samples1(Obuffer buffer) {\r
 842     final float[] vp = actual_v;\r
 843     // int inc = v_inc;\r
 844     final float[] tmpOut = _tmpOut;\r
 845     int dvp = 0;\r
 846 \r
 847     // fat chance of having this loop unroll\r
 848     for (int i = 0; i < 32; i++) {\r
 849       final float[] dp = d16[i];\r
 850       float pcm_sample;\r
 851 \r
 852       pcm_sample =\r
 853           (float) (((vp[1 + dvp] * dp[0]) + (vp[0 + dvp] * dp[1]) + (vp[15 + dvp] * dp[2])\r
 854               + (vp[14 + dvp] * dp[3]) + (vp[13 + dvp] * dp[4]) + (vp[12 + dvp] * dp[5])\r
 855               + (vp[11 + dvp] * dp[6]) + (vp[10 + dvp] * dp[7]) + (vp[9 + dvp] * dp[8])\r
 856               + (vp[8 + dvp] * dp[9]) + (vp[7 + dvp] * dp[10]) + (vp[6 + dvp] * dp[11])\r
 857               + (vp[5 + dvp] * dp[12]) + (vp[4 + dvp] * dp[13]) + (vp[3 + dvp] * dp[14]) + (vp[2 + dvp] * dp[15])) * scalefactor);\r
 858 \r
 859       tmpOut[i] = pcm_sample;\r
 860 \r
 861       dvp += 16;\r
 862     } // for\r
 863   }\r
 864 \r
 865   private void compute_pcm_samples2(Obuffer buffer) {\r
 866     final float[] vp = actual_v;\r
 867 \r
 868     // int inc = v_inc;\r
 869     final float[] tmpOut = _tmpOut;\r
 870     int dvp = 0;\r
 871 \r
 872     // fat chance of having this loop unroll\r
 873     for (int i = 0; i < 32; i++) {\r
 874       final float[] dp = d16[i];\r
 875       float pcm_sample;\r
 876 \r
 877       pcm_sample =\r
 878           (float) (((vp[2 + dvp] * dp[0]) + (vp[1 + dvp] * dp[1]) + (vp[0 + dvp] * dp[2])\r
 879               + (vp[15 + dvp] * dp[3]) + (vp[14 + dvp] * dp[4]) + (vp[13 + dvp] * dp[5])\r
 880               + (vp[12 + dvp] * dp[6]) + (vp[11 + dvp] * dp[7]) + (vp[10 + dvp] * dp[8])\r
 881               + (vp[9 + dvp] * dp[9]) + (vp[8 + dvp] * dp[10]) + (vp[7 + dvp] * dp[11])\r
 882               + (vp[6 + dvp] * dp[12]) + (vp[5 + dvp] * dp[13]) + (vp[4 + dvp] * dp[14]) + (vp[3 + dvp] * dp[15])) * scalefactor);\r
 883 \r
 884       tmpOut[i] = pcm_sample;\r
 885 \r
 886       dvp += 16;\r
 887     } // for\r
 888   }\r
 889 \r
 890   private void compute_pcm_samples3(Obuffer buffer) {\r
 891     final float[] vp = actual_v;\r
 892 \r
 893     int idx = 0;\r
 894     // int inc = v_inc;\r
 895     final float[] tmpOut = _tmpOut;\r
 896     int dvp = 0;\r
 897 \r
 898     // fat chance of having this loop unroll\r
 899     for (int i = 0; i < 32; i++) {\r
 900       final float[] dp = d16[i];\r
 901       float pcm_sample;\r
 902 \r
 903       pcm_sample =\r
 904           (float) (((vp[3 + dvp] * dp[0]) + (vp[2 + dvp] * dp[1]) + (vp[1 + dvp] * dp[2])\r
 905               + (vp[0 + dvp] * dp[3]) + (vp[15 + dvp] * dp[4]) + (vp[14 + dvp] * dp[5])\r
 906               + (vp[13 + dvp] * dp[6]) + (vp[12 + dvp] * dp[7]) + (vp[11 + dvp] * dp[8])\r
 907               + (vp[10 + dvp] * dp[9]) + (vp[9 + dvp] * dp[10]) + (vp[8 + dvp] * dp[11])\r
 908               + (vp[7 + dvp] * dp[12]) + (vp[6 + dvp] * dp[13]) + (vp[5 + dvp] * dp[14]) + (vp[4 + dvp] * dp[15])) * scalefactor);\r
 909 \r
 910       tmpOut[i] = pcm_sample;\r
 911 \r
 912       dvp += 16;\r
 913     } // for\r
 914   }\r
 915 \r
 916   private void compute_pcm_samples4(Obuffer buffer) {\r
 917     final float[] vp = actual_v;\r
 918 \r
 919     // int inc = v_inc;\r
 920     final float[] tmpOut = _tmpOut;\r
 921     int dvp = 0;\r
 922 \r
 923     // fat chance of having this loop unroll\r
 924     for (int i = 0; i < 32; i++) {\r
 925       final float[] dp = d16[i];\r
 926       float pcm_sample;\r
 927 \r
 928       pcm_sample =\r
 929           (float) (((vp[4 + dvp] * dp[0]) + (vp[3 + dvp] * dp[1]) + (vp[2 + dvp] * dp[2])\r
 930               + (vp[1 + dvp] * dp[3]) + (vp[0 + dvp] * dp[4]) + (vp[15 + dvp] * dp[5])\r
 931               + (vp[14 + dvp] * dp[6]) + (vp[13 + dvp] * dp[7]) + (vp[12 + dvp] * dp[8])\r
 932               + (vp[11 + dvp] * dp[9]) + (vp[10 + dvp] * dp[10]) + (vp[9 + dvp] * dp[11])\r
 933               + (vp[8 + dvp] * dp[12]) + (vp[7 + dvp] * dp[13]) + (vp[6 + dvp] * dp[14]) + (vp[5 + dvp] * dp[15])) * scalefactor);\r
 934 \r
 935       tmpOut[i] = pcm_sample;\r
 936 \r
 937       dvp += 16;\r
 938     } // for\r
 939   }\r
 940 \r
 941   private void compute_pcm_samples5(Obuffer buffer) {\r
 942     final float[] vp = actual_v;\r
 943 \r
 944     // int inc = v_inc;\r
 945     final float[] tmpOut = _tmpOut;\r
 946     int dvp = 0;\r
 947 \r
 948     // fat chance of having this loop unroll\r
 949     for (int i = 0; i < 32; i++) {\r
 950       final float[] dp = d16[i];\r
 951       float pcm_sample;\r
 952 \r
 953       pcm_sample =\r
 954           (float) (((vp[5 + dvp] * dp[0]) + (vp[4 + dvp] * dp[1]) + (vp[3 + dvp] * dp[2])\r
 955               + (vp[2 + dvp] * dp[3]) + (vp[1 + dvp] * dp[4]) + (vp[0 + dvp] * dp[5])\r
 956               + (vp[15 + dvp] * dp[6]) + (vp[14 + dvp] * dp[7]) + (vp[13 + dvp] * dp[8])\r
 957               + (vp[12 + dvp] * dp[9]) + (vp[11 + dvp] * dp[10]) + (vp[10 + dvp] * dp[11])\r
 958               + (vp[9 + dvp] * dp[12]) + (vp[8 + dvp] * dp[13]) + (vp[7 + dvp] * dp[14]) + (vp[6 + dvp] * dp[15])) * scalefactor);\r
 959 \r
 960       tmpOut[i] = pcm_sample;\r
 961 \r
 962       dvp += 16;\r
 963     } // for\r
 964   }\r
 965 \r
 966   private void compute_pcm_samples6(Obuffer buffer) {\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) (((vp[6 + dvp] * dp[0]) + (vp[5 + dvp] * dp[1]) + (vp[4 + dvp] * dp[2])\r
 979               + (vp[3 + dvp] * dp[3]) + (vp[2 + dvp] * dp[4]) + (vp[1 + dvp] * dp[5])\r
 980               + (vp[0 + dvp] * dp[6]) + (vp[15 + dvp] * dp[7]) + (vp[14 + dvp] * dp[8])\r
 981               + (vp[13 + dvp] * dp[9]) + (vp[12 + dvp] * dp[10]) + (vp[11 + dvp] * dp[11])\r
 982               + (vp[10 + dvp] * dp[12]) + (vp[9 + dvp] * dp[13]) + (vp[8 + dvp] * dp[14]) + (vp[7 + dvp] * dp[15])) * scalefactor);\r
 983 \r
 984       tmpOut[i] = pcm_sample;\r
 985 \r
 986       dvp += 16;\r
 987     } // for\r
 988   }\r
 989 \r
 990   private void compute_pcm_samples7(Obuffer buffer) {\r
 991     final float[] vp = actual_v;\r
 992 \r
 993     // int inc = v_inc;\r
 994     final float[] tmpOut = _tmpOut;\r
 995     int dvp = 0;\r
 996 \r
 997     // fat chance of having this loop unroll\r
 998     for (int i = 0; i < 32; i++) {\r
 999       final float[] dp = d16[i];\r
1000       float pcm_sample;\r
1001 \r
1002       pcm_sample =\r
1003           (float) (((vp[7 + dvp] * dp[0]) + (vp[6 + dvp] * dp[1]) + (vp[5 + dvp] * dp[2])\r
1004               + (vp[4 + dvp] * dp[3]) + (vp[3 + dvp] * dp[4]) + (vp[2 + dvp] * dp[5])\r
1005               + (vp[1 + dvp] * dp[6]) + (vp[0 + dvp] * dp[7]) + (vp[15 + dvp] * dp[8])\r
1006               + (vp[14 + dvp] * dp[9]) + (vp[13 + dvp] * dp[10]) + (vp[12 + dvp] * dp[11])\r
1007               + (vp[11 + dvp] * dp[12]) + (vp[10 + dvp] * dp[13]) + (vp[9 + dvp] * dp[14]) + (vp[8 + dvp] * dp[15])) * scalefactor);\r
1008 \r
1009       tmpOut[i] = pcm_sample;\r
1010 \r
1011       dvp += 16;\r
1012     } // for\r
1013   }\r
1014 \r
1015   private void compute_pcm_samples8(Obuffer buffer) {\r
1016     final float[] vp = actual_v;\r
1017 \r
1018     // int inc = v_inc;\r
1019     final float[] tmpOut = _tmpOut;\r
1020     int dvp = 0;\r
1021 \r
1022     // fat chance of having this loop unroll\r
1023     for (int i = 0; i < 32; i++) {\r
1024       final float[] dp = d16[i];\r
1025       float pcm_sample;\r
1026 \r
1027       pcm_sample =\r
1028           (float) (((vp[8 + dvp] * dp[0]) + (vp[7 + dvp] * dp[1]) + (vp[6 + dvp] * dp[2])\r
1029               + (vp[5 + dvp] * dp[3]) + (vp[4 + dvp] * dp[4]) + (vp[3 + dvp] * dp[5])\r
1030               + (vp[2 + dvp] * dp[6]) + (vp[1 + dvp] * dp[7]) + (vp[0 + dvp] * dp[8])\r
1031               + (vp[15 + dvp] * dp[9]) + (vp[14 + dvp] * dp[10]) + (vp[13 + dvp] * dp[11])\r
1032               + (vp[12 + dvp] * dp[12]) + (vp[11 + dvp] * dp[13]) + (vp[10 + dvp] * dp[14]) + (vp[9 + dvp] * dp[15])) * scalefactor);\r
1033 \r
1034       tmpOut[i] = pcm_sample;\r
1035 \r
1036       dvp += 16;\r
1037     } // for\r
1038   }\r
1039 \r
1040   private void compute_pcm_samples9(Obuffer buffer) {\r
1041     final float[] vp = actual_v;\r
1042 \r
1043     // int inc = v_inc;\r
1044     final float[] tmpOut = _tmpOut;\r
1045     int dvp = 0;\r
1046 \r
1047     // fat chance of having this loop unroll\r
1048     for (int i = 0; i < 32; i++) {\r
1049       final float[] dp = d16[i];\r
1050       float pcm_sample;\r
1051 \r
1052       pcm_sample =\r
1053           (float) (((vp[9 + dvp] * dp[0]) + (vp[8 + dvp] * dp[1]) + (vp[7 + dvp] * dp[2])\r
1054               + (vp[6 + dvp] * dp[3]) + (vp[5 + dvp] * dp[4]) + (vp[4 + dvp] * dp[5])\r
1055               + (vp[3 + dvp] * dp[6]) + (vp[2 + dvp] * dp[7]) + (vp[1 + dvp] * dp[8])\r
1056               + (vp[0 + dvp] * dp[9]) + (vp[15 + dvp] * dp[10]) + (vp[14 + dvp] * dp[11])\r
1057               + (vp[13 + dvp] * dp[12]) + (vp[12 + dvp] * dp[13]) + (vp[11 + dvp] * dp[14]) + (vp[10 + dvp] * dp[15])) * scalefactor);\r
1058 \r
1059       tmpOut[i] = pcm_sample;\r
1060 \r
1061       dvp += 16;\r
1062     } // for\r
1063   }\r
1064 \r
1065   private void compute_pcm_samples10(Obuffer buffer) {\r
1066     final float[] vp = actual_v;\r
1067     // int inc = v_inc;\r
1068     final float[] tmpOut = _tmpOut;\r
1069     int dvp = 0;\r
1070 \r
1071     // fat chance of having this loop unroll\r
1072     for (int i = 0; i < 32; i++) {\r
1073       final float[] dp = d16[i];\r
1074       float pcm_sample;\r
1075 \r
1076       pcm_sample =\r
1077           (float) (((vp[10 + dvp] * dp[0]) + (vp[9 + dvp] * dp[1]) + (vp[8 + dvp] * dp[2])\r
1078               + (vp[7 + dvp] * dp[3]) + (vp[6 + dvp] * dp[4]) + (vp[5 + dvp] * dp[5])\r
1079               + (vp[4 + dvp] * dp[6]) + (vp[3 + dvp] * dp[7]) + (vp[2 + dvp] * dp[8])\r
1080               + (vp[1 + dvp] * dp[9]) + (vp[0 + dvp] * dp[10]) + (vp[15 + dvp] * dp[11])\r
1081               + (vp[14 + dvp] * dp[12]) + (vp[13 + dvp] * dp[13]) + (vp[12 + dvp] * dp[14]) + (vp[11 + dvp] * dp[15])) * scalefactor);\r
1082 \r
1083       tmpOut[i] = pcm_sample;\r
1084 \r
1085       dvp += 16;\r
1086     } // for\r
1087   }\r
1088 \r
1089   private void compute_pcm_samples11(Obuffer buffer) {\r
1090     final float[] vp = actual_v;\r
1091 \r
1092     // int inc = v_inc;\r
1093     final float[] tmpOut = _tmpOut;\r
1094     int dvp = 0;\r
1095 \r
1096     // fat chance of having this loop unroll\r
1097     for (int i = 0; i < 32; i++) {\r
1098       final float[] dp = d16[i];\r
1099       float pcm_sample;\r
1100 \r
1101       pcm_sample =\r
1102           (float) (((vp[11 + dvp] * dp[0]) + (vp[10 + dvp] * dp[1]) + (vp[9 + dvp] * dp[2])\r
1103               + (vp[8 + dvp] * dp[3]) + (vp[7 + dvp] * dp[4]) + (vp[6 + dvp] * dp[5])\r
1104               + (vp[5 + dvp] * dp[6]) + (vp[4 + dvp] * dp[7]) + (vp[3 + dvp] * dp[8])\r
1105               + (vp[2 + dvp] * dp[9]) + (vp[1 + dvp] * dp[10]) + (vp[0 + dvp] * dp[11])\r
1106               + (vp[15 + dvp] * dp[12]) + (vp[14 + dvp] * dp[13]) + (vp[13 + dvp] * dp[14]) + (vp[12 + dvp] * dp[15])) * scalefactor);\r
1107 \r
1108       tmpOut[i] = pcm_sample;\r
1109 \r
1110       dvp += 16;\r
1111     } // for\r
1112   }\r
1113 \r
1114   private void compute_pcm_samples12(Obuffer buffer) {\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 =\r
1126           (float) (((vp[12 + dvp] * dp[0]) + (vp[11 + dvp] * dp[1]) + (vp[10 + dvp] * dp[2])\r
1127               + (vp[9 + dvp] * dp[3]) + (vp[8 + dvp] * dp[4]) + (vp[7 + dvp] * dp[5])\r
1128               + (vp[6 + dvp] * dp[6]) + (vp[5 + dvp] * dp[7]) + (vp[4 + dvp] * dp[8])\r
1129               + (vp[3 + dvp] * dp[9]) + (vp[2 + dvp] * dp[10]) + (vp[1 + dvp] * dp[11])\r
1130               + (vp[0 + dvp] * dp[12]) + (vp[15 + dvp] * dp[13]) + (vp[14 + dvp] * dp[14]) + (vp[13 + dvp] * dp[15])) * scalefactor);\r
1131 \r
1132       tmpOut[i] = pcm_sample;\r
1133 \r
1134       dvp += 16;\r
1135     } // for\r
1136   }\r
1137 \r
1138   private void compute_pcm_samples13(Obuffer buffer) {\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) (((vp[13 + dvp] * dp[0]) + (vp[12 + dvp] * dp[1]) + (vp[11 + dvp] * dp[2])\r
1152               + (vp[10 + dvp] * dp[3]) + (vp[9 + dvp] * dp[4]) + (vp[8 + dvp] * dp[5])\r
1153               + (vp[7 + dvp] * dp[6]) + (vp[6 + dvp] * dp[7]) + (vp[5 + dvp] * dp[8])\r
1154               + (vp[4 + dvp] * dp[9]) + (vp[3 + dvp] * dp[10]) + (vp[2 + dvp] * dp[11])\r
1155               + (vp[1 + dvp] * dp[12]) + (vp[0 + dvp] * dp[13]) + (vp[15 + dvp] * dp[14]) + (vp[14 + dvp] * dp[15])) * scalefactor);\r
1156 \r
1157       tmpOut[i] = pcm_sample;\r
1158 \r
1159       dvp += 16;\r
1160     } // for\r
1161   }\r
1162 \r
1163   private void compute_pcm_samples14(Obuffer buffer) {\r
1164     final float[] vp = actual_v;\r
1165 \r
1166     // int inc = v_inc;\r
1167     final float[] tmpOut = _tmpOut;\r
1168     int dvp = 0;\r
1169 \r
1170     // fat chance of having this loop unroll\r
1171     for (int i = 0; i < 32; i++) {\r
1172       final float[] dp = d16[i];\r
1173       float pcm_sample;\r
1174 \r
1175       pcm_sample =\r
1176           (float) (((vp[14 + dvp] * dp[0]) + (vp[13 + dvp] * dp[1]) + (vp[12 + dvp] * dp[2])\r
1177               + (vp[11 + dvp] * dp[3]) + (vp[10 + dvp] * dp[4]) + (vp[9 + dvp] * dp[5])\r
1178               + (vp[8 + dvp] * dp[6]) + (vp[7 + dvp] * dp[7]) + (vp[6 + dvp] * dp[8])\r
1179               + (vp[5 + dvp] * dp[9]) + (vp[4 + dvp] * dp[10]) + (vp[3 + dvp] * dp[11])\r
1180               + (vp[2 + dvp] * dp[12]) + (vp[1 + dvp] * dp[13]) + (vp[0 + dvp] * dp[14]) + (vp[15 + dvp] * dp[15])) * scalefactor);\r
1181 \r
1182       tmpOut[i] = pcm_sample;\r
1183 \r
1184       dvp += 16;\r
1185     } // for\r
1186   }\r
1187 \r
1188   private void compute_pcm_samples15(Obuffer buffer) {\r
1189     final float[] vp = actual_v;\r
1190 \r
1191     // int inc = v_inc;\r
1192     final float[] tmpOut = _tmpOut;\r
1193     int dvp = 0;\r
1194 \r
1195     // fat chance of having this loop unroll\r
1196     for (int i = 0; i < 32; i++) {\r
1197       float pcm_sample;\r
1198       final float dp[] = d16[i];\r
1199       pcm_sample =\r
1200           (float) (((vp[15 + dvp] * dp[0]) + (vp[14 + dvp] * dp[1]) + (vp[13 + dvp] * dp[2])\r
1201               + (vp[12 + dvp] * dp[3]) + (vp[11 + dvp] * dp[4]) + (vp[10 + dvp] * dp[5])\r
1202               + (vp[9 + dvp] * dp[6]) + (vp[8 + dvp] * dp[7]) + (vp[7 + dvp] * dp[8])\r
1203               + (vp[6 + dvp] * dp[9]) + (vp[5 + dvp] * dp[10]) + (vp[4 + dvp] * dp[11])\r
1204               + (vp[3 + dvp] * dp[12]) + (vp[2 + dvp] * dp[13]) + (vp[1 + dvp] * dp[14]) + (vp[0 + dvp] * dp[15])) * scalefactor);\r
1205 \r
1206       tmpOut[i] = pcm_sample;\r
1207       dvp += 16;\r
1208     } // for\r
1209   }\r
1210 \r
1211   private void compute_pcm_samples(Obuffer buffer) {\r
1212 \r
1213     switch (actual_write_pos) {\r
1214     case 0:\r
1215       compute_pcm_samples0(buffer);\r
1216       break;\r
1217     case 1:\r
1218       compute_pcm_samples1(buffer);\r
1219       break;\r
1220     case 2:\r
1221       compute_pcm_samples2(buffer);\r
1222       break;\r
1223     case 3:\r
1224       compute_pcm_samples3(buffer);\r
1225       break;\r
1226     case 4:\r
1227       compute_pcm_samples4(buffer);\r
1228       break;\r
1229     case 5:\r
1230       compute_pcm_samples5(buffer);\r
1231       break;\r
1232     case 6:\r
1233       compute_pcm_samples6(buffer);\r
1234       break;\r
1235     case 7:\r
1236       compute_pcm_samples7(buffer);\r
1237       break;\r
1238     case 8:\r
1239       compute_pcm_samples8(buffer);\r
1240       break;\r
1241     case 9:\r
1242       compute_pcm_samples9(buffer);\r
1243       break;\r
1244     case 10:\r
1245       compute_pcm_samples10(buffer);\r
1246       break;\r
1247     case 11:\r
1248       compute_pcm_samples11(buffer);\r
1249       break;\r
1250     case 12:\r
1251       compute_pcm_samples12(buffer);\r
1252       break;\r
1253     case 13:\r
1254       compute_pcm_samples13(buffer);\r
1255       break;\r
1256     case 14:\r
1257       compute_pcm_samples14(buffer);\r
1258       break;\r
1259     case 15:\r
1260       compute_pcm_samples15(buffer);\r
1261       break;\r
1262     }\r
1263 \r
1264     if (buffer != null) {\r
1265       buffer.appendSamples(channel, _tmpOut);\r
1266     }\r
1267 \r
1268     /*\r
1269      * // MDM: I was considering putting in quality control for // low-spec\r
1270      * CPUs, but the performance gain (about 10-15%) // did not justify the\r
1271      * considerable drop in audio quality. switch (inc) { case 16:\r
1272      * buffer.appendSamples(channel, tmpOut); break; case 32: for (int i=0;\r
1273      * i<16; i++) { buffer.append(channel, (short)tmpOut[i]);\r
1274      * buffer.append(channel, (short)tmpOut[i]); } break; case 64: for (int i=0;\r
1275      * i<8; i++) { buffer.append(channel, (short)tmpOut[i]);\r
1276      * buffer.append(channel, (short)tmpOut[i]); buffer.append(channel,\r
1277      * (short)tmpOut[i]); buffer.append(channel, (short)tmpOut[i]); } break;\r
1278      * \r
1279      * }\r
1280      */\r
1281   }\r
1282 \r
1283   /**\r
1284    * Calculate 32 PCM samples and put the into the Obuffer-object.\r
1285    */\r
1286 \r
1287   public void calculate_pcm_samples(Obuffer buffer) {\r
1288     compute_new_v();\r
1289     compute_pcm_samples(buffer);\r
1290 \r
1291     actual_write_pos = (actual_write_pos + 1) & 0xf;\r
1292     actual_v = (actual_v == v1) ? v2 : v1;\r
1293 \r
1294     // initialize samples[]:\r
1295     // for (register float *floatp = samples + 32; floatp > samples; )\r
1296     // *--floatp = 0.0f;\r
1297 \r
1298     // MDM: this may not be necessary. The Layer III decoder always\r
1299     // outputs 32 subband samples, but I haven't checked layer I & II.\r
1300     for (int p = 0; p < 32; p++)\r
1301       samples[p] = 0.0f;\r
1302   }\r
1303 \r
1304   private static final double MY_PI = 3.14159265358979323846;\r
1305   private static final float cos1_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI / 64.0)));\r
1306   private static final float cos3_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 3.0 / 64.0)));\r
1307   private static final float cos5_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 5.0 / 64.0)));\r
1308   private static final float cos7_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 7.0 / 64.0)));\r
1309   private static final float cos9_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 9.0 / 64.0)));\r
1310   private static final float cos11_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 11.0 / 64.0)));\r
1311   private static final float cos13_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 13.0 / 64.0)));\r
1312   private static final float cos15_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 15.0 / 64.0)));\r
1313   private static final float cos17_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 17.0 / 64.0)));\r
1314   private static final float cos19_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 19.0 / 64.0)));\r
1315   private static final float cos21_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 21.0 / 64.0)));\r
1316   private static final float cos23_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 23.0 / 64.0)));\r
1317   private static final float cos25_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 25.0 / 64.0)));\r
1318   private static final float cos27_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 27.0 / 64.0)));\r
1319   private static final float cos29_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 29.0 / 64.0)));\r
1320   private static final float cos31_64 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 31.0 / 64.0)));\r
1321   private static final float cos1_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI / 32.0)));\r
1322   private static final float cos3_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 3.0 / 32.0)));\r
1323   private static final float cos5_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 5.0 / 32.0)));\r
1324   private static final float cos7_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 7.0 / 32.0)));\r
1325   private static final float cos9_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 9.0 / 32.0)));\r
1326   private static final float cos11_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 11.0 / 32.0)));\r
1327   private static final float cos13_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 13.0 / 32.0)));\r
1328   private static final float cos15_32 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 15.0 / 32.0)));\r
1329   private static final float cos1_16 = (float) (1.0 / (2.0 * Math.cos(MY_PI / 16.0)));\r
1330   private static final float cos3_16 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 3.0 / 16.0)));\r
1331   private static final float cos5_16 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 5.0 / 16.0)));\r
1332   private static final float cos7_16 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 7.0 / 16.0)));\r
1333   private static final float cos1_8 = (float) (1.0 / (2.0 * Math.cos(MY_PI / 8.0)));\r
1334   private static final float cos3_8 = (float) (1.0 / (2.0 * Math.cos(MY_PI * 3.0 / 8.0)));\r
1335   private static final float cos1_4 = (float) (1.0 / (2.0 * Math.cos(MY_PI / 4.0)));\r
1336 \r
1337   // Note: These values are not in the same order\r
1338   // as in Annex 3-B.3 of the ISO/IEC DIS 11172-3\r
1339   // private float d[] = {0.000000000, -4.000442505};\r
1340 \r
1341   /**\r
1342    * d[] split into subarrays of length 16. This provides for more faster access\r
1343    * by allowing a block of 16 to be addressed with constant offset.\r
1344    **/\r
1345   private static float d16[][] = null;\r
1346 \r
1347   /**\r
1348    * Converts a 1D array into a number of smaller arrays. This is used to\r
1349    * achieve offset + constant indexing into an array. Each sub-array represents\r
1350    * a block of values of the original array.\r
1351    * \r
1352    * @param array\r
1353    *          The array to split up into blocks.\r
1354    * @param blockSize\r
1355    *          The size of the blocks to split the array into. This must be an\r
1356    *          exact divisor of the length of the array, or some data will be\r
1357    *          lost from the main array.\r
1358    * \r
1359    * @return An array of arrays in which each element in the returned array will\r
1360    *         be of length <code>blockSize</code>.\r
1361    */\r
1362   static private float[][] splitArray(final float[] array, final int blockSize) {\r
1363     int size = array.length / blockSize;\r
1364     float[][] split = new float[size][];\r
1365     for (int i = 0; i < size; i++) {\r
1366       split[i] = subArray(array, i * blockSize, blockSize);\r
1367     }\r
1368     return split;\r
1369   }\r
1370 \r
1371   /**\r
1372    * Returns a subarray of an existing array.\r
1373    * \r
1374    * @param array\r
1375    *          The array to retrieve a subarra from.\r
1376    * @param offs\r
1377    *          The offset in the array that corresponds to the first index of the\r
1378    *          subarray.\r
1379    * @param len\r
1380    *          The number of indeces in the subarray.\r
1381    * @return The subarray, which may be of length 0.\r
1382    */\r
1383   static private float[] subArray(final float[] array, final int offs, int len) {\r
1384     if (offs + len > array.length) {\r
1385       len = array.length - offs;\r
1386     }\r
1387 \r
1388     if (len < 0)\r
1389       len = 0;\r
1390 \r
1391     float[] subarray = new float[len];\r
1392     for (int i = 0; i < len; i++) {\r
1393       subarray[i] = array[offs + i];\r
1394     }\r
1395 \r
1396     return subarray;\r
1397   }\r
1398 \r
1399   // The original data for d[]. This data is loaded from a file\r
1400   // to reduce the overall package size and to improve performance.\r
1401 \r
1402   static final float d[] = { 0.000000000f, -0.000442505f, 0.003250122f, -0.007003784f,\r
1403       0.031082153f, -0.078628540f, 0.100311279f, -0.572036743f, 1.144989014f, 0.572036743f,\r
1404       0.100311279f, 0.078628540f, 0.031082153f, 0.007003784f, 0.003250122f, 0.000442505f,\r
1405       -0.000015259f, -0.000473022f, 0.003326416f, -0.007919312f, 0.030517578f, -0.084182739f,\r
1406       0.090927124f, -0.600219727f, 1.144287109f, 0.543823242f, 0.108856201f, 0.073059082f,\r
1407       0.031478882f, 0.006118774f, 0.003173828f, 0.000396729f, -0.000015259f, -0.000534058f,\r
1408       0.003387451f, -0.008865356f, 0.029785156f, -0.089706421f, 0.080688477f, -0.628295898f,\r
1409       1.142211914f, 0.515609741f, 0.116577148f, 0.067520142f, 0.031738281f, 0.005294800f,\r
1410       0.003082275f, 0.000366211f, -0.000015259f, -0.000579834f, 0.003433228f, -0.009841919f,\r
1411       0.028884888f, -0.095169067f, 0.069595337f, -0.656219482f, 1.138763428f, 0.487472534f,\r
1412       0.123474121f, 0.061996460f, 0.031845093f, 0.004486084f, 0.002990723f, 0.000320435f,\r
1413       -0.000015259f, -0.000625610f, 0.003463745f, -0.010848999f, 0.027801514f, -0.100540161f,\r
1414       0.057617188f, -0.683914185f, 1.133926392f, 0.459472656f, 0.129577637f, 0.056533813f,\r
1415       0.031814575f, 0.003723145f, 0.002899170f, 0.000289917f, -0.000015259f, -0.000686646f,\r
1416       0.003479004f, -0.011886597f, 0.026535034f, -0.105819702f, 0.044784546f, -0.711318970f,\r
1417       1.127746582f, 0.431655884f, 0.134887695f, 0.051132202f, 0.031661987f, 0.003005981f,\r
1418       0.002792358f, 0.000259399f, -0.000015259f, -0.000747681f, 0.003479004f, -0.012939453f,\r
1419       0.025085449f, -0.110946655f, 0.031082153f, -0.738372803f, 1.120223999f, 0.404083252f,\r
1420       0.139450073f, 0.045837402f, 0.031387329f, 0.002334595f, 0.002685547f, 0.000244141f,\r
1421       -0.000030518f, -0.000808716f, 0.003463745f, -0.014022827f, 0.023422241f, -0.115921021f,\r
1422       0.016510010f, -0.765029907f, 1.111373901f, 0.376800537f, 0.143264771f, 0.040634155f,\r
1423       0.031005859f, 0.001693726f, 0.002578735f, 0.000213623f, -0.000030518f, -0.000885010f,\r
1424       0.003417969f, -0.015121460f, 0.021575928f, -0.120697021f, 0.001068115f, -0.791213989f,\r
1425       1.101211548f, 0.349868774f, 0.146362305f, 0.035552979f, 0.030532837f, 0.001098633f,\r
1426       0.002456665f, 0.000198364f, -0.000030518f, -0.000961304f, 0.003372192f, -0.016235352f,\r
1427       0.019531250f, -0.125259399f, -0.015228271f, -0.816864014f, 1.089782715f, 0.323318481f,\r
1428       0.148773193f, 0.030609131f, 0.029937744f, 0.000549316f, 0.002349854f, 0.000167847f,\r
1429       -0.000030518f, -0.001037598f, 0.003280640f, -0.017349243f, 0.017257690f, -0.129562378f,\r
1430       -0.032379150f, -0.841949463f, 1.077117920f, 0.297210693f, 0.150497437f, 0.025817871f,\r
1431       0.029281616f, 0.000030518f, 0.002243042f, 0.000152588f, -0.000045776f, -0.001113892f,\r
1432       0.003173828f, -0.018463135f, 0.014801025f, -0.133590698f, -0.050354004f, -0.866363525f,\r
1433       1.063217163f, 0.271591187f, 0.151596069f, 0.021179199f, 0.028533936f, -0.000442505f,\r
1434       0.002120972f, 0.000137329f, -0.000045776f, -0.001205444f, 0.003051758f, -0.019577026f,\r
1435       0.012115479f, -0.137298584f, -0.069168091f, -0.890090942f, 1.048156738f, 0.246505737f,\r
1436       0.152069092f, 0.016708374f, 0.027725220f, -0.000869751f, 0.002014160f, 0.000122070f,\r
1437       -0.000061035f, -0.001296997f, 0.002883911f, -0.020690918f, 0.009231567f, -0.140670776f,\r
1438       -0.088775635f, -0.913055420f, 1.031936646f, 0.221984863f, 0.151962280f, 0.012420654f,\r
1439       0.026840210f, -0.001266479f, 0.001907349f, 0.000106812f, -0.000061035f, -0.001388550f,\r
1440       0.002700806f, -0.021789551f, 0.006134033f, -0.143676758f, -0.109161377f, -0.935195923f,\r
1441       1.014617920f, 0.198059082f, 0.151306152f, 0.008316040f, 0.025909424f, -0.001617432f,\r
1442       0.001785278f, 0.000106812f, -0.000076294f, -0.001480103f, 0.002487183f, -0.022857666f,\r
1443       0.002822876f, -0.146255493f, -0.130310059f, -0.956481934f, 0.996246338f, 0.174789429f,\r
1444       0.150115967f, 0.004394531f, 0.024932861f, -0.001937866f, 0.001693726f, 0.000091553f,\r
1445       -0.000076294f, -0.001586914f, 0.002227783f, -0.023910522f, -0.000686646f, -0.148422241f,\r
1446       -0.152206421f, -0.976852417f, 0.976852417f, 0.152206421f, 0.148422241f, 0.000686646f,\r
1447       0.023910522f, -0.002227783f, 0.001586914f, 0.000076294f, -0.000091553f, -0.001693726f,\r
1448       0.001937866f, -0.024932861f, -0.004394531f, -0.150115967f, -0.174789429f, -0.996246338f,\r
1449       0.956481934f, 0.130310059f, 0.146255493f, -0.002822876f, 0.022857666f, -0.002487183f,\r
1450       0.001480103f, 0.000076294f, -0.000106812f, -0.001785278f, 0.001617432f, -0.025909424f,\r
1451       -0.008316040f, -0.151306152f, -0.198059082f, -1.014617920f, 0.935195923f, 0.109161377f,\r
1452       0.143676758f, -0.006134033f, 0.021789551f, -0.002700806f, 0.001388550f, 0.000061035f,\r
1453       -0.000106812f, -0.001907349f, 0.001266479f, -0.026840210f, -0.012420654f, -0.151962280f,\r
1454       -0.221984863f, -1.031936646f, 0.913055420f, 0.088775635f, 0.140670776f, -0.009231567f,\r
1455       0.020690918f, -0.002883911f, 0.001296997f, 0.000061035f, -0.000122070f, -0.002014160f,\r
1456       0.000869751f, -0.027725220f, -0.016708374f, -0.152069092f, -0.246505737f, -1.048156738f,\r
1457       0.890090942f, 0.069168091f, 0.137298584f, -0.012115479f, 0.019577026f, -0.003051758f,\r
1458       0.001205444f, 0.000045776f, -0.000137329f, -0.002120972f, 0.000442505f, -0.028533936f,\r
1459       -0.021179199f, -0.151596069f, -0.271591187f, -1.063217163f, 0.866363525f, 0.050354004f,\r
1460       0.133590698f, -0.014801025f, 0.018463135f, -0.003173828f, 0.001113892f, 0.000045776f,\r
1461       -0.000152588f, -0.002243042f, -0.000030518f, -0.029281616f, -0.025817871f, -0.150497437f,\r
1462       -0.297210693f, -1.077117920f, 0.841949463f, 0.032379150f, 0.129562378f, -0.017257690f,\r
1463       0.017349243f, -0.003280640f, 0.001037598f, 0.000030518f, -0.000167847f, -0.002349854f,\r
1464       -0.000549316f, -0.029937744f, -0.030609131f, -0.148773193f, -0.323318481f, -1.089782715f,\r
1465       0.816864014f, 0.015228271f, 0.125259399f, -0.019531250f, 0.016235352f, -0.003372192f,\r
1466       0.000961304f, 0.000030518f, -0.000198364f, -0.002456665f, -0.001098633f, -0.030532837f,\r
1467       -0.035552979f, -0.146362305f, -0.349868774f, -1.101211548f, 0.791213989f, -0.001068115f,\r
1468       0.120697021f, -0.021575928f, 0.015121460f, -0.003417969f, 0.000885010f, 0.000030518f,\r
1469       -0.000213623f, -0.002578735f, -0.001693726f, -0.031005859f, -0.040634155f, -0.143264771f,\r
1470       -0.376800537f, -1.111373901f, 0.765029907f, -0.016510010f, 0.115921021f, -0.023422241f,\r
1471       0.014022827f, -0.003463745f, 0.000808716f, 0.000030518f, -0.000244141f, -0.002685547f,\r
1472       -0.002334595f, -0.031387329f, -0.045837402f, -0.139450073f, -0.404083252f, -1.120223999f,\r
1473       0.738372803f, -0.031082153f, 0.110946655f, -0.025085449f, 0.012939453f, -0.003479004f,\r
1474       0.000747681f, 0.000015259f, -0.000259399f, -0.002792358f, -0.003005981f, -0.031661987f,\r
1475       -0.051132202f, -0.134887695f, -0.431655884f, -1.127746582f, 0.711318970f, -0.044784546f,\r
1476       0.105819702f, -0.026535034f, 0.011886597f, -0.003479004f, 0.000686646f, 0.000015259f,\r
1477       -0.000289917f, -0.002899170f, -0.003723145f, -0.031814575f, -0.056533813f, -0.129577637f,\r
1478       -0.459472656f, -1.133926392f, 0.683914185f, -0.057617188f, 0.100540161f, -0.027801514f,\r
1479       0.010848999f, -0.003463745f, 0.000625610f, 0.000015259f, -0.000320435f, -0.002990723f,\r
1480       -0.004486084f, -0.031845093f, -0.061996460f, -0.123474121f, -0.487472534f, -1.138763428f,\r
1481       0.656219482f, -0.069595337f, 0.095169067f, -0.028884888f, 0.009841919f, -0.003433228f,\r
1482       0.000579834f, 0.000015259f, -0.000366211f, -0.003082275f, -0.005294800f, -0.031738281f,\r
1483       -0.067520142f, -0.116577148f, -0.515609741f, -1.142211914f, 0.628295898f, -0.080688477f,\r
1484       0.089706421f, -0.029785156f, 0.008865356f, -0.003387451f, 0.000534058f, 0.000015259f,\r
1485       -0.000396729f, -0.003173828f, -0.006118774f, -0.031478882f, -0.073059082f, -0.108856201f,\r
1486       -0.543823242f, -1.144287109f, 0.600219727f, -0.090927124f, 0.084182739f, -0.030517578f,\r
1487       0.007919312f, -0.003326416f, 0.000473022f, 0.000015259f };\r
1488 \r
1489 }\r