interpretor/faust-0.9.47mr3/architecture/music.lib

   1 /************************************************************************
   2  ************************************************************************
   3         FAUST library file
   4         Copyright (C) 2003-2011 GRAME, Centre National de Creation Musicale
   5     ---------------------------------------------------------------------
   6     This program is free software; you can redistribute it and/or modify
   7     it under the terms of the GNU Lesser General Public License as
   8         published by the Free Software Foundation; either version 2.1 of the
   9         License, or (at your option) any later version.
  10
  11     This program is distributed in the hope that it will be useful,
  12     but WITHOUT ANY WARRANTY; without even the implied warranty of
  13     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  14     GNU Lesser General Public License for more details.
  15
  16     You should have received a copy of the GNU Lesser General Public
  17         License along with the GNU C Library; if not, write to the Free
  18         Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
  19         02111-1307 USA.
  20  ************************************************************************
  21  ************************************************************************/
  22
  23 declare name "Music Library";
  24 declare author "GRAME";
  25 declare copyright "GRAME";
  26 declare version "1.0";
  27 declare license "LGPL";
  28
  29 import("math.lib");
  30
  31 //-----------------------------------------------
  32 //                                      DELAY LINE
  33 //-----------------------------------------------
  34 frac(n)                 = n-int(n);
  35 index(n)                = &(n-1) ~ +(1);                // n = 2**i
  36 //delay(n,d,x)  = rwtable(n, 0.0, index(n), x, (index(n)-int(d)) & (n-1));
  37 delay(n,d,x)    = x@(int(d)&(n-1));
  38 fdelay(n,d,x)   = delay(n,int(d),x)*(1 - frac(d)) + delay(n,int(d)+1,x)*frac(d);
  39
  40
  41 delay1s(d)              = delay(65536,d);
  42 delay2s(d)              = delay(131072,d);
  43 delay5s(d)              = delay(262144,d);
  44 delay10s(d)             = delay(524288,d);
  45 delay21s(d)             = delay(1048576,d);
  46 delay43s(d)             = delay(2097152,d);
  47
  48 fdelay1s(d)             = fdelay(65536,d);
  49 fdelay2s(d)             = fdelay(131072,d);
  50 fdelay5s(d)             = fdelay(262144,d);
  51 fdelay10s(d)    = fdelay(524288,d);
  52 fdelay21s(d)    = fdelay(1048576,d);
  53 fdelay43s(d)    = fdelay(2097152,d);
  54
  55 millisec        = SR/1000.0;
  56
  57 time1s  = hslider("time", 0, 0,  1000, 0.1)*millisec;
  58 time2s  = hslider("time", 0, 0,  2000, 0.1)*millisec;
  59 time5s  = hslider("time", 0, 0,  5000, 0.1)*millisec;
  60 time10s = hslider("time", 0, 0, 10000, 0.1)*millisec;
  61 time21s = hslider("time", 0, 0, 21000, 0.1)*millisec;
  62 time43s = hslider("time", 0, 0, 43000, 0.1)*millisec;
  63
  64
  65 echo1s  = vgroup("echo  1000", +~(delay(65536,   int(hslider("millisecond", 0, 0,       1000, 0.10)*millisec)-1) * (hslider("feedback", 0, 0,  100, 0.1)/100.0)));
  66 echo2s  = vgroup("echo  2000", +~(delay(131072,  int(hslider("millisecond", 0, 0,       2000, 0.25)*millisec)-1) * (hslider("feedback", 0, 0,  100, 0.1)/100.0)));
  67 echo5s  = vgroup("echo  5000", +~(delay(262144,  int(hslider("millisecond", 0, 0,       5000, 0.50)*millisec)-1) * (hslider("feedback", 0, 0,  100, 0.1)/100.0)));
  68 echo10s = vgroup("echo 10000", +~(delay(524288,  int(hslider("millisecond", 0, 0,  10000, 1.00)*millisec)-1) * (hslider("feedback", 0, 0,  100, 0.1)/100.0)));
  69 echo21s = vgroup("echo 21000", +~(delay(1048576, int(hslider("millisecond", 0, 0,  21000, 1.00)*millisec)-1) * (hslider("feedback", 0, 0,  100, 0.1)/100.0)));
  70 echo43s = vgroup("echo 43000", +~(delay(2097152, int(hslider("millisecond", 0, 0,  43000, 1.00)*millisec)-1) * (hslider("feedback", 0, 0,  100, 0.1)/100.0)));
  71
  72
  73 //--------------------------sdelay(N,it,dt)----------------------------
  74 // s(mooth)delay : a mono delay that doesn't click and doesn't
  75 // transpose when the delay time is changed. It takes 4 input signals
  76 // and produces a delayed output signal
  77 //
  78 // USAGE :      ... : sdelay(N,it,dt) : ...
  79 //
  80 // Where :
  81 //      <N>  = maximal delay in samples (must be a constant power of 2, for example 65536)
  82 //      <it> = interpolation time (in samples) for example 1024
  83 //      <dt> = delay time (in samples)
  84 //  <  > = input signal we want to delay
  85 //--------------------------------------------------------------------------
  86
  87 sdelay(N, it, dt) = ctrl(it,dt),_ : ddi(N)
  88
  89         with {
  90
  91                 //---------------------------ddi(N,i,d0,d1)-------------------------------
  92                 //      DDI (Double Delay with Interpolation) : the input signal is sent to two
  93                 //      delay lines. The outputs of these delay lines are crossfaded with
  94                 //      an interpolation stage. By acting on this interpolation value one
  95                 //      can move smoothly from one delay to another. When <i> is 0 we can
  96                 //      freely change the delay time <d1> of line 1, when it is 1 we can freely change
  97                 //      the delay time <d0> of line 0.
  98                 //
  99                 //      <N>  = maximal delay in samples (must be a power of 2, for example 65536)
 100                 //      <i>  = interpolation value between 0 and 1 used to crossfade the outputs of the
 101                 //                 two delay lines (0.0: first delay line, 1.0: second delay line)
 102                 //      <d0> = delay time of delay line 0 in samples between 0 and <N>-1
 103                 //      <d1> = delay time of delay line 1 in samples between 0 and <N>-1
 104                 //  <  > = the input signal we want to delay
 105                 //-------------------------------------------------------------------------
 106                 ddi(N, i, d0, d1) = _ <: delay(N,d0), delay(N,d1) : interpolate(i);
 107
 108
 109                 //----------------------------ctrl(it,dt)------------------------------------
 110                 //      Control logic for a Double Delay with Interpolation according to two
 111                 //
 112                 //      USAGE : ctrl(it,dt)
 113                 //  where :
 114                 //      <it> an interpolation time (in samples, for example 256)
 115                 //      <dt> a delay time (in samples)
 116                 //
 117                 //      ctrl produces 3 outputs : an interpolation value <i> and two delay
 118                 //      times <d0> and <d1>. These signals are used to control a ddi (Double Delay with Interpolation).
 119                 //      The principle is to detect changes in the input delay time dt, then to
 120                 //      change the delay time of the delay line currently unused and then by a
 121                 //      smooth crossfade to remove the first delay line and activate the second one.
 122                 //
 123                 //      The control logic has an internal state controlled by 4 elements
 124                 //      <v> : the interpolation variation (0, 1/it, -1/it)
 125                 //      <i> : the interpolation value (between 0 and 1)
 126                 //      <d0>: the delay time of line 0
 127                 //      <d1>: the delay time of line 1
 128                 //
 129                 //      Please note that the last stage (!,_,_,_) cut <v> because it is only
 130                 //      used internally.
 131                 //-------------------------------------------------------------------------
 132                 ctrl(it, dt) = \(v,ip,d0,d1).( (nv, nip, nd0, nd1)
 133                         with {
 134
 135                                 // interpolation variation
 136                                 nv = if (v!=0.0,                                                        // if variation we are interpolating
 137                                                 if( (ip>0.0) & (ip<1.0), v , 0),        //              should we continue or not ?
 138                                          if ((ip==0.0) & (dt!=d0),  1.0/it,             // if true xfade from dl0 to dl1
 139                                          if ((ip==1.0) & (dt!=d1), -1.0/it,             // if true xfade from dl1 to dl0
 140                                          0)));                                                                  // nothing to change
 141                                 // interpolation value
 142                                 nip = ip+nv : min(1.0) : max(0.0);
 143
 144                                 // update delay time of line 0 if needed
 145                                 nd0 = if ((ip >= 1.0) & (d1!=dt), dt, d0);
 146
 147                                 // update delay time of line 0 if needed
 148                                 nd1 = if ((ip <= 0.0) & (d0!=dt), dt, d1);
 149
 150                         } ) ~ (_,_,_,_) : (!,_,_,_);
 151         };
 152
 153
 154
 155
 156 //-----------------------------------------------
 157 //                      Tempo, beats and pulses
 158 //-----------------------------------------------
 159
 160 tempo(t)        = (60*SR)/t;                    // tempo(t) -> samples
 161
 162 period(p)       = %(int(p))~+(1);               // signal en dent de scie de periode p
 163 pulse(t)        = period(t)==0;                 // pulse (10000...) de periode p
 164 pulsen(n,t) = period(t)<n;                      // pulse (1110000...) de taille n et de periode p
 165 beat(t)         = pulse(tempo(t));              // pulse au tempo t
 166
 167
 168
 169 //-----------------------------------------------
 170 //      conversions between db and linear values
 171 //-----------------------------------------------
 172
 173 db2linear(x)    = pow(10, x/20.0);
 174 linear2db(x)    = 20*log10(x);
 175
 176
 177
 178
 179
 180 //===============================================
 181 //                      Random and Noise generators
 182 //===============================================
 183
 184
 185 //-----------------------------------------------
 186 //                      noise : Noise generator
 187 //-----------------------------------------------
 188
 189 random          = +(12345) ~ *(1103515245);
 190 RANDMAX         = 2147483647.0;
 191
 192 noise           = random / RANDMAX;
 193
 194
 195 //-----------------------------------------------
 196 // Generates multiple decorrelated random numbers
 197 // in parallel. Expects n>0.
 198 //-----------------------------------------------
 199
 200 multirandom(n) = randomize(n) ~_
 201 with {
 202         randomize (1)   = +(12345) : *(1103515245);
 203         randomize (n)   = randomize(1) <: randomize(n-1),_;
 204 };
 205
 206
 207 //-----------------------------------------------
 208 // Generates multiple decorrelated noises
 209 // in parallel. Expects n>0.
 210 //-----------------------------------------------
 211
 212 multinoise(n) = multirandom(n) : par(i,n,/(RANDMAX))
 213 with {
 214         RANDMAX = 2147483647.0;
 215 };
 216
 217
 218 //------------------------------------------------
 219
 220 noises(N,i) = multinoise(N) : selector(i,N);
 221
 222
 223 //-----------------------------------------------
 224 //                      osc(freq) : Sinusoidal Oscillator
 225 //-----------------------------------------------
 226
 227 tablesize       = 1 << 16;
 228 samplingfreq    = SR;
 229
 230 time            = (+(1)~_ ) - 1;                        // 0,1,2,3,...
 231 sinwaveform     = float(time)*(2.0*PI)/float(tablesize) : sin;
 232
 233 decimal(x)      = x - floor(x);
 234 phase(freq)     = freq/float(samplingfreq) : (+ : decimal) ~ _ : *(float(tablesize));
 235 osc(freq)       = rdtable(tablesize, sinwaveform, int(phase(freq)) );
 236 osci(freq)      = s1 + d * (s2 - s1)
 237                 with {
 238                         i = int(phase(freq));
 239                         d = decimal(phase(freq));
 240                         s1 = rdtable(tablesize+1,sinwaveform,i);
 241                         s2 = rdtable(tablesize+1,sinwaveform,i+1);};
 242
 243
 244 //-----------------------------------------------
 245 //                      ADSR envelop
 246 //-----------------------------------------------
 247
 248 // a,d,s,r = attack (#samples), decay (sec), sustain (percentage), release (sec)
 249 // t       = trigger signal
 250
 251 adsr(a,d,s,r,t) = env ~ (_,_) : (!,_) // the 2 'state' signals are fed back
 252 with {
 253     env (p2,y) =
 254         (t>0) & (p2|(y>=1)),          // p2 = decay-sustain phase
 255         (y + p1*u - (p2&(y>s))*v*y - p3*w*y)    // y  = envelop signal
 256         *((p3==0)|(y>=eps)) // cut off tails to prevent denormals
 257     with {
 258         p1 = (p2==0) & (t>0) & (y<1);         // p1 = attack phase
 259         p3 = (t<=0) & (y>0);                  // p3 = release phase
 260         // #samples in attack, decay, release, must be >0
 261         na = SR*a+(a==0.0); nd = SR*d+(d==0.0); nr = SR*r+(r==0.0);
 262         // correct zero sustain level
 263         z = s+(s==0.0)*db2linear(-60);
 264         // attack, decay and (-60dB) release rates
 265         u = 1/na; v = 1-pow(z, 1/nd); w = 1-1/pow(z*db2linear(60), 1/nr);
 266         // values below this threshold are considered zero in the release phase
 267         eps = db2linear(-120);
 268     };
 269 };
 270
 271
 272 //-----------------------------------------------
 273 //                      Spatialisation
 274 //-----------------------------------------------
 275
 276 panner(c) = _ <: *(1-c), *(c);
 277
 278 bus2 = _,_;
 279 bus3 = _,_,_;
 280 bus4 = _,_,_,_;
 281 bus5 = _,_,_,_,_;
 282 bus6 = _,_,_,_,_,_;
 283 bus7 = _,_,_,_,_,_,_;
 284 bus8 = _,_,_,_,_,_,_,_;
 285
 286 gain2(g) = *(g),*(g);
 287 gain3(g) = *(g),*(g),*(g);
 288 gain4(g) = *(g),*(g),*(g),*(g);
 289 gain5(g) = *(g),*(g),*(g),*(g),*(g);
 290 gain6(g) = *(g),*(g),*(g),*(g),*(g),*(g);
 291 gain7(g) = *(g),*(g),*(g),*(g),*(g),*(g),*(g);
 292 gain8(g) = *(g),*(g),*(g),*(g),*(g),*(g),*(g),*(g);
 293
 294
 295 //------------------------------------------------------
 296 //
 297 //                                          GMEM SPAT
 298 //      n-outputs spatializer
 299 //      implementation of L. Pottier
 300 //
 301 //------------------------------------------------------
 302 //
 303 //  n = number of outputs
 304 //      r = rotation (between 0 et 1)
 305 //      d = distance of the source (between 0 et 1)
 306 //
 307 //------------------------------------------------------
 308 spat(n,a,d)     = _ <: par(i, n, *( scaler(i, n, a, d) : smooth(0.9999) ))
 309         with {
 310                 scaler(i,n,a,d) = (d/2.0+0.5)
 311                                                 * sqrt( max(0.0, 1.0 - abs(fmod(a+0.5+float(n-i)/n, 1.0) - 0.5) * n * d) );
 312                 smooth(c) = *(1-c) : +~*(c);
 313         };
 314
 315
 316
 317 //--------------- Second Order Generic Transfert Function -------------------------
 318 // TF2(b0,b1,b2,a1,a2)
 319 //
 320 //---------------------------------------------------------------------------------
 321
 322 TF2(b0,b1,b2,a1,a2) = sub ~ conv2(a1,a2) : conv3(b0,b1,b2)
 323         with {
 324                 conv3(k0,k1,k2,x)       = k0*x + k1*x' + k2*x'';
 325                 conv2(k0,k1,x)          = k0*x + k1*x';
 326                 sub(x,y)                        = y-x;
 327         };
 328
 329
 330 /*************************** Break Point Functions ***************************
 331
 332 bpf is an environment (a group of related definitions) tha can be used to
 333 create break-point functions. It contains three functions :
 334   - start(x,y) to start a break-point function
 335   - end(x,y) to end a break-point function
 336   - point(x,y) to add intermediate points to a break-point function
 337
 338 A minimal break-point function must contain at least a start and an end point :
 339
 340   f = bpf.start(x0,y0) : bpf.end(x1,y1);
 341
 342 A more involved break-point function can contains any number of intermediate
 343 points
 344
 345   f = bpf.start(x0,y0) : bpf.point(x1,y1) : bpf.point(x2,y2) : bpf.end(x3,y3);
 346
 347 In any case the x_{i} must be in increasing order (for all i, x_{i} < x_{i+1})
 348
 349 For example the following definition :
 350
 351   f = bpf.start(x0,y0) : ... : bpf.point(xi,yi) : ... : bpf.end(xn,yn);
 352
 353 implements a break-point function f such that :
 354
 355   f(x) = y_{0} when x < x_{0}
 356   f(x) = y_{n} when x > x_{n}
 357   f(x) = y_{i} + (y_{i+1}-y_{i})*(x-x_{i})/(x_{i+1}-x_{i}) when x_{i} <= x and x < x_{i+1}
 358
 359 ******************************************************************************/
 360
 361 bpf = environment
 362 {
 363   // Start a break-point function
 364   start(x0,y0) = \(x).(x0,y0,x,y0);
 365
 366   // Add a break-point
 367   point(x1,y1) = \(x0,y0,x,y).(x1, y1, x , if (x < x0, y, if (x < x1, y0 + (x-x0)*(y1-y0)/(x1-x0), y1)));
 368
 369   // End a break-point function
 370   end  (x1,y1) = \(x0,y0,x,y).(if (x < x0, y, if (x < x1, y0 + (x-x0)*(y1-y0)/(x1-x0), y1)));
 371
 372   // definition of if
 373   if (c,t,e) = select2(c,e,t);
 374 };
 375
 376
 377