开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > 流媒体/Mpeg4/MP4 > 查看源码
mpeglayer2.cpp
资源名称：NETVIDEO.rar　[点击查看]
上传用户：sun1608
上传日期：2007-02-02
资源大小：6116k
文件大小：36k
源码类别：

流媒体/Mpeg4/MP4

开发平台：

Visual C++

mpeglayer2.cpp：源码内容
  1. /* MPEG/WAVE Sound library
  3.    (C) 1997 by Jung woo-jae */
  5. // Mpeglayer2.cc
  6. // It's for MPEG Layer 2
  8. #ifdef HAVE_CONFIG_H
  9. #include "config.h"
  10. #endif
  11. #if defined(_WIN32) && defined(_MSC_VER)
  12. // disable warnings about double to float conversions
  13. #pragma warning(disable: 4244 4305)
  14. #endif
  16. #include "MPEGaudio.h"
  18. #define MAXTABLE 2
  20. // Tables for layer 2
  21. static const int bitalloclengthtable[MAXTABLE][MAXSUBBAND]=
  22. {{4,4,3,3,3,3,3,3,3,3,3,3,3,3,3,3, 3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3},
  23.  {4,4,4,4,4,4,4,4,4,4,4,3,3,3,3,3, 3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2}};
  25. static const REAL group5bits[27*3]=
  26. {
  27.   -2.0/3.0, -2.0/3.0, -2.0/3.0,
  28.        0.0, -2.0/3.0, -2.0/3.0,
  29.    2.0/3.0, -2.0/3.0, -2.0/3.0,
  30.   -2.0/3.0,      0.0, -2.0/3.0,
  31.        0.0,      0.0, -2.0/3.0,
  32.    2.0/3.0,      0.0, -2.0/3.0,
  33.   -2.0/3.0,  2.0/3.0, -2.0/3.0,
  34.        0.0,  2.0/3.0, -2.0/3.0,
  35.    2.0/3.0,  2.0/3.0, -2.0/3.0,
  36.   -2.0/3.0, -2.0/3.0,      0.0,
  37.        0.0, -2.0/3.0,      0.0,
  38.    2.0/3.0, -2.0/3.0,      0.0,
  39.   -2.0/3.0,      0.0,      0.0,
  40.        0.0,      0.0,      0.0,
  41.    2.0/3.0,      0.0,      0.0,
  42.   -2.0/3.0,  2.0/3.0,      0.0,
  43.        0.0,  2.0/3.0,      0.0,
  44.    2.0/3.0,  2.0/3.0,      0.0,
  45.   -2.0/3.0, -2.0/3.0,  2.0/3.0,
  46.        0.0, -2.0/3.0,  2.0/3.0,
  47.    2.0/3.0, -2.0/3.0,  2.0/3.0,
  48.   -2.0/3.0,      0.0,  2.0/3.0,
  49.        0.0,      0.0,  2.0/3.0,
  50.    2.0/3.0,      0.0,  2.0/3.0,
  51.   -2.0/3.0,  2.0/3.0,  2.0/3.0,
  52.        0.0,  2.0/3.0,  2.0/3.0,
  53.    2.0/3.0,  2.0/3.0,  2.0/3.0
  54. };
  56. static const REAL group7bits[125*3]=
  57. {
  58.  -0.8,-0.8,-0.8, -0.4,-0.8,-0.8, 0.0,-0.8,-0.8, 0.4,-0.8,-0.8, 0.8,-0.8,-0.8,
  59.  -0.8,-0.4,-0.8, -0.4,-0.4,-0.8, 0.0,-0.4,-0.8, 0.4,-0.4,-0.8, 0.8,-0.4,-0.8,
  60.  -0.8, 0.0,-0.8, -0.4, 0.0,-0.8, 0.0, 0.0,-0.8, 0.4, 0.0,-0.8, 0.8, 0.0,-0.8,
  61.  -0.8, 0.4,-0.8, -0.4, 0.4,-0.8, 0.0, 0.4,-0.8, 0.4, 0.4,-0.8, 0.8, 0.4,-0.8,
  62.  -0.8, 0.8,-0.8, -0.4, 0.8,-0.8, 0.0, 0.8,-0.8, 0.4, 0.8,-0.8, 0.8, 0.8,-0.8,
  63.  -0.8,-0.8,-0.4, -0.4,-0.8,-0.4, 0.0,-0.8,-0.4, 0.4,-0.8,-0.4, 0.8,-0.8,-0.4,
  64.  -0.8,-0.4,-0.4, -0.4,-0.4,-0.4, 0.0,-0.4,-0.4, 0.4,-0.4,-0.4, 0.8,-0.4,-0.4,
  65.  -0.8, 0.0,-0.4, -0.4, 0.0,-0.4, 0.0, 0.0,-0.4, 0.4, 0.0,-0.4, 0.8, 0.0,-0.4,
  66.  -0.8, 0.4,-0.4, -0.4, 0.4,-0.4, 0.0, 0.4,-0.4, 0.4, 0.4,-0.4, 0.8, 0.4,-0.4,
  67.  -0.8, 0.8,-0.4, -0.4, 0.8,-0.4, 0.0, 0.8,-0.4, 0.4, 0.8,-0.4, 0.8, 0.8,-0.4,
  68.  -0.8,-0.8, 0.0, -0.4,-0.8, 0.0, 0.0,-0.8, 0.0, 0.4,-0.8, 0.0, 0.8,-0.8, 0.0,
  69.  -0.8,-0.4, 0.0, -0.4,-0.4, 0.0, 0.0,-0.4, 0.0, 0.4,-0.4, 0.0, 0.8,-0.4, 0.0,
  70.  -0.8, 0.0, 0.0, -0.4, 0.0, 0.0, 0.0, 0.0, 0.0, 0.4, 0.0, 0.0, 0.8, 0.0, 0.0,
  71.  -0.8, 0.4, 0.0, -0.4, 0.4, 0.0, 0.0, 0.4, 0.0, 0.4, 0.4, 0.0, 0.8, 0.4, 0.0,
  72.  -0.8, 0.8, 0.0, -0.4, 0.8, 0.0, 0.0, 0.8, 0.0, 0.4, 0.8, 0.0, 0.8, 0.8, 0.0,
  73.  -0.8,-0.8, 0.4, -0.4,-0.8, 0.4, 0.0,-0.8, 0.4, 0.4,-0.8, 0.4, 0.8,-0.8, 0.4,
  74.  -0.8,-0.4, 0.4, -0.4,-0.4, 0.4, 0.0,-0.4, 0.4, 0.4,-0.4, 0.4, 0.8,-0.4, 0.4,
  75.  -0.8, 0.0, 0.4, -0.4, 0.0, 0.4, 0.0, 0.0, 0.4, 0.4, 0.0, 0.4, 0.8, 0.0, 0.4,
  76.  -0.8, 0.4, 0.4, -0.4, 0.4, 0.4, 0.0, 0.4, 0.4, 0.4, 0.4, 0.4, 0.8, 0.4, 0.4,
  77.  -0.8, 0.8, 0.4, -0.4, 0.8, 0.4, 0.0, 0.8, 0.4, 0.4, 0.8, 0.4, 0.8, 0.8, 0.4,
  78.  -0.8,-0.8, 0.8, -0.4,-0.8, 0.8, 0.0,-0.8, 0.8, 0.4,-0.8, 0.8, 0.8,-0.8, 0.8,
  79.  -0.8,-0.4, 0.8, -0.4,-0.4, 0.8, 0.0,-0.4, 0.8, 0.4,-0.4, 0.8, 0.8,-0.4, 0.8,
  80.  -0.8, 0.0, 0.8, -0.4, 0.0, 0.8, 0.0, 0.0, 0.8, 0.4, 0.0, 0.8, 0.8, 0.0, 0.8,
  81.  -0.8, 0.4, 0.8, -0.4, 0.4, 0.8, 0.0, 0.4, 0.8, 0.4, 0.4, 0.8, 0.8, 0.4, 0.8,
  82.  -0.8, 0.8, 0.8, -0.4, 0.8, 0.8, 0.0, 0.8, 0.8, 0.4, 0.8, 0.8, 0.8, 0.8, 0.8
  83. };
  85. static const REAL group10bits[729*3]=
  86. {
  87.  -8.0/9.0,-8.0/9.0,-8.0/9.0, -6.0/9.0,-8.0/9.0,-8.0/9.0, -4.0/9.0,-8.0/9.0,-8.0/9.0,
  88.  -2.0/9.0,-8.0/9.0,-8.0/9.0,      0.0,-8.0/9.0,-8.0/9.0,  2.0/9.0,-8.0/9.0,-8.0/9.0,
  89.   4.0/9.0,-8.0/9.0,-8.0/9.0,  6.0/9.0,-8.0/9.0,-8.0/9.0,  8.0/9.0,-8.0/9.0,-8.0/9.0,
  90.  -8.0/9.0,-6.0/9.0,-8.0/9.0, -6.0/9.0,-6.0/9.0,-8.0/9.0, -4.0/9.0,-6.0/9.0,-8.0/9.0,
  91.  -2.0/9.0,-6.0/9.0,-8.0/9.0,      0.0,-6.0/9.0,-8.0/9.0,  2.0/9.0,-6.0/9.0,-8.0/9.0,
  92.   4.0/9.0,-6.0/9.0,-8.0/9.0,  6.0/9.0,-6.0/9.0,-8.0/9.0,  8.0/9.0,-6.0/9.0,-8.0/9.0,
  93.  -8.0/9.0,-4.0/9.0,-8.0/9.0, -6.0/9.0,-4.0/9.0,-8.0/9.0, -4.0/9.0,-4.0/9.0,-8.0/9.0,
  94.  -2.0/9.0,-4.0/9.0,-8.0/9.0,      0.0,-4.0/9.0,-8.0/9.0,  2.0/9.0,-4.0/9.0,-8.0/9.0,
  95.   4.0/9.0,-4.0/9.0,-8.0/9.0,  6.0/9.0,-4.0/9.0,-8.0/9.0,  8.0/9.0,-4.0/9.0,-8.0/9.0,
  96.  -8.0/9.0,-2.0/9.0,-8.0/9.0, -6.0/9.0,-2.0/9.0,-8.0/9.0, -4.0/9.0,-2.0/9.0,-8.0/9.0,
  97.  -2.0/9.0,-2.0/9.0,-8.0/9.0,      0.0,-2.0/9.0,-8.0/9.0,  2.0/9.0,-2.0/9.0,-8.0/9.0,
  98.   4.0/9.0,-2.0/9.0,-8.0/9.0,  6.0/9.0,-2.0/9.0,-8.0/9.0,  8.0/9.0,-2.0/9.0,-8.0/9.0,
  99.  -8.0/9.0,     0.0,-8.0/9.0, -6.0/9.0,     0.0,-8.0/9.0, -4.0/9.0,     0.0,-8.0/9.0,
  100.  -2.0/9.0,     0.0,-8.0/9.0,      0.0,     0.0,-8.0/9.0,  2.0/9.0,     0.0,-8.0/9.0,
  101.   4.0/9.0,     0.0,-8.0/9.0,  6.0/9.0,     0.0,-8.0/9.0,  8.0/9.0,     0.0,-8.0/9.0,
  102.  -8.0/9.0, 2.0/9.0,-8.0/9.0, -6.0/9.0, 2.0/9.0,-8.0/9.0, -4.0/9.0, 2.0/9.0,-8.0/9.0,
  103.  -2.0/9.0, 2.0/9.0,-8.0/9.0,      0.0, 2.0/9.0,-8.0/9.0,  2.0/9.0, 2.0/9.0,-8.0/9.0,
  104.   4.0/9.0, 2.0/9.0,-8.0/9.0,  6.0/9.0, 2.0/9.0,-8.0/9.0,  8.0/9.0, 2.0/9.0,-8.0/9.0,
  105.  -8.0/9.0, 4.0/9.0,-8.0/9.0, -6.0/9.0, 4.0/9.0,-8.0/9.0, -4.0/9.0, 4.0/9.0,-8.0/9.0,
  106.  -2.0/9.0, 4.0/9.0,-8.0/9.0,      0.0, 4.0/9.0,-8.0/9.0,  2.0/9.0, 4.0/9.0,-8.0/9.0,
  107.   4.0/9.0, 4.0/9.0,-8.0/9.0,  6.0/9.0, 4.0/9.0,-8.0/9.0,  8.0/9.0, 4.0/9.0,-8.0/9.0,
  108.  -8.0/9.0, 6.0/9.0,-8.0/9.0, -6.0/9.0, 6.0/9.0,-8.0/9.0, -4.0/9.0, 6.0/9.0,-8.0/9.0,
  109.  -2.0/9.0, 6.0/9.0,-8.0/9.0,      0.0, 6.0/9.0,-8.0/9.0,  2.0/9.0, 6.0/9.0,-8.0/9.0,
  110.   4.0/9.0, 6.0/9.0,-8.0/9.0,  6.0/9.0, 6.0/9.0,-8.0/9.0,  8.0/9.0, 6.0/9.0,-8.0/9.0,
  111.  -8.0/9.0, 8.0/9.0,-8.0/9.0, -6.0/9.0, 8.0/9.0,-8.0/9.0, -4.0/9.0, 8.0/9.0,-8.0/9.0,
  112.  -2.0/9.0, 8.0/9.0,-8.0/9.0,      0.0, 8.0/9.0,-8.0/9.0,  2.0/9.0, 8.0/9.0,-8.0/9.0,
  113.   4.0/9.0, 8.0/9.0,-8.0/9.0,  6.0/9.0, 8.0/9.0,-8.0/9.0,  8.0/9.0, 8.0/9.0,-8.0/9.0,
  114.  -8.0/9.0,-8.0/9.0,-6.0/9.0, -6.0/9.0,-8.0/9.0,-6.0/9.0, -4.0/9.0,-8.0/9.0,-6.0/9.0,
  115.  -2.0/9.0,-8.0/9.0,-6.0/9.0,      0.0,-8.0/9.0,-6.0/9.0,  2.0/9.0,-8.0/9.0,-6.0/9.0,
  116.   4.0/9.0,-8.0/9.0,-6.0/9.0,  6.0/9.0,-8.0/9.0,-6.0/9.0,  8.0/9.0,-8.0/9.0,-6.0/9.0,
  117.  -8.0/9.0,-6.0/9.0,-6.0/9.0, -6.0/9.0,-6.0/9.0,-6.0/9.0, -4.0/9.0,-6.0/9.0,-6.0/9.0,
  118.  -2.0/9.0,-6.0/9.0,-6.0/9.0,      0.0,-6.0/9.0,-6.0/9.0,  2.0/9.0,-6.0/9.0,-6.0/9.0,
  119.   4.0/9.0,-6.0/9.0,-6.0/9.0,  6.0/9.0,-6.0/9.0,-6.0/9.0,  8.0/9.0,-6.0/9.0,-6.0/9.0,
  120.  -8.0/9.0,-4.0/9.0,-6.0/9.0, -6.0/9.0,-4.0/9.0,-6.0/9.0, -4.0/9.0,-4.0/9.0,-6.0/9.0,
  121.  -2.0/9.0,-4.0/9.0,-6.0/9.0,      0.0,-4.0/9.0,-6.0/9.0,  2.0/9.0,-4.0/9.0,-6.0/9.0,
  122.   4.0/9.0,-4.0/9.0,-6.0/9.0,  6.0/9.0,-4.0/9.0,-6.0/9.0,  8.0/9.0,-4.0/9.0,-6.0/9.0,
  123.  -8.0/9.0,-2.0/9.0,-6.0/9.0, -6.0/9.0,-2.0/9.0,-6.0/9.0, -4.0/9.0,-2.0/9.0,-6.0/9.0,
  124.  -2.0/9.0,-2.0/9.0,-6.0/9.0,      0.0,-2.0/9.0,-6.0/9.0,  2.0/9.0,-2.0/9.0,-6.0/9.0,
  125.   4.0/9.0,-2.0/9.0,-6.0/9.0,  6.0/9.0,-2.0/9.0,-6.0/9.0,  8.0/9.0,-2.0/9.0,-6.0/9.0,
  126.  -8.0/9.0,     0.0,-6.0/9.0, -6.0/9.0,     0.0,-6.0/9.0, -4.0/9.0,     0.0,-6.0/9.0,
  127.  -2.0/9.0,     0.0,-6.0/9.0,      0.0,     0.0,-6.0/9.0,  2.0/9.0,     0.0,-6.0/9.0,
  128.   4.0/9.0,     0.0,-6.0/9.0,  6.0/9.0,     0.0,-6.0/9.0,  8.0/9.0,     0.0,-6.0/9.0,
  129.  -8.0/9.0, 2.0/9.0,-6.0/9.0, -6.0/9.0, 2.0/9.0,-6.0/9.0, -4.0/9.0, 2.0/9.0,-6.0/9.0,
  130.  -2.0/9.0, 2.0/9.0,-6.0/9.0,      0.0, 2.0/9.0,-6.0/9.0,  2.0/9.0, 2.0/9.0,-6.0/9.0,
  131.   4.0/9.0, 2.0/9.0,-6.0/9.0,  6.0/9.0, 2.0/9.0,-6.0/9.0,  8.0/9.0, 2.0/9.0,-6.0/9.0,
  132.  -8.0/9.0, 4.0/9.0,-6.0/9.0, -6.0/9.0, 4.0/9.0,-6.0/9.0, -4.0/9.0, 4.0/9.0,-6.0/9.0,
  133.  -2.0/9.0, 4.0/9.0,-6.0/9.0,      0.0, 4.0/9.0,-6.0/9.0,  2.0/9.0, 4.0/9.0,-6.0/9.0,
  134.   4.0/9.0, 4.0/9.0,-6.0/9.0,  6.0/9.0, 4.0/9.0,-6.0/9.0,  8.0/9.0, 4.0/9.0,-6.0/9.0,
  135.  -8.0/9.0, 6.0/9.0,-6.0/9.0, -6.0/9.0, 6.0/9.0,-6.0/9.0, -4.0/9.0, 6.0/9.0,-6.0/9.0,
  136.  -2.0/9.0, 6.0/9.0,-6.0/9.0,      0.0, 6.0/9.0,-6.0/9.0,  2.0/9.0, 6.0/9.0,-6.0/9.0,
  137.   4.0/9.0, 6.0/9.0,-6.0/9.0,  6.0/9.0, 6.0/9.0,-6.0/9.0,  8.0/9.0, 6.0/9.0,-6.0/9.0,
  138.  -8.0/9.0, 8.0/9.0,-6.0/9.0, -6.0/9.0, 8.0/9.0,-6.0/9.0, -4.0/9.0, 8.0/9.0,-6.0/9.0,
  139.  -2.0/9.0, 8.0/9.0,-6.0/9.0,      0.0, 8.0/9.0,-6.0/9.0,  2.0/9.0, 8.0/9.0,-6.0/9.0,
  140.   4.0/9.0, 8.0/9.0,-6.0/9.0,  6.0/9.0, 8.0/9.0,-6.0/9.0,  8.0/9.0, 8.0/9.0,-6.0/9.0,
  141.  -8.0/9.0,-8.0/9.0,-4.0/9.0, -6.0/9.0,-8.0/9.0,-4.0/9.0, -4.0/9.0,-8.0/9.0,-4.0/9.0,
  142.  -2.0/9.0,-8.0/9.0,-4.0/9.0,      0.0,-8.0/9.0,-4.0/9.0,  2.0/9.0,-8.0/9.0,-4.0/9.0,
  143.   4.0/9.0,-8.0/9.0,-4.0/9.0,  6.0/9.0,-8.0/9.0,-4.0/9.0,  8.0/9.0,-8.0/9.0,-4.0/9.0,
  144.  -8.0/9.0,-6.0/9.0,-4.0/9.0, -6.0/9.0,-6.0/9.0,-4.0/9.0, -4.0/9.0,-6.0/9.0,-4.0/9.0,
  145.  -2.0/9.0,-6.0/9.0,-4.0/9.0,      0.0,-6.0/9.0,-4.0/9.0,  2.0/9.0,-6.0/9.0,-4.0/9.0,
  146.   4.0/9.0,-6.0/9.0,-4.0/9.0,  6.0/9.0,-6.0/9.0,-4.0/9.0,  8.0/9.0,-6.0/9.0,-4.0/9.0,
  147.  -8.0/9.0,-4.0/9.0,-4.0/9.0, -6.0/9.0,-4.0/9.0,-4.0/9.0, -4.0/9.0,-4.0/9.0,-4.0/9.0,
  148.  -2.0/9.0,-4.0/9.0,-4.0/9.0,      0.0,-4.0/9.0,-4.0/9.0,  2.0/9.0,-4.0/9.0,-4.0/9.0,
  149.   4.0/9.0,-4.0/9.0,-4.0/9.0,  6.0/9.0,-4.0/9.0,-4.0/9.0,  8.0/9.0,-4.0/9.0,-4.0/9.0,
  150.  -8.0/9.0,-2.0/9.0,-4.0/9.0, -6.0/9.0,-2.0/9.0,-4.0/9.0, -4.0/9.0,-2.0/9.0,-4.0/9.0,
  151.  -2.0/9.0,-2.0/9.0,-4.0/9.0,      0.0,-2.0/9.0,-4.0/9.0,  2.0/9.0,-2.0/9.0,-4.0/9.0,
  152.   4.0/9.0,-2.0/9.0,-4.0/9.0,  6.0/9.0,-2.0/9.0,-4.0/9.0,  8.0/9.0,-2.0/9.0,-4.0/9.0,
  153.  -8.0/9.0,     0.0,-4.0/9.0, -6.0/9.0,     0.0,-4.0/9.0, -4.0/9.0,     0.0,-4.0/9.0,
  154.  -2.0/9.0,     0.0,-4.0/9.0,      0.0,     0.0,-4.0/9.0,  2.0/9.0,     0.0,-4.0/9.0,
  155.   4.0/9.0,     0.0,-4.0/9.0,  6.0/9.0,     0.0,-4.0/9.0,  8.0/9.0,     0.0,-4.0/9.0,
  156.  -8.0/9.0, 2.0/9.0,-4.0/9.0, -6.0/9.0, 2.0/9.0,-4.0/9.0, -4.0/9.0, 2.0/9.0,-4.0/9.0,
  157.  -2.0/9.0, 2.0/9.0,-4.0/9.0,      0.0, 2.0/9.0,-4.0/9.0,  2.0/9.0, 2.0/9.0,-4.0/9.0,
  158.   4.0/9.0, 2.0/9.0,-4.0/9.0,  6.0/9.0, 2.0/9.0,-4.0/9.0,  8.0/9.0, 2.0/9.0,-4.0/9.0,
  159.  -8.0/9.0, 4.0/9.0,-4.0/9.0, -6.0/9.0, 4.0/9.0,-4.0/9.0, -4.0/9.0, 4.0/9.0,-4.0/9.0,
  160.  -2.0/9.0, 4.0/9.0,-4.0/9.0,      0.0, 4.0/9.0,-4.0/9.0,  2.0/9.0, 4.0/9.0,-4.0/9.0,
  161.   4.0/9.0, 4.0/9.0,-4.0/9.0,  6.0/9.0, 4.0/9.0,-4.0/9.0,  8.0/9.0, 4.0/9.0,-4.0/9.0,
  162.  -8.0/9.0, 6.0/9.0,-4.0/9.0, -6.0/9.0, 6.0/9.0,-4.0/9.0, -4.0/9.0, 6.0/9.0,-4.0/9.0,
  163.  -2.0/9.0, 6.0/9.0,-4.0/9.0,      0.0, 6.0/9.0,-4.0/9.0,  2.0/9.0, 6.0/9.0,-4.0/9.0,
  164.   4.0/9.0, 6.0/9.0,-4.0/9.0,  6.0/9.0, 6.0/9.0,-4.0/9.0,  8.0/9.0, 6.0/9.0,-4.0/9.0,
  165.  -8.0/9.0, 8.0/9.0,-4.0/9.0, -6.0/9.0, 8.0/9.0,-4.0/9.0, -4.0/9.0, 8.0/9.0,-4.0/9.0,
  166.  -2.0/9.0, 8.0/9.0,-4.0/9.0,      0.0, 8.0/9.0,-4.0/9.0,  2.0/9.0, 8.0/9.0,-4.0/9.0,
  167.   4.0/9.0, 8.0/9.0,-4.0/9.0,  6.0/9.0, 8.0/9.0,-4.0/9.0,  8.0/9.0, 8.0/9.0,-4.0/9.0,
  168.  -8.0/9.0,-8.0/9.0,-2.0/9.0, -6.0/9.0,-8.0/9.0,-2.0/9.0, -4.0/9.0,-8.0/9.0,-2.0/9.0,
  169.  -2.0/9.0,-8.0/9.0,-2.0/9.0,      0.0,-8.0/9.0,-2.0/9.0,  2.0/9.0,-8.0/9.0,-2.0/9.0,
  170.   4.0/9.0,-8.0/9.0,-2.0/9.0,  6.0/9.0,-8.0/9.0,-2.0/9.0,  8.0/9.0,-8.0/9.0,-2.0/9.0,
  171.  -8.0/9.0,-6.0/9.0,-2.0/9.0, -6.0/9.0,-6.0/9.0,-2.0/9.0, -4.0/9.0,-6.0/9.0,-2.0/9.0,
  172.  -2.0/9.0,-6.0/9.0,-2.0/9.0,      0.0,-6.0/9.0,-2.0/9.0,  2.0/9.0,-6.0/9.0,-2.0/9.0,
  173.   4.0/9.0,-6.0/9.0,-2.0/9.0,  6.0/9.0,-6.0/9.0,-2.0/9.0,  8.0/9.0,-6.0/9.0,-2.0/9.0,
  174.  -8.0/9.0,-4.0/9.0,-2.0/9.0, -6.0/9.0,-4.0/9.0,-2.0/9.0, -4.0/9.0,-4.0/9.0,-2.0/9.0,
  175.  -2.0/9.0,-4.0/9.0,-2.0/9.0,      0.0,-4.0/9.0,-2.0/9.0,  2.0/9.0,-4.0/9.0,-2.0/9.0,
  176.   4.0/9.0,-4.0/9.0,-2.0/9.0,  6.0/9.0,-4.0/9.0,-2.0/9.0,  8.0/9.0,-4.0/9.0,-2.0/9.0,
  177.  -8.0/9.0,-2.0/9.0,-2.0/9.0, -6.0/9.0,-2.0/9.0,-2.0/9.0, -4.0/9.0,-2.0/9.0,-2.0/9.0,
  178.  -2.0/9.0,-2.0/9.0,-2.0/9.0,      0.0,-2.0/9.0,-2.0/9.0,  2.0/9.0,-2.0/9.0,-2.0/9.0,
  179.   4.0/9.0,-2.0/9.0,-2.0/9.0,  6.0/9.0,-2.0/9.0,-2.0/9.0,  8.0/9.0,-2.0/9.0,-2.0/9.0,
  180.  -8.0/9.0,     0.0,-2.0/9.0, -6.0/9.0,     0.0,-2.0/9.0, -4.0/9.0,     0.0,-2.0/9.0,
  181.  -2.0/9.0,     0.0,-2.0/9.0,      0.0,     0.0,-2.0/9.0,  2.0/9.0,     0.0,-2.0/9.0,
  182.   4.0/9.0,     0.0,-2.0/9.0,  6.0/9.0,     0.0,-2.0/9.0,  8.0/9.0,     0.0,-2.0/9.0,
  183.  -8.0/9.0, 2.0/9.0,-2.0/9.0, -6.0/9.0, 2.0/9.0,-2.0/9.0, -4.0/9.0, 2.0/9.0,-2.0/9.0,
  184.  -2.0/9.0, 2.0/9.0,-2.0/9.0,      0.0, 2.0/9.0,-2.0/9.0,  2.0/9.0, 2.0/9.0,-2.0/9.0,
  185.   4.0/9.0, 2.0/9.0,-2.0/9.0,  6.0/9.0, 2.0/9.0,-2.0/9.0,  8.0/9.0, 2.0/9.0,-2.0/9.0,
  186.  -8.0/9.0, 4.0/9.0,-2.0/9.0, -6.0/9.0, 4.0/9.0,-2.0/9.0, -4.0/9.0, 4.0/9.0,-2.0/9.0,
  187.  -2.0/9.0, 4.0/9.0,-2.0/9.0,      0.0, 4.0/9.0,-2.0/9.0,  2.0/9.0, 4.0/9.0,-2.0/9.0,
  188.   4.0/9.0, 4.0/9.0,-2.0/9.0,  6.0/9.0, 4.0/9.0,-2.0/9.0,  8.0/9.0, 4.0/9.0,-2.0/9.0,
  189.  -8.0/9.0, 6.0/9.0,-2.0/9.0, -6.0/9.0, 6.0/9.0,-2.0/9.0, -4.0/9.0, 6.0/9.0,-2.0/9.0,
  190.  -2.0/9.0, 6.0/9.0,-2.0/9.0,      0.0, 6.0/9.0,-2.0/9.0,  2.0/9.0, 6.0/9.0,-2.0/9.0,
  191.   4.0/9.0, 6.0/9.0,-2.0/9.0,  6.0/9.0, 6.0/9.0,-2.0/9.0,  8.0/9.0, 6.0/9.0,-2.0/9.0,
  192.  -8.0/9.0, 8.0/9.0,-2.0/9.0, -6.0/9.0, 8.0/9.0,-2.0/9.0, -4.0/9.0, 8.0/9.0,-2.0/9.0,
  193.  -2.0/9.0, 8.0/9.0,-2.0/9.0,      0.0, 8.0/9.0,-2.0/9.0,  2.0/9.0, 8.0/9.0,-2.0/9.0,
  194.   4.0/9.0, 8.0/9.0,-2.0/9.0,  6.0/9.0, 8.0/9.0,-2.0/9.0,  8.0/9.0, 8.0/9.0,-2.0/9.0,
  195.  -8.0/9.0,-8.0/9.0,     0.0, -6.0/9.0,-8.0/9.0,     0.0, -4.0/9.0,-8.0/9.0,     0.0,
  196.  -2.0/9.0,-8.0/9.0,     0.0,      0.0,-8.0/9.0,     0.0,  2.0/9.0,-8.0/9.0,     0.0,
  197.   4.0/9.0,-8.0/9.0,     0.0,  6.0/9.0,-8.0/9.0,     0.0,  8.0/9.0,-8.0/9.0,     0.0,
  198.  -8.0/9.0,-6.0/9.0,     0.0, -6.0/9.0,-6.0/9.0,     0.0, -4.0/9.0,-6.0/9.0,     0.0,
  199.  -2.0/9.0,-6.0/9.0,     0.0,      0.0,-6.0/9.0,     0.0,  2.0/9.0,-6.0/9.0,     0.0,
  200.   4.0/9.0,-6.0/9.0,     0.0,  6.0/9.0,-6.0/9.0,     0.0,  8.0/9.0,-6.0/9.0,     0.0,
  201.  -8.0/9.0,-4.0/9.0,     0.0, -6.0/9.0,-4.0/9.0,     0.0, -4.0/9.0,-4.0/9.0,     0.0,
  202.  -2.0/9.0,-4.0/9.0,     0.0,      0.0,-4.0/9.0,     0.0,  2.0/9.0,-4.0/9.0,     0.0,
  203.   4.0/9.0,-4.0/9.0,     0.0,  6.0/9.0,-4.0/9.0,     0.0,  8.0/9.0,-4.0/9.0,     0.0,
  204.  -8.0/9.0,-2.0/9.0,     0.0, -6.0/9.0,-2.0/9.0,     0.0, -4.0/9.0,-2.0/9.0,     0.0,
  205.  -2.0/9.0,-2.0/9.0,     0.0,      0.0,-2.0/9.0,     0.0,  2.0/9.0,-2.0/9.0,     0.0,
  206.   4.0/9.0,-2.0/9.0,     0.0,  6.0/9.0,-2.0/9.0,     0.0,  8.0/9.0,-2.0/9.0,     0.0,
  207.  -8.0/9.0,     0.0,     0.0, -6.0/9.0,     0.0,     0.0, -4.0/9.0,     0.0,     0.0,
  208.  -2.0/9.0,     0.0,     0.0,      0.0,     0.0,     0.0,  2.0/9.0,     0.0,     0.0,
  209.   4.0/9.0,     0.0,     0.0,  6.0/9.0,     0.0,     0.0,  8.0/9.0,     0.0,     0.0,
  210.  -8.0/9.0, 2.0/9.0,     0.0, -6.0/9.0, 2.0/9.0,     0.0, -4.0/9.0, 2.0/9.0,     0.0,
  211.  -2.0/9.0, 2.0/9.0,     0.0,      0.0, 2.0/9.0,     0.0,  2.0/9.0, 2.0/9.0,     0.0,
  212.   4.0/9.0, 2.0/9.0,     0.0,  6.0/9.0, 2.0/9.0,     0.0,  8.0/9.0, 2.0/9.0,     0.0,
  213.  -8.0/9.0, 4.0/9.0,     0.0, -6.0/9.0, 4.0/9.0,     0.0, -4.0/9.0, 4.0/9.0,     0.0,
  214.  -2.0/9.0, 4.0/9.0,     0.0,      0.0, 4.0/9.0,     0.0,  2.0/9.0, 4.0/9.0,     0.0,
  215.   4.0/9.0, 4.0/9.0,     0.0,  6.0/9.0, 4.0/9.0,     0.0,  8.0/9.0, 4.0/9.0,     0.0,
  216.  -8.0/9.0, 6.0/9.0,     0.0, -6.0/9.0, 6.0/9.0,     0.0, -4.0/9.0, 6.0/9.0,     0.0,
  217.  -2.0/9.0, 6.0/9.0,     0.0,      0.0, 6.0/9.0,     0.0,  2.0/9.0, 6.0/9.0,     0.0,
  218.   4.0/9.0, 6.0/9.0,     0.0,  6.0/9.0, 6.0/9.0,     0.0,  8.0/9.0, 6.0/9.0,     0.0,
  219.  -8.0/9.0, 8.0/9.0,     0.0, -6.0/9.0, 8.0/9.0,     0.0, -4.0/9.0, 8.0/9.0,     0.0,
  220.  -2.0/9.0, 8.0/9.0,     0.0,      0.0, 8.0/9.0,     0.0,  2.0/9.0, 8.0/9.0,     0.0,
  221.   4.0/9.0, 8.0/9.0,     0.0,  6.0/9.0, 8.0/9.0,     0.0,  8.0/9.0, 8.0/9.0,     0.0,
  222.  -8.0/9.0,-8.0/9.0, 2.0/9.0, -6.0/9.0,-8.0/9.0, 2.0/9.0, -4.0/9.0,-8.0/9.0, 2.0/9.0,
  223.  -2.0/9.0,-8.0/9.0, 2.0/9.0,      0.0,-8.0/9.0, 2.0/9.0,  2.0/9.0,-8.0/9.0, 2.0/9.0,
  224.   4.0/9.0,-8.0/9.0, 2.0/9.0,  6.0/9.0,-8.0/9.0, 2.0/9.0,  8.0/9.0,-8.0/9.0, 2.0/9.0,
  225.  -8.0/9.0,-6.0/9.0, 2.0/9.0, -6.0/9.0,-6.0/9.0, 2.0/9.0, -4.0/9.0,-6.0/9.0, 2.0/9.0,
  226.  -2.0/9.0,-6.0/9.0, 2.0/9.0,      0.0,-6.0/9.0, 2.0/9.0,  2.0/9.0,-6.0/9.0, 2.0/9.0,
  227.   4.0/9.0,-6.0/9.0, 2.0/9.0,  6.0/9.0,-6.0/9.0, 2.0/9.0,  8.0/9.0,-6.0/9.0, 2.0/9.0,
  228.  -8.0/9.0,-4.0/9.0, 2.0/9.0, -6.0/9.0,-4.0/9.0, 2.0/9.0, -4.0/9.0,-4.0/9.0, 2.0/9.0,
  229.  -2.0/9.0,-4.0/9.0, 2.0/9.0,      0.0,-4.0/9.0, 2.0/9.0,  2.0/9.0,-4.0/9.0, 2.0/9.0,
  230.   4.0/9.0,-4.0/9.0, 2.0/9.0,  6.0/9.0,-4.0/9.0, 2.0/9.0,  8.0/9.0,-4.0/9.0, 2.0/9.0,
  231.  -8.0/9.0,-2.0/9.0, 2.0/9.0, -6.0/9.0,-2.0/9.0, 2.0/9.0, -4.0/9.0,-2.0/9.0, 2.0/9.0,
  232.  -2.0/9.0,-2.0/9.0, 2.0/9.0,      0.0,-2.0/9.0, 2.0/9.0,  2.0/9.0,-2.0/9.0, 2.0/9.0,
  233.   4.0/9.0,-2.0/9.0, 2.0/9.0,  6.0/9.0,-2.0/9.0, 2.0/9.0,  8.0/9.0,-2.0/9.0, 2.0/9.0,
  234.  -8.0/9.0,     0.0, 2.0/9.0, -6.0/9.0,     0.0, 2.0/9.0, -4.0/9.0,     0.0, 2.0/9.0,
  235.  -2.0/9.0,     0.0, 2.0/9.0,      0.0,     0.0, 2.0/9.0,  2.0/9.0,     0.0, 2.0/9.0,
  236.   4.0/9.0,     0.0, 2.0/9.0,  6.0/9.0,     0.0, 2.0/9.0,  8.0/9.0,     0.0, 2.0/9.0,
  237.  -8.0/9.0, 2.0/9.0, 2.0/9.0, -6.0/9.0, 2.0/9.0, 2.0/9.0, -4.0/9.0, 2.0/9.0, 2.0/9.0,
  238.  -2.0/9.0, 2.0/9.0, 2.0/9.0,      0.0, 2.0/9.0, 2.0/9.0,  2.0/9.0, 2.0/9.0, 2.0/9.0,
  239.   4.0/9.0, 2.0/9.0, 2.0/9.0,  6.0/9.0, 2.0/9.0, 2.0/9.0,  8.0/9.0, 2.0/9.0, 2.0/9.0,
  240.  -8.0/9.0, 4.0/9.0, 2.0/9.0, -6.0/9.0, 4.0/9.0, 2.0/9.0, -4.0/9.0, 4.0/9.0, 2.0/9.0,
  241.  -2.0/9.0, 4.0/9.0, 2.0/9.0,      0.0, 4.0/9.0, 2.0/9.0,  2.0/9.0, 4.0/9.0, 2.0/9.0,
  242.   4.0/9.0, 4.0/9.0, 2.0/9.0,  6.0/9.0, 4.0/9.0, 2.0/9.0,  8.0/9.0, 4.0/9.0, 2.0/9.0,
  243.  -8.0/9.0, 6.0/9.0, 2.0/9.0, -6.0/9.0, 6.0/9.0, 2.0/9.0, -4.0/9.0, 6.0/9.0, 2.0/9.0,
  244.  -2.0/9.0, 6.0/9.0, 2.0/9.0,      0.0, 6.0/9.0, 2.0/9.0,  2.0/9.0, 6.0/9.0, 2.0/9.0,
  245.   4.0/9.0, 6.0/9.0, 2.0/9.0,  6.0/9.0, 6.0/9.0, 2.0/9.0,  8.0/9.0, 6.0/9.0, 2.0/9.0,
  246.  -8.0/9.0, 8.0/9.0, 2.0/9.0, -6.0/9.0, 8.0/9.0, 2.0/9.0, -4.0/9.0, 8.0/9.0, 2.0/9.0,
  247.  -2.0/9.0, 8.0/9.0, 2.0/9.0,      0.0, 8.0/9.0, 2.0/9.0,  2.0/9.0, 8.0/9.0, 2.0/9.0,
  248.   4.0/9.0, 8.0/9.0, 2.0/9.0,  6.0/9.0, 8.0/9.0, 2.0/9.0,  8.0/9.0, 8.0/9.0, 2.0/9.0,
  249.  -8.0/9.0,-8.0/9.0, 4.0/9.0, -6.0/9.0,-8.0/9.0, 4.0/9.0, -4.0/9.0,-8.0/9.0, 4.0/9.0,
  250.  -2.0/9.0,-8.0/9.0, 4.0/9.0,      0.0,-8.0/9.0, 4.0/9.0,  2.0/9.0,-8.0/9.0, 4.0/9.0,
  251.   4.0/9.0,-8.0/9.0, 4.0/9.0,  6.0/9.0,-8.0/9.0, 4.0/9.0,  8.0/9.0,-8.0/9.0, 4.0/9.0,
  252.  -8.0/9.0,-6.0/9.0, 4.0/9.0, -6.0/9.0,-6.0/9.0, 4.0/9.0, -4.0/9.0,-6.0/9.0, 4.0/9.0,
  253.  -2.0/9.0,-6.0/9.0, 4.0/9.0,      0.0,-6.0/9.0, 4.0/9.0,  2.0/9.0,-6.0/9.0, 4.0/9.0,
  254.   4.0/9.0,-6.0/9.0, 4.0/9.0,  6.0/9.0,-6.0/9.0, 4.0/9.0,  8.0/9.0,-6.0/9.0, 4.0/9.0,
  255.  -8.0/9.0,-4.0/9.0, 4.0/9.0, -6.0/9.0,-4.0/9.0, 4.0/9.0, -4.0/9.0,-4.0/9.0, 4.0/9.0,
  256.  -2.0/9.0,-4.0/9.0, 4.0/9.0,      0.0,-4.0/9.0, 4.0/9.0,  2.0/9.0,-4.0/9.0, 4.0/9.0,
  257.   4.0/9.0,-4.0/9.0, 4.0/9.0,  6.0/9.0,-4.0/9.0, 4.0/9.0,  8.0/9.0,-4.0/9.0, 4.0/9.0,
  258.  -8.0/9.0,-2.0/9.0, 4.0/9.0, -6.0/9.0,-2.0/9.0, 4.0/9.0, -4.0/9.0,-2.0/9.0, 4.0/9.0,
  259.  -2.0/9.0,-2.0/9.0, 4.0/9.0,      0.0,-2.0/9.0, 4.0/9.0,  2.0/9.0,-2.0/9.0, 4.0/9.0,
  260.   4.0/9.0,-2.0/9.0, 4.0/9.0,  6.0/9.0,-2.0/9.0, 4.0/9.0,  8.0/9.0,-2.0/9.0, 4.0/9.0,
  261.  -8.0/9.0,     0.0, 4.0/9.0, -6.0/9.0,     0.0, 4.0/9.0, -4.0/9.0,     0.0, 4.0/9.0,
  262.  -2.0/9.0,     0.0, 4.0/9.0,      0.0,     0.0, 4.0/9.0,  2.0/9.0,     0.0, 4.0/9.0,
  263.   4.0/9.0,     0.0, 4.0/9.0,  6.0/9.0,     0.0, 4.0/9.0,  8.0/9.0,     0.0, 4.0/9.0,
  264.  -8.0/9.0, 2.0/9.0, 4.0/9.0, -6.0/9.0, 2.0/9.0, 4.0/9.0, -4.0/9.0, 2.0/9.0, 4.0/9.0,
  265.  -2.0/9.0, 2.0/9.0, 4.0/9.0,      0.0, 2.0/9.0, 4.0/9.0,  2.0/9.0, 2.0/9.0, 4.0/9.0,
  266.   4.0/9.0, 2.0/9.0, 4.0/9.0,  6.0/9.0, 2.0/9.0, 4.0/9.0,  8.0/9.0, 2.0/9.0, 4.0/9.0,
  267.  -8.0/9.0, 4.0/9.0, 4.0/9.0, -6.0/9.0, 4.0/9.0, 4.0/9.0, -4.0/9.0, 4.0/9.0, 4.0/9.0,
  268.  -2.0/9.0, 4.0/9.0, 4.0/9.0,      0.0, 4.0/9.0, 4.0/9.0,  2.0/9.0, 4.0/9.0, 4.0/9.0,
  269.   4.0/9.0, 4.0/9.0, 4.0/9.0,  6.0/9.0, 4.0/9.0, 4.0/9.0,  8.0/9.0, 4.0/9.0, 4.0/9.0,
  270.  -8.0/9.0, 6.0/9.0, 4.0/9.0, -6.0/9.0, 6.0/9.0, 4.0/9.0, -4.0/9.0, 6.0/9.0, 4.0/9.0,
  271.  -2.0/9.0, 6.0/9.0, 4.0/9.0,      0.0, 6.0/9.0, 4.0/9.0,  2.0/9.0, 6.0/9.0, 4.0/9.0,
  272.   4.0/9.0, 6.0/9.0, 4.0/9.0,  6.0/9.0, 6.0/9.0, 4.0/9.0,  8.0/9.0, 6.0/9.0, 4.0/9.0,
  273.  -8.0/9.0, 8.0/9.0, 4.0/9.0, -6.0/9.0, 8.0/9.0, 4.0/9.0, -4.0/9.0, 8.0/9.0, 4.0/9.0,
  274.  -2.0/9.0, 8.0/9.0, 4.0/9.0,      0.0, 8.0/9.0, 4.0/9.0,  2.0/9.0, 8.0/9.0, 4.0/9.0,
  275.   4.0/9.0, 8.0/9.0, 4.0/9.0,  6.0/9.0, 8.0/9.0, 4.0/9.0,  8.0/9.0, 8.0/9.0, 4.0/9.0,
  276.  -8.0/9.0,-8.0/9.0, 6.0/9.0, -6.0/9.0,-8.0/9.0, 6.0/9.0, -4.0/9.0,-8.0/9.0, 6.0/9.0,
  277.  -2.0/9.0,-8.0/9.0, 6.0/9.0,      0.0,-8.0/9.0, 6.0/9.0,  2.0/9.0,-8.0/9.0, 6.0/9.0,
  278.   4.0/9.0,-8.0/9.0, 6.0/9.0,  6.0/9.0,-8.0/9.0, 6.0/9.0,  8.0/9.0,-8.0/9.0, 6.0/9.0,
  279.  -8.0/9.0,-6.0/9.0, 6.0/9.0, -6.0/9.0,-6.0/9.0, 6.0/9.0, -4.0/9.0,-6.0/9.0, 6.0/9.0,
  280.  -2.0/9.0,-6.0/9.0, 6.0/9.0,      0.0,-6.0/9.0, 6.0/9.0,  2.0/9.0,-6.0/9.0, 6.0/9.0,
  281.   4.0/9.0,-6.0/9.0, 6.0/9.0,  6.0/9.0,-6.0/9.0, 6.0/9.0,  8.0/9.0,-6.0/9.0, 6.0/9.0,
  282.  -8.0/9.0,-4.0/9.0, 6.0/9.0, -6.0/9.0,-4.0/9.0, 6.0/9.0, -4.0/9.0,-4.0/9.0, 6.0/9.0,
  283.  -2.0/9.0,-4.0/9.0, 6.0/9.0,      0.0,-4.0/9.0, 6.0/9.0,  2.0/9.0,-4.0/9.0, 6.0/9.0,
  284.   4.0/9.0,-4.0/9.0, 6.0/9.0,  6.0/9.0,-4.0/9.0, 6.0/9.0,  8.0/9.0,-4.0/9.0, 6.0/9.0,
  285.  -8.0/9.0,-2.0/9.0, 6.0/9.0, -6.0/9.0,-2.0/9.0, 6.0/9.0, -4.0/9.0,-2.0/9.0, 6.0/9.0,
  286.  -2.0/9.0,-2.0/9.0, 6.0/9.0,      0.0,-2.0/9.0, 6.0/9.0,  2.0/9.0,-2.0/9.0, 6.0/9.0,
  287.   4.0/9.0,-2.0/9.0, 6.0/9.0,  6.0/9.0,-2.0/9.0, 6.0/9.0,  8.0/9.0,-2.0/9.0, 6.0/9.0,
  288.  -8.0/9.0,     0.0, 6.0/9.0, -6.0/9.0,     0.0, 6.0/9.0, -4.0/9.0,     0.0, 6.0/9.0,
  289.  -2.0/9.0,     0.0, 6.0/9.0,      0.0,     0.0, 6.0/9.0,  2.0/9.0,     0.0, 6.0/9.0,
  290.   4.0/9.0,     0.0, 6.0/9.0,  6.0/9.0,     0.0, 6.0/9.0,  8.0/9.0,     0.0, 6.0/9.0,
  291.  -8.0/9.0, 2.0/9.0, 6.0/9.0, -6.0/9.0, 2.0/9.0, 6.0/9.0, -4.0/9.0, 2.0/9.0, 6.0/9.0,
  292.  -2.0/9.0, 2.0/9.0, 6.0/9.0,      0.0, 2.0/9.0, 6.0/9.0,  2.0/9.0, 2.0/9.0, 6.0/9.0,
  293.   4.0/9.0, 2.0/9.0, 6.0/9.0,  6.0/9.0, 2.0/9.0, 6.0/9.0,  8.0/9.0, 2.0/9.0, 6.0/9.0,
  294.  -8.0/9.0, 4.0/9.0, 6.0/9.0, -6.0/9.0, 4.0/9.0, 6.0/9.0, -4.0/9.0, 4.0/9.0, 6.0/9.0,
  295.  -2.0/9.0, 4.0/9.0, 6.0/9.0,      0.0, 4.0/9.0, 6.0/9.0,  2.0/9.0, 4.0/9.0, 6.0/9.0,
  296.   4.0/9.0, 4.0/9.0, 6.0/9.0,  6.0/9.0, 4.0/9.0, 6.0/9.0,  8.0/9.0, 4.0/9.0, 6.0/9.0,
  297.  -8.0/9.0, 6.0/9.0, 6.0/9.0, -6.0/9.0, 6.0/9.0, 6.0/9.0, -4.0/9.0, 6.0/9.0, 6.0/9.0,
  298.  -2.0/9.0, 6.0/9.0, 6.0/9.0,      0.0, 6.0/9.0, 6.0/9.0,  2.0/9.0, 6.0/9.0, 6.0/9.0,
  299.   4.0/9.0, 6.0/9.0, 6.0/9.0,  6.0/9.0, 6.0/9.0, 6.0/9.0,  8.0/9.0, 6.0/9.0, 6.0/9.0,
  300.  -8.0/9.0, 8.0/9.0, 6.0/9.0, -6.0/9.0, 8.0/9.0, 6.0/9.0, -4.0/9.0, 8.0/9.0, 6.0/9.0,
  301.  -2.0/9.0, 8.0/9.0, 6.0/9.0,      0.0, 8.0/9.0, 6.0/9.0,  2.0/9.0, 8.0/9.0, 6.0/9.0,
  302.   4.0/9.0, 8.0/9.0, 6.0/9.0,  6.0/9.0, 8.0/9.0, 6.0/9.0,  8.0/9.0, 8.0/9.0, 6.0/9.0,
  303.  -8.0/9.0,-8.0/9.0, 8.0/9.0, -6.0/9.0,-8.0/9.0, 8.0/9.0, -4.0/9.0,-8.0/9.0, 8.0/9.0,
  304.  -2.0/9.0,-8.0/9.0, 8.0/9.0,      0.0,-8.0/9.0, 8.0/9.0,  2.0/9.0,-8.0/9.0, 8.0/9.0,
  305.   4.0/9.0,-8.0/9.0, 8.0/9.0,  6.0/9.0,-8.0/9.0, 8.0/9.0,  8.0/9.0,-8.0/9.0, 8.0/9.0,
  306.  -8.0/9.0,-6.0/9.0, 8.0/9.0, -6.0/9.0,-6.0/9.0, 8.0/9.0, -4.0/9.0,-6.0/9.0, 8.0/9.0,
  307.  -2.0/9.0,-6.0/9.0, 8.0/9.0,      0.0,-6.0/9.0, 8.0/9.0,  2.0/9.0,-6.0/9.0, 8.0/9.0,
  308.   4.0/9.0,-6.0/9.0, 8.0/9.0,  6.0/9.0,-6.0/9.0, 8.0/9.0,  8.0/9.0,-6.0/9.0, 8.0/9.0,
  309.  -8.0/9.0,-4.0/9.0, 8.0/9.0, -6.0/9.0,-4.0/9.0, 8.0/9.0, -4.0/9.0,-4.0/9.0, 8.0/9.0,
  310.  -2.0/9.0,-4.0/9.0, 8.0/9.0,      0.0,-4.0/9.0, 8.0/9.0,  2.0/9.0,-4.0/9.0, 8.0/9.0,
  311.   4.0/9.0,-4.0/9.0, 8.0/9.0,  6.0/9.0,-4.0/9.0, 8.0/9.0,  8.0/9.0,-4.0/9.0, 8.0/9.0,
  312.  -8.0/9.0,-2.0/9.0, 8.0/9.0, -6.0/9.0,-2.0/9.0, 8.0/9.0, -4.0/9.0,-2.0/9.0, 8.0/9.0,
  313.  -2.0/9.0,-2.0/9.0, 8.0/9.0,      0.0,-2.0/9.0, 8.0/9.0,  2.0/9.0,-2.0/9.0, 8.0/9.0,
  314.   4.0/9.0,-2.0/9.0, 8.0/9.0,  6.0/9.0,-2.0/9.0, 8.0/9.0,  8.0/9.0,-2.0/9.0, 8.0/9.0,
  315.  -8.0/9.0,     0.0, 8.0/9.0, -6.0/9.0,     0.0, 8.0/9.0, -4.0/9.0,     0.0, 8.0/9.0,
  316.  -2.0/9.0,     0.0, 8.0/9.0,      0.0,     0.0, 8.0/9.0,  2.0/9.0,     0.0, 8.0/9.0,
  317.   4.0/9.0,     0.0, 8.0/9.0,  6.0/9.0,     0.0, 8.0/9.0,  8.0/9.0,     0.0, 8.0/9.0,
  318.  -8.0/9.0, 2.0/9.0, 8.0/9.0, -6.0/9.0, 2.0/9.0, 8.0/9.0, -4.0/9.0, 2.0/9.0, 8.0/9.0,
  319.  -2.0/9.0, 2.0/9.0, 8.0/9.0,      0.0, 2.0/9.0, 8.0/9.0,  2.0/9.0, 2.0/9.0, 8.0/9.0,
  320.   4.0/9.0, 2.0/9.0, 8.0/9.0,  6.0/9.0, 2.0/9.0, 8.0/9.0,  8.0/9.0, 2.0/9.0, 8.0/9.0,
  321.  -8.0/9.0, 4.0/9.0, 8.0/9.0, -6.0/9.0, 4.0/9.0, 8.0/9.0, -4.0/9.0, 4.0/9.0, 8.0/9.0,
  322.  -2.0/9.0, 4.0/9.0, 8.0/9.0,      0.0, 4.0/9.0, 8.0/9.0,  2.0/9.0, 4.0/9.0, 8.0/9.0,
  323.   4.0/9.0, 4.0/9.0, 8.0/9.0,  6.0/9.0, 4.0/9.0, 8.0/9.0,  8.0/9.0, 4.0/9.0, 8.0/9.0,
  324.  -8.0/9.0, 6.0/9.0, 8.0/9.0, -6.0/9.0, 6.0/9.0, 8.0/9.0, -4.0/9.0, 6.0/9.0, 8.0/9.0,
  325.  -2.0/9.0, 6.0/9.0, 8.0/9.0,      0.0, 6.0/9.0, 8.0/9.0,  2.0/9.0, 6.0/9.0, 8.0/9.0,
  326.   4.0/9.0, 6.0/9.0, 8.0/9.0,  6.0/9.0, 6.0/9.0, 8.0/9.0,  8.0/9.0, 6.0/9.0, 8.0/9.0,
  327.  -8.0/9.0, 8.0/9.0, 8.0/9.0, -6.0/9.0, 8.0/9.0, 8.0/9.0, -4.0/9.0, 8.0/9.0, 8.0/9.0,
  328.  -2.0/9.0, 8.0/9.0, 8.0/9.0,      0.0, 8.0/9.0, 8.0/9.0,  2.0/9.0, 8.0/9.0, 8.0/9.0,
  329.   4.0/9.0, 8.0/9.0, 8.0/9.0,  6.0/9.0, 8.0/9.0, 8.0/9.0,  8.0/9.0, 8.0/9.0, 8.0/9.0
  330. };
  333. static const REAL *grouptableA[16] =
  334. { 0,group5bits,group7bits,group10bits,0,0,0,0,0,0,0,0,0,0,0,0};
  335. static const REAL *grouptableB1[16] =
  336. { 0,group5bits,0,0,0,0,0,0,0,0,0,0,0,0,0,0 };
  337. static const REAL *grouptableB234[16] =
  338. { 0,group5bits,group7bits,0,group10bits,0,0,0,0,0,0,0,0,0,0,0};
  340. static const int codelengthtableA[16] =
  341. { 0, 5, 7, 10, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15 };
  342. static const int codelengthtableB1[16] =
  343. { 0, 5, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16 };
  344. static const int codelengthtableB2[16] =
  345. { 0, 5, 7, 3, 10, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 16 };
  346. static const int codelengthtableB3[8] = { 0, 5, 7, 3, 10, 4, 5, 16 };
  347. static const int codelengthtableB4[4] = { 0, 5, 7, 16 };
  350. static const REAL factortableA[16] =
  351. { 0.0,        1.0/2.0,    1.0/4.0,    1.0/8.0,
  352.   1.0/8.0,    1.0/16.0,   1.0/32.0,   1.0/64.0,
  353.   1.0/128.0,  1.0/256.0,  1.0/512.0,  1.0/1024.0,
  354.   1.0/2048.0, 1.0/4096.0, 1.0/8192.0, 1.0/16384.0 };
  355. static const REAL factortableB1[16] =
  356. { 0.0,        1.0/2.0,    1.0/4.0,     1.0/8.0, 
  357.   1.0/16.0,   1.0/32.0,   1.0/64.0,    1.0/128.0, 
  358.   1.0/256.0,  1.0/512.0,  1.0/1024.0,  1.0/2048.0,
  359.   1.0/4096.0, 1.0/8192.0, 1.0/16384.0, 1.0/32768.0 };
  360. static const REAL factortableB2[16] =
  361. { 0.0,        1.0/2.0,    1.0/4.0,    1.0/4.0,
  362.   1.0/8.0,    1.0/8.0,    1.0/16.0,   1.0/32.0,
  363.   1.0/64.0,   1.0/128.0,  1.0/256.0,  1.0/512.0,
  364.   1.0/1024.0, 1.0/2048.0, 1.0/4096.0, 1.0/32768.0 };
  365. static const REAL factortableB3[8] =
  366. { 0.0, 1.0/2.0, 1.0/4.0, 1.0/4.0, 1.0/8.0, 1.0/8.0, 1.0/16.0, 1.0/32768.0 };
  367. static const REAL factortableB4[4] = { 0.0, 1.0/2.0, 1.0/4.0, 1.0/32768.0 };
  370. static const REAL ctableA[16]=
  371. {           0.0, 1.33333333333, 1.60000000000, 1.77777777777,
  372.   1.06666666666, 1.03225806452, 1.01587301587, 1.00787401575,
  373.   1.00392156863, 1.00195694716, 1.00097751711, 1.00048851979,
  374.   1.00024420024, 1.00012208522, 1.00006103888, 1.00003051851};
  375. static const REAL ctableB1[16]=
  376. {           0.0, 1.33333333333, 1.14285714286, 1.06666666666,
  377.   1.03225806452, 1.01587301587, 1.00787401575, 1.00392156863,
  378.   1.00195694716, 1.00097751711, 1.00048851979, 1.00024420024,
  379.   1.00012208522, 1.00006103888, 1.00003051851, 1.00001525902};
  380. static const REAL ctableB2[16] =
  381. {           0.0, 1.33333333333, 1.60000000000, 1.14285714286,
  382.   1.77777777777, 1.06666666666, 1.03225806452, 1.01587301587,
  383.   1.00787401575, 1.00392156863, 1.00195694716, 1.00097751711,
  384.   1.00048851979, 1.00024420024, 1.00012208522, 1.00001525902};
  385. static const REAL ctableB3[8] =
  386. {           0.0, 1.33333333333, 1.60000000000, 1.14285714286,
  387.   1.77777777777, 1.06666666666, 1.03225806452, 1.00001525902 };
  388. static const REAL ctableB4[4] = 
  389. { 0.0, 1.33333333333, 1.60000000000, 1.00001525902 };
  392. static const REAL dtableA[16]=
  393. {           0.0, 0.50000000000, 0.50000000000, 0.50000000000,
  394.   0.12500000000, 0.06250000000, 0.03125000000, 0.01562500000,
  395.   0.00781250000, 0.00390625000, 0.00195312500, 0.00097656250,
  396.   0.00048828125, 0.00024414063, 0.00012207031, 0.00006103516};
  398. static const REAL dtableB1[16]=
  399. {           0.0, 0.50000000000, 0.25000000000, 0.12500000000,
  400.   0.06250000000, 0.03125000000, 0.01562500000, 0.00781250000,
  401.   0.00390625000, 0.00195312500, 0.00097656250, 0.00048828125,
  402.   0.00024414063, 0.00012207031, 0.00006103516, 0.00003051758};
  404. static const REAL dtableB2[16]=
  405. { 0.0,           0.50000000000, 0.50000000000, 0.25000000000,
  406.   0.50000000000, 0.12500000000, 0.06250000000, 0.03125000000,
  407.   0.01562500000, 0.00781250000, 0.00390625000, 0.00195312500,
  408.   0.00097656250, 0.00048828125, 0.00024414063, 0.00003051758};
  410. static const REAL dtableB3[8]=
  411. {           0.0, 0.50000000000, 0.50000000000, 0.25000000000,
  412.   0.50000000000, 0.12500000000, 0.06250000000, 0.00003051758};
  414. static const REAL dtableB4[4]=
  415. {0.0, 0.50000000000, 0.50000000000, 0.00003051758};
  418. // Mpeg layer 2
  419. void MPEGaudio::extractlayer2(void)
  420. {
  421.   REAL fraction[MAXCHANNEL][3][MAXSUBBAND];
  422.   unsigned int bitalloc[MAXCHANNEL][MAXSUBBAND],
  423.                scaleselector[MAXCHANNEL][MAXSUBBAND];
  424.   REAL scalefactor[2][3][MAXSUBBAND];
  426.   const REAL *group[MAXCHANNEL][MAXSUBBAND];
  427.   unsigned int codelength[MAXCHANNEL][MAXSUBBAND];
  428.   REAL factor[MAXCHANNEL][MAXSUBBAND];
  429.   REAL c[MAXCHANNEL][MAXSUBBAND],d[MAXCHANNEL][MAXSUBBAND];
  431.   int s=stereobound,n=subbandnumber;
  434. // Bitalloc
  435.   {
  436.     register int i;
  437.     register const int *t=bitalloclengthtable[tableindex];
  439.     for(i=0;i<s;i++,t++)
  440.     {
  441.       bitalloc[LS][i]=getbits(*t);
  442.       bitalloc[RS][i]=getbits(*t);
  443.     }
  444.     for(;i<n;i++,t++)
  445.       bitalloc[LS][i]=bitalloc[RS][i]=getbits(*t);
  446.   }
  449.   // Scale selector
  450.   if(inputstereo)
  451.     for(register int i=0;i<n;i++)
  452.     {
  453.       if(bitalloc[LS][i])scaleselector[LS][i]=getbits(2);
  454.       if(bitalloc[RS][i])scaleselector[RS][i]=getbits(2);
  455.     }
  456.   else
  457.     for(register int i=0;i<n;i++)
  458.       if(bitalloc[LS][i])scaleselector[LS][i]=getbits(2);
  460.   // Scale index
  461.   {
  462.     register int i,j;
  464.     for(i=0;i<n;i++)
  465.     {
  466.       if((j=bitalloc[LS][i]))
  467.       {
  468. if(!tableindex)
  469. {
  470.   group[LS][i]=grouptableA[j];
  471.   codelength[LS][i]=codelengthtableA[j];
  472.   factor[LS][i]=factortableA[j];
  473.   c[LS][i]=ctableA[j];
  474.   d[LS][i]=dtableA[j];
  475. }
  476. else
  477. {
  478.   if(i<=2)
  479.   {
  480.     group[LS][i]=grouptableB1[j];
  481.     codelength[LS][i]=codelengthtableB1[j];
  482.     factor[LS][i]=factortableB1[j];
  483.     c[LS][i]=ctableB1[j];
  484.     d[LS][i]=dtableB1[j];
  485.   }
  486.   else
  487.   {
  488.     group[LS][i]=grouptableB234[j];
  489.     if(i<=10)
  490.     {
  491.       codelength[LS][i]=codelengthtableB2[j];
  492.       factor[LS][i]=factortableB2[j];
  493.       c[LS][i]=ctableB2[j];
  494.       d[LS][i]=dtableB2[j];
  495.     }
  496.     else if(i<=22)
  497.     {
  498.       codelength[LS][i]=codelengthtableB3[j];
  499.       factor[LS][i]=factortableB3[j];
  500.       c[LS][i]=ctableB3[j];
  501.       d[LS][i]=dtableB3[j];
  502.     }
  503.     else
  504.     {
  505.       codelength[LS][i]=codelengthtableB4[j];
  506.       factor[LS][i]=factortableB4[j];
  507.       c[LS][i]=ctableB4[j];
  508.       d[LS][i]=dtableB4[j];
  509.     }
  510.   }
  511. }
  513. switch(scaleselector[LS][i])
  514. {
  515.   case 0:scalefactor[LS][0][i]=scalefactorstable[getbits(6)];
  516.  scalefactor[LS][1][i]=scalefactorstable[getbits(6)];
  517.  scalefactor[LS][2][i]=scalefactorstable[getbits(6)];
  518.  break;
  519.   case 1:scalefactor[LS][0][i]=
  520.  scalefactor[LS][1][i]=scalefactorstable[getbits(6)];
  521.  scalefactor[LS][2][i]=scalefactorstable[getbits(6)];
  522.  break;
  523.   case 2:scalefactor[LS][0][i]=
  524.  scalefactor[LS][1][i]=
  525.  scalefactor[LS][2][i]=scalefactorstable[getbits(6)];
  526.  break;
  527.   case 3:scalefactor[LS][0][i]=scalefactorstable[getbits(6)];
  528.  scalefactor[LS][1][i]=
  529.  scalefactor[LS][2][i]=scalefactorstable[getbits(6)];
  530.  break;
  531. }
  532.       }
  534.       if(inputstereo && (j=bitalloc[RS][i]))
  535.       {
  536. if(!tableindex)
  537. {
  538.   group[RS][i]=grouptableA[j];
  539.   codelength[RS][i]=codelengthtableA[j];
  540.   factor[RS][i]=factortableA[j];
  541.   c[RS][i]=ctableA[j];
  542.   d[RS][i]=dtableA[j];
  543. }
  544. else
  545. {
  546.   if(i<=2)
  547.   {
  548.     group[RS][i]=grouptableB1[j];
  549.     codelength[RS][i]=codelengthtableB1[j];
  550.     factor[RS][i]=factortableB1[j];
  551.     c[RS][i]=ctableB1[j];
  552.     d[RS][i]=dtableB1[j];
  553.   }
  554.   else
  555.   {
  556.     group[RS][i]=grouptableB234[j];
  557.     if(i<=10)
  558.     {
  559.       codelength[RS][i]=codelengthtableB2[j];
  560.       factor[RS][i]=factortableB2[j];
  561.       c[RS][i]=ctableB2[j];
  562.       d[RS][i]=dtableB2[j];
  563.     }
  564.     else if(i<=22)
  565.     {
  566.       codelength[RS][i]=codelengthtableB3[j];
  567.       factor[RS][i]=factortableB3[j];
  568.       c[RS][i]=ctableB3[j];
  569.       d[RS][i]=dtableB3[j];
  570.     }
  571.     else
  572.     {
  573.       codelength[RS][i]=codelengthtableB4[j];
  574.       factor[RS][i]=factortableB4[j];
  575.       c[RS][i]=ctableB4[j];
  576.       d[RS][i]=dtableB4[j];
  577.     }
  578.   }
  579. }
  581. switch(scaleselector[RS][i])
  582. {
  583.   case 0 : scalefactor[RS][0][i]=scalefactorstable[getbits(6)];
  584.    scalefactor[RS][1][i]=scalefactorstable[getbits(6)];
  585.    scalefactor[RS][2][i]=scalefactorstable[getbits(6)];
  586.    break;
  587.   case 1 : scalefactor[RS][0][i]=
  588.    scalefactor[RS][1][i]=scalefactorstable[getbits(6)];
  589.    scalefactor[RS][2][i]=scalefactorstable[getbits(6)];
  590.    break;
  591.   case 2 : scalefactor[RS][0][i]=
  592.    scalefactor[RS][1][i]=
  593.    scalefactor[RS][2][i]=scalefactorstable[getbits(6)];
  594.    break;
  595.   case 3 : scalefactor[RS][0][i]=scalefactorstable[getbits(6)];
  596.    scalefactor[RS][1][i]=
  597.    scalefactor[RS][2][i]=scalefactorstable[getbits(6)];
  598.    break;
  599. }
  600.       }
  601.     }
  602.   }
  605. // Read Sample
  606.   {
  607.     register int i;
  609.     for(int l=0;l<SCALEBLOCK;l++)
  610.     {
  611.       // Read Sample
  612.       for(i=0;i<s;i++)
  613.       {
  614. if(bitalloc[LS][i])
  615. {
  616.   if(group[LS][i])
  617.   {
  618.     register const REAL *s;
  619.     int code=getbits(codelength[LS][i]);
  621.     code+=code<<1;
  622.             if (code > 2184) {
  623. //printf("fraction LS OverFlow code %d -> 2184 (1)n", code);
  624.               code=2184;
  625.             }
  626.     s=group[LS][i]+code;
  628.     fraction[LS][0][i]=s[0];
  629.     fraction[LS][1][i]=s[1];
  630.     fraction[LS][2][i]=s[2];
  631.   }
  632.   else
  633.   {
  634.     fraction[LS][0][i]=
  635.       REAL(getbits(codelength[LS][i]))*factor[LS][i]-1.0;
  636.     fraction[LS][1][i]=
  637.       REAL(getbits(codelength[LS][i]))*factor[LS][i]-1.0;
  638.     fraction[LS][2][i]=
  639.       REAL(getbits(codelength[LS][i]))*factor[LS][i]-1.0;
  640.   }
  641. }
  642. else fraction[LS][0][i]=fraction[LS][1][i]=fraction[LS][2][i]=0.0;
  644. if(inputstereo && bitalloc[RS][i])
  645. {
  646.   if(group[RS][i])
  647.   {
  648.     const REAL *s;
  649.     int code=getbits(codelength[RS][i]);
  651.     code+=code<<1;
  652.             if (code > 2184) {
  653. //printf("fraction LS OverFlow code %d -> 2184 (2)n", code);
  654.               code=2184;
  655.             }
  656.     s=group[RS][i]+code;
  658.     fraction[RS][0][i]=s[0];
  659.     fraction[RS][1][i]=s[1];
  660.     fraction[RS][2][i]=s[2];
  661.   }
  662.   else
  663.   {
  664.     fraction[RS][0][i]=
  665.       REAL(getbits(codelength[RS][i]))*factor[RS][i]-1.0;
  666.     fraction[RS][1][i]=
  667.       REAL(getbits(codelength[RS][i]))*factor[RS][i]-1.0;
  668.     fraction[RS][2][i]=
  669.       REAL(getbits(codelength[RS][i]))*factor[RS][i]-1.0;
  670.   }
  671. }
  672. else fraction[RS][0][i]=fraction[RS][1][i]=fraction[RS][2][i]=0.0;
  673.       }
  675.       for(;i<n;i++)
  676.       {
  677. if(bitalloc[LS][i])
  678. {
  679.   if(group[LS][i])
  680.   {
  681.     register const REAL *s;
  682.     int code=getbits(codelength[LS][i]);
  684.     code+=code<<1;
  685.     s=group[LS][i]+code;
  687.     fraction[LS][0][i]=fraction[RS][0][i]=s[0];
  688.     fraction[LS][1][i]=fraction[RS][1][i]=s[1];
  689.     fraction[LS][2][i]=fraction[RS][2][i]=s[2];
  690.   }
  691.   else
  692.   {
  693.     fraction[LS][0][i]=fraction[RS][0][i]=
  694.       REAL(getbits(codelength[LS][i]))*factor[LS][i]-1.0;
  695.     fraction[LS][1][i]=fraction[RS][1][i]=
  696.       REAL(getbits(codelength[LS][i]))*factor[LS][i]-1.0;
  697.     fraction[LS][2][i]=fraction[RS][2][i]=
  698.       REAL(getbits(codelength[LS][i]))*factor[LS][i]-1.0;
  699.   }
  700. }
  701. else fraction[LS][0][i]=fraction[LS][1][i]=fraction[LS][2][i]=
  702.      fraction[RS][0][i]=fraction[RS][1][i]=fraction[RS][2][i]=0.0;
  703.       }
  707.       //Fraction
  708.       if(outputstereo)
  709. for(i=0;i<n;i++)
  710. {
  711.   if(bitalloc[LS][i])
  712.   {
  713.     if(!group[LS][i])
  714.     {
  715.       fraction[LS][0][i]=(fraction[LS][0][i]+d[LS][i])*c[LS][i];
  716.       fraction[LS][1][i]=(fraction[LS][1][i]+d[LS][i])*c[LS][i];
  717.       fraction[LS][2][i]=(fraction[LS][2][i]+d[LS][i])*c[LS][i];
  718.     }
  720.     register REAL t=scalefactor[LS][l>>2][i];
  721.     fraction[LS][0][i]*=t;
  722.     fraction[LS][1][i]*=t;
  723.     fraction[LS][2][i]*=t;
  724.   }
  726.   if(bitalloc[RS][i])
  727.   {
  728.     if(!group[RS][i])
  729.     {
  730.       fraction[RS][0][i]=(fraction[RS][0][i]+d[RS][i])*c[RS][i];
  731.       fraction[RS][1][i]=(fraction[RS][1][i]+d[RS][i])*c[RS][i];
  732.       fraction[RS][2][i]=(fraction[RS][2][i]+d[RS][i])*c[RS][i];
  733.     }
  735.     register REAL t=scalefactor[RS][l>>2][i];
  736.     fraction[RS][0][i]*=t;
  737.     fraction[RS][1][i]*=t;
  738.     fraction[RS][2][i]*=t;
  739.   }
  740. }
  741.       else
  742. for(i=0;i<n;i++)
  743.   if(bitalloc[LS][i])
  744.   {
  745.     if(!group[LS][i])
  746.     {
  747.       fraction[LS][0][i]=(fraction[LS][0][i]+d[LS][i])*c[LS][i];
  748.       fraction[LS][1][i]=(fraction[LS][1][i]+d[LS][i])*c[LS][i];
  749.       fraction[LS][2][i]=(fraction[LS][2][i]+d[LS][i])*c[LS][i];
  750.     }
  752.     register REAL t=scalefactor[LS][l>>2][i];
  753.     fraction[LS][0][i]*=t;
  754.     fraction[LS][1][i]*=t;
  755.     fraction[LS][2][i]*=t;
  756.   }
  758.       for(;i<MAXSUBBAND;i++)
  759. fraction[LS][0][i]=fraction[LS][1][i]=fraction[LS][2][i]=
  760. fraction[RS][0][i]=fraction[RS][1][i]=fraction[RS][2][i]=0.0;
  762.       for(i=0;i<3;i++)
  763. subbandsynthesis(fraction[LS][i],fraction[RS][i]);
  764.     }
  765.   }
  766. }