cubecurl
[SugarCubes.git] / DanHorwitz.pde
index b15bad9d1b6a801a8fd6ebe84610043f3248f8c0..d084609e620c2dde46d31ce2e02c09f674da5ffe 100755 (executable)
@@ -1,38 +1,47 @@
 //----------------------------------------------------------------------------------------------------------------------------------\r
 public class Pong extends DPat {\r
        SinLFO x,y,z,dx,dy,dz; \r
-       float cRad;     _DhP pSize;\r
-       Pick  pChoose;\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 = xdMax/15;\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, 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
+               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     , 3     );\r
+           pChoose = addPick   ("Animiation"   , 0, 2, new String[] {"Pong", "Ball", "Cone"}   );\r
        }\r
 \r
-       void    StartRun(double deltaMs)        { cRad = xdMax*pSize.Val()/6; }\r
+       void    StartRun(double deltaMs)        { cRad = mMax.x*pSize.Val()/6; }\r
        color   CalcPoint(xyz p)                {\r
-               xyz v = new xyz(x.getValuef(), y.getValuef(), z.getValuef());\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: 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
+               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
 //----------------------------------------------------------------------------------------------------------------------------------\r
 public class NDat {\r
-       float xz, yz, zz, hue, sat, speed, angle, den;\r
-       float xoff,yoff,zoff;\r
-       NDat (float _hue, float _sat, float _xz, float _yz, float _zz, float _den, float _speed, float _angle) {\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
@@ -40,64 +49,87 @@ public class NDat {
 \r
 public class Noise extends DPat\r
 {\r
-       int             CurAnim = -1, numAnims = 6;\r
-       float           zTime   = random(10000), zTheta=0;\r
-       float           rtime   = 0, ttime      = 0, transAdd=0;\r
+       int                     CurAnim, iSymm;\r
        int             XSym=1,YSym=2,RadSym=3;\r
-       int                     iSymm;\r
-       ArrayList       noises  = new ArrayList();\r
-       _DhP            pSpeed , pDensity, pRotZ;\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
-               pRotZ   = addParam("RotZ"       , .5 ); pSpeed          = addParam("Fast", .55);\r
-               pDensity= addParam("Dens"       , .5);\r
-               pSymm   = addPick("Symmetry", 0, 4);    pChoose         = addPick("Animation", 1, 6);\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*(pRotZ .Val()-.5)*.01        ;\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
-                       noises.clear(); CurAnim = pChoose.Cur(); ttime = rtime;\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: noises.add(new NDat(0  ,0  ,75 ,75 ,150,45 ,3  ,0  )); pSharp.Set(1 ); break;   // drip\r
-                       case 1: noises.add(new NDat(0  ,0  ,100,100,200,45 ,3  ,180)); pSharp.Set(0 ); break;   // clouds\r
-                       case 2: noises.add(new NDat(0  ,0  ,2  ,400,2  ,20 ,3  ,0  )); pSharp.Set(.5); break;   // rain\r
-                       case 3: noises.add(new NDat(40 ,1  ,100,100,200,10 ,1  ,180)); \r
-                                       noises.add(new NDat(0  ,1  ,100,100,200,10 ,5  ,180)); pSharp.Set(0 ); break;   // fire 1\r
-                       case 4: noises.add(new NDat(0  ,1  ,40 ,40 ,40 ,15 ,2.5,180));\r
-                                       noises.add(new NDat(20 ,1  ,40 ,40 ,40 ,15 ,4  ,0  ));\r
-                                       noises.add(new NDat(40 ,1  ,40 ,40 ,40 ,15 ,2  ,90 ));\r
-                                       noises.add(new NDat(60 ,1  ,40 ,40 ,40 ,15 ,3  ,-90)); pSharp.Set(.5); break;   // machine\r
-                       case 5: noises.add(new NDat(0  ,1  ,400,100,2  ,15 ,3  ,90 ));\r
-                                       noises.add(new NDat(20 ,1  ,400,100,2  ,15 ,2.5,0  ));\r
-                                       noises.add(new NDat(40 ,1  ,100,100,2  ,15 ,2  ,180));\r
-                                       noises.add(new NDat(60 ,1  ,100,100,2  ,15 ,1.5,270)); pSharp.Set(.5); break;   // spark\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
 \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
-               xyz v   = P.RotateZ(xyzMid,zTheta);\r
-               if (iSymm == XSym && v.x > xdMax/2) v.x = xdMax-v.x;\r
-               if (iSymm == YSym && v.y > ydMax/2) v.y = ydMax-v.y;\r
-\r
-               for (int i=0;i<noises.size(); i++) {\r
-                       NDat  n     = (NDat) noises.get(i);\r
-                       float zx    = zTime * n.speed * sin(radians(n.angle)),\r
-                                 zy    = zTime * n.speed * cos(radians(n.angle));\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(v,xyzMid)/n.xz)\r
-                                                                                : noise(v.x/n.xz+zx+n.xoff,v.y/n.yz+zy+n.yoff,v.z/n.zz+n.zoff))\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
@@ -110,106 +142,429 @@ public class Noise extends DPat
 //----------------------------------------------------------------------------------------------------------------------------------\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
-       _DhP    pAmp, pRad;\r
-       _DhP    pRotX, pRotY, pRotZ;\r
-       xyz             Theta = new xyz();\r
-       Pick    pTimePattern, pTempoMult, pShape, pForm;\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
-           pAmp                = addParam("Amp" , .2);\r
-               pRotX           = addParam("RotX", .5);\r
-               pRotY           = addParam("RotY", .5);\r
-               pRotZ           = addParam("RotZ", .5);\r
-           pRad                = addParam("Rad" , .1           );\r
-               pTimePattern= addPick ("TimePattern", 0 , 5             );\r
-               pTempoMult      = addPick ("TimeMult"   , 0 , 6         );\r
-               pShape          = addPick ("Shape"              , 0 , 16        );\r
-               pForm           = addPick ("Form"               , 0 , 3         );\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;\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
-               a = pAmp.Val();\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.x = pRotX.ZeroOrOne() ? 0 :(pRotX.Val()-.5)*10;\r
-               Theta.y = pRotY.ZeroOrOne() ? 0 :(pRotY.Val()-.5)*10;\r
-               Theta.z = pRotZ.ZeroOrOne() ? 0 :(pRotZ.Val()-.5)*10;\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) { CurRandTempo = int(random(4)); } LastMeasure = t;\r
-               \r
-               switch (pTempoMult.Cur()) {\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
-                       case 5:         t = (t*pow(2,CurRandTempo))%1.;         break;\r
                }\r
 \r
-               if (t<LastBeat) { cPrev = cCur; cCur = cCur.setRand(); } LastBeat = t;\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
-               switch (pTimePattern.Cur()) {\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
+               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
-               cMid            = cPrev.interpolate(t,cCur);\r
-               cMidNorm        = cMid.setNorm();\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
-               xyz V           = new xyz();\r
-               xyz P           = Px.setNorm();\r
-               if (Theta.x>0)  P = P.RotateZ(xyzHalf, Theta.x);\r
-               if (Theta.y>0)  P = P.RotateZ(xyzHalf, Theta.y);\r
-               if (Theta.z>0)  P = P.RotateZ(xyzHalf, Theta.z);\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(P.x, P.z);\r
-               float yt        = map(t,0,1,.5-a/2,.5+a/2);\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 = new xyz(P.x, yt                                                             , P.z);                                         break;  // bouncing line\r
-                       case 1:         V = new xyz(P.x, map(cos(PI*t * P.x),-1,1,0,1)  , P.z);                                         break;  // top tap\r
-                       case 2:         V = new xyz(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 = new xyz(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
+               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
-                       case 4:         V = new xyz(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 = new xyz(P.x, a*map((P.x-.5)*(P.x-.5),0,.25,0,t-.5)+.5, P.z);                        break;  // wings\r
-                       case 6:         V = new xyz(P.x, a*map((mp -.5)*(mp -.5),0,.25,0,t-.5)+.5, P.z);                        break;  // wings\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
-                       case 7:         V = new xyz(cMid.x,cMid.y,cMid.z);\r
-                                               return color(0,0,c1c(1 - (V.distance(Px) > (pRad.getValuef()+.1)*150?1:0)) );           // sphere\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
-                       case 8:         V = new xyz(cMid.x,cMid.y,cMid.z);\r
-                                               return color(0,0,c1c(1 - CalcCone(Px,V,xyzMid) * 0.02 > .5?1:0));                                       // cone\r
+dLattice DL_;\r
+//----------------------------------------------------------------------------------------------------------------------------------\r
+class Worms extends SCPattern {\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
+       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
-               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
+       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
-                       default:        return color(0,0,c1c(P.y < V.y ?1:0));\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