home *** CD-ROM | disk | FTP | other *** search
/ PC World 2007 January / PCWorld_2007-01_cd.bin / temacd / reaper / reaper147-install.exe / Effects / SStillwell / majortom < prev    next >
Text File  |  2006-10-16  |  5KB  |  153 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.   dcoffset = 0.0000000000000000000000001;
  41.   i=0;
  42.   loop(
  43.     120,
  44.     attimes[i] = ((0.08924 / i) + (0.60755 / (i ^ 2)) - 0.00006)+dcoffset;
  45.     i+=1;
  46.   );
  47.   attime=0.010;
  48.   reltime=0.100;
  49.   maxover=0;
  50.   ratio=0;
  51.   cratio=0;
  52.   rundb=0;
  53.   overdb=0;
  54.   maxover=0;
  55.   atcoef=exp(-1/(attime * srate));
  56.   relcoef=exp(-1/(reltime * srate));
  57.   fbacoef=exp(-1000/(2 * srate)); // 2 msec. opto attack for feedback detection
  58.   fbrcoef=exp(-1000/(200 * srate)); // 200 msec. opto release for feedback detection
  59.   sidechain = 0;
  60.   automakeup = 0;
  61.  
  62. @slider
  63.   thresh = slider1;
  64.   threshv = exp(thresh * db2log);
  65.   ratio = slider2;
  66.   softknee = slider4;
  67.   cthresh = (softknee ? (thresh -3) : thresh);
  68.   cthreshv = exp(cthresh * db2log);
  69.   sidechain = slider5;
  70.   automakeup = slider6;
  71.   automakeup ? (
  72.     autogain = (abs(thresh) - (abs(thresh)/max(1,ratio-1)))/2;
  73.   ) : (
  74.     autogain = 0;
  75.   );
  76.   makeup = slider3;
  77.   makeupv = exp((makeup+autogain) * db2log);
  78.   RMSdet = slider7;
  79.   RMSdet ? (
  80.     rmscoef=exp(-1000/(10 * srate));       // 10 ms RMS window
  81.   ) : (
  82.     rmscoef=exp(-1000/(0.0025 * srate));  // 2.5 us Peak detector
  83.   );
  84.   opto = slider8;
  85.  
  86. @sample
  87.   sidechain ? (
  88.     aspl0 = abs(spl2);
  89.     aspl1 = abs(spl3);
  90.   ) : (
  91.     opto ? (
  92.       ospl = ospl0 * ospl0 + ospl1 * ospl1;
  93.       ospl > runospl ? (
  94.         runospl = ospl + atcoef * (runospl - ospl);
  95.       ) : (
  96.         runospl = ospl + relcoef * (runospl - ospl);
  97.       );
  98.       ospl = sqrt(max(0,runospl));
  99.  
  100. //      ospl = 2.08136898 * log(ospl)*log2db;
  101. //      ospl = ospl - 6;            // reduce feedback to -6 dB
  102. //      ospl = exp(ospl * db2log);
  103.       ospl *= 0.5;
  104.  
  105.       aspl0 = abs(ospl);
  106.       aspl1 = abs(ospl);
  107.     ) : (
  108.       aspl0 = abs(spl0);
  109.       aspl1 = abs(spl1);
  110.     );
  111.   );
  112.  
  113.   RMSDet ? (  
  114.     ave = (aspl0 * aspl0) + (aspl1 * aspl1);
  115.     runave = ave + rmscoef * (runave - ave);
  116.     det = sqrt(max(0,runave));
  117.   ) : (
  118.     maxspl = max(aspl0, aspl1);
  119.     maxspl = maxspl * maxspl;
  120.     runave = maxspl + rmscoef * (runave - maxspl);
  121.     det = sqrt(max(0,runave));
  122.   );
  123.   overdb = 2.08136898 * log(det/cthreshv) * log2db;
  124.   overdb > maxover ? (
  125.     maxover = overdb;
  126.     attime = attimes[max(0,floor(abs(overdb)))];   // attack time per formula
  127.     atcoef = exp(-1/(attime * srate));
  128.     reltime = overdb / 125;                        // release at constant 125 dB/sec.
  129.     relcoef = exp(-1/(reltime * srate));
  130.   );
  131.   overdb = max(0,overdb);
  132.  
  133.   overdb > rundb ? (
  134.     rundb = overdb + atcoef * (rundb - overdb);
  135.   ) : (
  136.     rundb = overdb + relcoef * (rundb - overdb);
  137.   );
  138.   overdb = rundb;
  139.  
  140.   cratio = (softknee ? (1 + (ratio-1) * min(overdb, 6) / 6) : ratio);
  141.   
  142.   gr = -overdb * (cratio-1)/cratio;
  143.   grv = exp(gr * db2log);
  144.   
  145.   runmax = maxover + relcoef * (runmax - maxover);  // highest peak for setting att/rel decays in reltime
  146.   maxover = runmax;
  147.  
  148.   spl0 *= grv * makeupv;
  149.   spl1 *= grv * makeupv;  
  150.  
  151.   ospl0 = spl0;
  152.   ospl1 = spl1;
  153.