home *** CD-ROM | disk | FTP | other *** search
/ PC World 2007 March / PCWorld_2007-03_cd.bin / audio-video / reaper / reaper174-install.exe / Effects / SStillwell / rbj1073 < prev    next >
Text File  |  2007-01-15  |  8KB  |  295 lines
  1. // Copyright 2006, Thomas Scott Stillwell
  2. // All rights reserved.
  3. //
  4. //Redistribution and use in source and binary forms, with or without modification, are permitted 
  5. //provided that the following conditions are met:
  6. //
  7. //Redistributions of source code must retain the above copyright notice, this list of conditions 
  8. //and the following disclaimer. 
  9. //
  10. //Redistributions in binary form must reproduce the above copyright notice, this list of conditions 
  11. //and the following disclaimer in the documentation and/or other materials provided with the distribution. 
  12. //
  13. //The name of Thomas Scott Stillwell may not be used to endorse or 
  14. //promote products derived from this software without specific prior written permission. 
  15. //
  16. //THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR 
  17. //IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND 
  18. //FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS 
  19. //BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES 
  20. //(INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR 
  21. //PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, 
  22. //STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF 
  23. //THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  24.  
  25. desc:1073-ish EQ, based on RBJ Filter Cookbook
  26.  
  27. slider1:0<0,4,1{Off,50,80,160,300}>HPF
  28. slider2:0<0,4,1{Off,35,60,110,220}>Low Shelf Freq
  29. slider3:0<-20,20,0.1>Low Boost/Cut
  30. slider4:0<0,5,1{360,700,1.6k,3.2k,4.8k,7.2k}>Mid Freq
  31. slider5:0<-20,20,0.1>Mid Boost/Cut
  32. slider6:0<-20,20,0.1>Hi Shelf (12k) Boost/Cut
  33. slider7:0<-20,10,0.1>Gain
  34.  
  35. @init
  36.   hpf = 0;
  37.   gain1 = 0;
  38.   freq1 = 50;
  39.   a1 = 1;
  40.   s1 = 1;
  41.   q1 = 1 / (sqrt((a1 + 1/a1)*(1/s1 - 1) + 2));
  42.   w01 = 2 * $pi * freq1/srate;
  43.   cosw01 = cos(w01);
  44.   sinw01 = sin(w01);
  45.   alpha1 = sinw01 / (2 * q1);
  46.  
  47.   b01 = (1 + cosw01)/2;
  48.   b11 = -(1 + cosw01);
  49.   b21 = (1 + cosw01)/2;
  50.   a01 = 1 + alpha1;
  51.   a11 = -2 * cosw01;
  52.   a21 = 1 - alpha1;
  53.   b01 /= a01;
  54.   b11 /= a01;
  55.   b21 /= a01;
  56.   a11 /= a01;
  57.   a21 /= a01;
  58.  
  59.   lshelf = 0;
  60.   gain3 = 0;
  61.   freq3 = 35;
  62.   a3 = 10^(gain3/40);
  63.   s3 = 2;
  64.   q3 = 1 / (sqrt((a3 + 1/a3)*(1/s3 - 1) + 2));
  65.   w03 = 2 * $pi * freq3/srate;
  66.   cosw03 = cos(w03);
  67.   sinw03 = sin(w03);
  68.   alpha3 = sinw03 / (2 * q3);
  69.  
  70.   b03 = a3 * ((a3+1) - (a3-1)*cosw03 + 2*sqrt(a3)*alpha3);
  71.   b13 = 2 * a3 * ((a3-1) - (a3+1)*cosw03);
  72.   b23 = a3 * ((a3+1) - (a3-1)*cosw03 - 2*sqrt(a3)*alpha3);
  73.   a03 = (a3+1) + (a3-1)*cosw03 + 2*sqrt(a3)*alpha3;
  74.   a13 = -2 * ((a3-1) + (a3+1)*cosw03);
  75.   a23 = (a3+1)+(a3-1)*cosw03-2*sqrt(a3)*alpha3;
  76.   b03 /= a03;
  77.   b13 /= a03;
  78.   b23 /= a03;
  79.   a13 /= a03;
  80.   a23 /= a03;
  81.  
  82.   gain5 = 0;
  83.   freq5 = 360;
  84.   a5 = 10^(gain5/20);
  85.   q5 = 1.4;
  86.   w05 = 2 * $pi * freq5/srate;
  87.   cosw05 = cos(w05);
  88.   sinw05 = sin(w05);
  89.   alpha5 = sinw05 / (2 * q5);
  90.  
  91.   b05 = 1 + alpha5 * a5;
  92.   b15 = -2 * cosw05;
  93.   b25 = 1 - alpha5 * a5;
  94.   a05 = 1 + alpha5 / a5;
  95.   a15 = -2 * cosw05;
  96.   a25 = 1 - alpha5 / a5;
  97.   b05 /= a05;
  98.   b15 /= a05;
  99.   b25 /= a05;
  100.   a15 /= a05;
  101.   a25 /= a05;
  102.  
  103.   gain7 = 0;
  104.   freq7 = 12000;
  105.   a7 = 10^(gain7/40);
  106.   s7 = 0.3;
  107.   q7 = 1 / (sqrt((a7 + 1/a7)*(1/s7 - 1) + 2));
  108.   w07 = 2 * $pi * freq7/srate;
  109.   cosw07 = cos(w07);
  110.   sinw07 = sin(w07);
  111.   alpha7 = sinw07 / (2 * q7);
  112.  
  113.   b07 = a7 * ((a7+1) + (a7-1)*cosw07 + 2*sqrt(a7)*alpha7);
  114.   b17 = -2*a7*((a7-1) + (a7+1)*cosw07);
  115.   b27 = a7*((a7+1) + (a7-1)*cosw07 - 2*sqrt(a7)*alpha7);
  116.   a07 = (a7+1) - (a7-1)*cosw07 + 2*sqrt(a7)*alpha7;
  117.   a17 = 2*((a7-1) - (a7+1)*cosw07);
  118.   a27 = (a7+1)-(a7-1)*cosw07 - 2*sqrt(a7)*alpha7;
  119.   b07 /= a07;
  120.   b17 /= a07;
  121.   b27 /= a07;
  122.   a17 /= a07;
  123.   a27 /= a07;
  124.  
  125.   gain = 1;
  126.  
  127. cDenorm = 10^-30;
  128.  
  129. @slider
  130.   freq1 = (slider1 == 0 ? 50 : (slider1 == 1 ? 50 : (slider1 == 2 ? 80 : (slider1 == 3 ? 160 : 300))));
  131.   freq3 = (slider2 == 0 ? 35 : (slider2 == 1 ? 35 : (slider2 == 2 ? 60 : (slider2 == 3 ? 110 : 220))));
  132.   gain3 = slider3;
  133.   freq5 = (slider4 == 0 ? 360 : (slider4 == 1 ? 700 : (slider4 == 2 ? 1600 : (slider4 == 3 ? 3200 : (slider4 == 4 ? 4800 : 7200)))));
  134.   gain5 = slider5;
  135.   gain7 = slider6;
  136.  
  137.   gain = 10^(slider7/20);
  138.  
  139.   slider1 == 0 ? hpf = 0 : hpf = 1;
  140.   slider2 == 0 ? lshelf = 0 : lshelf = 1;
  141.  
  142.   a1 = 1;
  143.   s1 = 1;
  144.   q1 = 1 / (sqrt((a1 + 1/a1)*(1/s1 - 1) + 2));
  145.   w01 = 2 * $pi * freq1/srate;
  146.   cosw01 = cos(w01);
  147.   sinw01 = sin(w01);
  148.   alpha1 = sinw01 / (2 * q1);
  149.  
  150.   b01 = (1 + cosw01)/2;
  151.   b11 = -(1 + cosw01);
  152.   b21 = (1 + cosw01)/2;
  153.   a01 = 1 + alpha1;
  154.   a11 = -2 * cosw01;
  155.   a21 = 1 - alpha1;
  156.   b01 /= a01;
  157.   b11 /= a01;
  158.   b21 /= a01;
  159.   a11 /= a01;
  160.   a21 /= a01;
  161.  
  162.   a3 = 10^(gain3/40);
  163.   s3 = 2;
  164.   q3 = 1 / (sqrt((a3 + 1/a3)*(1/s3 - 1) + 2));
  165.   w03 = 2 * $pi * freq3/srate;
  166.   cosw03 = cos(w03);
  167.   sinw03 = sin(w03);
  168.   alpha3 = sinw03 / (2 * q3);
  169.  
  170.   b03 = a3 * ((a3+1) - (a3-1)*cosw03 + 2*sqrt(a3)*alpha3);
  171.   b13 = 2 * a3 * ((a3-1) - (a3+1)*cosw03);
  172.   b23 = a3 * ((a3+1) - (a3-1)*cosw03 - 2*sqrt(a3)*alpha3);
  173.   a03 = (a3+1) + (a3-1)*cosw03 + 2*sqrt(a3)*alpha3;
  174.   a13 = -2 * ((a3-1) + (a3+1)*cosw03);
  175.   a23 = (a3+1)+(a3-1)*cosw03-2*sqrt(a3)*alpha3;
  176.   b03 /= a03;
  177.   b13 /= a03;
  178.   b23 /= a03;
  179.   a13 /= a03;
  180.   a23 /= a03;
  181.  
  182.   a5 = 10^(gain5/20);
  183.   q5 = 1.4;
  184.   w05 = 2 * $pi * freq5/srate;
  185.   cosw05 = cos(w05);
  186.   sinw05 = sin(w05);
  187.   alpha5 = sinw05 / (2 * q5);
  188.  
  189.   b05 = 1 + alpha5 * a5;
  190.   b15 = -2 * cosw05;
  191.   b25 = 1 - alpha5 * a5;
  192.   a05 = 1 + alpha5 / a5;
  193.   a15 = -2 * cosw05;
  194.   a25 = 1 - alpha5 / a5;
  195.   b05 /= a05;
  196.   b15 /= a05;
  197.   b25 /= a05;
  198.   a15 /= a05;
  199.   a25 /= a05;
  200.  
  201.   a7 = 10^(gain7/40);
  202.   freq7 = 12000;
  203.   s7 = 0.3;
  204.   q7 = 1 / (sqrt((a7 + 1/a7)*(1/s7 - 1) + 2));
  205.   w07 = 2 * $pi * freq7/srate;
  206.   cosw07 = cos(w07);
  207.   sinw07 = sin(w07);
  208.   alpha7 = sinw07 / (2 * q7);
  209.  
  210.   b07 = a7 * ((a7+1) + (a7-1)*cosw07 + 2*sqrt(a7)*alpha7);
  211.   b17 = -2*a7*((a7-1) + (a7+1)*cosw07);
  212.   b27 = a7*((a7+1) + (a7-1)*cosw07 - 2*sqrt(a7)*alpha7);
  213.   a07 = (a7+1) - (a7-1)*cosw07 + 2*sqrt(a7)*alpha7;
  214.   a17 = 2*((a7-1) - (a7+1)*cosw07);
  215.   a27 = (a7+1)-(a7-1)*cosw07 - 2*sqrt(a7)*alpha7;
  216.   b07 /= a07;
  217.   b17 /= a07;
  218.   b27 /= a07;
  219.   a17 /= a07;
  220.   a27 /= a07;
  221.  
  222. @sample
  223.  
  224.   hpf != 0 ? (
  225.   ospl0 = spl0;
  226.   spl0 = b01 * spl0 + b11 * xl11 + b21 * xl21 - a11 * yl11 - a21 * yl21;
  227.   xl21 = xl11;
  228.   xl11 = ospl0;
  229.   yl21 = yl11;
  230.   yl11 = abs(spl0) < cDenorm ? 0 : spl0 ;
  231.  
  232.   ospl1 = spl1;
  233.   spl1 = b01 * spl1 + b11 * xr11 + b21 * xr21 - a11 * yr11 - a21 * yr21;
  234.   xr21 = xr11;
  235.   xr11 = ospl1;
  236.   yr21 = yr11;
  237.   yr11 = abs(spl1) < cDenorm ? 0 : spl1 ;
  238.   );
  239.  
  240. spl0 += cDenorm;
  241. spl1 += cDenorm;
  242.  
  243.   lshelf != 0 && gain3 != 0 ? (
  244.   ospl0 = spl0;
  245.   spl0 = b03 * spl0 + b13 * xl13 + b23 * xl23 - a13 * yl13 - a23 * yl23;
  246.   xl23 = xl13;
  247.   xl13 = ospl0;
  248.   yl23 = yl13;
  249.   yl13 = spl0;
  250.  
  251.   ospl1 = spl1;
  252.   spl1 = b03 * spl1 + b13 * xr13 + b23 * xr23 - a13 * yr13 - a23 * yr23;
  253.   xr23 = xr13;
  254.   xr13 = ospl1;
  255.   yr23 = yr13;
  256.   yr13 = spl1;
  257.   );
  258.  
  259.   gain5 != 0 ? (
  260.   ospl0 = spl0;
  261.   spl0 = b05 * spl0 + b15 * xl15 + b25 * xl25 - a15 * yl15 - a25 * yl25;
  262.   xl25 = xl15;
  263.   xl15 = ospl0;
  264.   yl25 = yl15;
  265.   yl15 = spl0;
  266.  
  267.   ospl1 = spl1;
  268.   spl1 = b05 * spl1 + b15 * xr15 + b25 * xr25 - a15 * yr15 - a25 * yr25;
  269.   xr25 = xr15;
  270.   xr15 = ospl1;
  271.   yr25 = yr15;
  272.   yr15 = spl1;
  273.   );
  274.  
  275.   gain7 != 0 ? (
  276.   ospl0 = spl0;
  277.   spl0 = b07 * spl0 + b17 * xl17 + b27 * xl27 - a17 * yl17 - a27 * yl27;
  278.   xl27 = xl17;
  279.   xl17 = ospl0;
  280.   yl27 = yl17;
  281.   yl17 = spl0;
  282.  
  283.   ospl1 = spl1;
  284.   spl1 = b07 * spl1 + b17 * xr17 + b27 * xr27 - a17 * yr17 - a27 * yr27;
  285.   xr27 = xr17;
  286.   xr17 = ospl1;
  287.   yr27 = yr17;
  288.   yr17 = spl1;
  289.   );
  290.  
  291.   spl0 *= gain;
  292.   spl1 *= gain;
  293.  
  294.  
  295.