cubecenter, not quite right
[SugarCubes.git] / DanHorwitz.pde
old mode 100644 (file)
new mode 100755 (executable)
index 64991c2..d084609
-public class _P extends BasicParameter {
-       _P(String label, double value)                                  { super(label,value);                                                   }
-       public void     updateValue     (double value)          { super.updateValue(value);                                             }
-       public float    Val                     ()                                      { return getValuef();                                                   }
-       public int              Pick            (int b)                         { return int(Val())==1 ? b-1 : int(b*Val());    } 
-}
-
-float  c1c                     (float a)       { return 100*constrain(a,0,1); }
-float  CalcCone        (float x1, float y1, float z1, float x2, float y2, float z2) {
-                                               return degrees( acos ( ( x1*x2 + y1*y2 + z1*z2 ) / (sqrt(x1*x1+y1*y1+z1*z1) * sqrt(x2*x2+y2*y2+z2*z2)) ) );}
-float          xMax,yMax,zMax;
-float          zTime   = random(10000);
-_P                     pChoose = new _P("ANIM", 0);
-//----------------------------------------------------------------------------------------------------------------------------------
-class Pong extends SCPattern {
-       SinLFO x,y,z,dx,dy,dz; 
-       float cRad;     _P size;
-
-       Pong(GLucose glucose) {
-               super(glucose);
-               xMax = model.xMax; yMax = model.yMax; zMax = model.zMax;
-               cRad = xMax/15;
-               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, xMax - cRad, 0)).trigger();  x.modulateDurationBy(dx);
-               addModulator(y  = new SinLFO(cRad, yMax - cRad, 0)).trigger();  y.modulateDurationBy(dy);
-               addModulator(z  = new SinLFO(cRad, zMax - cRad, 0)).trigger();  z.modulateDurationBy(dz);
-           addParameter(pChoose);
-           addParameter(size = new _P("SIZE", 0.4));
-       }
-
-       color Calc(float px, float py, float pz, float vx, float vy, float vz) {
-               switch(pChoose.Pick(3)) {
-               /*spot*/        case 0: return color(0,0,c1c(1 - CalcCone(vx-xMax/2, vy, vz-zMax/2, px-xMax/2, py, pz-zMax/2)*max(.02,.45-size.Val())));
-               /*ball*/        case 1: return color(0,0,c1c(1 - dist(vx,vy,vz,px,py,pz)*.5/cRad));
-               /*ballz*/       default:return color(0,0,c1c(1 - min(   dist(vx,vy,vz,px,py,pz),
-                                                                                                                       dist(vx,vy,vz,xMax-px,yMax-py,zMax-pz))*.5/cRad)); 
-               }               
-       }
-
-       public void run(int deltaMs) {
-               cRad = xMax*size.Val()/6;
-               for (Point p : model.points) {
-                       colors[p.index] = Calc(p.fx, p.fy, p.fz, x.getValuef(), y.getValuef(), z.getValuef());
-               }
-       }
-}
-//----------------------------------------------------------------------------------------------------------------------------------
-class NDat {
-       float xz, yz, zz, hue, sat, speed, dir, den, contrast;
-       float xoff,yoff,zoff;
-       NDat (float _hue, float _sat, float _xz, float _yz, float _zz, float _con, float _den, float _speed, float _dir) {
-               hue=_hue; sat=_sat; xz=_xz; yz=_yz; zz =_zz; contrast=_con; den=_den; speed=_speed; dir=_dir;
-               xoff = random(100e3); yoff = random(100e3); zoff = random(100e3);
-       }
-}
-
-class Noise extends SCPattern
-{
-       int     CurAnim =       -1      ;
-       ArrayList noises = new ArrayList();
-       _P pSpeed, pMir, pBright, pContrast, pDensity, pDir;
-
-       Noise(GLucose glucose) {
-               super(glucose);                                                                  addParameter(pChoose   );
-               addParameter(pSpeed             = new _P("MPH"  , .55)); addParameter(pMir              = new _P("MIR"  , 0     ));
-               addParameter(pBright    = new _P("BRTE" ,  1 )); addParameter(pDir              = new _P("DIR"  , 0 ));
-               addParameter(pContrast  = new _P("CNTR" ,  .5)); addParameter(pDensity  = new _P("DENS" , .5));
-       }
-
-       public void run(int deltaMs) {
-               zTime += deltaMs*(pSpeed.Val()-.5)*.002;
-               int anim = pChoose.Pick(6);
-               if (anim != CurAnim) {
-                       noises.clear(); CurAnim = anim; switch(anim) {
-                       //_hue, _sat, _xz, _yz, _zz, _con, _den, _speed, _dir
-                       case 0: noises.add(new NDat(0  ,0  ,100,100,200,0  ,40 ,3  ,2)); break; // clouds
-                       case 1: noises.add(new NDat(0  ,0  ,75 ,75 ,150,1  ,45 ,3  ,0)); break; // drip
-                       case 2: noises.add(new NDat(0  ,0  ,2  ,400,2  ,1  ,40 ,3  ,0)); break; // rain
-                       case 3: noises.add(new NDat(40 ,100,100,100,200,0  ,20 ,1  ,2)); 
-                                       noises.add(new NDat(0  ,100,100,100,200,0  ,20 ,5   ,2)); break;        // fire 1
-                       case 4: noises.add(new NDat(0  ,100,40 ,40 ,40 ,.5 ,25 ,2.5 ,2));
-                                       noises.add(new NDat(20 ,100,40 ,40 ,40 ,.5 ,25 ,4   ,0));
-                                       noises.add(new NDat(40 ,100,40 ,40 ,40 ,.5 ,25 ,2   ,1));
-                                       noises.add(new NDat(60 ,100,40 ,40 ,40 ,.5 ,25 ,3   ,3)); break;        // machine
-                       case 5: noises.add(new NDat(0  ,100,400,100,2  ,1  ,20 ,3   ,2));
-                                       noises.add(new NDat(20 ,100,400,100,2  ,1  ,20 ,2.5 ,0));
-                                       noises.add(new NDat(40 ,400,100,100,2  ,1  ,20 ,2   ,1));
-                                       noises.add(new NDat(60 ,400,100,100,2  ,1  ,20 ,1.5 ,3)); break;        // spark
-                       default: break;
-                       }
-               }
-
-               for (Point p : model.points) {
-                       color c = color(0,0,0);
-                       int mir = pMir.Pick(5);
-
-                       for (int i=0;i<noises.size(); i++) { NDat n = (NDat) noises.get(i);
-                               float vx=p.fx, vy=p.fy, vz=p.fz;
-//                             if (vx > xMax/2 && (mir == 1 || mir == 3 || mir == 4)) { vx = xMax-vx; }
-                               if (vy > yMax/2 && (mir == 2 || mir == 3 || mir == 4)) { vy = yMax-vy; }
-//                             if (vz > zMax/2 && (                                            mir == 4)) { vz = zMax-vz; }
-                               float deg = radians(90*(n.dir + pDir.Pick(4)));
-                               float zx  = zTime * n.speed * sin(deg);
-                               float zy  = zTime * n.speed * cos(deg);
-                               float b   = noise(vx/n.xz+zx+n.xoff,vy/n.yz+zy+n.yoff,vz/n.zz+n.zoff)*1.6-.3 + n.den/100 + pDensity.Val() -1;
-                               float con = 1/constrain(2-n.contrast - 2*pContrast.Val(),0,1);
-                               b = b < .5 ? pow(b,con) : 1-pow(1-b,con);
-                               c = blendColor(c,color(n.hue,n.sat,c1c(b * pBright.Val())),ADD);
-                       }
-                       colors[p.index] = c;
-               }
-       }
-}
-//----------------------------------------------------------------------------------------------------------------------------------
+//----------------------------------------------------------------------------------------------------------------------------------\r
+public class Pong extends DPat {\r
+       SinLFO x,y,z,dx,dy,dz; \r
+       float cRad;     DParam pSize;\r
+       Pick    pChoose;\r
+       xyz     v = new xyz(), vMir =  new xyz();\r
+\r
+       Pong(GLucose glucose) {\r
+               super(glucose);\r
+               cRad = mMax.x/10;\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, mMax.x - cRad, 0)).trigger();        x.modulateDurationBy(dx);\r
+               addModulator(y  = new SinLFO(cRad, mMax.y - cRad, 0)).trigger();        y.modulateDurationBy(dy);\r
+               addModulator(z  = new SinLFO(cRad, mMax.z - cRad, 0)).trigger();        z.modulateDurationBy(dz);\r
+           pSize       = addParam      ("Size"                 , 0.4   );\r
+           pChoose = addPick   ("Animiation"   , 0, 2, new String[] {"Pong", "Ball", "Cone"}   );\r
+       }\r
+\r
+       void    StartRun(double deltaMs)        { cRad = mMax.x*pSize.Val()/6; }\r
+       color   CalcPoint(xyz p)                {\r
+               v.set(x.getValuef(), y.getValuef(), z.getValuef());\r
+               v.z=0;p.z=0;// ignore z dimension\r
+               switch(pChoose.Cur()) {\r
+               case 0: vMir.set(mMax); vMir.subtract(p);\r
+                               return color(0,0,c1c(1 - min(v.distance(p), v.distance(vMir))*.5/cRad));        // balls\r
+               case 1: return color(0,0,c1c(1 - v.distance(p)*.5/cRad));                                                       // ball\r
+               case 2: vMir.set(mMax.x/2,0,mMax.z/2);\r
+                               return color(0,0,c1c(1 - CalcCone(p,v,vMir) * 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 , zTheta=0, zSin, zCos, rtime, ttime, transAdd;\r
+       DParam          pSpeed , pDensity;\r
+       Pick            pChoose, pSymm;\r
+       int                     _ND = 4;\r
+       NDat            N[] = new NDat[_ND];\r
+\r
+       Noise(GLucose glucose) {\r
+               super(glucose);\r
+               pSpeed          = addParam("Fast"       , .55);\r
+               pDensity        = addParam("Dens"        , .5);\r
+               pSymm           = addPick("Symmetry" , 0, 3, new String[] {"None", "X", "Y", "Radial"}  );\r
+               pChoose         = addPick("Animation", 6, 7, new String[] {"Drip", "Cloud", "Rain", "Fire", "Machine", "Spark","VWave", "Wave"} );\r
+               for (int i=0; i<_ND; i++) N[i] = new NDat();\r
+       }\r
+\r
+       void StartPattern() { zTime = random(500); zTheta=0; rtime = 0; ttime = 0; transAdd=0; }\r
+       void StartRun(double deltaMs) {\r
+               zTime   += deltaMs*(pSpeed.Val()-.5)*.002       ;\r
+               zTheta  += deltaMs*(pSpin .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
+                       pSpin           .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(mCtr, zSin, zCos);\r
+               \r
+               if (CurAnim == 6 || CurAnim == 7) {\r
+                       P.setNorm();\r
+                       return color(0,0, 100 * (\r
+                                                       constrain(1-50*(1-pDensity.Val())*abs(P.y-sin(zTime*10  + P.x*(300))*.5 - .5),0,1) + \r
+                       (CurAnim == 7 ? constrain(1-50*(1-pDensity.Val())*abs(P.x-sin(zTime*10  + P.y*(300))*.5 - .5),0,1) : 0))\r
+                       );\r
+               }                       \r
+                       \r
+               if (iSymm == XSym && P.x > mMax.x/2) P.x = mMax.x-P.x;\r
+               if (iSymm == YSym && P.y > mMax.y/2) P.y = mMax.y-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,mCtr)/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
+       public class rAngle {\r
+               float   prvA, dstA, c;\r
+               float   prvR, dstR, r;          \r
+               float   _cos, _sin, x, y;\r
+               float   fixAngle        (float a, float b) { return a<b ?\r
+                                                                               (abs(a-b) > abs(a+2*PI-b) ? a : a+2*PI) :\r
+                                                                               (abs(a-b) > abs(a-2*PI-b) ? a : a-2*PI) ; }\r
+               float   getX(float r)   {       return mCtr.x + _cos*r; }\r
+               float   getY(float r)   {       return mCtr.y + _sin*r; }\r
+               void    move()                  {       c               = interp(t,prvA,dstA); \r
+                                                                       r               = interp(t,prvR,dstR);\r
+                                                                       _cos    = cos(c);       _sin    = sin(c);\r
+                                      x                = getX(r);      y               = getY(r);              }               \r
+               void    set()                   {       prvA    = dstA;         dstA    = random(2*PI);         prvA = fixAngle(prvA, dstA);\r
+                                                                       prvR    = dstR;         dstR    = random(mCtr.y);                                                                       }\r
+       }\r
+\r
+       int             nBeats  =       0;\r
+       DParam  pAmp, pRadius, pBounce;\r
+\r
+       float   t,amp,rad,bnc;\r
+       float   zTheta=0;\r
+       ArrayList<rWave> waves = new ArrayList<rWave>(10);\r
+\r
+       rAngle  a1 = new rAngle(), a2 = new rAngle(),\r
+                       a3 = new rAngle(), a4 = new rAngle();\r
+       xyz             cPrev   = new xyz(), cRand      = new xyz(),\r
+                       cMid    = new xyz(), V          = new xyz(),\r
+                       Theta   = new xyz(), TSin       = new xyz(),\r
+                       TCos    = new xyz(), cMidNorm = new xyz(),\r
+                       Pn              = new xyz();\r
+       float   LastBeat=3, LastMeasure=3;\r
+       int             CurRandTempo = 1, CurRandTPat = 1;\r
+\r
+       Pick    pTimePattern, pTempoMult, pShape;\r
+       int             RandCube;\r
+\r
+       Play(GLucose glucose) {\r
+               super(glucose);\r
+           pRadius             = addParam("Rad"        , .1    );\r
+               pBounce         = addParam("Bnc"        , .2    );\r
+           pAmp                = addParam("Amp"        , .2    );\r
+               pTempoMult      = addPick ("TMult"      , 0 , 5         , new String[] {"1x", "2x", "4x", "8x", "16x", "Rand"   }       );\r
+               pTimePattern= addPick ("TPat"   , 6 , 7         , new String[] {"Bounce", "Sin", "Roll", "Quant", "Accel", "Deccel", "Slide", "Rand"}   );\r
+               pShape          = addPick ("Shape"      , 3 , 15        , new String[] {"Line", "Tap", "V", "RandV",\r
+                                                                                                                                       "Pyramid", "Wings", "W2", "Clock",\r
+                                                                                                                                       "Triangle", "Quad", "Sphere", "Cone",\r
+                                                                                                                                       "Noise", "Wave", "?", "?"}                                              );\r
+       }\r
+\r
+       public class rWave {\r
+               float v0, a0, x0, t,damp,a;\r
+               boolean bDone=false;\r
+               final float len=8;\r
+               rWave(float _x0, float _a0, float _v0, float _damp) { x0=_x0*len; a0=_a0; v0=_v0; t=0; damp = _damp; }\r
+               void move(double deltaMs) {\r
+                       t += deltaMs*.001;\r
+                       if (t>4) bDone=true;\r
+               }\r
+               float val(float _x) {\r
+                       _x*=len;\r
+                       float dist = t*v0 - abs(_x-x0);\r
+                       if (dist<0) { a=1; return 0; }\r
+                       a  = a0*exp(-dist*damp) * exp(-abs(_x-x0)/(.2*len)); // * max(0,1-t/dur)\r
+                       return  -a*sin(dist);\r
+               }\r
+       }\r
+\r
+       void StartPattern() { zTheta=0; }\r
+       void StartRun(double deltaMs) {\r
+               t       = lx.tempo.rampf();\r
+               amp = pAmp.Val();\r
+               rad     = pRadius.getValuef();\r
+               bnc     = pBounce.getValuef();          \r
+               zTheta  += deltaMs*(pSpin .Val()-.5)*.01;\r
+\r
+               Theta   .set(pRotX.Val()*PI*2, pRotY.Val()*PI*2, pRotZ.Val()*PI*2 + zTheta);\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(3) < 1) { CurRandTempo = int(random(4)); if (CurRandTempo == 3) CurRandTempo = int(random(4));       }\r
+                       if (random(3) < 1) { CurRandTPat  = pShape.Cur() > 6 ? 2+int(random(5)) : int(random(7));                                       }\r
+               } LastMeasure = t;\r
+                       \r
+               int nTempo = pTempoMult  .Cur(); if (nTempo == 5) nTempo = CurRandTempo;\r
+               int nTPat  = pTimePattern.Cur(); if (nTPat  == 7) 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
+               int i=0; while (i< waves.size()) {\r
+                       rWave w = waves.get(i);\r
+                       w.move(deltaMs); if (w.bDone) waves.remove(i); else i++;\r
+               }\r
+\r
+               if ((t<LastBeat && !pKey.b) || DG.KeyPressed>-1) {\r
+                       waves.add(new rWave(\r
+                                               pKey.b ? map(DG.KeyPressed,0,7,0,1) : random(1),                // location\r
+                                               bnc*10,                 // bounciness\r
+                                               7,                              // velocity\r
+                                               2*(1-amp)));    // dampiness\r
+                       DG.KeyPressed=-1;\r
+                       if (waves.size() > 5) waves.remove(0);\r
+               }\r
+               \r
+               if (t<LastBeat) {\r
+                       cPrev.set(cRand); cRand.setRand();\r
+                       a1.set(); a2.set(); a3.set(); a4.set();\r
+               } LastBeat = t;\r
+\r
+               switch (nTPat) {\r
+                       case 0:         t = sin(PI*t);                                                  break;  // bounce\r
+                       case 1:         t = norm(sin(2*PI*(t+PI/2)),-1,1);              break;  // sin\r
+                       case 2:         t = t;                                                                  break;  // roll\r
+                       case 3:         t = constrain(int(t*8)/7.,0,1);                 break;  // quant\r
+                       case 4:         t = t*t*t;                                                              break;  // accel\r
+                       case 5:         t = sin(PI*t*.5);                                               break;  // deccel\r
+                       case 6:         t = .5*(1-cos(PI*t));                                   break;  // slide\r
+               }\r
+               \r
+               cMid.set                (cPrev);        cMid.interpolate        (t,cRand);\r
+               cMidNorm.set    (cMid);         cMidNorm.setNorm();\r
+               a1.move(); a2.move(); a3.move(); a4.move();\r
+       }\r
+\r
+       color CalcPoint(xyz Px) {\r
+               if (Theta.x != 0) Px.RotateX(mCtr, TSin.x, TCos.x);\r
+               if (Theta.y != 0) Px.RotateY(mCtr, TSin.y, TCos.y);\r
+               if (Theta.z != 0) Px.RotateZ(mCtr, TSin.z, TCos.z);\r
+               \r
+               Pn.set(Px); Pn.setNorm();\r
+\r
+               float mp        = min(Pn.x, Pn.z);\r
+               float yt        = map(t,0,1,.5-bnc/2,.5+bnc/2);\r
+               float r,d;\r
+\r
+               switch (pShape.Cur()) {\r
+               case 0:         V.set(Pn.x, yt                                                          , Pn.z);                                                        break;  // bouncing line\r
+               case 1:         V.set(Pn.x, map(cos(PI*t * Pn.x),-1,1,0,1)  , Pn.z);                                                    break;  // top tap\r
+               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\r
+               case 3:         V.set(Pn.x, Pn.x < cMidNorm.x ? map(Pn.x,0,cMidNorm.x, .5,yt) :\r
+                                                                                               map(Pn.x,cMidNorm.x,1, yt,.5), Pn.z);                           break;  //  Random V shape\r
+\r
+               case 4:         V.set(Pn.x,     .5*(Pn.x < cMidNorm.x ?         map(Pn.x,0,cMidNorm.x, .5,yt) :\r
+                                                                                                               map(Pn.x,cMidNorm.x,1, yt,.5)) +\r
+                                                       .5*(Pn.z < cMidNorm.z ?         map(Pn.z,0,cMidNorm.z, .5,yt) :\r
+                                                                                                               map(Pn.z,cMidNorm.z,1, yt,.5)), Pn.z);          break;  //  Random Pyramid shape\r
+                                                                                                       \r
+               case 5:         V.set(Pn.x, bnc*map((Pn.x-.5)*(Pn.x-.5),0,.25,0,t-.5)+.5, Pn.z);                                break;  // wings\r
+               case 6:         V.set(Pn.x, bnc*map((mp  -.5)*(mp  -.5),0,.25,0,t-.5)+.5, Pn.z);                                break;  // wings\r
+\r
+               case 7:         d = min(\r
+                                               distToSeg(Px.x, Px.y, a1.getX(70),a1.getY(70), mCtr.x, mCtr.y),\r
+                                               distToSeg(Px.x, Px.y, a2.getX(40),a2.getY(40), mCtr.x, mCtr.y));\r
+                                       d = constrain(30*(rad*40-d),0,100);\r
+                                       return color(0,max(0,150-d), d); // clock\r
+\r
+               case 8:         r = amp*200 * map(bnc,0,1,1,sin(PI*t));\r
+                                       d = min(\r
+                                               distToSeg(Px.x, Px.y, a1.getX(r),a1.getY(r), a2.getX(r),a2.getY(r)),\r
+                                               distToSeg(Px.x, Px.y, a2.getX(r),a2.getY(r), a3.getX(r),a3.getY(r)),\r
+                                               distToSeg(Px.x, Px.y, a3.getX(r),a3.getY(r), a1.getX(r),a1.getY(r))                             // triangle\r
+                                               );\r
+                                       d = constrain(30*(rad*40-d),0,100);\r
+                                       return color(0,max(0,150-d), d); // clock\r
+\r
+               case 9:         r = amp*200 * map(bnc,0,1,1,sin(PI*t));\r
+                                       d = min(\r
+                                               distToSeg(Px.x, Px.y, a1.getX(r),a1.getY(r), a2.getX(r),a2.getY(r)),\r
+                                               distToSeg(Px.x, Px.y, a2.getX(r),a2.getY(r), a3.getX(r),a3.getY(r)),\r
+                                               distToSeg(Px.x, Px.y, a3.getX(r),a3.getY(r), a4.getX(r),a4.getY(r)),\r
+                                               distToSeg(Px.x, Px.y, a4.getX(r),a4.getY(r), a1.getX(r),a1.getY(r))                             // quad\r
+                                       );\r
+                                       d = constrain(30*(rad*40-d),0,100);\r
+                                       return color(0,max(0,150-d), d); // clock\r
+\r
+               case 10:\r
+                                       r = map(bnc,0,1,a1.r,amp*200*sin(PI*t));\r
+                                       return color(0,0,c1c(.9+2*rad - dist(Px.x,Px.y,a1.getX(r),a1.getY(r))*.03) );           // sphere\r
+\r
+               case 11:\r
+                                       Px.z=mCtr.z; cMid.z=mCtr.z;\r
+                                       return color(0,0,c1c(1 - CalcCone(Px,cMid,mCtr) * 0.02 > .5?1:0));                              // cone\r
+\r
+               case 12:        return color(100 + noise(Pn.x,Pn.y,Pn.z + (NoiseMove+50000)/1000.)*200,\r
+                                               85,c1c(Pn.y < noise(Pn.x + NoiseMove/2000.,Pn.z)*(1+amp)-amp/2.-.1 ? 1 : 0));   // noise\r
+\r
+               case 13:        float y=0; for (rWave w : waves) y += .5*w.val(Pn.x);\r
+                                       V.set(Pn.x, .7+y, Pn.z);\r
+                                       break;\r
+\r
+               default:        return color(0,0,0);\r
+               }\r
+\r
+               return color(0,\r
+                               150-c1c(1 - V.distance(Pn)/rad),\r
+                               c1c(1 - V.distance(Pn)/rad));\r
+       }\r
+}\r
+//----------------------------------------------------------------------------------------------------------------------------------\r
+// 0 - TLB, L (b), BLB, B (l)          // Fwd , Down, Back, Up\r
+// 4 - TLF, F (l), BLF, L (f)          // Fwd , Down, Back, Up\r
+// 8 - TRF, R (f), BRF, F (r)          // Back, Down, Fwd , Up\r
+// 12- TRB, B (r), BRB, R (b)          // Back, Down, Fwd , Up\r
+// 1->7, 15->9\r
+\r
+class dBolt {\r
+       dStrip v, h;\r
+       int    vpos, hpos;\r
+       dBolt(dStrip _v, dStrip _h, int _vpos, int _hpos) {\r
+               v=_v; h=_h; vpos=_vpos; hpos=_hpos;\r
+               if (v.b0 == null) { v.b0=this; h.b0=this; } \r
+               else                      { v.b1=this; h.b1=this; }\r
+       }\r
+}\r
+\r
+class dVertex {\r
+       dStrip  s1      , s2    ;\r
+       int     dir1, dir2      ;\r
+       dVertex(dStrip s, int d) {\r
+               int _a = (s.iS%4==1)? (d==1? 5: 3) :\r
+                                (s.iS%4==3)? (d==1? 9:11) :\r
+                                (d==1)         ? (s.Top()?4:12)   : (s.Top()?12:4);\r
+               dir1 = d * (s.isVert() ? -1 : 1);\r
+               dir2 = d;\r
+               s1       = DL_.DS[s.iCube() + ((s.iS+_a) % 16)];\r
+               s2       = DL_.DS[d == 1 ?      (s.idx == s.iFace()+3 ? s.idx-3 : s.idx+1):\r
+                                                               (s.idx == s.iFace()   ? s.idx+3 : s.idx-1)];\r
+               swapout(1 , 6);\r
+               swapout(15,-6);\r
+       }\r
+       void swapout(int a, int b) {\r
+               if (s1.iS == a) { s1 = DL_.DS[s1.idx + b]; dir1 = -dir1; }\r
+               if (s2.iS == a) { s2 = DL_.DS[s2.idx + b]; dir2 = -dir2; }\r
+       }\r
+\r
+}\r
+\r
+class dStrip  { // THIS WAS SUCH A PAIN!\r
+       int row, col, ci, idx, iS, axis;        // 1-y, 2-left, 3-right\r
+       Strip s; \r
+       boolean bTop;   // direction: top ccw, bottom cw.\r
+\r
+       boolean Top   (){       return axis!=1 &&  bTop ;               }\r
+       boolean Bottom(){       return axis!=1 && !bTop;                }\r
+       boolean isVert(){       return axis==1;                                 }\r
+       boolean isHorz(){       return axis!=1;                                 }\r
+       int     iCube (){       return 16*floor(idx/16);                }\r
+       int     iFace (){       return iCube() + 4*floor(iS/4); }\r
+\r
+       void    init(Strip _s, int _i, int _row, int _col)  {\r
+               idx = _i; row = _row; col = _col; s = _s;\r
+               iS      = idx%16; bTop = (iS%4==0);\r
+               ci  = s.points.get(0).index;\r
+               DL_.DQ[col][row] = iCube();\r
+               switch (iS) {\r
+                       case 4: case 6 : case 12: case 14:      axis=2; break;\r
+                       case 0: case 2 : case 8 : case 10:      axis=3; break;\r
+                       default:                                                        axis=1; break;\r
+               }\r
+       }\r
+\r
+       void addBolts() {\r
+               v0 = new dVertex(this, 1);\r
+               v1 = new dVertex(this,-1);\r
+\r
+\r
+               if (iS == 7 && col != 0 && row != 0)                                                                    // left bottom\r
+                       new dBolt(this, DL_.GetStrip(row-1,col-1,(col % 2 == 1) ? 8 : 12),  \r
+                                                       4, (col % 2 == 1) ? 6 : 9);     \r
+\r
+               if (iS == 7 && col != 0 && row < MaxCubeHeight*2-2)                                     // left top\r
+                       new dBolt(this, DL_.GetStrip(row+1,col-1,(col % 2 == 1) ? 10 : 14), \r
+                                                       11, (col % 2 == 1) ? 9 : 6);\r
+\r
+               if (iS == 9 && col < NumBackTowers-1 && row < MaxCubeHeight*2-2)                // right top\r
+                       new dBolt(this, DL_.GetStrip(row+1,col+1,(col % 2 == 1) ? 6 : 2), \r
+                                                       4, (col % 2 == 1) ? 6 : 9);\r
+\r
+               if (iS == 9 && col < NumBackTowers-1 && row != 0)                                               // right bottom\r
+                       new dBolt(this, DL_.GetStrip(row-1,col+1,(col % 2 == 1) ? 4 : 0), \r
+                                                       11, (col % 2 == 1) ? 9 : 6);\r
+       }\r
+\r
+       dBolt   b0, b1;\r
+       dVertex v0, v1;\r
+}\r
+\r
+class dCursor {\r
+       dStrip  s, sNext;\r
+       int     nLast,pos,posNext,end;  // 0 - 65535\r
+       int     dir;                    // 1 or -1\r
+       color   clr;\r
+       \r
+       dCursor(color _c) { clr=_c;}\r
+\r
+       boolean isDone()                                        { return pos==end; }\r
+       void    set(dStrip _s, int _dir)        { \r
+                       s=_s; dir=_dir; pos = 0; end=65536; nLast=-1; sNext=null;\r
+       }\r
+\r
+       boolean MakeTurn(dBolt b) {\r
+               int nEnd=       (s.isVert() ? b.vpos : b.hpos) <<12;\r
+                       nEnd=   (dir==1 ? nEnd : 65536-nEnd);\r
+               if (nEnd < pos) return false;\r
+               if (s.isVert()) { sNext = b.h; posNext = b.hpos<<12; end = nEnd; }\r
+               else                    { sNext = b.v; posNext = b.vpos<<12; end = nEnd; }\r
+               return true;\r
+       }\r
+\r
+       void    PickNext()      {\r
+               if (sNext != null) {\r
+                       if (end == 65536) exit();\r
+                       end                     = 65536;        \r
+                       pos                     = posNext;\r
+                       dir                     = randDir(); \r
+                       if (dir<0) pos = end-pos;\r
+                       s                       = sNext; sNext = null;\r
+                       nLast           = -1;\r
+                       return;// could switch again!!\r
+               } else {\r
+                       dVertex v = (dir == 1 ? s.v0 : s.v1);\r
+                       int r = floor(random(2));\r
+                       set(r==0 ? v.s1 : v.s2,r==0 ? v.dir1 : v.dir2);\r
+               }\r
+\r
+               // plan to turn the corner\r
+               if (random(6)<1 && s.b0 != null && MakeTurn(s.b0)) return;\r
+               if (random(6)<1 && s.b1 != null && MakeTurn(s.b1)) return;\r
+       }\r
+}\r
+\r
+int randDir() { return round(random(1))*2-1; }\r
+\r
+class dLattice {\r
+       int             iTowerStrips=-1;\r
+       dStrip[]        DS = new dStrip[glucose.model.strips.size()];\r
+       int[][]         DQ = new int[NumBackTowers][MaxCubeHeight*2];\r
+\r
+       int             nStrips() { return iTowerStrips; }\r
+       dStrip GetStrip (int row, int col, int off) { return DS[DQ[col][row]+off]; }\r
+       dLattice() {\r
+               DL_=this;\r
+               int   col = 0, row = -2, i=-1;\r
+               for (Strip strip : glucose.model.strips  ) { i++; \r
+                       if (i % 16 == 0) row+=2;\r
+                       if (row >= MaxCubeHeight*2-1) { col++; row = (col%2==1)?1:0; }\r
+                       if (col >= NumBackTowers) continue;\r
+                       iTowerStrips++  ;\r
+                       dStrip s = DS[iTowerStrips] = new dStrip();\r
+                       s.init(strip, iTowerStrips, row, col);\r
+               }\r
+\r
+               for (int j=0; j<iTowerStrips; j++) DS[j].addBolts();\r
+       }\r
+       dStrip  rand()                          { return DS[floor(random(iTowerStrips))];               }\r
+       void    setRand(dCursor c)      { c.set(rand(),randDir());                      }\r
+}\r
+\r
+dLattice DL_;\r
+//----------------------------------------------------------------------------------------------------------------------------------\r
+class Worms extends SCPattern {\r
+\r
+       int draw(dCursor c, int nAmount) {\r
+               int nFrom       = max(c.nLast+1,c.pos >> 12);\r
+               int     nTo     = min(15,(c.pos+nAmount) >> 12); c.nLast=nTo;\r
+               int     nMv     = min(nAmount, c.end-c.pos);\r
+                       c.pos   += nMv;\r
+               for (int i = nFrom; i <= nTo; i++) {\r
+                       int n = c.s.ci + (c.dir>0 ? i : 15-i);\r
+                       colors[n] = c.clr;\r
+               }\r
+               return nAmount - nMv;\r
+       }\r
+\r
+       float   StripsPerSec    = 6;\r
+       float   TrailTime               = 1500;\r
+       int     Cursors                 = 40;\r
+       dCursor cur[] = new dCursor[Cursors];\r
+\r
+       Worms(GLucose glucose) { \r
+               super(glucose); \r
+               if (DL_ == null) DL_ = new dLattice();\r
+               for (int i=0; i<Cursors; i++) { cur[i] = new dCursor(color(random(360),50+random(50),50+random(50))); DL_.setRand(cur[i]); }\r
+       }\r
+\r
+       void run(double deltaMs) {\r
+                       //Test Joints\r
+               if (false) {\r
+                       for (int j=0; j<DL_.nStrips(); j++) {\r
+                                               dStrip s =DL_.DS[j]; dBolt d = s.b0;\r
+                                               if (d != null) {for (int i=0;i<16;i++) {\r
+                                                       if (s == d.v && i <= d.vpos) colors[d.v.ci+i]   = color(0,0,30);\r
+                                                       if (s == d.h && i <= d.hpos) colors[d.h.ci+i]   = color(0,0,30);\r
+                                               }}\r
+\r
+                                               d = s.b1; \r
+                                               if (d != null) {for (int i=0;i<16;i++) {\r
+                                                       if (s == d.v && i >= d.vpos) colors[d.v.ci+i]   = color(0,0,30);\r
+                                                       if (s == d.h && i >= d.hpos) colors[d.h.ci+i]   = color(0,0,30);\r
+                                               }}\r
+                       }\r
+               } else {\r
+                       for (int i=0; i<Cursors; i++) {\r
+                               int nLeft = floor((float)deltaMs*.001*StripsPerSec * 65536);\r
+                               while(nLeft > 0) {      \r
+                                       nLeft = draw(cur[i], nLeft);\r
+                                       if (cur[i].isDone()) cur[i].PickNext(); \r
+                               }\r
+                       }\r
+\r
+                       for (int i=0,s=model.points.size(); i<s; i++) {\r
+                               float b = brightness(colors[i]); \r
+                               color c = colors[i];\r
+                               if (b>0) colors[i] = color(hue(c), saturation(c), \r
+                                               (float)(b-100*deltaMs/TrailTime));\r
+                       }\r
+               }\r
+       }\r
+}\r
+//----------------------------------------------------------------------------------------------------------------------------------\r