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

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

流媒体/Mpeg4/MP4

开发平台：

Visual C++

filter.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(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.   }
  37.   q0=p0+pf;q1=p1+pe;q2=p2+pd;q3=p3+pc;
  38.   q4=p4+pb;q5=p5+pa;q6=p6+p9;q7=p7+p8;
  39.   q8=hcos_32[0]*(p0-pf);q9=hcos_32[1]*(p1-pe);
  40.   qa=hcos_32[2]*(p2-pd);qb=hcos_32[3]*(p3-pc);
  41.   qc=hcos_32[4]*(p4-pb);qd=hcos_32[5]*(p5-pa);
  42.   qe=hcos_32[6]*(p6-p9);qf=hcos_32[7]*(p7-p8);
  44.   p0=q0+q7;p1=q1+q6;p2=q2+q5;p3=q3+q4;
  45.   p4=hcos_16[0]*(q0-q7);p5=hcos_16[1]*(q1-q6);
  46.   p6=hcos_16[2]*(q2-q5);p7=hcos_16[3]*(q3-q4);
  47.   p8=q8+qf;p9=q9+qe;pa=qa+qd;pb=qb+qc;
  48.   pc=hcos_16[0]*(q8-qf);pd=hcos_16[1]*(q9-qe);
  49.   pe=hcos_16[2]*(qa-qd);pf=hcos_16[3]*(qb-qc);
  51.   q0=p0+p3;q1=p1+p2;q2=hcos_8[0]*(p0-p3);q3=hcos_8[1]*(p1-p2);
  52.   q4=p4+p7;q5=p5+p6;q6=hcos_8[0]*(p4-p7);q7=hcos_8[1]*(p5-p6);
  53.   q8=p8+pb;q9=p9+pa;qa=hcos_8[0]*(p8-pb);qb=hcos_8[1]*(p9-pa);
  54.   qc=pc+pf;qd=pd+pe;qe=hcos_8[0]*(pc-pf);qf=hcos_8[1]*(pd-pe);
  56.   p0=q0+q1;p1=hcos_4*(q0-q1);p2=q2+q3;p3=hcos_4*(q2-q3);
  57.   p4=q4+q5;p5=hcos_4*(q4-q5);p6=q6+q7;p7=hcos_4*(q6-q7);
  58.   p8=q8+q9;p9=hcos_4*(q8-q9);pa=qa+qb;pb=hcos_4*(qa-qb);
  59.   pc=qc+qd;pd=hcos_4*(qc-qd);pe=qe+qf;pf=hcos_4*(qe-qf);
  61.   {
  62.     register REAL tmp;
  64.     tmp=p6+p7;
  65.     OUT2(36)=-(p5+tmp);
  66.     OUT2(44)=-(p4+tmp);
  67.     tmp=pb+pf;
  68.     OUT1(10,tmp);
  69.     OUT1(6,pd+tmp);
  70.     tmp=pe+pf;
  71.     OUT2(46)=-(p8+pc+tmp);
  72.     OUT2(34)=-(p9+pd+tmp);
  73.     tmp+=pa+pb;
  74.     OUT2(38)=-(pd+tmp);
  75.     OUT2(42)=-(pc+tmp);
  76.     OUT1(2,p9+pd+pf);
  77.     OUT1(4,p5+p7);
  78.     OUT2(48)=-p0;
  79.     out2[0]=-(out1[0]=p1);
  80.     OUT1( 8,p3);
  81.     OUT1(12,p7);
  82.     OUT1(14,pf);
  83.     OUT2(40)=-(p2+p3);
  84.   }
  86.   {
  87.     register REAL *x=fraction;
  89.     p0=hcos_64[ 0]*(x[ 0]-x[31]);p1=hcos_64[ 1]*(x[ 1]-x[30]);
  90.     p2=hcos_64[ 2]*(x[ 2]-x[29]);p3=hcos_64[ 3]*(x[ 3]-x[28]);
  91.     p4=hcos_64[ 4]*(x[ 4]-x[27]);p5=hcos_64[ 5]*(x[ 5]-x[26]);
  92.     p6=hcos_64[ 6]*(x[ 6]-x[25]);p7=hcos_64[ 7]*(x[ 7]-x[24]);
  93.     p8=hcos_64[ 8]*(x[ 8]-x[23]);p9=hcos_64[ 9]*(x[ 9]-x[22]);
  94.     pa=hcos_64[10]*(x[10]-x[21]);pb=hcos_64[11]*(x[11]-x[20]);
  95.     pc=hcos_64[12]*(x[12]-x[19]);pd=hcos_64[13]*(x[13]-x[18]);
  96.     pe=hcos_64[14]*(x[14]-x[17]);pf=hcos_64[15]*(x[15]-x[16]);
  97.   }
  99.   q0=p0+pf;q1=p1+pe;q2=p2+pd;q3=p3+pc;
  100.   q4=p4+pb;q5=p5+pa;q6=p6+p9;q7=p7+p8;
  101.   q8=hcos_32[0]*(p0-pf);q9=hcos_32[1]*(p1-pe);
  102.   qa=hcos_32[2]*(p2-pd);qb=hcos_32[3]*(p3-pc);
  103.   qc=hcos_32[4]*(p4-pb);qd=hcos_32[5]*(p5-pa);
  104.   qe=hcos_32[6]*(p6-p9);qf=hcos_32[7]*(p7-p8);
  106.   p0=q0+q7;p1=q1+q6;p2=q2+q5;p3=q3+q4;
  107.   p4=hcos_16[0]*(q0-q7);p5=hcos_16[1]*(q1-q6);
  108.   p6=hcos_16[2]*(q2-q5);p7=hcos_16[3]*(q3-q4);
  109.   p8=q8+qf;p9=q9+qe;pa=qa+qd;pb=qb+qc;
  110.   pc=hcos_16[0]*(q8-qf);pd=hcos_16[1]*(q9-qe);
  111.   pe=hcos_16[2]*(qa-qd);pf=hcos_16[3]*(qb-qc);
  113.   q0=p0+p3;q1=p1+p2;q2=hcos_8[0]*(p0-p3);q3=hcos_8[1]*(p1-p2);
  114.   q4=p4+p7;q5=p5+p6;q6=hcos_8[0]*(p4-p7);q7=hcos_8[1]*(p5-p6);
  115.   q8=p8+pb;q9=p9+pa;qa=hcos_8[0]*(p8-pb);qb=hcos_8[1]*(p9-pa);
  116.   qc=pc+pf;qd=pd+pe;qe=hcos_8[0]*(pc-pf);qf=hcos_8[1]*(pd-pe);
  118.   p0=q0+q1;p1=hcos_4*(q0-q1);
  119.   p2=q2+q3;p3=hcos_4*(q2-q3);
  120.   p4=q4+q5;p5=hcos_4*(q4-q5);
  121.   p6=q6+q7;p7=hcos_4*(q6-q7);
  122.   p8=q8+q9;p9=hcos_4*(q8-q9);
  123.   pa=qa+qb;pb=hcos_4*(qa-qb);
  124.   pc=qc+qd;pd=hcos_4*(qc-qd);
  125.   pe=qe+qf;pf=hcos_4*(qe-qf);
  127.   {
  128.     REAL tmp;
  130.     tmp=pd+pf;
  131.     OUT1(5,p5+p7+pb+tmp);
  132.     tmp+=p9;
  133.     OUT1(1,p1+tmp);
  134.     OUT2(33)=-(p1+pe+tmp);
  135.     tmp+=p5+p7;
  136.     OUT1(3,tmp);
  137.     OUT2(35)=-(p6+pe+tmp);
  138.     tmp=pa+pb+pc+pd+pe+pf;
  139.     OUT2(39)=-(p2+p3+tmp-pc);
  140.     OUT2(43)=-(p4+p6+p7+tmp-pd);
  141.     OUT2(37)=-(p5+p6+p7+tmp-pc);
  142.     OUT2(41)=-(p2+p3+tmp-pd);
  143.     tmp=p8+pc+pe+pf;
  144.     OUT2(47)=-(p0+tmp);
  145.     OUT2(45)=-(p4+p6+p7+tmp);
  146.     tmp=pb+pf;
  147.     OUT1(11,p7+tmp);
  148.     tmp+=p3;
  149.     OUT1( 9,tmp);
  150.     OUT1( 7,pd+tmp);
  151.     OUT1(13,p7+pf);
  152.     OUT1(15,pf);
  153.   }
  154. }
  157. #define SAVE 
  158.         raw=(int)(r*scalefactor); 
  159.         if(raw>MAXSCALE)raw=MAXSCALE;else if(raw<MINSCALE)raw=MINSCALE; 
  160. putraw(raw);
  161. #define OS  r=*vp * *dp++
  162. #define XX  vp+=15;r+=*vp * *dp++
  163. #define OP  r+=*--vp * *dp++
  165. void MPEGaudio::generatesingle(void)
  166. {
  167.   int i;
  168.   register REAL r, *vp;
  169.   register const REAL *dp;
  170.   int raw;
  172.   i=32;
  173.   dp=filter;
  174.   vp=calcbufferL[currentcalcbuffer]+calcbufferoffset; 
  175. // actual_v+actual_write_pos;
  177.   switch (calcbufferoffset)
  178.   {
  179.     case  0:for(;i;i--,vp+=15){
  180.               OS;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  181.       SAVE;}break;
  182.     case  1:for(;i;i--,vp+=15){
  183.       OS;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  184.       SAVE;}break;
  185.     case  2:for(;i;i--,vp+=15){
  186.               OS;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  187.       SAVE;}break;
  188.     case  3:for(;i;i--,vp+=15){
  189.       OS;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  190.       SAVE;}break;
  191.     case  4:for(;i;i--,vp+=15){
  192.               OS;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  193.       SAVE;}break;
  194.     case  5:for(;i;i--,vp+=15){
  195.       OS;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  196.       SAVE;}break;
  197.     case  6:for(;i;i--,vp+=15){
  198.       OS;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;
  199.       SAVE;}break;
  200.     case  7:for(;i;i--,vp+=15){
  201.       OS;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;
  202.       SAVE;}break;
  203.     case  8:for(;i;i--,vp+=15){
  204.               OS;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;
  205.       SAVE;}break;
  206.     case  9:for(;i;i--,vp+=15){
  207.       OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;
  208.       SAVE;}break;
  209.     case 10:for(;i;i--,vp+=15){
  210.               OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;
  211.       SAVE;}break;
  212.     case 11:for(;i;i--,vp+=15){
  213.               OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;
  214.       SAVE;}break;
  215.     case 12:for(;i;i--,vp+=15){
  216.       OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;
  217.       SAVE;}break;
  218.     case 13:for(;i;i--,vp+=15){
  219.       OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;
  220.       SAVE;}break;
  221.     case 14:for(;i;i--,vp+=15){
  222.       OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;
  223.       SAVE;}break;
  224.     case 15:for(;i;i--,vp+=31){
  225.               OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  226.       SAVE;}break;
  227.   }
  228. }
  230. #undef OS
  231. #undef XX
  232. #undef OP
  233. #undef SAVE
  235. #define SAVE 
  236.         raw=(int)(r1*scalefactor);  
  237.         if(raw>MAXSCALE)raw=MAXSCALE;else if(raw<MINSCALE)raw=MINSCALE; 
  238. putraw(raw);  
  239.         raw=(int)(r2*scalefactor);  
  240.         if(raw>MAXSCALE)raw=MAXSCALE;else if(raw<MINSCALE)raw=MINSCALE; 
  241. putraw(raw);
  242. #define OS r1=*vp1 * *dp; 
  243.            r2=*vp2 * *dp++ 
  244. #define XX vp1+=15;r1+=*vp1 * *dp; 
  245.    vp2+=15;r2+=*vp2 * *dp++
  246. #define OP r1+=*--vp1 * *dp; 
  247.    r2+=*--vp2 * *dp++
  250. void MPEGaudio::generate(void)
  251. {
  252.   int i;
  253.   REAL r1,r2;
  254.   register REAL *vp1,*vp2;
  255.   register const REAL *dp;
  256.   int raw;
  258.   dp=filter;
  259.   vp1=calcbufferL[currentcalcbuffer]+calcbufferoffset;
  260.   vp2=calcbufferR[currentcalcbuffer]+calcbufferoffset;
  261. // actual_v+actual_write_pos;
  263.   i=32;
  264.   switch (calcbufferoffset)
  265.   {
  266.     case  0:for(;i;i--,vp1+=15,vp2+=15){
  267.               OS;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  268.       SAVE;}break;
  269.     case  1:for(;i;i--,vp1+=15,vp2+=15){
  270.       OS;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  271.       SAVE;}break;
  272.     case  2:for(;i;i--,vp1+=15,vp2+=15){
  273.               OS;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  274.       SAVE;}break;
  275.     case  3:for(;i;i--,vp1+=15,vp2+=15){
  276.       OS;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  277.       SAVE;}break;
  278.     case  4:for(;i;i--,vp1+=15,vp2+=15){
  279.               OS;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  280.       SAVE;}break;
  281.     case  5:for(;i;i--,vp1+=15,vp2+=15){
  282.       OS;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  283.       SAVE;}break;
  284.     case  6:for(;i;i--,vp1+=15,vp2+=15){
  285.       OS;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;OP;
  286.       SAVE;}break;
  287.     case  7:for(;i;i--,vp1+=15,vp2+=15){
  288.       OS;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;OP;
  289.       SAVE;}break;
  290.     case  8:for(;i;i--,vp1+=15,vp2+=15){
  291.               OS;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;OP;
  292.       SAVE;}break;
  293.     case  9:for(;i;i--,vp1+=15,vp2+=15){
  294.       OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;OP;
  295.       SAVE;}break;
  296.     case 10:for(;i;i--,vp1+=15,vp2+=15){
  297.               OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;OP;
  298.       SAVE;}break;
  299.     case 11:for(;i;i--,vp1+=15,vp2+=15){
  300.               OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;OP;
  301.       SAVE;}break;
  302.     case 12:for(;i;i--,vp1+=15,vp2+=15){
  303.       OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;OP;
  304.       SAVE;}break;
  305.     case 13:for(;i;i--,vp1+=15,vp2+=15){
  306.       OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;OP;
  307.       SAVE;}break;
  308.     case 14:for(;i;i--,vp1+=15,vp2+=15){
  309.       OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;XX;
  310.       SAVE;}break;
  311.     case 15:for(;i;i--,vp1+=31,vp2+=31){
  312.               OS;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;OP;
  313.       SAVE;}break;
  314.   }
  315. }
  318. void MPEGaudio::subbandsynthesis(REAL *fractionL,REAL *fractionR)
  319. {
  320.   if(downfrequency)
  321.   {
  322.     subbandsynthesis_2(fractionL,fractionR);
  323.     return;
  324.   }
  326.   computebuffer(fractionL,calcbufferL);
  327.   if(!outputstereo)generatesingle();
  328.   else
  329.   {
  330.     computebuffer(fractionR,calcbufferR);
  331.     generate();
  332.   }
  334.   if(calcbufferoffset<15)calcbufferoffset++;
  335.   else calcbufferoffset=0;
  337.   currentcalcbuffer^=1;
  338. }