// Copyright 2006, Thomas Scott Stillwell
// All rights reserved.
//
//Redistribution and use in source and binary forms, with or without modification, are permitted 
//provided that the following conditions are met:
//
//Redistributions of source code must retain the above copyright notice, this list of conditions 
//and the following disclaimer. 
//
//Redistributions in binary form must reproduce the above copyright notice, this list of conditions 
//and the following disclaimer in the documentation and/or other materials provided with the distribution. 
//
//The name of Thomas Scott Stillwell may not be used to endorse or 
//promote products derived from this software without specific prior written permission. 
//
//THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR 
//IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND 
//FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS 
//BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES 
//(INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR 
//PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, 
//STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF 
//THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
desc:peak-eating limiter with program-dependent release

slider1:0.0<-30.0,0.0,0.1>Threshold
slider2:-0.1<-20.0,0.0,0.1>Ceiling
slider3:0<0,1200,1>Release (ms), 0=auto

@init
  pi = 3.1415926535;
  log2db = 8.6858896380650365530225783783321; // 20 / ln(10)
  db2log = 0.11512925464970228420089957273422; // ln(10) / 20 
  attime=0.004;
  reltime=0.200;
  atcoef=exp(-1/(attime * srate));
  relcoef=exp(-1/(reltime * srate));

@slider
  thresh=exp(slider1 * db2log);
  threshdb=slider1;
  ceiling=exp(slider2 * db2log);
  ceildb=slider2;
  makeup=exp((ceildb-threshdb) * db2log);
  makeupdb=ceildb-threshdb;
  peakdb = ceildb+25;
  peaklvl = exp(peakdb * db2log);
  release=slider3/1000;
  release==0 ? (reltime = overdb / 125) : (reltime = release);
  relcoef=exp(-1/(reltime * srate));


@sample
  spl0=spl0*makeup;
  spl1=spl1*makeup;

  abs0=abs(spl0);
  abs1=abs(spl1);
  overdbv = max(abs0, abs1);
  overdb = 2.08136898 * log(overdbv) * log2db - ceilingdb;

  overdb > rundb ? (                                  // calculate overshoot (and thus GR)
    rundb = overdb + atcoef * (rundb - overdb);       //    use attack and release to 
  ) : (                                               //    smoothly scale gain reduction
    rundb = overdb + relcoef * (rundb - overdb);
  );
  overdb = max(0,rundb);

  release == 0 ? (                                    // Auto-release?  Compute release time
    overdb > maxover ? (
      maxover = overdb;
      reltime = overdb / 125;                         // release at constant 125 dB/sec.
      relcoef = exp(-1/(reltime * srate));
    );
    runmax = maxover + relcoef * (runmax - maxover);  // highest peak for setting att/rel decays in reltime
    maxover = runmax;
  );
  
  gr = exp(overdb * db2log);

  spl0*=gr;
  spl1*=gr;
  
  spl0=min(ceiling,abs(spl0))*sign(spl0);             // If it gets past the limiter, whack it.
  spl1=min(ceiling,abs(spl1))*sign(spl1);             // since we don't lookahead, some stuff will
