blob: f1a6dda77825b190845ffe6dc152e23b691cd3c7 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
|
/************************************************************************************
*
* Wavechild670 v0.1
*
* tubemodel.cpp
*
* By Peter Raffensperger 10 July 2012
*
* Reference:
* Toward a Wave Digital Filter Model of the Fairchild 670 Limiter, Raffensperger, P. A., (2012).
* Proc. of the 15th International Conference on Digital Audio Effects (DAFx-12),
* York, UK, September 17-21, 2012.
*
* Note:
* Fairchild (R) a registered trademark of Avid Technology, Inc., which is in no way associated or
* affiliated with the author.
*
* License:
* Wavechild670 is licensed under the GNU GPL v2 license. If you use this
* software in an academic context, we would appreciate it if you referenced the original
* paper.
*
************************************************************************************/
#include "tubemodel.h"
const Real TriodeRemoteCutoff6386::fa = -0.1961135820501605;
const Real TriodeRemoteCutoff6386::aa = 3.980508168e-08;
const Real TriodeRemoteCutoff6386::ab = 2.3830020303;
const Real TriodeRemoteCutoff6386::e = -1.8;
const Real TriodeRemoteCutoff6386::d = 0.5;
const Real TriodeRemoteCutoff6386::c = 10.0;
const Real TriodeRemoteCutoff6386::g = 5.0;
const Real TriodeRemoteCutoff6386::h = 0.5;
Real WDFTubeInterface::iterateNewtonRaphson(Real x, Real dx){
/*
x(n+1) = x(n) - Fn(x)/Fn'(x)
Fn'(x) = dFn(x)/dx = (Fn(x + dx) - Fn(x))/(dx) #Finite difference approximation
x(n+1) = x(n) - dx*Fn(x)/(Fn(x+dx) - Fn(x))
*/
Real F = evaluateImplicitEquation(x);
Real xNew = x - dx*F/(evaluateImplicitEquation(x + dx) - F);
Assert(isfinite(xNew));
return xNew;
}
|
