diff options
Diffstat (limited to 'apps/plugins/pdbox/PDa/extra/notch.c')
-rw-r--r-- | apps/plugins/pdbox/PDa/extra/notch.c | 95 |
1 files changed, 1 insertions, 94 deletions
diff --git a/apps/plugins/pdbox/PDa/extra/notch.c b/apps/plugins/pdbox/PDa/extra/notch.c index 49be753ff2..0ab4332166 100644 --- a/apps/plugins/pdbox/PDa/extra/notch.c +++ b/apps/plugins/pdbox/PDa/extra/notch.c | |||
@@ -82,97 +82,4 @@ void notch_setup(void) | |||
82 | class_addbang(notch_class,notch_bang); | 82 | class_addbang(notch_class,notch_bang); |
83 | class_addfloat(notch_class,notch_float); | 83 | class_addfloat(notch_class,notch_float); |
84 | } | 84 | } |
85 | 85 | ||
86 | |||
87 | |||
88 | |||
89 | |||
90 | /* (C) Guenter Geiger <geiger@epy.co.at> */ | ||
91 | |||
92 | |||
93 | /* | ||
94 | |||
95 | These filter coefficients computations are taken from | ||
96 | http://www.harmony-central.com/Computer/Programming/Audio-EQ-Cookbook.txt | ||
97 | |||
98 | written by Robert Bristow-Johnson | ||
99 | |||
100 | */ | ||
101 | |||
102 | #include "m_pd.h" | ||
103 | #ifdef NT | ||
104 | #pragma warning( disable : 4244 ) | ||
105 | #pragma warning( disable : 4305 ) | ||
106 | #endif | ||
107 | #include <math.h> | ||
108 | #include "filters.h" | ||
109 | |||
110 | |||
111 | |||
112 | /* ------------------- notch ----------------------------*/ | ||
113 | |||
114 | static t_class *notch_class; | ||
115 | |||
116 | void notch_bang(t_rbjfilter *x) | ||
117 | { | ||
118 | t_atom at[5]; | ||
119 | t_float omega = e_omega(x->x_freq,x->x_rate); | ||
120 | t_float alpha = e_alpha(x->x_bw* 0.01,omega); | ||
121 | t_float b1 = -2.*cos(omega); | ||
122 | t_float b0 = 1; | ||
123 | t_float b2 = b0; | ||
124 | t_float a0 = 1 + alpha; | ||
125 | t_float a1 = -2.*cos(omega); | ||
126 | t_float a2 = 1 - alpha; | ||
127 | |||
128 | /* post("bang %f %f %f",x->x_freq, x->x_gain, x->x_bw); */ | ||
129 | |||
130 | if (!check_stability(-a1/a0,-a2/a0,b0/a0,b1/a0,b2/a0)) { | ||
131 | post("notch: filter unstable -> resetting"); | ||
132 | a0=1.;a1=0.;a2=0.; | ||
133 | b0=1.;b1=0.;b2=0.; | ||
134 | } | ||
135 | |||
136 | SETFLOAT(at,-a1/a0); | ||
137 | SETFLOAT(at+1,-a2/a0); | ||
138 | SETFLOAT(at+2,b0/a0); | ||
139 | SETFLOAT(at+3,b1/a0); | ||
140 | SETFLOAT(at+4,b2/a0); | ||
141 | |||
142 | outlet_list(x->x_obj.ob_outlet,&s_list,5,at); | ||
143 | } | ||
144 | |||
145 | |||
146 | void notch_float(t_rbjfilter *x,t_floatarg f) | ||
147 | { | ||
148 | x->x_freq = f; | ||
149 | notch_bang(x); | ||
150 | } | ||
151 | |||
152 | |||
153 | static void *notch_new(t_floatarg f,t_floatarg bw) | ||
154 | { | ||
155 | t_rbjfilter *x = (t_rbjfilter *)pd_new(notch_class); | ||
156 | |||
157 | x->x_rate = 44100.0; | ||
158 | outlet_new(&x->x_obj,&s_float); | ||
159 | /* floatinlet_new(&x->x_obj, &x->x_gain); */ | ||
160 | floatinlet_new(&x->x_obj, &x->x_bw); | ||
161 | if (f > 0.) x->x_freq = f; | ||
162 | if (bw > 0.) x->x_bw = bw; | ||
163 | return (x); | ||
164 | } | ||
165 | |||
166 | |||
167 | void notch_setup(void) | ||
168 | { | ||
169 | notch_class = class_new(gensym("notch"), (t_newmethod)notch_new, 0, | ||
170 | sizeof(t_rbjfilter), 0,A_DEFFLOAT,A_DEFFLOAT,0); | ||
171 | class_addbang(notch_class,notch_bang); | ||
172 | class_addfloat(notch_class,notch_float); | ||
173 | } | ||
174 | |||
175 | |||
176 | |||
177 | |||
178 | |||