home *** CD-ROM | disk | FTP | other *** search
/ PC World 2006 October / PCWorld_2006-10_cd.bin / audio-video / reaper / reaper105-install.exe / Effects / Compressors / majortom
Text File  |  2006-09-02  |  5KB  |  151 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. //THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR 
  18. //IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND 
  19. //FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS 
  20. //BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES 
  21. //(INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR 
  22. //PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, 
  23. //STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF 
  24. //THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  25. desc:program dependent compressor
  26.  
  27. slider1:0<-60,0,0.1>Threshold (dB)
  28. slider2:1<1,20,0.1>Ratio
  29. slider3:0<-20,20,0.1>Gain
  30. slider4:0<0,1,1{No,Yes}>Soft knee?
  31. slider5:0<0,1,1{Normal,Sidechain}>Detector input
  32. slider6:0<0,1,1{No,Yes}>Auto-makeup?
  33. slider7:0<0,1,1{Peak,RMS}>Detection
  34. slider8:0<0,1,1{Feedforward,Feedback}>Detection source
  35.  
  36. @init
  37.   log2db = 8.6858896380650365530225783783321; // 20 / ln(10)
  38.   db2log = 0.11512925464970228420089957273422; // ln(10) / 20 
  39.   dcoffset = 0.0000000000000000000000001;
  40.   i=0;
  41.   loop(
  42.     120,
  43.     attimes[i] = ((0.08924 / i) + (0.60755 / (i ^ 2)) - 0.00006)+dcoffset;
  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 = log(ospl)*log2db;
  100.       ospl = ospl - 6;            // reduce feedback to -6 dB
  101.       ospl = exp(ospl * db2log);
  102.  
  103.       aspl0 = abs(ospl);
  104.       aspl1 = abs(ospl);
  105.     ) : (
  106.       aspl0 = abs(spl0);
  107.       aspl1 = abs(spl1);
  108.     );
  109.   );
  110.  
  111.   RMSDet ? (  
  112.     ave = (aspl0 * aspl0) + (aspl1 * aspl1);
  113.     runave = ave + rmscoef * (runave - ave);
  114.     det = sqrt(max(0,runave));
  115.   ) : (
  116.     maxspl = max(aspl0, aspl1);
  117.     maxspl = maxspl * maxspl;
  118.     runave = maxspl + rmscoef * (runave - maxspl);
  119.     det = sqrt(max(0,runave));
  120.   );
  121.   overdb = (log(det) * log2db) - (log(cthreshv) * log2db);
  122.   overdb > maxover ? (
  123.     maxover = overdb;
  124.     attime = attimes[max(0,floor(abs(overdb)))];   // attack time per formula
  125.     attcoef = exp(-1/(attime * srate));
  126.     reltime = overdb / 125;                        // release at constant 125 dB/sec.
  127.     relcoef = exp(-1/(reltime * srate));
  128.   );
  129.   overdb = max(0,overdb);
  130.  
  131.   overdb > rundb ? (
  132.     rundb = overdb + atcoef * (rundb - overdb);
  133.   ) : (
  134.     rundb = overdb + relcoef * (rundb - overdb);
  135.   );
  136.   overdb = rundb;
  137.  
  138.   cratio = (softknee ? (1 + (ratio-1) * min(overdb, 6) / 6) : ratio);
  139.   
  140.   gr = -overdb * (cratio-1)/cratio;
  141.   grv = exp(gr * db2log);
  142.   
  143.   runmax = maxover + relcoef * (runmax - maxover);  // highest peak for setting att/rel decays in reltime
  144.   maxover = runmax;
  145.  
  146.   spl0 *= grv * makeupv;
  147.   spl1 *= grv * makeupv;  
  148.  
  149.   ospl0 = spl0;
  150.   ospl1 = spl1;
  151.