Apat and added acos to spherycolor, not included in color yet but working
[SugarCubes.git] / DanHorwitz.pde
old mode 100755 (executable)
new mode 100644 (file)
index d73b6c4..ab23fa5
-//----------------------------------------------------------------------------------------------------------------------------------\r
-public class Pong extends DPat {\r
-       SinLFO x,y,z,dx,dy,dz; \r
-       float cRad;     DParam pSize;\r
-       Pick  pChoose;\r
-\r
-       Pong(GLucose glucose) {\r
-               super(glucose);\r
-               cRad = xdMax/15;\r
-               addModulator(dx = new SinLFO(6000,  500, 30000  )).trigger();\r
-               addModulator(dy = new SinLFO(3000,  500, 22472  )).trigger();\r
-               addModulator(dz = new SinLFO(1000,  500, 18420  )).trigger();\r
-               addModulator(x  = new SinLFO(cRad, xdMax - cRad, 0)).trigger(); x.modulateDurationBy(dx);\r
-               addModulator(y  = new SinLFO(cRad, ydMax - cRad, 0)).trigger(); y.modulateDurationBy(dy);\r
-               addModulator(z  = new SinLFO(cRad, zdMax - cRad, 0)).trigger(); z.modulateDurationBy(dz);\r
-           pSize       = addParam      ("Size"                 , 0.4   );\r
-           pChoose = addPick   ("Animiation"   , 0     , 3, new String[] {"Pong", "Ball", "Cone"}      );\r
-       }\r
-\r
-       void    StartRun(double deltaMs)        { cRad = xdMax*pSize.Val()/6; }\r
-       color   CalcPoint(xyz p)                {\r
-               xyz v = new xyz(x.getValuef(), y.getValuef(), z.getValuef());\r
-               switch(pChoose.Cur()) {\r
-               case 0: return color(0,0,c1c(1 - min(v.distance(p), v.distance(xyzdMax.minus(p)))*.5/cRad));                    // balls\r
-               case 1: return color(0,0,c1c(1 - v.distance(p)*.5/cRad));                                                                                               // ball\r
-               case 2: return color(0,0,c1c(1 - CalcCone(p,v,new xyz(xdMax/2,0,zdMax/2)) * max(.02,.45-pSize.Val())));  // spot\r
-               }\r
-               return color(0,0,0);\r
-       }               \r
-}\r
-//----------------------------------------------------------------------------------------------------------------------------------\r
-public class NDat {\r
-       float   xz, yz, zz, hue, sat, speed, angle, den;\r
-       float   xoff,yoff,zoff;\r
-       float   sinAngle, cosAngle;\r
-       boolean isActive;\r
-       NDat              () { isActive=false; }\r
-       boolean Active() { return isActive; }\r
-       void    set     (float _hue, float _sat, float _xz, float _yz, float _zz, float _den, float _speed, float _angle) {\r
-               isActive = true;\r
-               hue=_hue; sat=_sat; xz=_xz; yz=_yz; zz =_zz; den=_den; speed=_speed; angle=_angle;\r
-               xoff = random(100e3); yoff = random(100e3); zoff = random(100e3);\r
-       }\r
-}\r
-\r
-public class Noise extends DPat\r
-{\r
-       int                     CurAnim, iSymm;\r
-       int             XSym=1,YSym=2,RadSym=3;\r
-       float           zTime   = random(10000), zTheta=0, zSin, zCos;\r
-       float           rtime   = 0, ttime      = 0, transAdd=0;\r
-       DParam          pSpeed , pDensity, pRotZ;\r
-       Pick            pChoose, pSymm;\r
-       int                     _ND = 4;\r
-       NDat            N[] = new NDat[_ND];\r
-\r
-       Noise(GLucose glucose) {\r
-               super(glucose);\r
-               pRotZ   = addParam("RotZ"        , .5 );        pSpeed          = addParam("Fast", .55);\r
-               pDensity= addParam("Dens"        , .5);\r
-               pSymm   = addPick("Symmetry" , 0, 4, new String[] {"None", "X", "Y", "Radial"}  );\r
-               pChoose = addPick("Animation", 1, 6, new String[] {"Drip", "Cloud", "Rain", "Fire", "Machine", "Spark"} );\r
-               for (int i=0; i<_ND; i++) N[i] = new NDat();\r
-       }\r
-\r
-       void StartRun(double deltaMs) {\r
-               zTime   += deltaMs*(pSpeed.Val()-.5)*.002       ;\r
-               zTheta  += deltaMs*(pRotZ .Val()-.5)*.01        ;\r
-               rtime   += deltaMs;\r
-               iSymm    = pSymm.Cur();\r
-               transAdd = 1*(1 - constrain(rtime - ttime,0,1000)/1000);\r
-               zSin    = sin(zTheta);\r
-               zCos    = cos(zTheta);\r
-\r
-               if (pChoose.Cur() != CurAnim) {\r
-                       CurAnim = pChoose.Cur(); ttime = rtime;\r
-                       pRotZ           .reset();       zTheta          = 0;\r
-                       pDensity        .reset();       pSpeed          .reset();       \r
-                       for (int i=0; i<_ND; i++) { N[i].isActive = false; }\r
-                       \r
-                       switch(CurAnim) {\r
-                       //                          hue sat xz  yz  zz  den mph angle\r
-                       case 0: N[0].set(0  ,0  ,75 ,75 ,150,45 ,3  ,0  ); pSharp.Set(1 ); break;       // drip\r
-                       case 1: N[0].set(0  ,0  ,100,100,200,45 ,3  ,180); pSharp.Set(0 ); break;       // clouds\r
-                       case 2: N[0].set(0  ,0  ,2  ,400,2  ,20 ,3  ,0  ); pSharp.Set(.5); break;       // rain\r
-                       case 3: N[0].set(40 ,1  ,100,100,200,10 ,1  ,180); \r
-                                       N[1].set(0  ,1  ,100,100,200,10 ,5  ,180); pSharp.Set(0 ); break;       // fire 1\r
-                       case 4: N[0].set(0  ,1  ,40 ,40 ,40 ,15 ,2.5,180);\r
-                                       N[1].set(20 ,1  ,40 ,40 ,40 ,15 ,4  ,0  );\r
-                                       N[2].set(40 ,1  ,40 ,40 ,40 ,15 ,2  ,90 );\r
-                                       N[3].set(60 ,1  ,40 ,40 ,40 ,15 ,3  ,-90); pSharp.Set(.5); break; // machine\r
-                       case 5: N[0].set(0  ,1  ,400,100,2  ,15 ,3  ,90 );\r
-                                       N[1].set(20 ,1  ,400,100,2  ,15 ,2.5,0  );\r
-                                       N[2].set(40 ,1  ,100,100,2  ,15 ,2  ,180);\r
-                                       N[3].set(60 ,1  ,100,100,2  ,15 ,1.5,270); pSharp.Set(.5); break; // spark\r
-                       }\r
-\r
-                       DG.UpdateLights();\r
-               }\r
-               \r
-               for (int i=0; i<_ND; i++) if (N[i].Active()) {\r
-                       N[i].sinAngle = sin(radians(N[i].angle));\r
-                       N[i].cosAngle = cos(radians(N[i].angle));\r
-               }\r
-       }\r
-\r
-       color CalcPoint(xyz P) {\r
-               color c = 0;\r
-               P.RotateZ(xyzMid, zSin, zCos);\r
-               \r
-               if (iSymm == XSym && P.x > xdMax/2) P.x = xdMax-P.x;\r
-               if (iSymm == YSym && P.y > ydMax/2) P.y = ydMax-P.y;\r
-\r
-               for (int i=0;i<_ND; i++) if (N[i].Active()) {\r
-                       NDat  n     = N[i];\r
-                       float zx    = zTime * n.speed * n.sinAngle,\r
-                                 zy    = zTime * n.speed * n.cosAngle;\r
-                       \r
-                       float b     = (iSymm==RadSym ? noise(zTime*n.speed+n.xoff-Dist(P,xyzMid)/n.xz)\r
-                                                                                : noise(P.x/n.xz+zx+n.xoff,P.y/n.yz+zy+n.yoff,P.z/n.zz+n.zoff))\r
-                                                       *1.8;\r
-\r
-                       b +=    n.den/100 -.4 + pDensity.Val() -1;\r
-                       b +=    transAdd;\r
-                       c =     blendColor(c,color(n.hue,100*n.sat,c1c(b)),ADD);\r
-               }\r
-               return c;\r
-       }\r
-}\r
-//----------------------------------------------------------------------------------------------------------------------------------\r
-public class Play extends DPat\r
-{\r
-       int             nBeats  =       0;\r
-       DParam  pAmp, pRad;\r
-       DParam  pRotX, pRotY, pRotZ;\r
-       xyz             Theta   = new xyz();\r
-       xyz             TSin    = new xyz();\r
-       xyz             TCos    = new xyz();\r
-       \r
-       Pick    pTimePattern, pTempoMult, pShape, pForm;\r
-       int             RandCube;\r
-\r
-       Play(GLucose glucose) {\r
-               super(glucose);\r
-               pRotX           = addParam("RotX", .5);\r
-               pRotY           = addParam("RotY", .5);\r
-               pRotZ           = addParam("RotZ", .5);\r
-           pAmp                = addParam("Amp" , .2);\r
-           pRad                = addParam("Rad"        , .1            );\r
-               pTempoMult      = addPick ("TMult"      , 5 , 6         , new String[] {"1x", "2x", "4x", "8x", "16x", "Rand"   }       );\r
-               pTimePattern= addPick ("TPat"   , 5 , 6         , new String[] {"Bounce", "Sin", "Roll", "Quant", "Accel", "Rand"       }       );\r
-               pShape          = addPick ("Shape"      , 0 , 10        , new String[] {"Line", "Tap", "V", "RandV", "Pyramid",\r
-                                                                                                                                       "Wings", "W2", "Sphere", "Cone", "Noise" }      );\r
-               pForm           = addPick ("Form"       , 0 , 3         , new String[] {"Bar", "Volume", "Fade"                                 }       );\r
-       }\r
-\r
-       float   t,a;\r
-       xyz             cPrev = new xyz(), cCur = new xyz(), cMid = new xyz(), cMidNorm;\r
-       float   LastBeat=3, LastMeasure=3;\r
-       int             CurRandTempo = 1, CurRandTPat = 1;\r
-\r
-       void StartRun(double deltaMs) {\r
-               t = lx.tempo.rampf();\r
-               a = pAmp.Val();\r
-\r
-               Theta   .set(pRotX.Val()*PI*2, pRotY.Val()*PI*2, pRotZ.Val()*PI*2);\r
-               TSin    .set(sin(Theta.x), sin(Theta.y), sin(Theta.z));\r
-               TCos    .set(cos(Theta.x), cos(Theta.y), cos(Theta.z));\r
-\r
-               if (t<LastMeasure) {\r
-                       if (random(2) < 1) CurRandTempo = int(random(4));\r
-                       if (random(2) < 1) CurRandTPat  = int(random(5));\r
-               } LastMeasure = t;\r
-                       \r
-               int nTempo = pTempoMult  .Cur(); if (nTempo == 5) nTempo = CurRandTempo;\r
-               int nTPat  = pTimePattern.Cur(); if (nTPat  == 5) nTPat  = CurRandTPat ;\r
-\r
-               switch (nTempo) {\r
-                       case 0:         t = t;                                                          break;\r
-                       case 1:         t = (t*2. )%1.;                                         break;\r
-                       case 2:         t = (t*4. )%1.;                                         break;\r
-                       case 3:         t = (t*8. )%1.;                                         break;\r
-                       case 4:         t = (t*16.)%1.;                                         break;\r
-               }\r
-\r
-               if (t<LastBeat) { cPrev = cCur; cCur = cCur.setRand(); } LastBeat = t;\r
-\r
-               switch (nTPat) {\r
-                       case 0:         t = sin(PI*t);                                          break;\r
-                       case 1:         t = norm(sin(2*PI*(t+PI/2)),-1,1);      break;\r
-                       case 2:         t = t;                                                          break;\r
-                       case 3:         t = constrain(int(t*8)/7.,0,1);         break;\r
-                       case 4:         t = t*t*t;                                                      break;\r
-               }\r
-\r
-               \r
-               cMid                            = cPrev.interpolate(t,cCur);\r
-               cMidNorm                        = cMid.setNorm();\r
-       }\r
-\r
-       color CalcPoint(xyz Px) {\r
-               xyz V           =                       new xyz();\r
-               xyz P           =                       Px.setNorm();\r
-               if (Theta.x != 0)               P.RotateX(xyzHalf, TSin.x, TCos.x);\r
-               if (Theta.y != 0)               P.RotateY(xyzHalf, TSin.y, TCos.y);\r
-               if (Theta.z != 0)               P.RotateZ(xyzHalf, TSin.z, TCos.z);\r
-\r
-               float mp        = min(P.x, P.z);\r
-               float yt        = map(t,0,1,.5-a/2,.5+a/2);\r
-\r
-               switch (pShape.Cur()) {\r
-                       case 0:         V.set(P.x, yt                                                     , P.z);                                                       break;  // bouncing line\r
-                       case 1:         V.set(P.x, map(cos(PI*t * P.x),-1,1,0,1)  , P.z);                                                       break;  // top tap\r
-                       case 2:         V.set(P.x, a*map(P.x<.5?P.x:1-P.x,0,.5 ,0,t-.5)+.5, P.z);                                       break;  // V shape\r
-                       case 3:         V.set(P.x, P.x < cMidNorm.x ? map(P.x,0,cMidNorm.x, .5,yt) :\r
-                                                                                                         map(P.x,cMidNorm.x,1, yt,.5), P.z);                   break;  //  Random V shape\r
-\r
-                       case 4:         V.set(P.x,      .5*(P.x < cMidNorm.x ?  map(P.x,0,cMidNorm.x, .5,yt) :\r
-                                                                                                                       map(P.x,cMidNorm.x,1, yt,.5)) +\r
-                                                                       .5*(P.z < cMidNorm.z ?  map(P.z,0,cMidNorm.z, .5,yt) :\r
-                                                                                                                       map(P.z,cMidNorm.z,1, yt,.5)), P.z);    break;  //  Random Pyramid shape\r
-                                                                                                                       \r
-                       case 5:         V.set(P.x, a*map((P.x-.5)*(P.x-.5),0,.25,0,t-.5)+.5, P.z);                                      break;  // wings\r
-                       case 6:         V.set(P.x, a*map((mp -.5)*(mp -.5),0,.25,0,t-.5)+.5, P.z);                                      break;  // wings\r
-\r
-                       case 7:         V.set(cMid.x,cMid.y,cMid.z);\r
-                                               return color(0,0,c1c(1 - (V.distance(Px) > (pRad.getValuef()+.1)*150?1:0)) );           // sphere\r
-\r
-                       case 8:         V.set(cMid.x,cMid.y,cMid.z);\r
-                                               return color(0,0,c1c(1 - CalcCone(Px,V,xyzMid) * 0.02 > .5?1:0));                                       // cone\r
-\r
-                       case 9:         return color(100 + noise(P.x,P.y,P.z + (NoiseMove+50000)/1000.)*200,\r
-                                                       85,c1c(P.y < noise(P.x + NoiseMove/2000.,P.z)*(1+a)-a/2.-.1 ? 1 : 0));                  //\r
-               }\r
-\r
-\r
-               switch (pForm.Cur()) {\r
-                       case 0:         return color(0,0,c1c(1 - V.distance(P)/pRad.getValuef() > .5?1:0));\r
-                       case 1:         return color(0,0,c1c(P.y < V.y ?1:0));\r
-                       case 2:         return color(0,0,c1c(1 - V.distance(P)/pRad.getValuef()));\r
-\r
-                       default:        return color(0,0,c1c(P.y < V.y ?1:0));\r
-               }\r
-       }\r
-}\r
-//----------------------------------------------------------------------------------------------------------------------------------\r
+//----------------------------------------------------------------------------------------------------------------------------------
+public class Pong extends DPat {
+       SinLFO x,y,z,dx,dy,dz;
+       float cRad;     BasicParameter pSize;
+       Pick    pChoose;
+       PVector v = new PVector(), vMir =  new PVector();
+
+       Pong(GLucose glucose) {
+               super(glucose);
+               cRad = mMax.x/10;
+               addModulator(dx = new SinLFO(6000,  500, 30000  )).trigger();
+               addModulator(dy = new SinLFO(3000,  500, 22472  )).trigger();
+               addModulator(dz = new SinLFO(1000,  500, 18420  )).trigger();
+               addModulator(x  = new SinLFO(cRad, mMax.x - cRad, 0)).trigger();        x.modulateDurationBy(dx);
+               addModulator(y  = new SinLFO(cRad, mMax.y - cRad, 0)).trigger();        y.modulateDurationBy(dy);
+               addModulator(z  = new SinLFO(cRad, mMax.z - cRad, 0)).trigger();        z.modulateDurationBy(dz);
+           pSize       = addParam      ("Size"                 , 0.4   );
+           pChoose = addPick   ("Animiation"   , 2, 2, new String[] {"Pong", "Ball", "Cone"}   );
+       }
+
+       void    StartRun(double deltaMs)        { cRad = mMax.x*val(pSize)/6; }
+       color   CalcPoint(PVector p)            {
+               v.set(x.getValuef(), y.getValuef(), z.getValuef());
+               v.z=0;p.z=0;// ignore z dimension
+               switch(pChoose.Cur()) {
+               case 0: vMir.set(mMax); vMir.sub(p);
+                               return lx.hsb(lxh(),100,c1c(1 - min(v.dist(p), v.dist(vMir))*.5/cRad));         // balls
+               case 1: return lx.hsb(lxh(),100,c1c(1 - v.dist(p)*.5/cRad));                                                    // ball
+               case 2: vMir.set(mMax.x/2,0,mMax.z/2);
+                               return lx.hsb(lxh(),100,c1c(1 - calcCone(p,v,vMir) * max(.02,.45-val(pSize))));         // spot
+               }
+               return lx.hsb(0,0,0);
+       }
+}
+//----------------------------------------------------------------------------------------------------------------------------------
+public class NDat {
+       float   xz, yz, zz, hue, speed, angle, den;
+       float   xoff,yoff,zoff;
+       float   sinAngle, cosAngle;
+       boolean isActive;
+       NDat              () { isActive=false; }
+       boolean Active() { return isActive; }
+       void    set     (float _hue, float _xz, float _yz, float _zz, float _den, float _speed, float _angle) {
+               isActive = true;
+               hue=_hue; xz=_xz; yz=_yz; zz =_zz; den=_den; speed=_speed; angle=_angle;
+               xoff = random(100e3); yoff = random(100e3); zoff = random(100e3);
+       }
+}
+
+public class Noise extends DPat
+{
+       int                             CurAnim, iSymm;
+       int                     XSym=1,YSym=2,RadSym=3;
+       float                   zTime , zTheta=0, zSin, zCos, rtime, ttime;
+       BasicParameter  pSpeed , pDensity, pSharp;
+       Pick                    pChoose, pSymm;
+       int                             _ND = 4;
+       NDat                    N[] = new NDat[_ND];
+
+       Noise(GLucose glucose) {
+               super(glucose);
+               pSpeed          = addParam("Fast"       , .55);
+               pDensity        = addParam("Dens"        , .5);
+               pSharp          = addParam("Shrp"        ,  0);
+               pSymm           = addPick("Symmetry" , 0, 3, new String[] {"None", "X", "Y", "Radial"}  );
+               pChoose         = addPick("Animation", 6, 7, new String[] {"Drip", "Cloud", "Rain", "Fire", "Machine", "Spark","VWave", "Wave"} );
+               for (int i=0; i<_ND; i++) N[i] = new NDat();
+       }
+
+       void onActive() { zTime = random(500); zTheta=0; rtime = 0; ttime = 0; }
+
+       void StartRun(double deltaMs) {
+               zTime   += deltaMs*(val(pSpeed)-.5)*.002        ;
+               zTheta  += deltaMs*(spin()-.5)*.01      ;
+               rtime   += deltaMs;
+               iSymm    = pSymm.Cur();
+               zSin    = sin(zTheta);
+               zCos    = cos(zTheta);
+
+               if (pChoose.Cur() != CurAnim) {
+                       CurAnim = pChoose.Cur(); ttime = rtime;
+                       pSpin           .reset();       zTheta          = 0;
+                       pDensity        .reset();       pSpeed          .reset();
+                       for (int i=0; i<_ND; i++) { N[i].isActive = false; }
+                       
+                       switch(CurAnim) {
+                       //                          hue xz yz zz den mph angle
+                       case 0: N[0].set(0  ,75 ,75 ,150,45 ,3  ,0  ); pSharp.setValue(1 ); break;      // drip
+                       case 1: N[0].set(0  ,100,100,200,45 ,3  ,180); pSharp.setValue(0 ); break;      // clouds
+                       case 2: N[0].set(0  ,2  ,400,2  ,20 ,3  ,0  ); pSharp.setValue(.5); break;      // rain
+                       case 3: N[0].set(40 ,100,100,200,10 ,1  ,180); 
+                                       N[1].set(0  ,100,100,200,10 ,5  ,180); pSharp.setValue(0 ); break;      // fire 1
+                       case 4: N[0].set(0  ,40 ,40 ,40 ,15 ,2.5,180);
+                                       N[1].set(20 ,40 ,40 ,40 ,15 ,4  ,0  );
+                                       N[2].set(40 ,40 ,40 ,40 ,15 ,2  ,90 );
+                                       N[3].set(60 ,40 ,40 ,40 ,15 ,3  ,-90); pSharp.setValue(.5); break; // machine
+                       case 5: N[0].set(0  ,400,100,2  ,15 ,3  ,90 );
+                                       N[1].set(20 ,400,100,2  ,15 ,2.5,0  );
+                                       N[2].set(40 ,100,100,2  ,15 ,2  ,180);
+                                       N[3].set(60 ,100,100,2  ,15 ,1.5,270); pSharp.setValue(.5); break; // spark
+                       }
+               }
+               
+               for (int i=0; i<_ND; i++) if (N[i].Active()) {
+                       N[i].sinAngle = sin(radians(N[i].angle));
+                       N[i].cosAngle = cos(radians(N[i].angle));
+               }
+       }
+
+       color CalcPoint(PVector p) {
+               color c = 0;
+               rotateZ(p, mCtr, zSin, zCos);
+
+               if (CurAnim == 6 || CurAnim == 7) {
+                       setNorm(p);
+                       return lx.hsb(lxh(),100, 100 * (
+                                                       constrain(1-50*(1-val(pDensity))*abs(p.y-sin(zTime*10  + p.x*(300))*.5 - .5),0,1) + 
+                       (CurAnim == 7 ? constrain(1-50*(1-val(pDensity))*abs(p.x-sin(zTime*10  + p.y*(300))*.5 - .5),0,1) : 0))
+                       );
+               }                       
+
+               if (iSymm == XSym && p.x > mMax.x/2) p.x = mMax.x-p.x;
+               if (iSymm == YSym && p.y > mMax.y/2) p.y = mMax.y-p.y;
+
+               for (int i=0;i<_ND; i++) if (N[i].Active()) {
+                       NDat  n     = N[i];
+                       float zx    = zTime * n.speed * n.sinAngle,
+                                 zy    = zTime * n.speed * n.cosAngle;
+
+                       float b     = (iSymm==RadSym ? noise(zTime*n.speed+n.xoff-p.dist(mCtr)/n.xz)
+                                                                                : noise(p.x/n.xz+zx+n.xoff,p.y/n.yz+zy+n.yoff,p.z/n.zz+n.zoff))
+                                                       *1.8;
+
+                       b +=    n.den/100 -.4 + val(pDensity) -1;
+                       c =     blendColor(c,lx.hsb(lxh()+n.hue,100,c1c(b)),ADD);
+               }
+               return c;
+       }
+}
+//----------------------------------------------------------------------------------------------------------------------------------
+public class Play extends DPat
+{
+       public class rAngle {
+               float   prvA, dstA, c;
+               float   prvR, dstR, r;          
+               float   _cos, _sin, x, y;
+               float   fixAngle        (float a, float b) { return a<b ?
+                                                                               (abs(a-b) > abs(a+2*PI-b) ? a : a+2*PI) :
+                                                                               (abs(a-b) > abs(a-2*PI-b) ? a : a-2*PI) ; }
+               float   getX(float r)   {       return mCtr.x + _cos*r; }
+               float   getY(float r)   {       return mCtr.y + _sin*r; }
+               void    move()                  {       c               = interp(t,prvA,dstA); 
+                                                                       r               = interp(t,prvR,dstR);
+                                                                       _cos    = cos(c);       _sin    = sin(c);
+                                      x                = getX(r);      y               = getY(r);              }               
+               void    set()                   {       prvA    = dstA;         dstA    = random(2*PI);         prvA = fixAngle(prvA, dstA);
+                                                                       prvR    = dstR;         dstR    = random(mCtr.y);                                                                       }
+       }
+
+       BasicParameter  pAmp, pRadius, pBounce;
+       Pick                    pTimePattern, pTempoMult, pShape;
+
+       ArrayList<rWave> waves = new ArrayList<rWave>(10);
+
+       int             nBeats  =       0;
+       float   t,amp,rad,bnc,zTheta=0;
+
+       rAngle  a1              = new rAngle(), a2                      = new rAngle(),
+                       a3              = new rAngle(), a4                      = new rAngle();
+       PVector cPrev   = new PVector(), cRand          = new PVector(),
+                       cMid    = new PVector(), V                      = new PVector(),
+                       theta   = new PVector(), tSin           = new PVector(),
+                       tCos    = new PVector(), cMidNorm       = new PVector(),
+                       Pn              = new PVector();
+       float   LastBeat=3, LastMeasure=3;
+       int             curRandTempo = 1, curRandTPat = 1;
+
+       Play(GLucose glucose) {
+               super(glucose);
+           pRadius             = addParam("Rad"        , .1    );
+               pBounce         = addParam("Bnc"        , .2    );
+           pAmp                = addParam("Amp"        , .2    );
+               pTempoMult      = addPick ("TMult"      , 5 , 5         , new String[] {"1x", "2x", "4x", "8x", "16x", "Rand"   }       );
+               pTimePattern= addPick ("TPat"   , 7 , 7         , new String[] {"Bounce", "Sin", "Roll", "Quant", "Accel", "Deccel", "Slide", "Rand"}   );
+               pShape          = addPick ("Shape"      , 7 , 15        , new String[] {"Line", "Tap", "V", "RandV",
+                                                                                                                                       "Pyramid", "Wings", "W2", "Clock",
+                                                                                                                                       "Triangle", "Quad", "Sphere", "Cone",
+                                                                                                                                       "Noise", "Wave", "?", "?"}                                              );
+       }
+
+       public class rWave {
+               float v0, a0, x0, t,damp,a;
+               boolean bDone=false;
+               final float len=8;
+               rWave(float _x0, float _a0, float _v0, float _damp) { x0=_x0*len; a0=_a0; v0=_v0; t=0; damp = _damp; }
+               void move(double deltaMs) {
+                       t += deltaMs*.001;
+                       if (t>4) bDone=true;
+               }
+               float val(float _x) {
+                       _x*=len;
+                       float dist = t*v0 - abs(_x-x0);
+                       if (dist<0) { a=1; return 0; }
+                       a  = a0*exp(-dist*damp) * exp(-abs(_x-x0)/(.2*len)); // * max(0,1-t/dur)
+                       return  -a*sin(dist);
+               }
+       }
+
+       void onReset()  { zTheta=0; super.onReset(); }
+       void onActive() { 
+               zTheta=0; 
+               while (lx.tempo.bpm() > 40) lx.tempo.setBpm(lx.tempo.bpm()/2);
+       }
+
+       int KeyPressed = -1;
+       boolean noteOn(Note note) {
+               int row = note.getPitch(), col = note.getChannel();
+               if (row == 57) {KeyPressed = col; return true; }
+               return super.noteOn(note);
+       }
+
+       void StartRun(double deltaMs) {
+               t       = lx.tempo.rampf();
+               amp = pAmp              .getValuef();
+               rad     = pRadius       .getValuef();
+               bnc     = pBounce       .getValuef();           
+               zTheta  += deltaMs*(val(pSpin)-.5)*.01;
+
+               theta   .set(val(pRotX)*PI*2, val(pRotY)*PI*2, val(pRotZ)*PI*2 + zTheta);
+               tSin    .set(sin(theta.x), sin(theta.y), sin(theta.z));
+               tCos    .set(cos(theta.x), cos(theta.y), cos(theta.z));
+
+               if (t<LastMeasure) {
+                       if (random(3) < 1) { curRandTempo = int(random(4)); if (curRandTempo == 3) curRandTempo = int(random(4));       }
+                       if (random(3) < 1) { curRandTPat  = pShape.Cur() > 6 ? 2+int(random(5)) : int(random(7));                                       }
+               } LastMeasure = t;
+                       
+               int nTempo = pTempoMult  .Cur(); if (nTempo == 5) nTempo = curRandTempo;
+               int nTPat  = pTimePattern.Cur(); if (nTPat  == 7) nTPat  = curRandTPat ;
+
+               switch (nTempo) {
+                       case 0:         t = t;                                                          break;
+                       case 1:         t = (t*2. )%1.;                                         break;
+                       case 2:         t = (t*4. )%1.;                                         break;
+                       case 3:         t = (t*8. )%1.;                                         break;
+                       case 4:         t = (t*16.)%1.;                                         break;
+               }
+
+               int i=0; while (i< waves.size()) {
+                       rWave w = waves.get(i);
+                       w.move(deltaMs); if (w.bDone) waves.remove(i); else i++;
+               }
+
+               if ((t<LastBeat && pShape.Cur()!=14) || KeyPressed>-1) {
+                       waves.add(new rWave(
+                                               KeyPressed>-1 ? map(KeyPressed,0,7,0,1) : random(1),            // location
+                                               bnc*10,                 // bounciness
+                                               7,                              // velocity
+                                               2*(1-amp)));    // dampiness
+                       KeyPressed=-1;
+                       if (waves.size() > 5) waves.remove(0);
+               }
+               
+               if (t<LastBeat) {
+                       cPrev.set(cRand); setRand(cRand);
+                       a1.set(); a2.set(); a3.set(); a4.set();
+               } LastBeat = t;
+
+               switch (nTPat) {
+                       case 0:         t = sin(PI*t);                                                  break;  // bounce
+                       case 1:         t = norm(sin(2*PI*(t+PI/2)),-1,1);              break;  // sin
+                       case 2:         t = t;                                                                  break;  // roll
+                       case 3:         t = constrain(int(t*8)/7.,0,1);                 break;  // quant
+                       case 4:         t = t*t*t;                                                              break;  // accel
+                       case 5:         t = sin(PI*t*.5);                                               break;  // deccel
+                       case 6:         t = .5*(1-cos(PI*t));                                   break;  // slide
+               }
+               
+               cMid.set                (cPrev);        interpolate(t,cMid,cRand);
+               cMidNorm.set    (cMid);         setNorm(cMidNorm);
+               a1.move(); a2.move(); a3.move(); a4.move();
+       }
+
+       color CalcPoint(PVector Px) {
+               if (theta.x != 0) rotateX(Px, mCtr, tSin.x, tCos.x);
+               if (theta.y != 0) rotateY(Px, mCtr, tSin.y, tCos.y);
+               if (theta.z != 0) rotateZ(Px, mCtr, tSin.z, tCos.z);
+               
+               Pn.set(Px); setNorm(Pn);
+
+               float mp        = min(Pn.x, Pn.z);
+               float yt        = map(t,0,1,.5-bnc/2,.5+bnc/2);
+               float r,d;
+
+               switch (pShape.Cur()) {
+               case 0:         V.set(Pn.x, yt                                                          , Pn.z);                                                        break;  // bouncing line
+               case 1:         V.set(Pn.x, map(cos(PI*t * Pn.x),-1,1,0,1)  , Pn.z);                                                    break;  // top tap
+               case 2:         V.set(Pn.x, bnc*map(Pn.x<.5?Pn.x:1-Pn.x,0,.5 ,0,t-.5)+.5, Pn.z);                                break;  // V shape
+               case 3:         V.set(Pn.x, Pn.x < cMidNorm.x ? map(Pn.x,0,cMidNorm.x, .5,yt) :
+                                                                                               map(Pn.x,cMidNorm.x,1, yt,.5), Pn.z);                           break;  //  Random V shape
+
+               case 4:         V.set(Pn.x,     .5*(Pn.x < cMidNorm.x ?         map(Pn.x,0,cMidNorm.x, .5,yt) :
+                                                                                                               map(Pn.x,cMidNorm.x,1, yt,.5)) +
+                                                       .5*(Pn.z < cMidNorm.z ?         map(Pn.z,0,cMidNorm.z, .5,yt) :
+                                                                                                               map(Pn.z,cMidNorm.z,1, yt,.5)), Pn.z);          break;  //  Random Pyramid shape
+                                                                                                       
+               case 5:         V.set(Pn.x, bnc*map((Pn.x-.5)*(Pn.x-.5),0,.25,0,t-.5)+.5, Pn.z);                                break;  // wings
+               case 6:         V.set(Pn.x, bnc*map((mp  -.5)*(mp  -.5),0,.25,0,t-.5)+.5, Pn.z);                                break;  // wings
+
+               case 7:         d = min(
+                                               distToSeg(Px.x, Px.y, a1.getX(70),a1.getY(70), mCtr.x, mCtr.y),
+                                               distToSeg(Px.x, Px.y, a2.getX(40),a2.getY(40), mCtr.x, mCtr.y));
+                                       d = constrain(30*(rad*40-d),0,100);
+                                       return lx.hsb(lxh(),100, d); // clock
+
+               case 8:         r = amp*200 * map(bnc,0,1,1,sin(PI*t));
+                                       d = min(
+                                               distToSeg(Px.x, Px.y, a1.getX(r),a1.getY(r), a2.getX(r),a2.getY(r)),
+                                               distToSeg(Px.x, Px.y, a2.getX(r),a2.getY(r), a3.getX(r),a3.getY(r)),
+                                               distToSeg(Px.x, Px.y, a3.getX(r),a3.getY(r), a1.getX(r),a1.getY(r))                             // triangle
+                                               );
+                                       d = constrain(30*(rad*40-d),0,100);
+                                       return lx.hsb(lxh(),100, d); // clock
+
+               case 9:         r = amp*200 * map(bnc,0,1,1,sin(PI*t));
+                                       d = min(
+                                               distToSeg(Px.x, Px.y, a1.getX(r),a1.getY(r), a2.getX(r),a2.getY(r)),
+                                               distToSeg(Px.x, Px.y, a2.getX(r),a2.getY(r), a3.getX(r),a3.getY(r)),
+                                               distToSeg(Px.x, Px.y, a3.getX(r),a3.getY(r), a4.getX(r),a4.getY(r)),
+                                               distToSeg(Px.x, Px.y, a4.getX(r),a4.getY(r), a1.getX(r),a1.getY(r))                             // quad
+                                       );
+                                       d = constrain(30*(rad*40-d),0,100);
+                                       return lx.hsb(lxh(),100, d); // clock
+
+               case 10:
+                                       r = map(bnc,0,1,a1.r,amp*200*sin(PI*t));
+                                       return lx.hsb(lxh(),100,c1c(.9+2*rad - dist(Px.x,Px.y,a1.getX(r),a1.getY(r))*.03) );            // sphere
+
+               case 11:
+                                       Px.z=mCtr.z; cMid.z=mCtr.z;
+                                       return lx.hsb(lxh(),100,c1c(1 - calcCone(Px,cMid,mCtr) * 0.02 > .5?1:0));                               // cone
+
+               case 12:        return lx.hsb(lxh() + noise(Pn.x,Pn.y,Pn.z + (NoiseMove+50000)/1000.)*200,
+                                               85,c1c(Pn.y < noise(Pn.x + NoiseMove/2000.,Pn.z)*(1+amp)-amp/2.-.1 ? 1 : 0));   // noise
+
+               case 13:        
+               case 14:        float y=0; for (rWave w : waves) y += .5*w.val(Pn.x);   // wave
+                                       V.set(Pn.x, .7+y, Pn.z);
+                                       break;
+
+               default:        return lx.hsb(0,0,0);
+               }
+
+               return lx.hsb(lxh(), 100, c1c(1 - V.dist(Pn)/rad));
+       }
+}
+//----------------------------------------------------------------------------------------------------------------------------------
+boolean dDebug = false;
+class dCursor {
+       dVertex vCur, vNext, vDest;
+       float   destSpeed;
+       int     posStop, pos,posNext;   // 0 - 65535
+       color   clr;
+
+       dCursor() {}
+
+       boolean isDone  ()                                                                      { return pos==posStop;                                                                           }
+       boolean atDest  ()                                                                      { return vCur.s==vDest.s || 
+                                                                                                                                xyDist(vCur.getPoint(0), vDest.getPoint(0)) < 12 || 
+                                                                                                                                xyDist(vCur.getPoint(0), vDest.getPoint(15))< 12;}
+       void    setCur  (dVertex _v, int _p)                            { p2=null; vCur=_v; pos=_p; pickNext();                                          }
+       void    setCur  (dPixel  _p)                                            { setCur(_p.v, _p.pos);                                                                          }
+       void    setNext (dVertex _v, int _p, int _s)            { vNext = _v; posNext = _p<<12; posStop = _s<<12;                        }
+       void    setDest (dVertex _v, float _speed)                      { vDest = _v; destSpeed = _speed;                                                        }
+       void    onDone  ()                                                                      { setCur(vNext, posNext); pickNext();                                            }
+
+       float   minDist;
+       int     nTurns;
+       boolean bRandEval;
+
+       void    evaluate(dVertex v, int p, int s) {
+               if (v == null) return; ++nTurns;
+               if (bRandEval) {
+                       if (random(nTurns) < 1) setNext(v,p,s); return; }
+               else {
+                       float d = xyDist(v.getPoint(15), vDest.getPoint(0));
+                       if (d <  minDist)                                       { minDist=d; setNext(v,p,s); }
+                       if (d == minDist && random(2)<1)        { minDist=d; setNext(v,p,s); }
+               }
+       }
+
+       void    evalTurn(dTurn t) { 
+               if (t == null || t.pos0<<12 <= pos) return; 
+               evaluate(t.v    ,    t.pos1, t.pos0);
+               evaluate(t.v.opp, 16-t.pos1, t.pos0);
+       }
+
+       void    pickNext()      {
+               bRandEval = random(.05+destSpeed) < .05; minDist=500; nTurns=0;
+               evaluate(vCur.c0, 0, 16);       evaluate(vCur.c1, 0, 16);
+               evaluate(vCur.c2, 0, 16);       evaluate(vCur.c3, 0, 16);
+               evalTurn(vCur.t0);                      evalTurn(vCur.t1);
+               evalTurn(vCur.t2);                      evalTurn(vCur.t3);
+       }
+
+       LXPoint         p1, p2; int i2;
+
+       int draw(int nAmount, SCPattern pat) {
+               int nFrom       = (pos    ) >> 12;
+               int     nMv     = min(nAmount, posStop-pos);
+               int     nTo     = min(15,(pos+nMv) >> 12);
+               dVertex v       = vCur;
+
+               if (dDebug) {   p1 = v.getPoint(nFrom); float d = (p2 == null ? 0 : pointDist(p1,p2)); if (d>5) { println("too wide! quitting: " + d); exit(); }}
+                                                               for (int i = nFrom; i <= nTo; i++) { pat.getColors()[v.ci          + v.dir*i     ] = clr; }
+               if (v.same != null)             for (int i = nFrom; i <= nTo; i++) { pat.getColors()[v.same.ci + v.same.dir*i] = clr; }
+
+               if (dDebug) {   p2 = v.getPoint(nTo); i2 = nTo; }
+
+               pos += nMv; return nAmount - nMv;
+                       }       
+}
+
+//----------------------------------------------------------------------------------------------------------------------------------
+class Worms extends SCPattern {
+       float   StripsPerSec    = 10;
+       float   TrailTime               = 3000;
+       int     numCursors              = 50;
+       ArrayList<dCursor> cur  = new ArrayList<dCursor>(30);
+
+       private GraphicEQ eq = null;
+
+       private BasicParameter pBeat      = new BasicParameter("BEAT",  0);
+       private BasicParameter pSpeed     = new BasicParameter("FAST", .2);
+       private BasicParameter pBlur      = new BasicParameter("BLUR", .3);
+       private BasicParameter pWorms     = new BasicParameter("WRMS", .3);
+       private BasicParameter pConfusion = new BasicParameter("CONF", .1);
+       private BasicParameter pEQ        = new BasicParameter("EQ"  ,  0);
+       private BasicParameter pSpawn     = new BasicParameter("DIR" ,  0);
+       private BasicParameter pColor     = new BasicParameter("CLR" ,  .1);
+
+       float   zMidLat = 82.;
+       float   nConfusion;
+       private final Click moveChase = new Click(1000);
+
+       PVector middle;
+       int     AnimNum() { return floor(pSpawn.getValuef()*(4-.01));   }
+       float   randX() { return random(model.xMax-model.xMin)+model.xMin; }
+       float   randY() { return random(model.yMax-model.yMin)+model.yMin; }
+       PVector randEdge() { 
+               return random(2) < 1 ?  new PVector(random(2)<1 ? model.xMin:model.xMax, randY(), zMidLat)      :
+                                                               new PVector(randX(), random(2)<1 ? model.yMin:model.yMax, zMidLat)      ;
+       }
+
+       Worms(GLucose glucose) {
+               super(glucose); 
+           addModulator(moveChase).start();
+           addParameter(pBeat);    addParameter(pSpeed);
+           addParameter(pBlur);    addParameter(pWorms);
+           addParameter(pEQ);      addParameter(pConfusion);
+               addParameter(pSpawn);   addParameter(pColor);
+
+           middle = new PVector(1.5*model.cx, 1.5*model.cy, 71);
+               if (lattice == null) lattice = new dLattice();
+               for (int i=0; i<numCursors; i++) { dCursor c = new dCursor(); reset(c); cur.add(c); }
+               onParameterChanged(pEQ); setNewDest();
+       }
+
+       void onParameterChanged(LXParameter parameter) {
+               super.onParameterChanged(parameter);
+               nConfusion = 1-pConfusion.getValuef();
+               for (int i=0; i<numCursors; i++) {
+                       if (parameter==pSpawn) reset(cur.get(i));
+                       cur.get(i).destSpeed = nConfusion;
+               }
+       }
+
+       float getClr() { return lx.getBaseHuef() + random(pColor.getValuef()*300); }
+       void reset(dCursor c) {
+               switch(AnimNum()) {
+                       case 0: c.clr = lx.hsb(getClr(),100,100);                       // middle to edges
+                                       c.setDest(lattice.getClosest(randEdge()).v, nConfusion);
+                                       c.setCur (lattice.getClosest(middle));
+                                       break;
+
+                       case 1: c.clr = lx.hsb(getClr(),100,100);                               // top to bottom
+                                       float xLin = randX();
+                                       c.setDest(lattice.getClosest(new PVector(xLin, 0         , zMidLat)).v, nConfusion);
+                                       c.setCur (lattice.getClosest(new PVector(xLin, model.yMax, zMidLat)));
+                                       break;
+
+                       case 2: c.clr = lx.hsb(getClr(),100,100); break;                // chase a point around
+
+                       case 3: boolean bLeft = random(2)<1;
+                                       c.clr = lx.hsb(getClr()+random(120),100,100);                           // sideways
+                                       float yLin = randX();
+                                       c.setDest(lattice.getClosest(new PVector(bLeft ? 0 : model.xMax,yLin,zMidLat)).v, nConfusion);
+                                       c.setCur (lattice.getClosest(new PVector(bLeft ? model.xMax : 0,yLin,zMidLat)));
+                                       break;
+               }
+               if (pBlur.getValuef() == 1 && random(2)<1) c.clr = lx.hsb(0,0,0);
+       }
+
+       void setNewDest() {
+               if (AnimNum() != 2) return;
+               PVector dest = new PVector(randX(), randY(), zMidLat);
+               for (int i=0; i<numCursors; i++) {
+                       cur.get(i).setDest(lattice.getClosest(dest).v, nConfusion);
+                       cur.get(i).clr = lx.hsb(getClr()+75,100,100);   // chase a point around
+               }
+       }
+
+       void run(double deltaMs) { 
+               if (deltaMs > 100) return;
+           if (moveChase.click()) setNewDest();
+
+           float fBass=0, fTreble=0;
+           if (pEQ.getValuef()>0) {            // EQ
+                   eq.run(deltaMs);
+                   fBass       = eq.getAverageLevel(0, 4);
+                   fTreble = eq.getAverageLevel(eq.numBands-7, 7);
+               }
+
+               if (pBlur.getValuef() < 1) {    // trails
+                       for (int i=0,s=model.points.size(); i<s; i++) {
+                               color c = colors[i]; float b = lx.b(c); 
+                               if (b>0) colors[i] = lx.hsb(lx.h(c), lx.s(c), constrain((float)(b-100*deltaMs/(pBlur.getValuef()*TrailTime)),0,100));
+                       }
+               }
+
+               int nWorms = floor(pWorms.getValuef() * numCursors * 
+                                        map(pEQ.getValuef(),0,1,1,constrain(2*fTreble,0,1)));
+
+               for (int i=0; i<nWorms; i++) {
+                       dCursor c = cur.get(i);
+                       int nLeft = floor((float)deltaMs*.001*StripsPerSec * 65536 * (5*pSpeed.getValuef()));
+                       nLeft *= (1 - lx.tempo.rampf()*pBeat.getValuef());
+                       while(nLeft > 0) { 
+                               nLeft = c.draw(nLeft,this); if (!c.isDone()) continue;
+                               c.onDone(); if (c.atDest()) reset(c);
+                       }
+               }
+       }
+
+
+       public void onActive() { if (eq == null) {
+               eq = new GraphicEQ(lx, 16);             eq.slope.setValue(0.6);
+               eq.level.setValue(0.65);                eq.range.setValue(0.35);
+               eq.release.setValue(0.4);
+       }}
+}
+//----------------------------------------------------------------------------------------------------------------------------------