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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > DVD > 查看源码
synthesis_filter.cc
资源名称：dvd-munitions.tar.gz　[点击查看]
上传用户：aoeyumen
上传日期：2007-01-06
资源大小：3329k
文件大小：28k
源码类别：

DVD

开发平台：

Unix_Linux

synthesis_filter.cc：源码内容
  1. /*
  2.    File: synthesis_filter.cc
  3. */
  5. #include <stdio.h>
  6. #include <stdlib.h>
  7. #include <fstream.h>
  8. #ifdef IRIX
  9. #include <dmedia/audio.h>
  10. #endif
  11. #ifdef SOLARIS
  12. #include <sys/audioio.h>
  13. #endif
  15. #include "athread.hh"
  17. #include "error.hh"
  18. #include "debug.hh"
  19. #include "util.hh"
  21. #include "sync.hh"
  22. #include "mpeg2const.hh"
  23. #include "mpeg2buff.hh"
  25. #include "astream.hh"
  26. #include "crc.hh"
  27. #include "header.hh"
  28. #include "obuffer.hh"
  29. #include "synthesis_filter.hh"
  32. const real SynthesisFilter::n[32][31] = {
  33.   {  0.707106781,  0.671558955,  0.634393284,  0.595699304,
  34.      0.555570233,  0.514102744,  0.471396737,  0.427555093,
  35.      0.382683432,  0.336889853,  0.290284677,  0.242980180,
  36.      0.195090322,  0.146730474,  0.098017140,  0.049067674,
  37.     -0.740951126, -0.773010454, -0.803207532, -0.831469612,
  38.     -0.857728610, -0.881921265, -0.903989293, -0.923879533,
  39.     -0.941544065, -0.956940336, -0.970031253, -0.980785280,
  40.     -0.989176510, -0.995184727, -0.998795456 },
  41.   { -0.707106781, -0.803207532, -0.881921265, -0.941544065,
  42.     -0.980785280, -0.998795456, -0.995184727, -0.970031253,
  43.     -0.923879532, -0.857728610, -0.773010453, -0.671558954,
  44.     -0.555570233, -0.427555093, -0.290284677, -0.146730474,
  45.      0.595699304,  0.471396737,  0.336889853,  0.195090322,
  46.      0.049067674, -0.098017141, -0.242980180, -0.382683433,
  47.     -0.514102744, -0.634393284, -0.740951126, -0.831469612,
  48.     -0.903989293, -0.956940336, -0.989176510 },
  49.   { -0.707106781, -0.514102744, -0.290284677, -0.049067674,
  50.      0.195090323,  0.427555094,  0.634393285,  0.803207532,
  51.      0.923879533,  0.989176510,  0.995184727,  0.941544065,
  52.      0.831469612,  0.671558955,  0.471396737,  0.242980180,
  53.      0.857728610,  0.956940336,  0.998795456,  0.980785280,
  54.      0.903989293,  0.773010453,  0.595699304,  0.382683432,
  55.      0.146730474, -0.098017141, -0.336889854, -0.555570233,
  56.     -0.740951126, -0.881921265, -0.970031253 },
  57.   {  0.707106782,  0.903989293,  0.995184727,  0.970031253,
  58.      0.831469612,  0.595699304,  0.290284677, -0.049067675,
  59.     -0.382683433, -0.671558955, -0.881921265, -0.989176510,
  60.     -0.980785280, -0.857728610, -0.634393284, -0.336889853,
  61.     -0.427555093, -0.098017140,  0.242980181,  0.555570234,
  62.      0.803207532,  0.956940336,  0.998795456,  0.923879532,
  63.      0.740951125,  0.471396737,  0.146730474, -0.195090322,
  64.     -0.514102744, -0.773010454, -0.941544065 },
  65.   {  0.707106781,  0.336889853, -0.098017141, -0.514102744,
  66.     -0.831469612, -0.989176510, -0.956940336, -0.740951125,
  67.     -0.382683432,  0.049067675,  0.471396737,  0.803207532,
  68.      0.980785281,  0.970031253,  0.773010453,  0.427555093,
  69.     -0.941544065, -0.995184727, -0.857728610, -0.555570233,
  70.     -0.146730474,  0.290284678,  0.671558955,  0.923879533,
  71.      0.998795456,  0.881921264,  0.595699304,  0.195090322,
  72.     -0.242980180, -0.634393284, -0.903989293 },
  73.   { -0.707106781, -0.970031253, -0.956940336, -0.671558954,
  74.     -0.195090321,  0.336889854,  0.773010454,  0.989176510,
  75.      0.923879532,  0.595699304,  0.098017140, -0.427555094,
  76.     -0.831469612, -0.998795456, -0.881921264, -0.514102744,
  77.      0.242980180, -0.290284677, -0.740951126, -0.980785280,
  78.     -0.941544065, -0.634393284, -0.146730474,  0.382683433,
  79.      0.803207532,  0.995184727,  0.903989293,  0.555570233,
  80.      0.049067674, -0.471396737, -0.857728610 },
  81.   { -0.707106781, -0.146730474,  0.471396737,  0.903989293,
  82.      0.980785280,  0.671558955,  0.098017140, -0.514102744,
  83.     -0.923879533, -0.970031253, -0.634393284, -0.049067674,
  84.      0.555570234,  0.941544065,  0.956940336,  0.595699304,
  85.      0.989176510,  0.881921264,  0.427555093, -0.195090322,
  86.     -0.740951126, -0.995184727, -0.857728610, -0.382683432,
  87.      0.242980181,  0.773010454,  0.998795456,  0.831469612,
  88.      0.336889853, -0.290284677, -0.803207532 },
  89.   {  0.707106782,  0.998795456,  0.773010453,  0.146730474,
  90.     -0.555570233, -0.970031253, -0.881921264, -0.336889853,
  91.      0.382683433,  0.903989293,  0.956940336,  0.514102744,
  92.     -0.195090322, -0.803207532, -0.995184727, -0.671558954,
  93.     -0.049067674,  0.634393285,  0.989176510,  0.831469612,
  94.      0.242980180, -0.471396737, -0.941544065, -0.923879532,
  95.     -0.427555093,  0.290284678,  0.857728610,  0.980785280,
  96.      0.595699304, -0.098017141, -0.740951126 },
  97.   {  0.707106781, -0.049067675, -0.773010454, -0.989176510,
  98.     -0.555570233,  0.242980181,  0.881921265,  0.941544065,
  99.      0.382683432, -0.427555094, -0.956940336, -0.857728610,
  100.     -0.195090321,  0.595699305,  0.995184727,  0.740951125,
  101.     -0.998795456, -0.634393284,  0.146730475,  0.831469613,
  102.      0.970031253,  0.471396737, -0.336889854, -0.923879533,
  103.     -0.903989293, -0.290284677,  0.514102745,  0.980785281,
  104.      0.803207531,  0.098017140, -0.671558955 },
  105.   { -0.707106781, -0.989176510, -0.471396736,  0.427555094,
  106.      0.980785281,  0.740951125, -0.098017141, -0.857728610,
  107.     -0.923879532, -0.242980179,  0.634393285,  0.998795456,
  108.      0.555570233, -0.336889854, -0.956940336, -0.803207531,
  109.     -0.146730475, -0.881921265, -0.903989293, -0.195090321,
  110.      0.671558955,  0.995184727,  0.514102744, -0.382683433,
  111.     -0.970031253, -0.773010453,  0.049067675,  0.831469613,
  112.      0.941544065,  0.290284677, -0.595699305 },
  113.   { -0.707106781,  0.242980181,  0.956940336,  0.740951125,
  114.     -0.195090322, -0.941544065, -0.773010453,  0.146730475,
  115.      0.923879533,  0.803207531, -0.098017141, -0.903989293,
  116.     -0.831469612,  0.049067675,  0.881921265,  0.857728610,
  117.      0.970031253,  0.290284677, -0.671558955, -0.980785280,
  118.     -0.336889853,  0.634393285,  0.989176510,  0.382683432,
  119.     -0.595699305, -0.995184727, -0.427555093,  0.555570234,
  120.      0.998795456,  0.471396737, -0.514102744 },
  121.   {  0.707106782,  0.941544065,  0.098017140, -0.857728610,
  122.     -0.831469612,  0.146730475,  0.956940336,  0.671558955,
  123.     -0.382683433, -0.998795456, -0.471396736,  0.595699305,
  124.      0.980785280,  0.242980180, -0.773010454, -0.903989293,
  125.      0.336889854,  0.995184727,  0.514102744, -0.555570233,
  126.     -0.989176510, -0.290284677,  0.740951126,  0.923879532,
  127.      0.049067674, -0.881921265, -0.803207531,  0.195090323,
  128.      0.970031253,  0.634393284, -0.427555094 },
  129.   {  0.707106781, -0.427555094, -0.995184727, -0.242980179,
  130.      0.831469613,  0.803207531, -0.290284677, -0.998795456,
  131.     -0.382683432,  0.740951126,  0.881921264, -0.146730475,
  132.     -0.980785280, -0.514102744,  0.634393285,  0.941544065,
  133.     -0.903989293,  0.098017141,  0.970031253,  0.555570233,
  134.     -0.595699305, -0.956940336, -0.049067674,  0.923879533,
  135.      0.671558955, -0.471396737, -0.989176510, -0.195090321,
  136.      0.857728610,  0.773010453, -0.336889854 },
  137.   { -0.707106781, -0.857728610,  0.290284678,  0.998795456,
  138.      0.195090322, -0.903989293, -0.634393284,  0.595699305,
  139.      0.923879532, -0.146730475, -0.995184727, -0.336889853,
  140.      0.831469613,  0.740951125, -0.471396737, -0.970031253,
  141.     -0.514102744, -0.956940336,  0.049067675,  0.980785281,
  142.      0.427555093, -0.773010454, -0.803207531,  0.382683433,
  143.      0.989176510,  0.098017140, -0.941544065, -0.555570233,
  144.      0.671558955,  0.881921264, -0.242980180 },
  145.   { -0.707106781,  0.595699305,  0.881921264, -0.336889854,
  146.     -0.980785280,  0.049067675,  0.995184727,  0.242980180,
  147.     -0.923879533, -0.514102744,  0.773010454,  0.740951125,
  148.     -0.555570233, -0.903989293,  0.290284678,  0.989176510,
  149.      0.803207531, -0.471396737, -0.941544065,  0.195090323,
  150.      0.998795456,  0.098017140, -0.970031253, -0.382683432,
  151.      0.857728610,  0.634393284, -0.671558955, -0.831469612,
  152.      0.427555094,  0.956940336, -0.146730475 },
  153.   {  0.707106782,  0.740951125, -0.634393284, -0.803207531,
  154.      0.555570234,  0.857728610, -0.471396737, -0.903989293,
  155.      0.382683433,  0.941544065, -0.290284677, -0.970031253,
  156.      0.195090323,  0.989176510, -0.098017141, -0.998795456,
  157.      0.671558955,  0.773010453, -0.595699305, -0.831469612,
  158.      0.514102745,  0.881921264, -0.427555094, -0.923879532,
  159.      0.336889854,  0.956940336, -0.242980180, -0.980785280,
  160.      0.146730475,  0.995184727, -0.049067675 },
  161.   {  0.707106781, -0.740951126, -0.634393284,  0.803207532,
  162.      0.555570233, -0.857728610, -0.471396736,  0.903989293,
  163.      0.382683432, -0.941544065, -0.290284677,  0.970031253,
  164.      0.195090322, -0.989176510, -0.098017140,  0.998795456,
  165.     -0.671558954,  0.773010454,  0.595699304, -0.831469612,
  166.     -0.514102744,  0.881921265,  0.427555093, -0.923879533,
  167.     -0.336889853,  0.956940336,  0.242980180, -0.980785280,
  168.     -0.146730474,  0.995184727,  0.049067674 },
  169.   { -0.707106781, -0.595699304,  0.881921265,  0.336889853,
  170.     -0.980785280, -0.049067674,  0.995184727, -0.242980180,
  171.     -0.923879532,  0.514102745,  0.773010453, -0.740951126,
  172.     -0.555570233,  0.903989293,  0.290284677, -0.989176510,
  173.     -0.803207532, -0.471396736,  0.941544065,  0.195090322,
  174.     -0.998795456,  0.098017141,  0.970031253, -0.382683433,
  175.     -0.857728610,  0.634393285,  0.671558955, -0.831469612,
  176.     -0.427555093,  0.956940336,  0.146730474 },
  177.   { -0.707106781,  0.857728610,  0.290284677, -0.998795456,
  178.      0.195090323,  0.903989293, -0.634393284, -0.595699304,
  179.      0.923879533,  0.146730474, -0.995184727,  0.336889854,
  180.      0.831469612, -0.740951126, -0.471396736,  0.970031253,
  181.      0.514102744, -0.956940336, -0.049067674,  0.980785280,
  182.     -0.427555094, -0.773010453,  0.803207532,  0.382683432,
  183.     -0.989176510,  0.098017141,  0.941544065, -0.555570233,
  184.     -0.671558954,  0.881921265,  0.242980180 },
  185.   {  0.707106782,  0.427555093, -0.995184727,  0.242980181,
  186.      0.831469612, -0.803207532, -0.290284677,  0.998795456,
  187.     -0.382683433, -0.740951125,  0.881921265,  0.146730474,
  188.     -0.980785280,  0.514102745,  0.634393284, -0.941544065,
  189.      0.903989293,  0.098017140, -0.970031253,  0.555570234,
  190.      0.595699304, -0.956940336,  0.049067675,  0.923879532,
  191.     -0.671558955, -0.471396736,  0.989176510, -0.195090322,
  192.     -0.857728610,  0.773010454,  0.336889853 },
  193.   {  0.707106781, -0.941544065,  0.098017141,  0.857728610,
  194.     -0.831469612, -0.146730474,  0.956940336, -0.671558955,
  195.     -0.382683432,  0.998795456, -0.471396737, -0.595699304,
  196.      0.980785281, -0.242980180, -0.773010453,  0.903989293,
  197.     -0.336889853,  0.995184727, -0.514102744, -0.555570233,
  198.      0.989176510, -0.290284677, -0.740951125,  0.923879533,
  199.     -0.049067675, -0.881921264,  0.803207532,  0.195090322,
  200.     -0.970031253,  0.634393285,  0.427555093 },
  201.   { -0.707106781, -0.242980179,  0.956940336, -0.740951126,
  202.     -0.195090321,  0.941544065, -0.773010454, -0.146730474,
  203.      0.923879532, -0.803207532, -0.098017140,  0.903989293,
  204.     -0.831469612, -0.049067674,  0.881921264, -0.857728610,
  205.     -0.970031253,  0.290284678,  0.671558955, -0.980785280,
  206.      0.336889854,  0.634393284, -0.989176510,  0.382683433,
  207.      0.595699304, -0.995184727,  0.427555094,  0.555570233,
  208.     -0.998795456,  0.471396737,  0.514102744 },
  209.   { -0.707106781,  0.989176510, -0.471396737, -0.427555093,
  210.      0.980785280, -0.740951126, -0.098017140,  0.857728610,
  211.     -0.923879533,  0.242980181,  0.634393284, -0.998795456,
  212.      0.555570234,  0.336889853, -0.956940336,  0.803207532,
  213.      0.146730474, -0.881921264,  0.903989293, -0.195090322,
  214.     -0.671558954,  0.995184727, -0.514102744, -0.382683432,
  215.      0.970031253, -0.773010454, -0.049067674,  0.831469612,
  216.     -0.941544065,  0.290284678,  0.595699304 },
  217.   {  0.707106782,  0.049067674, -0.773010453,  0.989176510,
  218.     -0.555570233, -0.242980179,  0.881921264, -0.941544065,
  219.      0.382683433,  0.427555093, -0.956940336,  0.857728610,
  220.     -0.195090322, -0.595699304,  0.995184727, -0.740951126,
  221.      0.998795456, -0.634393284, -0.146730474,  0.831469612,
  222.     -0.970031253,  0.471396737,  0.336889853, -0.923879532,
  223.      0.903989293, -0.290284677, -0.514102744,  0.980785280,
  224.     -0.803207532,  0.098017141,  0.671558955 },
  225.   {  0.707106781, -0.998795456,  0.773010454, -0.146730475,
  226.     -0.555570233,  0.970031253, -0.881921265,  0.336889854,
  227.      0.382683432, -0.903989293,  0.956940336, -0.514102744,
  228.     -0.195090321,  0.803207531, -0.995184727,  0.671558955,
  229.      0.049067675,  0.634393284, -0.989176510,  0.831469613,
  230.     -0.242980180, -0.471396736,  0.941544065, -0.923879533,
  231.      0.427555094,  0.290284677, -0.857728610,  0.980785281,
  232.     -0.595699305, -0.098017140,  0.740951125 },
  233.   { -0.707106781,  0.146730475,  0.471396737, -0.903989293,
  234.      0.980785281, -0.671558955,  0.098017141,  0.514102744,
  235.     -0.923879532,  0.970031253, -0.634393284,  0.049067675,
  236.      0.555570233, -0.941544065,  0.956940336, -0.595699305,
  237.     -0.989176510,  0.881921265, -0.427555094, -0.195090321,
  238.      0.740951125, -0.995184727,  0.857728610, -0.382683433,
  239.     -0.242980179,  0.773010453, -0.998795456,  0.831469613,
  240.     -0.336889854, -0.290284677,  0.803207531 },
  241.   { -0.707106781,  0.970031253, -0.956940336,  0.671558955,
  242.     -0.195090322, -0.336889853,  0.773010453, -0.989176510,
  243.      0.923879533, -0.595699305,  0.098017141,  0.427555093,
  244.     -0.831469612,  0.998795456, -0.881921265,  0.514102745,
  245.     -0.242980180, -0.290284677,  0.740951125, -0.980785280,
  246.      0.941544065, -0.634393284,  0.146730475,  0.382683432,
  247.     -0.803207531,  0.995184727, -0.903989293,  0.555570234,
  248.     -0.049067675, -0.471396736,  0.857728610 },
  249.   {  0.707106782, -0.336889854, -0.098017140,  0.514102744,
  250.     -0.831469612,  0.989176510, -0.956940336,  0.740951126,
  251.     -0.382683433, -0.049067674,  0.471396737, -0.803207531,
  252.      0.980785280, -0.970031253,  0.773010454, -0.427555094,
  253.      0.941544065, -0.995184727,  0.857728610, -0.555570233,
  254.      0.146730475,  0.290284677, -0.671558954,  0.923879532,
  255.     -0.998795456,  0.881921265, -0.595699305,  0.195090323,
  256.      0.242980180, -0.634393284,  0.903989293 },
  257.   {  0.707106781, -0.903989293,  0.995184727, -0.970031253,
  258.      0.831469613, -0.595699305,  0.290284678,  0.049067674,
  259.     -0.382683432,  0.671558955, -0.881921264,  0.989176510,
  260.     -0.980785280,  0.857728610, -0.634393284,  0.336889854,
  261.      0.427555094, -0.098017141, -0.242980179,  0.555570233,
  262.     -0.803207531,  0.956940336, -0.998795456,  0.923879533,
  263.     -0.740951126,  0.471396737, -0.146730475, -0.195090321,
  264.      0.514102744, -0.773010453,  0.941544065 },
  265.   { -0.707106781,  0.514102745, -0.290284677,  0.049067675,
  266.      0.195090322, -0.427555093,  0.634393284, -0.803207531,
  267.      0.923879532, -0.989176510,  0.995184727, -0.941544065,
  268.      0.831469613, -0.671558955,  0.471396737, -0.242980180,
  269.     -0.857728610,  0.956940336, -0.998795456,  0.980785281,
  270.     -0.903989293,  0.773010454, -0.595699305,  0.382683433,
  271.     -0.146730475, -0.098017140,  0.336889853, -0.555570233,
  272.      0.740951125, -0.881921264,  0.970031253 },
  273.   { -0.707106781,  0.803207531, -0.881921264,  0.941544065,
  274.     -0.980785280,  0.998795456, -0.995184727,  0.970031253,
  275.     -0.923879533,  0.857728610, -0.773010454,  0.671558955,
  276.     -0.555570233,  0.427555094, -0.290284677,  0.146730475,
  277.     -0.595699305,  0.471396737, -0.336889854,  0.195090323,
  278.     -0.049067675, -0.098017140,  0.242980180, -0.382683432,
  279.      0.514102744, -0.634393284,  0.740951125, -0.831469612,
  280.      0.903989293, -0.956940336,  0.989176510 },
  281.   {  0.707106782, -0.671558955,  0.634393285, -0.595699305,
  282.      0.555570234, -0.514102744,  0.471396737, -0.427555094,
  283.      0.382683433, -0.336889854,  0.290284678, -0.242980180,
  284.      0.195090323, -0.146730475,  0.098017141, -0.049067675,
  285.      0.740951125, -0.773010453,  0.803207531, -0.831469612,
  286.      0.857728610, -0.881921264,  0.903989293, -0.923879532,
  287.      0.941544065, -0.956940336,  0.970031253, -0.980785280,
  288.      0.989176510, -0.995184727,  0.998795456 }
  289. };
  292. const real SynthesisFilter::d[512] = {
  293.   // Warning: These values are not in the same order
  294.   // as in Annex 3-B.3 of the ISO/IEC DIS 11172-3
  295.    0.000000000, -0.000442505,  0.003250122, -0.007003784,
  296.    0.031082153, -0.078628540,  0.100311279, -0.572036743,
  297.    1.144989014,  0.572036743,  0.100311279,  0.078628540,
  298.    0.031082153,  0.007003784,  0.003250122,  0.000442505,
  299.   -0.000015259, -0.000473022,  0.003326416, -0.007919312,
  300.    0.030517578, -0.084182739,  0.090927124, -0.600219727,
  301.    1.144287109,  0.543823242,  0.108856201,  0.073059082,
  302.    0.031478882,  0.006118774,  0.003173828,  0.000396729,
  303.   -0.000015259, -0.000534058,  0.003387451, -0.008865356,
  304.    0.029785156, -0.089706421,  0.080688477, -0.628295898,
  305.    1.142211914,  0.515609741,  0.116577148,  0.067520142,
  306.    0.031738281,  0.005294800,  0.003082275,  0.000366211,
  307.   -0.000015259, -0.000579834,  0.003433228, -0.009841919,
  308.    0.028884888, -0.095169067,  0.069595337, -0.656219482,
  309.    1.138763428,  0.487472534,  0.123474121,  0.061996460,
  310.    0.031845093,  0.004486084,  0.002990723,  0.000320435,
  311.   -0.000015259, -0.000625610,  0.003463745, -0.010848999,
  312.    0.027801514, -0.100540161,  0.057617188, -0.683914185,
  313.    1.133926392,  0.459472656,  0.129577637,  0.056533813,
  314.    0.031814575,  0.003723145,  0.002899170,  0.000289917,
  315.   -0.000015259, -0.000686646,  0.003479004, -0.011886597,
  316.    0.026535034, -0.105819702,  0.044784546, -0.711318970,
  317.    1.127746582,  0.431655884,  0.134887695,  0.051132202,
  318.    0.031661987,  0.003005981,  0.002792358,  0.000259399,
  319.   -0.000015259, -0.000747681,  0.003479004, -0.012939453,
  320.    0.025085449, -0.110946655,  0.031082153, -0.738372803,
  321.    1.120223999,  0.404083252,  0.139450073,  0.045837402,
  322.    0.031387329,  0.002334595,  0.002685547,  0.000244141,
  323.   -0.000030518, -0.000808716,  0.003463745, -0.014022827,
  324.    0.023422241, -0.115921021,  0.016510010, -0.765029907,
  325.    1.111373901,  0.376800537,  0.143264771,  0.040634155,
  326.    0.031005859,  0.001693726,  0.002578735,  0.000213623,
  327.   -0.000030518, -0.000885010,  0.003417969, -0.015121460,
  328.    0.021575928, -0.120697021,  0.001068115, -0.791213989,
  329.    1.101211548,  0.349868774,  0.146362305,  0.035552979,
  330.    0.030532837,  0.001098633,  0.002456665,  0.000198364,
  331.   -0.000030518, -0.000961304,  0.003372192, -0.016235352,
  332.    0.019531250, -0.125259399, -0.015228271, -0.816864014,
  333.    1.089782715,  0.323318481,  0.148773193,  0.030609131,
  334.    0.029937744,  0.000549316,  0.002349854,  0.000167847,
  335.   -0.000030518, -0.001037598,  0.003280640, -0.017349243,
  336.    0.017257690, -0.129562378, -0.032379150, -0.841949463,
  337.    1.077117920,  0.297210693,  0.150497437,  0.025817871,
  338.    0.029281616,  0.000030518,  0.002243042,  0.000152588,
  339.   -0.000045776, -0.001113892,  0.003173828, -0.018463135,
  340.    0.014801025, -0.133590698, -0.050354004, -0.866363525,
  341.    1.063217163,  0.271591187,  0.151596069,  0.021179199,
  342.    0.028533936, -0.000442505,  0.002120972,  0.000137329,
  343.   -0.000045776, -0.001205444,  0.003051758, -0.019577026,
  344.    0.012115479, -0.137298584, -0.069168091, -0.890090942,
  345.    1.048156738,  0.246505737,  0.152069092,  0.016708374,
  346.    0.027725220, -0.000869751,  0.002014160,  0.000122070,
  347.   -0.000061035, -0.001296997,  0.002883911, -0.020690918,
  348.    0.009231567, -0.140670776, -0.088775635, -0.913055420,
  349.    1.031936646,  0.221984863,  0.151962280,  0.012420654,
  350.    0.026840210, -0.001266479,  0.001907349,  0.000106812,
  351.   -0.000061035, -0.001388550,  0.002700806, -0.021789551,
  352.    0.006134033, -0.143676758, -0.109161377, -0.935195923,
  353.    1.014617920,  0.198059082,  0.151306152,  0.008316040,
  354.    0.025909424, -0.001617432,  0.001785278,  0.000106812,
  355.   -0.000076294, -0.001480103,  0.002487183, -0.022857666,
  356.    0.002822876, -0.146255493, -0.130310059, -0.956481934,
  357.    0.996246338,  0.174789429,  0.150115967,  0.004394531,
  358.    0.024932861, -0.001937866,  0.001693726,  0.000091553,
  359.   -0.000076294, -0.001586914,  0.002227783, -0.023910522,
  360.   -0.000686646, -0.148422241, -0.152206421, -0.976852417,
  361.    0.976852417,  0.152206421,  0.148422241,  0.000686646,
  362.    0.023910522, -0.002227783,  0.001586914,  0.000076294,
  363.   -0.000091553, -0.001693726,  0.001937866, -0.024932861,
  364.   -0.004394531, -0.150115967, -0.174789429, -0.996246338,
  365.    0.956481934,  0.130310059,  0.146255493, -0.002822876,
  366.    0.022857666, -0.002487183,  0.001480103,  0.000076294,
  367.   -0.000106812, -0.001785278,  0.001617432, -0.025909424,
  368.   -0.008316040, -0.151306152, -0.198059082, -1.014617920,
  369.    0.935195923,  0.109161377,  0.143676758, -0.006134033,
  370.    0.021789551, -0.002700806,  0.001388550,  0.000061035,
  371.   -0.000106812, -0.001907349,  0.001266479, -0.026840210,
  372.   -0.012420654, -0.151962280, -0.221984863, -1.031936646,
  373.    0.913055420,  0.088775635,  0.140670776, -0.009231567,
  374.    0.020690918, -0.002883911,  0.001296997,  0.000061035,
  375.   -0.000122070, -0.002014160,  0.000869751, -0.027725220,
  376.   -0.016708374, -0.152069092, -0.246505737, -1.048156738,
  377.    0.890090942,  0.069168091,  0.137298584, -0.012115479,
  378.    0.019577026, -0.003051758,  0.001205444,  0.000045776,
  379.   -0.000137329, -0.002120972,  0.000442505, -0.028533936,
  380.   -0.021179199, -0.151596069, -0.271591187, -1.063217163,
  381.    0.866363525,  0.050354004,  0.133590698, -0.014801025,
  382.    0.018463135, -0.003173828,  0.001113892,  0.000045776,
  383.   -0.000152588, -0.002243042, -0.000030518, -0.029281616,
  384.   -0.025817871, -0.150497437, -0.297210693, -1.077117920,
  385.    0.841949463,  0.032379150,  0.129562378, -0.017257690,
  386.    0.017349243, -0.003280640,  0.001037598,  0.000030518,
  387.   -0.000167847, -0.002349854, -0.000549316, -0.029937744,
  388.   -0.030609131, -0.148773193, -0.323318481, -1.089782715,
  389.    0.816864014,  0.015228271,  0.125259399, -0.019531250,
  390.    0.016235352, -0.003372192,  0.000961304,  0.000030518,
  391.   -0.000198364, -0.002456665, -0.001098633, -0.030532837,
  392.   -0.035552979, -0.146362305, -0.349868774, -1.101211548,
  393.    0.791213989, -0.001068115,  0.120697021, -0.021575928,
  394.    0.015121460, -0.003417969,  0.000885010,  0.000030518,
  395.   -0.000213623, -0.002578735, -0.001693726, -0.031005859,
  396.   -0.040634155, -0.143264771, -0.376800537, -1.111373901,
  397.    0.765029907, -0.016510010,  0.115921021, -0.023422241,
  398.    0.014022827, -0.003463745,  0.000808716,  0.000030518,
  399.   -0.000244141, -0.002685547, -0.002334595, -0.031387329,
  400.   -0.045837402, -0.139450073, -0.404083252, -1.120223999,
  401.    0.738372803, -0.031082153,  0.110946655, -0.025085449,
  402.    0.012939453, -0.003479004,  0.000747681,  0.000015259,
  403.   -0.000259399, -0.002792358, -0.003005981, -0.031661987,
  404.   -0.051132202, -0.134887695, -0.431655884, -1.127746582,
  405.    0.711318970, -0.044784546,  0.105819702, -0.026535034,
  406.    0.011886597, -0.003479004,  0.000686646,  0.000015259,
  407.   -0.000289917, -0.002899170, -0.003723145, -0.031814575,
  408.   -0.056533813, -0.129577637, -0.459472656, -1.133926392,
  409.    0.683914185, -0.057617188,  0.100540161, -0.027801514,
  410.    0.010848999, -0.003463745,  0.000625610,  0.000015259,
  411.   -0.000320435, -0.002990723, -0.004486084, -0.031845093,
  412.   -0.061996460, -0.123474121, -0.487472534, -1.138763428,
  413.    0.656219482, -0.069595337,  0.095169067, -0.028884888,
  414.    0.009841919, -0.003433228,  0.000579834,  0.000015259,
  415.   -0.000366211, -0.003082275, -0.005294800, -0.031738281,
  416.   -0.067520142, -0.116577148, -0.515609741, -1.142211914,
  417.    0.628295898, -0.080688477,  0.089706421, -0.029785156,
  418.    0.008865356, -0.003387451,  0.000534058,  0.000015259,
  419.   -0.000396729, -0.003173828, -0.006118774, -0.031478882,
  420.   -0.073059082, -0.108856201, -0.543823242, -1.144287109,
  421.    0.600219727, -0.090927124,  0.084182739, -0.030517578,
  422.    0.007919312, -0.003326416,  0.000473022,  0.000015259
  423. };
  427. SynthesisFilter::SynthesisFilter (uint32 channelnumber, real factor){
  428.   register real *floatp, *floatp2;
  430.   // initialize v1[] and v2[]:
  431.   for (floatp = v1 + 512, floatp2 = v2 + 512; floatp > v1; )
  432.     *--floatp = *--floatp2 = 0.0;
  434.   channel = channelnumber;
  435.   number_of_samples = range_violations = written_samples = 0;
  436.   actual_v = v1;
  437.   actual_write_pos = 15;
  438.   scalefactor = factor;
  439. }
  442. void SynthesisFilter::compute_new_v(real *new_v){
  443.   {
  444.     unsigned int i;
  445.     register real sample = 0.0, *v = new_v;
  446.     register const real *cooefs;
  448.     // initialize new_v[]:
  449.     *v++ = sample; *v++ = sample; *v++ = sample; *v++ = sample;
  450.     *v++ = sample; *v++ = sample; *v++ = sample; *v++ = sample;
  451.     *v++ = sample; *v++ = sample; *v++ = sample; *v++ = sample;
  452.     *v++ = sample; *v++ = sample; *v++ = sample; *v++ = sample;
  453.     v = new_v + 33;
  454.     *v++ = sample; *v++ = sample; *v++ = sample; *v++ = sample;
  455.     *v++ = sample; *v++ = sample; *v++ = sample; *v++ = sample;
  456.     *v++ = sample; *v++ = sample; *v++ = sample; *v++ = sample;
  457.     *v++ = sample; *v++ = sample; *v++ = sample; *v++ = sample;
  459.     // compute new values:
  460.     for (i = 0; i < number_of_samples; ++i){
  461.       sample = samples[i];
  462.       cooefs = n[subbandnumbers[i]];
  463.       v = new_v;
  464.       // add to new_v[0]-new_v[15]: (new_v[16] remains always 0.0)
  465.       *v++ += sample * *cooefs++; *v++ += sample * *cooefs++;
  466.       *v++ += sample * *cooefs++; *v++ += sample * *cooefs++;
  467.       *v++ += sample * *cooefs++; *v++ += sample * *cooefs++;
  468.       *v++ += sample * *cooefs++; *v++ += sample * *cooefs++;
  469.       *v++ += sample * *cooefs++; *v++ += sample * *cooefs++;
  470.       *v++ += sample * *cooefs++; *v++ += sample * *cooefs++;
  471.       *v++ += sample * *cooefs++; *v++ += sample * *cooefs++;
  472.       *v++ += sample * *cooefs++; *v++ += sample * *cooefs++;
  473.       v = new_v + 33;
  474.       // add to new_v[33]-new_v[47]:
  475.       *v++ += sample * *cooefs++; *v++ += sample * *cooefs++;
  476.       *v++ += sample * *cooefs++; *v++ += sample * *cooefs++;
  477.       *v++ += sample * *cooefs++; *v++ += sample * *cooefs++;
  478.       *v++ += sample * *cooefs++; *v++ += sample * *cooefs++;
  479.       *v++ += sample * *cooefs++; *v++ += sample * *cooefs++;
  480.       *v++ += sample * *cooefs++; *v++ += sample * *cooefs++;
  481.       *v++ += sample * *cooefs++; *v++ += sample * *cooefs++;
  482.       *v++ += sample * *cooefs++;
  483.       *v   -= sample; // new_v[48]
  484.     }
  485.   }
  487.   // fold new values to there boundaries: (a few streams need this)
  488.   {
  489.     register real *v;
  490.     for (v = new_v + 15; v >= new_v; --v)
  491.       if (*v < -1.0)
  492. *v = -1.0;
  493.       else if (*v > 1.0)
  494. *v = 1.0;
  495.     for (v = new_v + 48; v >= new_v + 33; --v)
  496.       if (*v < -1.0)
  497. *v = -1.0;
  498.       else if (*v > 1.0)
  499. *v = 1.0;
  500.   }
  502.   {
  503.     // copy values -new_v[0-15] to new_v[32-17]:
  504.     register real *v1 = new_v, *v2 = new_v + 33;
  505.     *--v2 = -*v1++; *--v2 = -*v1++; *--v2 = -*v1++; *--v2 = -*v1++;
  506.     *--v2 = -*v1++; *--v2 = -*v1++; *--v2 = -*v1++; *--v2 = -*v1++;
  507.     *--v2 = -*v1++; *--v2 = -*v1++; *--v2 = -*v1++; *--v2 = -*v1++;
  508.     *--v2 = -*v1++; *--v2 = -*v1++; *--v2 = -*v1++; *--v2 = -*v1++;
  510.     // copy values new_v[33-47] to new_v[63-49]:
  511.     v1 = new_v + 33;
  512.     v2 = new_v + 64;
  513.     *--v2 = *v1++; *--v2 = *v1++; *--v2 = *v1++; *--v2 = *v1++;
  514.     *--v2 = *v1++; *--v2 = *v1++; *--v2 = *v1++; *--v2 = *v1++;
  515.     *--v2 = *v1++; *--v2 = *v1++; *--v2 = *v1++; *--v2 = *v1++;
  516.     *--v2 = *v1++; *--v2 = *v1++; *--v2 = *v1++;
  517.   }
  519.   // insert new values into v1[] and v2[]:
  520.   {
  521.     register int i;
  522.     register real *nnv, *nv;
  524.     // insert new_v[0-31] into actual v:
  525.     nnv = new_v;
  526.     nv = actual_v + actual_write_pos;
  527.     for (i = 0; i < 32; ++i, nv += 16)
  528.       *nv = *nnv++;
  530.     // insert new_v[32-63] into other v:
  531.     nv = (actual_v == v1 ? v2 : v1) + actual_write_pos;
  532.     for (; i < 64; ++i, nv += 16)
  533.       *nv = *nnv++;
  534.   }
  535. }
  538. void SynthesisFilter::compute_pcm_samples (Obuffer *buffer){
  539.   int i;
  540.   register real floatreg, *lower_bound, *vp;
  541.   register const real *dp;
  542.   int pcm_sample;
  544.   dp = d;
  545.   lower_bound = actual_v;
  546.   if (actual_write_pos != 15)
  547.     for (i = 0; i < 32; ++i)
  548.     {
  549.       floatreg = 0;
  550.       for (vp = lower_bound + actual_write_pos + 1; vp > lower_bound; )
  551.         floatreg += *--vp * *dp++;
  552.       for (vp += 16; vp > lower_bound + actual_write_pos + 1; )
  553.         floatreg += *--vp * *dp++;
  555.       pcm_sample = (int)(floatreg * scalefactor);
  556.       if (pcm_sample > 32767)
  557.         {
  558.           ++range_violations;
  559.           if (floatreg > max_violation)
  560.             max_violation = floatreg;
  561.           pcm_sample = 32767;
  562.         }
  563.       else if (pcm_sample < -32768)
  564.         {
  565.           ++range_violations;
  566.           if (-floatreg > max_violation)
  567.             max_violation = -floatreg;
  568.           pcm_sample = -32768;
  569.         }
  570.       buffer->append (channel, (int16)pcm_sample);
  571.       lower_bound += 16;
  572.     }
  573.   else
  574.     for (i = 0; i < 32; ++i)
  575.       {
  576.         floatreg = 0;
  577.         for (vp = lower_bound + 16; vp > lower_bound; )
  578.           floatreg += *--vp * *dp++;
  579.         
  580.         pcm_sample = (int)(floatreg * scalefactor);
  581.         if (pcm_sample > 32767)
  582.           {
  583.             ++range_violations;
  584.             if (floatreg > max_violation)
  585.               max_violation = floatreg;
  586.             pcm_sample = 32767;
  587.           }
  588.         else if (pcm_sample < -32768)
  589.           {
  590.             ++range_violations;
  591.             if (-floatreg > max_violation)
  592.               max_violation = -floatreg;
  593.             pcm_sample = -32768;
  594.           }
  595.         buffer->append (channel, (int16)pcm_sample);
  596.         lower_bound += 16;
  597.       }
  598. }
  601. void SynthesisFilter::calculate_pcm_samples (Obuffer *buffer){
  602.   static real new_v[64];
  604.   compute_new_v (new_v);
  605.   compute_pcm_samples (buffer);
  607.   written_samples += 32;
  608.   if (actual_write_pos < 15)
  609.     ++actual_write_pos;
  610.   else
  611.     actual_write_pos = 0;
  612.   actual_v = (actual_v == v1 ? v2 : v1);
  613.   number_of_samples = 0;
  614. }