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

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

流媒体/Mpeg4/MP4

开发平台：

Visual C++

filter_2.cpp：源码内容
  1. /* MPEG/WAVE Sound library
  3.    (C) 1997 by Jung woo-jae */
  5. // Filter.cc
  6. // Subbandsynthesis routines from maplay 1.2 for Linux
  7. // I've modified some macros for reducing source code.
  9. #ifdef HAVE_CONFIG_H
  10. #include "config.h"
  11. #endif
  13. #include "MPEGaudio.h"
  15. void MPEGaudio::computebuffer_2(REAL *fraction,REAL buffer[2][CALCBUFFERSIZE])
  16. {
  17.   REAL p0,p1,p2,p3,p4,p5,p6,p7,p8,p9,pa,pb,pc,pd,pe,pf;
  18.   REAL q0,q1,q2,q3,q4,q5,q6,q7,q8,q9,qa,qb,qc,qd,qe,qf;
  19.   REAL *out1,*out2;
  21.   out1=buffer[currentcalcbuffer]+calcbufferoffset;
  22.   out2=buffer[currentcalcbuffer^1]+calcbufferoffset;
  23. #define OUT1(v,t) out1[(32-(v))*16]   =(-(out1[(v)*16]=t))
  24. #define OUT2(v)   out2[(96-(v)-32)*16]=out2[((v)-32)*16]
  26.   // compute new values via a fast cosine transform:
  27.   /*  {
  28.     register REAL *x=fraction;
  30.     p0=x[ 0]+x[31];p1=x[ 1]+x[30];p2=x[ 2]+x[29];p3=x[ 3]+x[28];
  31.     p4=x[ 4]+x[27];p5=x[ 5]+x[26];p6=x[ 6]+x[25];p7=x[ 7]+x[24];
  32.     p8=x[ 8]+x[23];p9=x[ 9]+x[22];pa=x[10]+x[21];pb=x[11]+x[20];
  33.     pc=x[12]+x[19];pd=x[13]+x[18];pe=x[14]+x[17];pf=x[15]+x[16];
  34.   }
  36.   q0=p0+pf;q1=p1+pe;q2=p2+pd;q3=p3+pc;
  37.   q4=p4+pb;q5=p5+pa;q6=p6+p9;q7=p7+p8;
  38.   q8=hcos_32[0]*(p0-pf);q9=hcos_32[1]*(p1-pe);
  39.   qa=hcos_32[2]*(p2-pd);qb=hcos_32[3]*(p3-pc);
  40.   qc=hcos_32[4]*(p4-pb);qd=hcos_32[5]*(p5-pa);
  41.   qe=hcos_32[6]*(p6-p9);qf=hcos_32[7]*(p7-p8); */
  43.   {
  44.     register REAL *x=fraction;
  46.     q0=x[ 0]+x[15];q1=x[ 1]+x[14];q2=x[ 2]+x[13];q3=x[ 3]+x[12];
  47.     q4=x[ 4]+x[11];q5=x[ 5]+x[10];q6=x[ 6]+x[ 9];q7=x[ 7]+x[ 8];
  49.     q8=hcos_32[0]*(x[ 0]-x[15]);q9=hcos_32[1]*(x[ 1]-x[14]);
  50.     qa=hcos_32[2]*(x[ 2]-x[13]);qb=hcos_32[3]*(x[ 3]-x[12]);
  51.     qc=hcos_32[4]*(x[ 4]-x[11]);qd=hcos_32[5]*(x[ 5]-x[10]);
  52.     qe=hcos_32[6]*(x[ 6]-x[ 9]);qf=hcos_32[7]*(x[ 7]-x[ 8]);
  53.   }
  55.   p0=q0+q7;p1=q1+q6;p2=q2+q5;p3=q3+q4;
  56.   p4=hcos_16[0]*(q0-q7);p5=hcos_16[1]*(q1-q6);
  57.   p6=hcos_16[2]*(q2-q5);p7=hcos_16[3]*(q3-q4);
  58.   p8=q8+qf;p9=q9+qe;pa=qa+qd;pb=qb+qc;
  59.   pc=hcos_16[0]*(q8-qf);pd=hcos_16[1]*(q9-qe);
  60.   pe=hcos_16[2]*(qa-qd);pf=hcos_16[3]*(qb-qc);
  62.   q0=p0+p3;q1=p1+p2;q2=hcos_8[0]*(p0-p3);q3=hcos_8[1]*(p1-p2);
  63.   q4=p4+p7;q5=p5+p6;q6=hcos_8[0]*(p4-p7);q7=hcos_8[1]*(p5-p6);
  64.   q8=p8+pb;q9=p9+pa;qa=hcos_8[0]*(p8-pb);qb=hcos_8[1]*(p9-pa);
  65.   qc=pc+pf;qd=pd+pe;qe=hcos_8[0]*(pc-pf);qf=hcos_8[1]*(pd-pe);
  67.   p0=q0+q1;p1=hcos_4*(q0-q1);p2=q2+q3;p3=hcos_4*(q2-q3);
  68.   p4=q4+q5;p5=hcos_4*(q4-q5);p6=q6+q7;p7=hcos_4*(q6-q7);
  69.   p8=q8+q9;p9=hcos_4*(q8-q9);pa=qa+qb;pb=hcos_4*(qa-qb);
  70.   pc=qc+qd;pd=hcos_4*(qc-qd);pe=qe+qf;pf=hcos_4*(qe-qf);
  72.   {
  73.     register REAL tmp;
  75.     tmp=p6+p7;
  76.     OUT2(36)=-(p5+tmp);
  77.     OUT2(44)=-(p4+tmp);
  78.     tmp=pb+pf;
  79.     OUT1(10,tmp);
  80.     OUT1(6,pd+tmp);
  81.     tmp=pe+pf;
  82.     OUT2(46)=-(p8+pc+tmp);
  83.     OUT2(34)=-(p9+pd+tmp);
  84.     tmp+=pa+pb;
  85.     OUT2(38)=-(pd+tmp);
  86.     OUT2(42)=-(pc+tmp);
  87.     OUT1(2,p9+pd+pf);
  88.     OUT1(4,p5+p7);
  89.     OUT2(48)=-p0;
  90.     out2[0]=-(out1[0]=p1);
  91.     OUT1( 8,p3);
  92.     OUT1(12,p7);
  93.     OUT1(14,pf);
  94.     OUT2(40)=-(p2+p3);
  95.   }
  97.   {
  98.     register REAL *x=fraction;
  100.     /*    p0=hcos_64[ 0]*(x[ 0]-x[31]);p1=hcos_64[ 1]*(x[ 1]-x[30]);
  101.     p2=hcos_64[ 2]*(x[ 2]-x[29]);p3=hcos_64[ 3]*(x[ 3]-x[28]);
  102.     p4=hcos_64[ 4]*(x[ 4]-x[27]);p5=hcos_64[ 5]*(x[ 5]-x[26]);
  103.     p6=hcos_64[ 6]*(x[ 6]-x[25]);p7=hcos_64[ 7]*(x[ 7]-x[24]);
  104.     p8=hcos_64[ 8]*(x[ 8]-x[23]);p9=hcos_64[ 9]*(x[ 9]-x[22]);
  105.     pa=hcos_64[10]*(x[10]-x[21]);pb=hcos_64[11]*(x[11]-x[20]);
  106.     pc=hcos_64[12]*(x[12]-x[19]);pd=hcos_64[13]*(x[13]-x[18]);
  107.     pe=hcos_64[14]*(x[14]-x[17]);pf=hcos_64[15]*(x[15]-x[16]); */
  109.     p0=hcos_64[ 0]*x[ 0];p1=hcos_64[ 1]*x[ 1];
  110.     p2=hcos_64[ 2]*x[ 2];p3=hcos_64[ 3]*x[ 3];
  111.     p4=hcos_64[ 4]*x[ 4];p5=hcos_64[ 5]*x[ 5];
  112.     p6=hcos_64[ 6]*x[ 6];p7=hcos_64[ 7]*x[ 7];
  113.     p8=hcos_64[ 8]*x[ 8];p9=hcos_64[ 9]*x[ 9];
  114.     pa=hcos_64[10]*x[10];pb=hcos_64[11]*x[11];
  115.     pc=hcos_64[12]*x[12];pd=hcos_64[13]*x[13];
  116.     pe=hcos_64[14]*x[14];pf=hcos_64[15]*x[15];
  117.   }
  119.   q0=p0+pf;q1=p1+pe;q2=p2+pd;q3=p3+pc;
  120.   q4=p4+pb;q5=p5+pa;q6=p6+p9;q7=p7+p8;
  121.   q8=hcos_32[0]*(p0-pf);q9=hcos_32[1]*(p1-pe);
  122.   qa=hcos_32[2]*(p2-pd);qb=hcos_32[3]*(p3-pc);
  123.   qc=hcos_32[4]*(p4-pb);qd=hcos_32[5]*(p5-pa);
  124.   qe=hcos_32[6]*(p6-p9);qf=hcos_32[7]*(p7-p8);
  126.   p0=q0+q7;p1=q1+q6;p2=q2+q5;p3=q3+q4;
  127.   p4=hcos_16[0]*(q0-q7);p5=hcos_16[1]*(q1-q6);
  128.   p6=hcos_16[2]*(q2-q5);p7=hcos_16[3]*(q3-q4);
  129.   p8=q8+qf;p9=q9+qe;pa=qa+qd;pb=qb+qc;
  130.   pc=hcos_16[0]*(q8-qf);pd=hcos_16[1]*(q9-qe);
  131.   pe=hcos_16[2]*(qa-qd);pf=hcos_16[3]*(qb-qc);
  133.   q0=p0+p3;q1=p1+p2;q2=hcos_8[0]*(p0-p3);q3=hcos_8[1]*(p1-p2);
  134.   q4=p4+p7;q5=p5+p6;q6=hcos_8[0]*(p4-p7);q7=hcos_8[1]*(p5-p6);
  135.   q8=p8+pb;q9=p9+pa;qa=hcos_8[0]*(p8-pb);qb=hcos_8[1]*(p9-pa);
  136.   qc=pc+pf;qd=pd+pe;qe=hcos_8[0]*(pc-pf);qf=hcos_8[1]*(pd-pe);
  138.   p0=q0+q1;p1=hcos_4*(q0-q1);
  139.   p2=q2+q3;p3=hcos_4*(q2-q3);
  140.   p4=q4+q5;p5=hcos_4*(q4-q5);
  141.   p6=q6+q7;p7=hcos_4*(q6-q7);
  142.   p8=q8+q9;p9=hcos_4*(q8-q9);
  143.   pa=qa+qb;pb=hcos_4*(qa-qb);
  144.   pc=qc+qd;pd=hcos_4*(qc-qd);
  145.   pe=qe+qf;pf=hcos_4*(qe-qf);
  147.   {
  148.     REAL tmp;
  150.     tmp=pd+pf;
  151.     OUT1(5,p5+p7+pb+tmp);
  152.     tmp+=p9;
  153.     OUT1(1,p1+tmp);
  154.     OUT2(33)=-(p1+pe+tmp);
  155.     tmp+=p5+p7;
  156.     OUT1(3,tmp);
  157.     OUT2(35)=-(p6+pe+tmp);
  158.     tmp=pa+pb+pc+pd+pe+pf;
  159.     OUT2(39)=-(p2+p3+tmp-pc);
  160.     OUT2(43)=-(p4+p6+p7+tmp-pd);
  161.     OUT2(37)=-(p5+p6+p7+tmp-pc);
  162.     OUT2(41)=-(p2+p3+tmp-pd);
  163.     tmp=p8+pc+pe+pf;
  164.     OUT2(47)=-(p0+tmp);
  165.     OUT2(45)=-(p4+p6+p7+tmp);
  166.     tmp=pb+pf;
  167.     OUT1(11,p7+tmp);
  168.     tmp+=p3;
  169.     OUT1( 9,tmp);
  170.     OUT1( 7,pd+tmp);
  171.     OUT1(13,p7+pf);
  172.     OUT1(15,pf);
  173.   }
  174. }
  177. #define SAVE 
  178.         raw=(int)(r*scalefactor); 
  179.         if(raw>MAXSCALE)raw=MAXSCALE;else if(raw<MINSCALE)raw=MINSCALE; 
  180. putraw(raw); 
  181.         dp+=16;vp+=15+(15-14)
  182. #define OS   r=*vp * *dp++
  183. #define XX   vp+=15;r+=*vp * *dp++
  184. #define OP   r+=*--vp * *dp++
  186. void MPEGaudio::generatesingle_2(void)
  187. {
  188.   int i;
  189.   register REAL r, *vp;
  190.   register const REAL *dp;
  191.   int raw;
  193.   i=32/2;
  194.   dp=filter;
  195.   vp=calcbufferL[currentcalcbuffer]+calcbufferoffset; 
  196. // actual_v+actual_write_pos;
  198.   switch (calcbufferoffset)
  199.   {
  200.     case  0:for(;i;i--,vp+=15){
  201.               OS;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  202.       SAVE;}break;
  203.     case  1:for(;i;i--,vp+=15){
  204.       OS;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  205.       SAVE;}break;
  206.     case  2:for(;i;i--,vp+=15){
  207.               OS;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  208.       SAVE;}break;
  209.     case  3:for(;i;i--,vp+=15){
  210.       OS;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  211.       SAVE;}break;
  212.     case  4:for(;i;i--,vp+=15){
  213.               OS;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  214.       SAVE;}break;
  215.     case  5:for(;i;i--,vp+=15){
  216.       OS;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  217.       SAVE;}break;
  218.     case  6:for(;i;i--,vp+=15){
  219.       OS;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;
  220.       SAVE;}break;
  221.     case  7:for(;i;i--,vp+=15){
  222.       OS;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;
  223.       SAVE;}break;
  224.     case  8:for(;i;i--,vp+=15){
  225.               OS;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;
  226.       SAVE;}break;
  227.     case  9:for(;i;i--,vp+=15){
  228.       OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;
  229.       SAVE;}break;
  230.     case 10:for(;i;i--,vp+=15){
  231.               OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;
  232.       SAVE;}break;
  233.     case 11:for(;i;i--,vp+=15){
  234.               OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;
  235.       SAVE;}break;
  236.     case 12:for(;i;i--,vp+=15){
  237.       OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;
  238.       SAVE;}break;
  239.     case 13:for(;i;i--,vp+=15){
  240.       OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;
  241.       SAVE;}break;
  242.     case 14:for(;i;i--,vp+=15){
  243.       OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;
  244.       SAVE;}break;
  245.     case 15:for(;i;i--,vp+=31){
  246.               OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  247.       SAVE;}break;
  248.   }
  249. }
  251. #undef OS
  252. #undef XX
  253. #undef OP
  254. #undef SAVE
  256. #define SAVE 
  257.         raw=(int)(r1*scalefactor);  
  258.         if(raw>MAXSCALE)raw=MAXSCALE;else if(raw<MINSCALE)raw=MINSCALE; 
  259. putraw(raw);  
  260.         raw=(int)(r2*scalefactor);  
  261.         if(raw>MAXSCALE)raw=MAXSCALE;else if(raw<MINSCALE)raw=MINSCALE; 
  262. putraw(raw); 
  263.         dp+=16;vp1+=15+(15-14);vp2+=15+(15-14)
  264. #define OS r1=*vp1 * *dp; 
  265.            r2=*vp2 * *dp++ 
  266. #define XX vp1+=15;r1+=*vp1 * *dp; 
  267.    vp2+=15;r2+=*vp2 * *dp++
  268. #define OP r1+=*--vp1 * *dp; 
  269.    r2+=*--vp2 * *dp++
  272. void MPEGaudio::generate_2(void)
  273. {
  274.   int i;
  275.   REAL r1,r2;
  276.   register REAL *vp1,*vp2;
  277.   register const REAL *dp;
  278.   int raw;
  280.   dp=filter;
  281.   vp1=calcbufferL[currentcalcbuffer]+calcbufferoffset;
  282.   vp2=calcbufferR[currentcalcbuffer]+calcbufferoffset;
  283. // actual_v+actual_write_pos;
  285.   i=32/2;
  286.   switch (calcbufferoffset)
  287.   {
  288.     case  0:for(;i;i--,vp1+=15,vp2+=15){
  289.               OS;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  290.       SAVE;}break;
  291.     case  1:for(;i;i--,vp1+=15,vp2+=15){
  292.       OS;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  293.       SAVE;}break;
  294.     case  2:for(;i;i--,vp1+=15,vp2+=15){
  295.               OS;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  296.       SAVE;}break;
  297.     case  3:for(;i;i--,vp1+=15,vp2+=15){
  298.       OS;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  299.       SAVE;}break;
  300.     case  4:for(;i;i--,vp1+=15,vp2+=15){
  301.               OS;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  302.       SAVE;}break;
  303.     case  5:for(;i;i--,vp1+=15,vp2+=15){
  304.       OS;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  305.       SAVE;}break;
  306.     case  6:for(;i;i--,vp1+=15,vp2+=15){
  307.       OS;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;
  308.       SAVE;}break;
  309.     case  7:for(;i;i--,vp1+=15,vp2+=15){
  310.       OS;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;
  311.       SAVE;}break;
  312.     case  8:for(;i;i--,vp1+=15,vp2+=15){
  313.               OS;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;
  314.       SAVE;}break;
  315.     case  9:for(;i;i--,vp1+=15,vp2+=15){
  316.       OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;
  317.       SAVE;}break;
  318.     case 10:for(;i;i--,vp1+=15,vp2+=15){
  319.               OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;
  320.       SAVE;}break;
  321.     case 11:for(;i;i--,vp1+=15,vp2+=15){
  322.               OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;
  323.       SAVE;}break;
  324.     case 12:for(;i;i--,vp1+=15,vp2+=15){
  325.       OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;
  326.       SAVE;}break;
  327.     case 13:for(;i;i--,vp1+=15,vp2+=15){
  328.       OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;
  329.       SAVE;}break;
  330.     case 14:for(;i;i--,vp1+=15,vp2+=15){
  331.       OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;
  332.       SAVE;}break;
  333.     case 15:for(;i;i--,vp1+=31,vp2+=31){
  334.               OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  335.       SAVE;}break;
  336.   }
  337. }
  340. void MPEGaudio::subbandsynthesis_2(REAL *fractionL,REAL *fractionR)
  341. {
  342.   computebuffer_2(fractionL,calcbufferL);
  343.   if(!outputstereo)generatesingle_2();
  344.   else
  345.   {
  346.     computebuffer_2(fractionR,calcbufferR);
  347.     generate_2();
  348.   }
  350.   if(calcbufferoffset<15)calcbufferoffset++;
  351.   else calcbufferoffset=0;
  353.   currentcalcbuffer^=1;
  354. }