home *** CD-ROM | disk | FTP | other *** search
/ PC World 2008 February / PCWorld_2008-02_cd.bin / audio-video / reaper / reaper2028-install.exe / Effects / SStillwell / majortom < prev    next >
Text File  |  2007-12-03  |  5KB  |  152 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.  
  17.  
  18. //THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR 
  19. //IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND 
  20. //FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS 
  21. //BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES 
  22. //(INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR 
  23. //PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, 
  24. //STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF 
  25. //THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  26. desc:program dependent compressor
  27.  
  28. slider1:0<-60,0,0.1>Threshold (dB)
  29. slider2:1<1,20,0.1>Ratio
  30. slider3:0<-20,20,0.1>Gain
  31. slider4:0<0,1,1{No,Yes}>Soft knee?
  32. slider5:0<0,1,1{Normal,Sidechain}>Detector input
  33. slider6:0<0,1,1{No,Yes}>Auto-makeup?
  34. slider7:0<0,1,1{Peak,RMS}>Detection
  35. slider8:0<0,1,1{Feedforward,Feedback}>Detection source
  36.  
  37. @init
  38.   log2db = 8.6858896380650365530225783783321; // 20 / ln(10)
  39.   db2log = 0.11512925464970228420089957273422; // ln(10) / 20 
  40.   i=0;
  41.   loop(
  42.     120,
  43.     attimes[i] = ((0.08924 / i) + (0.60755 / (i ^ 2)) - 0.00006);
  44.     i+=1;
  45.   );
  46.   attime=0.010;
  47.   reltime=0.100;
  48.   maxover=0;
  49.   ratio=0;
  50.   cratio=0;
  51.   rundb=0;
  52.   overdb=0;
  53.   maxover=0;
  54.   atcoef=exp(-1/(attime * srate));
  55.   relcoef=exp(-1/(reltime * srate));
  56.   fbacoef=exp(-1000/(2 * srate)); // 2 msec. opto attack for feedback detection
  57.   fbrcoef=exp(-1000/(200 * srate)); // 200 msec. opto release for feedback detection
  58.   sidechain = 0;
  59.   automakeup = 0;
  60.  
  61. @slider
  62.   thresh = slider1;
  63.   threshv = exp(thresh * db2log);
  64.   ratio = slider2;
  65.   softknee = slider4;
  66.   cthresh = (softknee ? (thresh -3) : thresh);
  67.   cthreshv = exp(cthresh * db2log);
  68.   sidechain = slider5;
  69.   automakeup = slider6;
  70.   automakeup ? (
  71.     autogain = (abs(thresh) - (abs(thresh)/max(1,ratio-1)))/2;
  72.   ) : (
  73.     autogain = 0;
  74.   );
  75.   makeup = slider3;
  76.   makeupv = exp((makeup+autogain) * db2log);
  77.   RMSdet = slider7;
  78.   RMSdet ? (
  79.     rmscoef=exp(-1000/(10 * srate));       // 10 ms RMS window
  80.   ) : (
  81.     rmscoef=exp(-1000/(0.0025 * srate));  // 2.5 us Peak detector
  82.   );
  83.   opto = slider8;
  84.  
  85. @sample
  86.   sidechain ? (
  87.     aspl0 = abs(spl2);
  88.     aspl1 = abs(spl3);
  89.   ) : (
  90.     opto ? (
  91.       ospl = ospl0 * ospl0 + ospl1 * ospl1;
  92.       ospl > runospl ? (
  93.         runospl = ospl + atcoef * (runospl - ospl);
  94.       ) : (
  95.         runospl = ospl + relcoef * (runospl - ospl);
  96.       );
  97.       ospl = sqrt(max(0,runospl));
  98.  
  99. //      ospl = 2.08136898 * log(ospl)*log2db;
  100. //      ospl = ospl - 6;            // reduce feedback to -6 dB
  101. //      ospl = exp(ospl * db2log);
  102.       ospl *= 0.5;
  103.  
  104.       aspl0 = abs(ospl);
  105.       aspl1 = abs(ospl);
  106.     ) : (
  107.       aspl0 = abs(spl0);
  108.       aspl1 = abs(spl1);
  109.     );
  110.   );
  111.  
  112.   RMSDet ? (  
  113.     ave = (aspl0 * aspl0) + (aspl1 * aspl1);
  114.     runave = ave + rmscoef * (runave - ave);
  115.     det = sqrt(max(0,runave));
  116.   ) : (
  117.     maxspl = max(aspl0, aspl1);
  118.     maxspl = maxspl * maxspl;
  119.     runave = maxspl + rmscoef * (runave - maxspl);
  120.     det = sqrt(max(0,runave));
  121.   );
  122.   overdb = 2.08136898 * log(det/cthreshv) * log2db;
  123.   overdb > maxover ? (
  124.     maxover = overdb;
  125.     attime = attimes[max(0,floor(abs(overdb)))];   // attack time per formula
  126.     atcoef = exp(-1/(attime * srate));
  127.     reltime = overdb / 125;                        // release at constant 125 dB/sec.
  128.     relcoef = exp(-1/(reltime * srate));
  129.   );
  130.   overdb = max(0,overdb);
  131.  
  132.   overdb > rundb ? (
  133.     rundb = overdb + atcoef * (rundb - overdb);
  134.   ) : (
  135.     rundb = overdb + relcoef * (rundb - overdb);
  136.   );
  137.   overdb = rundb;
  138.  
  139.   cratio = (softknee ? (1 + (ratio-1) * min(overdb, 6) / 6) : ratio);
  140.   
  141.   gr = -overdb * (cratio-1)/cratio;
  142.   grv = exp(gr * db2log);
  143.   
  144.   runmax = maxover + relcoef * (runmax - maxover);  // highest peak for setting att/rel decays in reltime
  145.   maxover = runmax;
  146.  
  147.   spl0 *= grv * makeupv;
  148.   spl1 *= grv * makeupv;  
  149.  
  150.   ospl0 = spl0;
  151.   ospl1 = spl1;
  152.