//SLUGens released under the GNU GPL as extensions for SuperCollider 3, by Nick Collins, http://composerprogrammer.com/index.html

DoubleWell2 Forced DoubleWell Oscillator

DoubleWell2.ar(reset,ratex,ratey,f,w,delta,initx,inity, mul, add)

Improved Euler ODE solver implementation of the chaotic Forced Double Well Oscillator (see Strogatz, Steven H. (1994) Nonlinear Dynamics and Chaos. Addison-Wesley, Reading, MA. pp441-7). 

D2x +delta*Dx-x+x^3 = F*cos(w*t)

y = Dx

All inputs can have .kr rate UGens plugged in.

reset- restart with new initial conditions sampled from initx, inity

ratex- update rate for x

ratey- update rate for y 

f- equation constant

w- equation constant

delta- equation constant

initx- reset value for x

inity- reset value for y

note for me- add a t rate control? or arbitary oscillator input rather than Fcos(wt)

//Can be very noisy, high pitched and aliases badly- ie great fun, but watch out for your ears, filter, try out odd param settings. I've turned all the amplitudes down in the following

//defaults, high pitched! WARNING

{Out.ar(0,Pan2.ar(0.1*DoubleWell.ar(Impulse.kr(0),MouseX.kr(0.01,1,'exponential'),MouseY.kr(0.01,1,'exponential'),0.9,0.0001,0.25,LFNoise0.kr(1,0.4,0.5),LFNoise0.kr(1,0.4,0.5)),0.0))}.play

//2 per second

(

{Out.ar(0,Pan2.ar(

CombN.ar(

Resonz.ar(DoubleWell2.ar(Impulse.kr(2),MouseX.kr(0.01,0.1,'exponential'),MouseY.kr(0.01,0.1,'exponential'),3.9,MouseX.kr(0.0001,0.01,'exponential'),0.25,LFNoise0.kr(1,0.4,0.5),LFNoise0.kr(1,10.4,10.5)),1000,0.1),

0.05,0.05,2

)

,0.0))}.play

)

//plotting check of first samples

{DoubleWell2.ar(0,0.1,0.1,0.4,1,0.25,0,0.0)}.plot(0.05, minval:-1, maxval:1)

//gradually changing force

{DoubleWell2.ar(0,0.1,0.1,Line.kr(0.01,0.8,10),1,0.25,0,0.0)}.play

//controlling force- guitar feedback!

{DoubleWell2.ar(0,0.1,0.1,MouseX.kr(0.0,10.0),1,MouseY.kr(0.0,1.0),0,0.0)}.play

//triggering

{DoubleWell2.ar(Impulse.kr(MouseY.kr(1,10, 'exponential')),0.1,0.1,MouseX.kr(0.0,10.0),1,MouseY.kr(0.0,1.0),LFNoise0.kr(2.3,0.2,0.8),LFNoise1.kr(2.37,0.2,0.8))}.play

//overload patch

{DoubleWell2.ar(Dust.kr(LFNoise0.kr(0.5,0.2,0.4)),0.1,0.1,MouseX.kr(0.0,50.0),MouseY.kr(0.0,1.0),0.25, LFNoise0.kr(2.3,0.4,0.5),LFNoise1.kr(2.37,0.4,0.2))}.play