/*{
  "DESCRIPTION": "A 90s bedroom computer desk at golden hour. A beige CRT sits on a laminate hutch desk between two little beige speakers, next to a beige tower with a CD-ROM, a floppy drive and a blinking hard-drive light. The screen cycles through four programs: a green-hill desktop where a browser window opens and a page trickles in over dial-up, a gem-matching game whose board cascades, the starfield screensaver, and Mystify. Sunset light comes through venetian blinds and lays stripes across the desk, the monitor and the wall, and dust drifts in the beams. Around it are a ball mouse on a blue pad, a keyboard, a can of cola, an 'I love my computer' mug, floppies, a CD throwing a rainbow, a lava lamp, a shelf of jewel cases and big-box games, a little TV playing a platformer, a stack of tapes, posters, and an iron bed under a plaid blanket. SCREEN picks the program or cycles them. ARRANGE deals the desk things and the posters into one of twelve layouts. MIRROR flips the desk so the tower sits on the left, away from the window. CLUTTER runs from a tidy desk and half-empty shelves to cans, a stack of jewel cases, a magazine on the bed and paper balls under the desk. WOOD takes the desk from birch to walnut, WALL HUE repaints the walls, and YELLOWING ages the beige plastic. SPEED drives the slow camera push. SUN, SUN HEIGHT and GOLDEN set the light, BLINDS tilts the slats, and BEAMS and DUST fill the air. LAVA and GLOW are the lamp and the screens, SCANLINES is the CRT, FILM fades everything into a drugstore print, DATE STAMP adds the orange date in the corner, and EXPOSURE sets the overall brightness.",
  "CREDIT": "CC0",
  "ISFVSN": "2",
  "CATEGORIES": ["generator", "retro", "interior", "nostalgia", "raymarch"],
  "INPUTS": [
    { "NAME": "SCREEN", "TYPE": "long", "DEFAULT": 0, "VALUES": [0, 1, 2, 3, 4], "LABELS": ["cycle", "desktop", "gem game", "starfield", "mystify"] },
    { "NAME": "ARRANGE", "TYPE": "long", "DEFAULT": 0, "MIN": 0, "MAX": 11 },
    { "NAME": "MIRROR", "TYPE": "bool", "DEFAULT": false },
    { "NAME": "CLUTTER", "TYPE": "float", "DEFAULT": 0.6, "MIN": 0.0, "MAX": 1.0 },
    { "NAME": "WOOD", "TYPE": "float", "DEFAULT": 0.5, "MIN": 0.0, "MAX": 1.0 },
    { "NAME": "WALL_HUE", "TYPE": "float", "DEFAULT": 0.0, "MIN": 0.0, "MAX": 1.0, "LABEL": "wall hue" },
    { "NAME": "YELLOWING", "TYPE": "float", "DEFAULT": 0.5, "MIN": 0.0, "MAX": 1.0 },
    { "NAME": "SPEED", "TYPE": "float", "DEFAULT": 1.0, "MIN": 0.0, "MAX": 3.0 },
    { "NAME": "SUN", "TYPE": "float", "DEFAULT": 1.0, "MIN": 0.0, "MAX": 2.5 },
    { "NAME": "SUN_HEIGHT", "TYPE": "float", "DEFAULT": 0.4, "MIN": 0.0, "MAX": 1.0, "LABEL": "sun height" },
    { "NAME": "GOLDEN", "TYPE": "float", "DEFAULT": 0.7, "MIN": 0.0, "MAX": 1.0 },
    { "NAME": "BLINDS", "TYPE": "float", "DEFAULT": 0.3, "MIN": 0.0, "MAX": 1.0 },
    { "NAME": "BEAMS", "TYPE": "float", "DEFAULT": 1.0, "MIN": 0.0, "MAX": 3.0 },
    { "NAME": "DUST", "TYPE": "float", "DEFAULT": 1.0, "MIN": 0.0, "MAX": 3.0 },
    { "NAME": "LAVA", "TYPE": "float", "DEFAULT": 1.0, "MIN": 0.0, "MAX": 2.0 },
    { "NAME": "GLOW", "TYPE": "float", "DEFAULT": 1.0, "MIN": 0.0, "MAX": 2.0 },
    { "NAME": "SCANLINES", "TYPE": "float", "DEFAULT": 0.6, "MIN": 0.0, "MAX": 1.0 },
    { "NAME": "FILM", "TYPE": "float", "DEFAULT": 0.35, "MIN": 0.0, "MAX": 1.0 },
    { "NAME": "DATE_STAMP", "TYPE": "bool", "DEFAULT": true, "LABEL": "date stamp" },
    { "NAME": "EXPOSURE", "TYPE": "float", "DEFAULT": 1.0, "MIN": 0.3, "MAX": 3.0 }
  ]
}*/
// ---------------------------------------------------------------------------
// Layout, metres, y up. The camera looks down +z at the back wall (z = 0.95).
// The right wall is the plane x = 1.25 with a window cut in it, and venetian
// blinds hang just inside it at x = BLX. The desk is a hutch 1.6 m wide
// (x in [-0.8, 0.8]), with its top at y = 0.75 and a shelf at y = 1.30.
//
// Sunlight is never marched through the blinds. sunMask() projects a point
// along the sun direction onto the window and onto the blind plane and works
// out analytically how much of each slat's shadow it is in. So the stripes
// are exact, they soften with distance from the window, and the volumetric
// beams cost no map() calls at all. Objects in the room cast their shadows
// with one ordinary soft-shadow march.
//
// One loop does the march, the four normal taps, four AO taps and the shadow
// ray, and map() has a single call site. That keeps the compiled shader to one
// inlined copy of the scene.
// ---------------------------------------------------------------------------
#define PI 3.14159265
#define TAU 6.28318531
#define MX (-0.12)     // monitor centre x
#define SCY 1.015      // screen centre y
#define LX (-0.66)     // lava lamp
#define LZ 0.75
#define BLX 1.215      // plane the blinds hang in
#define SP 0.024       // slat pitch
#define SW 0.026       // slat width
vec3 gL; // toward the sun
vec3 gSun; // sun colour * intensity
float gCa; // slat angle cos/sin
float gSa;
float gMode; // 0 desktop, 1 gems, 2 starfield
float gMT; // time inside the current program
float gFade; // program switch blackout
vec3 gScr; // average colour the monitor throws into the room
float gS; // 1, or -1 with the desk mirrored (MIRROR)
float gSeed; // ARRANGE
vec2 gMug; // where the small things landed, desk frame xz
vec2 gFlp;
vec2 gCD;
vec2 gCan;
vec2 gJc; // the spare spots: jewel-case stack and a second can
vec2 gCan2;
vec4 gOn; // mug, floppies, CD, can on the desk at this CLUTTER
float gB; // nearest bounding volume (map)
float gR; // nearest real surface (map)
float hash11(float p) { p=fract(p*0.1031); p*=p+33.33; p*=p+p; return fract(p); }
float hash21(vec2 p) {
 vec3 p3=fract(vec3(p.xyx)*0.1031);
 p3+=dot(p3,p3.yzx+33.33);
 return fract((p3.x+p3.y)*p3.z);
}
vec2 hash22(vec2 p) {
 vec3 p3=fract(vec3(p.xyx)*vec3(0.1031,0.1030,0.0973));
 p3+=dot(p3,p3.yzx+33.33);
 return fract((p3.xx+p3.yz)*p3.zy);
}
float noise2(vec2 p) {
 vec2 i=floor(p),f=fract(p);
 vec2 u=f*f*(3.0-2.0*f);
 return mix(mix(hash21(i),hash21(i+vec2(1.0,0.0)),u.x),
   mix(hash21(i+vec2(0.0,1.0)),hash21(i+vec2(1.0,1.0)),u.x),u.y);
}
float fbm2(vec2 p) {
 float a=0.5,s=0.0;
 for (int i=0; i<4; i++) { s+=a*noise2(p); p=mat2(1.6,1.2,-1.2,1.6)*p+3.1; a*=0.5; }
 return s;
}
// ---- SDF primitives ---------------------------------------------------------
float sdBox(vec3 p,vec3 b) {
 vec3 q=abs(p)-b;
 return length(max(q,0.0))+min(max(q.x,max(q.y,q.z)),0.0);
}
float sdRBox(vec3 p,vec3 b,float r) { return sdBox(p,b-r)-r; }
float sdBox2(vec2 p,vec2 b) {
 vec2 q=abs(p)-b;
 return length(max(q,0.0))+min(max(q.x,q.y),0.0);
}
float sdCylY(vec3 p,float r,float h) {
 vec2 d=abs(vec2(length(p.xz),p.y))-vec2(r,h);
 return min(max(d.x,d.y),0.0)+length(max(d,0.0));
}
float sdEll(vec3 p,vec3 r) {
 float k0=length(p/r);
 float k1=length(p/(r*r));
 return k0*(k0-1.0)/k1;
}
float sdCap(vec3 p,vec3 a,vec3 b,float r) {
 vec3 pa=p-a,ba=b-a;
 float h=clamp(dot(pa,ba)/dot(ba,ba),0.0,1.0);
 return length(pa-ba*h)-r;
}
vec2 opU(vec2 a,vec2 b) { return a.x<b.x?a:b; }
// three-stop ramp: a at 0, b at 0.5, c at 1
vec3 ramp3(vec3 a,vec3 b,vec3 c,float x) { return x<0.5?mix(a,b,x*2.0):mix(b,c,x*2.0-1.0); }
vec3 hueRot(vec3 c,float a) {
 vec3 k=vec3(0.57735);
 return c*cos(a)+cross(k,c)*sin(a)+k*dot(k,c)*(1.0-cos(a));
}
// six spots on the desk the small things get dealt to (desk frame, xz)
vec2 slot(float k) {
 k=mod(k,6.0);
 return k<0.5?vec2(-0.62,0.36):k<1.5?vec2(-0.44,0.29):k<2.5?vec2(0.48,0.27)
  :k<3.5?vec2(0.66,0.30):k<4.5?vec2(-0.64,0.53):vec2(0.73,0.43);
}
mat2 rot(float a) { float c=cos(a),s=sin(a); return mat2(c,s,-s,c); }
// Materials
//  1 laminate   2 wall       3 carpet     4 beige plastic  5 CRT glass
//  6 jewel case 7 black      8 cola can   9 mug           10 lava glass
// 11 metal     12 slats     13 white paint 14 blanket     15 pillow
// 16 iron      17 game box  18 mousepad   19 floppy       20 keys
// 21 TV glass  22 CD        23 VHS
// Keycaps on top of the keyboard, q in the keyboard's frame.
float keysSDF(vec3 q) {
 float pt=0.0185;
 // main block, 15 x 6 (the back row is the F-keys)
 vec2 c=vec2(clamp(floor((q.x+0.218)/pt),0.0,14.0),
    clamp(floor((q.z+0.066)/pt),0.0,5.0));
 float kz=-0.066+(c.y+0.5)*pt+(c.y>4.5?0.003:0.0);
 float kx=-0.218+(c.x+0.5)*pt;
 float hx=0.0072;
 if (c.y<0.5&&c.x>3.5&&c.x<10.5) { kx=-0.218+7.5*pt; hx=3.5*pt-0.0021; }
 float d=sdRBox(vec3(q.x-kx,q.y-0.016,q.z-kz),vec3(hx,0.005,0.0072),0.002);
 // arrows, nav cluster and numpad; column 3 is the gap
 vec2 c2=vec2(clamp(floor((q.x-0.066)/pt),0.0,7.0),c.y);
 vec3 k2=vec3(q.x-0.066-(c2.x+0.5)*pt,q.y-0.016,q.z-kz);
 bool on=c2.x>3.5?c2.y<4.5
   :(c2.x<2.5&&(c2.y>2.5||c2.y<0.5||(c2.y<1.5&&c2.x>0.5&&c2.x<1.5)));
 float d2=on?sdRBox(k2,vec3(0.0072,0.005,0.0072),0.002)
    :max(max(-sdBox2(k2.xz,vec2(0.5*pt)),0.0)+0.0015,q.y-0.022);
 return min(d,d2);
}
vec2 map(vec3 p) {
 gB=1e5;
 // ---- room shell
 vec2 res=vec2(max(p.y,p.x-1.45),3.0);
 res=opU(res,vec2(sdBox(p-vec3(-0.6,1.3,1.05),vec3(2.05,1.4,0.1)),2.0));
 res=opU(res,vec2(p.x+2.6,2.0));
 res=opU(res,vec2(max(2.6-p.y,p.x-1.45),2.0));
 float rw=sdBox(p-vec3(1.35,1.3,-1.0),vec3(0.1,1.4,2.6));
 rw=max(rw,-sdBox(p-vec3(1.35,1.5,-0.05),vec3(0.2,0.55,0.6)));
 res=opU(res,vec2(rw,2.0));
 // ---- window trim, sill and the blinds
 float bw=sdBox(p-vec3(1.22,1.5,-0.05),vec3(0.07,0.62,0.72));
 if (bw<0.05) {
  vec3 wq=p-vec3(1.245,1.5,-0.05);
  float trim=max(sdBox(wq,vec3(0.012,0.6,0.65)),-sdBox(wq,vec3(0.1,0.55,0.6)));
  trim=min(trim,sdBox(p-vec3(1.215,0.935,-0.05),vec3(0.04,0.012,0.67)));
  trim=min(trim,sdBox(p-vec3(BLX,2.02,-0.05),vec3(0.022,0.022,0.6)));
  trim=min(trim,sdBox(p-vec3(BLX,1.0,-0.05),vec3(0.016,0.007,0.59)));
  res=opU(res,vec2(trim,13.0));
  float u=p.x-BLX;
  float v=mod(p.y,SP)-0.5*SP;
  vec2 sl=vec2(gCa*u+gSa*v,-gSa*u+gCa*v);
  float ds=sdBox2(sl,vec2(0.5*SW,0.0008));
  ds=max(ds,max(abs(p.z+0.05)-0.585,abs(p.y-1.5)-0.49));
  res=opU(res,vec2(ds,12.0));
 } else gB=min(gB,bw);
 // ---- the bed: plaid blanket, pillow, iron headboard against the wall
 float bb=sdBox(p-vec3(-1.65,0.6,-0.35),vec3(0.72,0.6,1.33));
 if (bb<0.05) {
  vec3 bq=p-vec3(-1.65,0.42,-0.35);
  float wr=0.006*sin(bq.z*23.0+2.0*sin(bq.x*9.0))*smoothstep(-0.1,0.1,bq.y);
  res=opU(res,vec2((sdRBox(bq,vec3(0.66,0.14,1.28),0.07)+wr)*0.9,14.0));
  res=opU(res,vec2(sdBox(p-vec3(-1.65,0.15,-0.35),vec3(0.62,0.13,1.25)),7.0));
  res=opU(res,vec2(sdEll(p-vec3(-1.65,0.62,0.70),vec3(0.34,0.075,0.17)),15.0));
  vec3 hq3=p-vec3(-1.65,0.0,0.905);
  float bx=hq3.x-0.1*clamp(floor(hq3.x/0.1+0.5),-6.0,6.0);
  float iron=max(length(vec2(bx,hq3.z))-0.010,abs(p.y-0.78)-0.27);
  iron=min(iron,sdCap(p,vec3(-2.31,1.05,0.905),vec3(-0.99,1.05,0.905),0.016));
  vec3 pq=vec3(abs(hq3.x)-0.66,p.y,hq3.z);
  iron=min(iron,sdCap(pq,vec3(0.0),vec3(0.0,1.08,0.0),0.02));
  iron=min(iron,length(pq-vec3(0.0,1.11,0.0))-0.032);
  res=opU(res,vec2(iron,16.0));
 } else gB=min(gB,bb);
 // ---- past CLUTTER 0.7: a magazine left on the bed, paper balls under the desk
 if (CLUTTER>0.7) {
  vec3 mq=p-vec3(-1.3,0.569,0.3);
  mq.xz=rot(0.5)*mq.xz;
  res=opU(res,vec2(sdBox(mq,vec3(0.105,0.003,0.14)),17.0));
  if (CLUTTER>0.82) {
   float pb=length(p-vec3(-0.42,0.034,0.62))-0.034;
   if (CLUTTER>0.94) pb=min(pb,length(p-vec3(0.15,0.03,0.7))-0.03);
   pb+=0.005*sin(p.x*90.0)*sin(p.y*80.0+1.0)*sin(p.z*85.0);
   res=opU(res,vec2(pb*0.7,15.0));
  }
 }
 // ---- the desk and everything on it (last: MIRROR flips p for it)
 float bd=sdBox(p-vec3(0.0,0.9,0.575),vec3(0.82,0.9,0.39));
 if (bd<0.05) {
  p.x*=gS;
  vec3 q=vec3(abs(p.x),p.y,p.z);
  float d=sdBox(p-vec3(0.0,0.735,0.575),vec3(0.80,0.015,0.375));
  d=min(d,sdBox(q-vec3(0.785,0.36,0.575),vec3(0.015,0.36,0.375)));
  d=min(d,sdBox(p-vec3(0.0,0.42,0.925),vec3(0.8,0.3,0.012)));
  d=min(d,sdBox(q-vec3(0.785,1.25,0.80),vec3(0.015,0.5,0.15)));
  d=min(d,sdBox(p-vec3(0.0,1.30,0.80),vec3(0.80,0.012,0.15)));
  d=min(d,sdBox(p-vec3(0.0,1.75,0.795),vec3(0.81,0.015,0.16)));
  d=min(d,sdBox(p-vec3(0.42,1.525,0.80),vec3(0.012,0.213,0.15)));
  d=min(d,sdBox(p-vec3(0.0,1.25,0.94),vec3(0.8,0.5,0.008)));
  res=opU(res,vec2(d,1.0));
  // things standing on the desktop
  float bi=sdBox(p-vec3(0.0,0.975,0.575),vec3(0.8,0.23,0.38));
  if (bi<0.03) {
   // CRT monitor: bezel with the tube recess, tapering bell, swivel base
   float bm=sdBox(p-vec3(MX,0.99,0.70),vec3(0.22,0.24,0.24));
   if (bm<0.03) {
    float bez=sdRBox(p-vec3(MX,1.0,0.515),vec3(0.205,0.19,0.035),0.016);
    vec3 rq=p-vec3(MX,SCY,0.48);
    bez=max(bez,-sdRBox(rq,vec3(0.168,0.128,0.03),0.01));
    float zz=clamp((p.z-0.55)/0.37,0.0,1.0);
    float s=1.0-0.38*zz;
    vec3 hq=p-vec3(MX,1.0-0.03*zz,0.735);
    hq.xy/=s;
    bez=min(bez,(sdBox(hq,vec3(0.178,0.166,0.173))-0.012)*s*0.9);
    bez=min(bez,sdCylY(p-vec3(MX,0.762,0.69),0.128,0.01)-0.002);
    bez=min(bez,sdRBox(p-vec3(MX,0.795,0.69),vec3(0.08,0.03,0.09),0.01));
    res=opU(res,vec2(bez,4.0));
    float gl=length(p-vec3(MX,SCY,0.488+1.4))-1.4;
    gl=max(gl,sdBox(rq-vec3(0.0,0.0,0.02),vec3(0.17,0.13,0.04)));
    res=opU(res,vec2(gl,5.0));
   } else gB=min(gB,bm);
   // tower: body, front bezel, drive bays
   float bt=sdBox(p-vec3(0.58,0.96,0.665),vec3(0.11,0.225,0.23));
   if (bt<0.03) {
    float tw=sdRBox(p-vec3(0.58,0.96,0.67),vec3(0.095,0.21,0.21),0.006);
    tw=min(tw,sdRBox(p-vec3(0.58,0.96,0.458),vec3(0.097,0.211,0.008),0.005));
    tw=max(tw,-sdBox(p-vec3(0.58,1.105,0.45),vec3(0.074,0.021,0.004)));
    tw=max(tw,-sdBox(p-vec3(0.58,1.01,0.45),vec3(0.052,0.012,0.004)));
    res=opU(res,vec2(tw,4.0));
   } else gB=min(gB,bt);
   // speakers either side of the monitor
   vec3 sq=p-vec3(p.x<-0.13?-0.47:0.21,0.865,0.60);
   res=opU(res,vec2(sdRBox(sq,vec3(0.055,0.115,0.07),0.008),4.0));
   // keyboard, tilted up at the back, with its keycaps
   float bk=sdBox(p-vec3(-0.10,0.77,0.315),vec3(0.24,0.035,0.09));
   if (bk<0.02) {
    vec3 kq=p-vec3(-0.10,0.762,0.315);
    kq.y-=0.08*kq.z;
  kq.x*=gS;
    res=opU(res,vec2(sdRBox(kq,vec3(0.232,0.011,0.078),0.006),4.0));
    res=opU(res,vec2(keysSDF(kq),20.0));
   } else gB=min(gB,bk);
   // ball mouse on a mousepad
   float bmo=sdBox(p-vec3(0.30,0.765,0.31),vec3(0.115,0.03,0.095));
   if (bmo<0.02) {
    vec3 pq=p-vec3(0.30,0.7515,0.31);
    float pad=max(sdBox2(pq.xz,vec2(0.10,0.08))-0.01,abs(pq.y)-0.0015);
    res=opU(res,vec2(pad,18.0));
    vec3 mo=p-vec3(0.31,0.753,0.30);
    mo.xz=rot(0.25)*mo.xz;
    float ms=max(sdEll(mo-vec3(0.0,0.004,0.0),vec3(0.031,0.021,0.056)),-mo.y);
    res=opU(res,vec2(ms,4.0));
   } else gB=min(gB,bmo);
   // mouse and keyboard leads
   float cab=sdCap(p,vec3(0.322,0.757,0.352),vec3(0.345,0.7535,0.44),0.0022);
   cab=min(cab,sdCap(p,vec3(0.345,0.7535,0.44),vec3(0.40,0.7535,0.62),0.0022));
   cab=min(cab,sdCap(p,vec3(-0.10,0.774,0.394),vec3(-0.085,0.7535,0.46),0.0024));
   cab=min(cab,sdCap(p,vec3(-0.085,0.7535,0.46),vec3(-0.04,0.7535,0.6),0.0024));
   res=opU(res,vec2(cab,4.0));
   // cola can
   if (gOn.w>0.5) res=opU(res,vec2(sdCylY(p-vec3(gCan.x,0.8115,gCan.y),0.031,0.0595)-0.002,8.0));
   // coffee mug with a handle
   vec3 uq=p-vec3(gMug.x,0.797,gMug.y);
   float mug=gOn.x>0.5?sdCylY(uq,0.04,0.047):1.0;
   mug=max(mug,-sdCylY(uq-vec3(0.0,0.035,0.0),0.035,0.047));
   if (gOn.x>0.5) mug=min(mug,sdCylY(uq-vec3(0.0,0.028,0.0),0.036,0.002));
   vec3 hq2=uq-vec3(0.032,0.0,-0.03);
   hq2.xz=rot(-0.7)*hq2.xz;
   if (gOn.x>0.5) mug=min(mug,length(vec2(length(hq2.xy)-0.024,hq2.z))-0.0065);
   res=opU(res,vec2(mug,9.0));
   // lava lamp: a revolved profile, base / glass / cap
   vec3 lq=p-vec3(LX,0.75,LZ);
   float lr=length(lq.xz);
   if (lr<0.1) {
    float y=lq.y;
    float r=y<0.12?mix(0.058,0.036,y/0.12)
      :y<0.19?mix(0.036,0.05,(y-0.12)/0.07)
      :y<0.35?mix(0.05,0.028,(y-0.19)/0.16)
      :mix(0.028,0.017,(y-0.35)/0.06);
    float lav=max((lr-r)*0.85,abs(y-0.205)-0.205);
    res=opU(res,vec2(lav,(y>0.12&&y<0.35)?10.0:11.0));
   } else gB=min(gB,lr-0.09);
   // two floppies and a CD lying about
   if (gOn.y>0.5) {
    vec3 fq=p-vec3(gFlp.x,0.7517,gFlp.y-0.005);
    fq.xz=rot(0.35)*fq.xz;
    float fl=sdRBox(fq,vec3(0.045,0.0016,0.047),0.001);
    vec3 fq2=p-vec3(gFlp.x-0.01,0.7551,gFlp.y+0.01);
    fq2.xz=rot(-0.2)*fq2.xz;
    fl=min(fl,sdRBox(fq2,vec3(0.045,0.0016,0.047),0.001));
    res=opU(res,vec2(fl,19.0));
   }
   if (CLUTTER>0.75) {
    vec3 sq2=p-vec3(gJc.x,0.7558,gJc.y);
    float jcs=1e3;
    for (int i=0; i<3; i++) {
     vec3 cq=sq2-vec3(0.0,0.0102*float(i),0.0);
     cq.xz=rot(0.3*sin(float(i)*2.4+gSeed))*cq.xz;
     jcs=min(jcs,sdBox(cq,vec3(0.0705,0.0049,0.0625)));
    }
    res=opU(res,vec2(jcs,6.0));
   }
   if (CLUTTER>0.88) res=opU(res,vec2(sdCylY(p-vec3(gCan2.x,0.8115,gCan2.y),0.031,0.0595)-0.002,8.0));
   vec3 dq=p-vec3(gCD.x,0.7508,gCD.y);
   if (gOn.z>0.5) res=opU(res,vec2(max(sdCylY(dq,0.06,0.0008),0.0075-length(dq.xz)),22.0));
  } else gB=min(gB,bi);
  // shelf: jewel cases, big-box games, a little TV, tapes in the cubby
  float bs=sdBox(p-vec3(0.0,1.46,0.80),vec3(0.79,0.16,0.15));
  if (bs<0.03) {
   float jc=clamp(floor((p.x+0.77)/0.0104),0.0,47.0);
   vec3 jq=p-vec3(-0.77+(jc+0.5)*0.0104,1.3745,0.80);
   float jd=sdBox(jq,vec3(0.0049,0.0625,0.0705));
   // a missing case: distance to the nearest neighbour's slab (lower bound)
   if (hash11(jc*3.1+gSeed*7.3)>0.4+CLUTTER)
    jd=max(abs(jq.x)<0.0052?0.0055-abs(jq.x):abs(jq.x)+0.0055,length(max(abs(jq.yz)-vec2(0.0625,0.0705),0.0)));
   res=opU(res,vec2(jd,6.0));
   float gc=clamp(floor((p.x+0.255)/0.05),0.0,5.0);
   float gh=0.098+0.024*hash11(gc*3.7+1.0+gSeed);
   vec3 gq=p-vec3(-0.255+(gc+0.5)*0.05,1.312+gh,0.85);
   float gd=sdBox(gq,vec3(0.0235,gh,0.088));
   if (hash11(gc*5.3+gSeed*3.1)>0.4+CLUTTER)
    gd=max(abs(gq.x)<0.025?0.0265-abs(gq.x):abs(gq.x)+0.0265,length(max(abs(vec2(p.y-1.434,gq.z))-vec2(0.122,0.088),0.0)));
   res=opU(res,vec2(gd,17.0));
   vec3 tsq=p-vec3(0.215,1.47,0.66);
   float tv=sdRBox(p-vec3(0.24,1.455,0.80),vec3(0.155,0.142,0.14),0.014);
   tv=max(tv,-sdRBox(tsq,vec3(0.105,0.08,0.02),0.008));
   res=opU(res,vec2(tv,7.0));
   float tg=length(p-vec3(0.215,1.47,0.668+0.8))-0.8;
   tg=max(tg,sdBox(tsq-vec3(0.0,0.0,0.02),vec3(0.106,0.081,0.03)));
   res=opU(res,vec2(tg,21.0));
   float vc=clamp(floor((p.y-1.312)/0.026),0.0,min(floor(1.5+6.0*CLUTTER),5.0));
   vec3 vq=p-vec3(0.60+0.012*(hash11(vc*5.1+gSeed)-0.5),1.312+(vc+0.5)*0.026,0.86);
   res=opU(res,vec2(sdBox(vq,vec3(0.094,0.0125,0.0515)),23.0));
  } else gB=min(gB,bs);
 } else gB=min(gB,bd);
 // gR is the nearest real surface; the march steps by the bounds as well,
 // but shadows and AO must never mistake a bound for a surface
 gR=res.x;
 res.x=min(res.x,gB);
 return res;
}
// ---- sunlight through the window and the blinds -----------------------------
// Fraction of direct sun reaching p. The aperture is the window opening seen
// along the sun direction; the blinds shadow it with stripes whose width
// follows from the slat angle (half-width w/2 |sin a - cos a tan(elevation)|)
// and whose edges blur with distance from the slats.
float sunMask(vec3 p,bool slats) {
 float s=(1.25-p.x)/gL.x;
 if (s<0.0) return 0.0;
 vec3 q=p+gL*s;
 float soft=0.003+0.012*s;
 float ap=smoothstep(-soft,soft,0.6-abs(q.z+0.05))*smoothstep(-soft,soft,0.55-abs(q.y-1.5));
 if (!slats) return ap;
 float s2=max((BLX-p.x)/gL.x,0.0);
 float y2=p.y+gL.y*s2;
 // the gap between two slat shadows is a band of half-width o round the
 // gap centre; blur it as box-filtered coverage so a wide blur tends to the
 // mean open fraction and a closed blind stays closed
 float e=abs(mod(y2+0.5*SP,SP)-0.5*SP);
 float o=max(0.5*SP-0.5*SW*abs(gSa-gCa*gL.y/gL.x),0.0);
 float sf=min(0.0008+0.011*s2,SP);
 float hb=max(o,0.5*sf);
 return ap*max(smoothstep(hb+0.5*sf,hb-0.5*sf,e)*o/hb,0.025);
}
// open length of the slat pattern below height y (cumulative)
float slatCum(float y,float hw) {
 float o=0.5*SP-hw;
 float u=mod(y,SP);
 return floor(y/SP)*2.0*o+min(u,o)+max(u-SP+o,0.0);
}
// Length of the view ray that lies in direct sun. A point's projection along
// the sun onto the window is linear in distance along the ray, so the lit
// stretch is a two-slab intersection, and the slat stripes (also linear in
// distance) integrate in closed form. Exact, and no sampling noise.
float beamLen(vec3 ro,vec3 rd,float tmax,float g) {
 float kx=1.0/gL.x;
 float s0=0.0,s1=tmax;
 if (rd.x>0.0) s1=min(s1,(1.25-ro.x)/rd.x);
 float ay=ro.y+gL.y*(1.25-ro.x)*kx-1.5,by=rd.y-gL.y*rd.x*kx;
 float az=ro.z+gL.z*(1.25-ro.x)*kx+0.05,bz=rd.z-gL.z*rd.x*kx;
 float hy=0.55+g,hz=0.6+g;
 if (abs(by)<1e-6) { if (abs(ay)>hy) return 0.0; }
 else { float t0=(-hy-ay)/by,t1=(hy-ay)/by; s0=max(s0,min(t0,t1)); s1=min(s1,max(t0,t1)); }
 if (abs(bz)<1e-6) { if (abs(az)>hz) return 0.0; }
 else { float t0=(-hz-az)/bz,t1=(hz-az)/bz; s0=max(s0,min(t0,t1)); s1=min(s1,max(t0,t1)); }
 if (s1<=s0) return 0.0;
 float hw=min(0.5*SW*abs(gSa-gCa*gL.y*kx),0.5*SP);
 float y0=ro.y+gL.y*(BLX-ro.x)*kx;
 float ya=y0+by*s0,yb=y0+by*s1;
 float open=abs(yb-ya)<1e-4?(s1-s0)*step(hw,abs(mod(ya,SP)-0.5*SP))
         :(slatCum(yb,hw)-slatCum(ya,hw))/by;
 return open+0.025*(s1-s0);
}
// light seen looking out through the blinds (for reflections in the glass)
vec3 windowSeen(vec3 ro,vec3 rd) {
 if (rd.x<=0.0) return vec3(0.0);
 float s=(BLX-ro.x)/rd.x;
 vec3 q=ro+rd*s;
 float w=step(abs(q.z+0.05),0.58)*step(abs(q.y-1.5),0.5);
 float v=abs(mod(q.y,SP)-0.5*SP)/(0.5*SP);
 return gSun*w*(0.18+0.5*smoothstep(0.55,0.9,v));
}
vec3 skyOut(vec3 rd) {
 vec3 c=mix(vec3(1.0,0.66,0.36),vec3(0.62,0.52,0.72),smoothstep(0.0,0.6,rd.y));
 c=mix(c,vec3(0.3,0.2,0.18),smoothstep(0.0,-0.12,rd.y));
 float sd=max(dot(rd,gL),0.0);
 return c*(0.5+0.5*SUN)*1.4+gSun*(0.4*pow(sd,40.0)+0.08*pow(sd,4.0));
}
// ---- screen programs --------------------------------------------------------
// pseudo text: word-shaped strokes. x in character cells, y 0..1 in the line.
// blur (0..1) fades the strokes to their average when they are sub-pixel.
float ptext(vec2 p,float seed,float blur) {
 float c=floor(p.x);
 float on=step(0.2,hash21(vec2(c,seed)));
 float tall=step(0.72,hash21(vec2(c,seed+7.0)));
 float body=step(0.22,p.y)*step(p.y,0.68+0.22*tall)*step(0.0,p.y)*step(p.y,1.0);
 float g=step(0.18,fract(p.x))*step(fract(p.x),0.82);
 return on*mix(body*g,0.28*step(0.15,p.y)*step(p.y,0.85),blur);
}
float segMask(float d) {
 // seven-segment masks, bit 0 = top, then clockwise, bit 6 = middle
 return d<0.5?63.0:d<1.5?6.0:d<2.5?91.0:d<3.5?79.0:d<4.5?102.0
  :d<5.5?109.0:d<6.5?125.0:d<7.5?7.0:d<8.5?127.0:111.0;
}
// signed distance to a seven-segment digit; x in [-0.5,0.5], y in [-1,1]
float seg7(vec2 p,float mask) {
 p.x-=0.12*p.y;
 float d=1e3;
 for (int i=0; i<7; i++) {
  float fi=float(i);
  if (mod(floor(mask/exp2(fi)),2.0)<0.5) continue;
  vec2 c=fi<0.5?vec2(0.0,1.0):fi<1.5?vec2(0.5,0.5):fi<2.5?vec2(0.5,-0.5)
   :fi<3.5?vec2(0.0,-1.0):fi<4.5?vec2(-0.5,-0.5):fi<5.5?vec2(-0.5,0.5):vec2(0.0);
  bool hz=fi<0.5||(fi>2.5&&fi<3.5)||fi>5.5;
  d=min(d,sdBox2(p-c,hz?vec2(0.34,0.075):vec2(0.075,0.34)));
 }
 return d;
}
// Desktop: rolling green hill, icons, taskbar, a browser window whose page
// loads top to bottom over a modem, and the arrow cursor wandering about.
vec3 desktop(vec2 s,float t,float fp) {
 t=mod(t,14.0);
 vec2 a=vec2(s.x*1.3333,s.y);
 vec3 col=mix(vec3(0.50,0.72,0.98),vec3(0.10,0.34,0.85),smoothstep(0.25,1.0,s.y));
 float cl=fbm2(vec2(a.x*2.6+0.3,a.y*6.0-1.0));
 col=mix(col,vec3(1.0),smoothstep(0.55,0.78,cl)*smoothstep(0.38,0.75,s.y)*0.95);
 float hy=0.30+0.20*exp(-pow((s.x-0.33)/0.40,2.0))-0.06*s.x;
 float hy2=0.33+0.07*exp(-pow((s.x-0.95)/0.22,2.0));
 if (s.y<max(hy,hy2)) {
  vec3 g=mix(vec3(0.17,0.42,0.05),vec3(0.45,0.74,0.15),clamp(s.y/hy,0.0,1.0)*0.8+0.2*(1.0-s.x));
  g=mix(g*0.78+vec3(0.0,0.0,0.04),g,step(s.y,hy));
  g*=0.88+0.22*noise2(a*vec2(30.0,60.0));
  col=g;
 }
 float blur=clamp(fp*150.0-0.5,0.0,1.0);
 // icons down the left
 float row=floor((0.97-s.y)/0.13);
 if (s.x<0.09&&row>=0.0&&row<4.0) {
  vec2 c=vec2(s.x-0.045,0.97-s.y-row*0.13-0.04);
  vec3 ic=row<0.5?vec3(0.85,0.88,0.95):row<1.5?vec3(0.95,0.8,0.25)
    :row<2.5?vec3(0.2,0.55,0.95):vec3(0.3,0.75,0.35);
  col=mix(col,ic*(1.0-5.0*c.y),step(sdBox2(c*vec2(1.333,1.0),vec2(0.025,0.025)),0.0));
  float lab=step(abs(c.y-0.055),0.009)*step(abs(c.x),0.035);
  col=mix(col,vec3(1.0),lab*mix(step(0.4,fract(s.x*120.0+row)),0.6,blur));
 }
 // taskbar
 if (s.y<0.055) {
  float v=s.y/0.055;
  col=mix(vec3(0.10,0.27,0.80),vec3(0.24,0.46,0.96),smoothstep(0.0,0.9,v))+0.15*smoothstep(0.85,1.0,v);
  if (s.x<0.12) {
   col=mix(vec3(0.18,0.52,0.16),vec3(0.36,0.74,0.30),v);
   col=mix(col,vec3(1.0),ptext(vec2((s.x-0.03)*70.0,(v-0.2)/0.6),3.0,blur)*step(s.x,0.1));
  }
  if (s.x>0.14&&s.x<0.36) col=mix(col,vec3(0.30,0.52,0.98),0.7*step(0.12,v)*step(v,0.88));
  if (s.x>0.85) col=mix(vec3(0.08,0.52,0.92),vec3(0.24,0.68,0.98),v);
 }
 // browser window
 float op=smoothstep(0.8,1.3,t);
 vec2 wh=vec2(0.32,0.36)*mix(0.15,1.0,op);
 vec2 w=s-vec2(0.52,0.55);
 if (op>0.01&&abs(w.x)<wh.x&&abs(w.y)<wh.y) {
  vec2 u=(w+wh)/(2.0*wh);
  float ty=1.0-u.y;
  float load=clamp((t-1.6)/8.0,0.0,1.0);
  vec3 wc=vec3(0.93,0.92,0.87);
  if (ty<0.06) {
   wc=mix(vec3(0.0,0.30,0.86),vec3(0.24,0.55,1.0),0.6*(1.0-ty/0.06)+0.2*u.x);
   wc=mix(wc,vec3(1.0),ptext(vec2((u.x-0.05)*60.0,(0.05-ty)/0.04),9.0,blur)*step(u.x,0.4)*step(0.05,u.x));
   float bx=u.x-0.84;
   if (bx>0.0&&abs(ty-0.03)<0.022&&mod(bx,0.055)<0.045)
    wc=bx>0.11?vec3(0.85,0.25,0.15):vec3(0.25,0.5,1.0);
  } else if (ty<0.12) {
   wc*=0.97;
   float ib=mod(u.x,0.07);
   if (ib<0.045&&u.x<0.6&&abs(ty-0.09)<0.02) wc=mix(vec3(0.3,0.55,0.3),vec3(0.3,0.4,0.75),step(0.3,u.x));
  } else if (ty<0.16) {
   if (u.x>0.08&&u.x<0.96&&abs(ty-0.14)<0.014) {
    wc=vec3(1.0);
    wc=mix(wc,vec3(0.1,0.1,0.4),ptext(vec2((u.x-0.1)*90.0,(0.15-ty)/0.018+0.5),4.0,blur)*step(u.x,0.45));
   }
  } else if (ty>0.955) {
   float seg=floor(u.x*60.0);
   if (u.x>0.62&&u.x<0.95&&abs(ty-0.977)<0.012&&mod(seg,1.0)==0.0)
    wc=(u.x-0.62)/0.33<load?vec3(0.1,0.2,0.75):wc*0.9;
  } else {
   float pv=(ty-0.16)/0.795;
   wc=vec3(1.0);
   if (pv<load) {
    if (pv<0.17) {
     wc=mix(vec3(0.08,0.06,0.4),vec3(0.55,0.1,0.55),u.x);
     wc=mix(wc,vec3(1.0,0.85,0.2),ptext(vec2((u.x-0.25)*22.0,(0.13-pv)/0.1),11.0,blur)*step(0.25,u.x)*step(u.x,0.7));
    } else if (u.x<0.2) {
     wc=vec3(0.86,0.9,1.0);
     float lr=fract((pv-0.2)*22.0);
     wc=mix(wc,vec3(0.0,0.0,0.8),ptext(vec2((u.x-0.02)*50.0,lr*1.6-0.3),floor((pv-0.2)*22.0),blur)*step(0.2,pv));
    } else {
     // an image that comes in interlaced, then lines of text
     vec2 ib=vec2(u.x-0.62,pv-0.36);
     if (abs(ib.x)<0.33&&abs(ib.y)<0.15) {
      float rows=floor(mod(t*6.0,1.0)*1.0);
      float il=step(0.5,fract(pv*60.0))+step(0.55,load);
      vec3 img=mix(vec3(0.95,0.5,0.2),vec3(0.2,0.3,0.7),ib.y/0.3+0.5);
      img=mix(img,vec3(0.1,0.15,0.1),step(ib.y+0.15,0.08+0.05*sin(ib.x*20.0)));
      wc=mix(vec3(0.75),img,clamp(il+rows,0.0,1.0));
     } else {
      float lr=fract((pv-0.53)*26.0);
      wc=mix(wc,vec3(0.12),ptext(vec2((u.x-0.24)*70.0,lr*1.5-0.25),floor(pv*26.0)+40.0,blur)
        *step(0.53,pv)*step(u.x,0.95-0.3*step(0.8,hash11(floor(pv*26.0)))));
     }
    }
   }
  }
  // bevel
  float e=min(min(u.x,1.0-u.x),min(u.y,1.0-u.y));
  col=mix(vec3(0.0,0.25,0.75),wc,step(0.006,e));
 }
 // cursor
 vec2 cp=vec2(0.55+0.3*sin(t*0.45),0.45+0.28*sin(t*0.31+1.0));
 vec2 q=(s-cp)*vec2(1.333,-1.0)/0.04;
 float inA=step(0.0,q.x)*step(q.x,0.62*q.y)*step(q.y,1.0-0.25*q.x);
 vec2 q2=q*1.18+vec2(0.07,0.11);
 float inB=step(0.0,q2.x)*step(q2.x,0.62*q2.y)*step(q2.y,1.0-0.25*q2.x);
 col=mix(col,vec3(0.0),inA);
 col=mix(col,vec3(1.0),inB);
 return col;
}
// gem board event j: x first column, y row, z width, w height of the match
vec4 gemEvent(float j) {
 float h1=hash11(j*0.6180339+0.37);
 float h2=hash11(j*1.3247+4.1);
 float h3=hash11(j*2.2134+7.7);
 bool hz=h1<0.55;
 float w=hz?3.0:1.0,hh=hz?1.0:3.0;
 return vec4(floor(h2*(9.0-w)),floor(h3*h3*(9.0-hh)),w,hh);
}
// one faceted gem, q in [-0.5,0.5]^2 of its cell
vec3 gem(vec2 q,float type,vec3 bg,float flash,float aa) {
 float n=type<0.5?4.0:type<1.5?6.0:type<2.5?4.0:type<3.5?8.0:type<4.5?12.0:type<5.5?3.0:10.0;
 float rt=type<1.5?0.0:type<2.5?PI*0.25:type<4.5?PI/8.0:type<5.5?-PI*0.5:0.0;
 vec3 gc=type<0.5?vec3(1.0,0.12,0.15):type<1.5?vec3(1.0,0.5,0.08):type<2.5?vec3(1.0,0.88,0.12)
   :type<3.5?vec3(0.15,0.95,0.3):type<4.5?vec3(0.15,0.45,1.0):type<5.5?vec3(0.85,0.25,1.0):vec3(0.92);
 vec2 p=q;
 if (type>1.5&&type<2.5) p.x/=0.78;
 if (type>4.5&&type<5.5) p.y+=0.06;
 float R=0.39*(1.0-0.5*flash);
 float sec=TAU/n;
 float an=atan(p.y,p.x)+rt;
 float k=floor(an/sec+0.5);
 float al=an-k*sec;
 float rin=R*cos(PI/n);
 float rr=length(p);
 float d=(type>3.5&&type<4.5)||type>5.5?rr-rin:rr*cos(al)-rin;
 float fa=k*sec-rt;
 float sh=0.55+0.45*dot(vec2(cos(fa),sin(fa)),vec2(-0.6,0.8));
 float table=step(rr*cos(al),rin*0.5);
 vec3 c=gc*mix(sh,1.08,table);
 c+=vec3(1.0)*0.8*smoothstep(0.06,0.0,length(p-vec2(-0.1,0.12)*R/0.39));
 c=mix(c,gc*0.35,smoothstep(-0.05,0.0,d));
 c=mix(c,vec3(1.0),flash*0.8);
 vec3 o=mix(bg,c,smoothstep(aa,-aa,d));
 return o+gc*flash*0.6*exp(-max(d,0.0)*12.0);
}
// Gem-matching game: a board of gems where matches of three flash and the
// column above drops in. Every twelfth event the board clears and refills.
vec3 gems(vec2 s,float t,float fp) {
 vec2 a=s*vec2(4.0,3.0);
 float nb=fbm2(a*1.2+vec2(t*0.03,0.0));
 vec3 col=mix(vec3(0.05,0.02,0.12),vec3(0.42,0.14,0.5),nb*nb*1.6);
 col+=vec3(0.15,0.25,0.6)*pow(fbm2(a*2.7-3.0),3.0);
 float EV=1.1;
 float uu=t/EV;
 float jn=floor(uu);
 float e=(uu-jn)*EV;
 float E=floor(jn/12.0);
 float j0=E*12.0;
 // left panel: score, the big level orb
 vec2 lp=a-vec2(0.62,1.5);
 float pnl=sdBox2(lp,vec2(0.48,1.3))-0.04;
 if (pnl<0.0) {
  col=mix(vec3(0.08,0.08,0.2),vec3(0.22,0.2,0.42),s.y);
  vec2 sb=lp-vec2(0.0,0.95);
  if (sdBox2(sb,vec2(0.4,0.13))<0.0) {
   col=vec3(0.02,0.01,0.03);
   float score=(jn-E*2.0)*30.0+150.0*E;
   float dc=floor((sb.x+0.4)/0.135);
   float dig=mod(floor(score/pow(10.0,5.0-dc)),10.0);
   float sd=seg7(vec2(mod(sb.x+0.4,0.135)-0.0675,sb.y)/vec2(0.135*0.8,0.1*0.9),segMask(dig));
   col=mix(col,vec3(1.0,0.7,0.2),smoothstep(0.08,0.0,sd)*step(dc,5.5));
  }
  float lvl=clamp((jn-j0+e/EV)/12.0,0.0,1.0);
  vec2 oq=lp-vec2(0.0,-0.65);
  float od=length(oq)-0.36;
  if (od<0.0) {
   vec3 oc=mix(vec3(0.25,0.05,0.4),vec3(0.85,0.3,1.0),step(oq.y/0.36*0.5+0.5,lvl));
   oc*=0.6+0.6*(1.0-length(oq-vec2(-0.1,0.12))/0.45);
   col=oc+vec3(1.0)*smoothstep(0.08,0.0,length(oq-vec2(-0.12,0.16)));
  }
  col=mix(col,vec3(0.7,0.75,0.95),smoothstep(0.025,0.0,abs(od)));
 }
 col=mix(col,vec3(0.75,0.78,0.95),smoothstep(0.03,0.0,abs(pnl)));
 vec2 bo=vec2(1.24,0.16);
 float cs=0.335;
 vec2 bp=(a-bo)/cs;
 float frame=sdBox2(bp-4.0,vec2(4.0));
 col=mix(col,vec3(0.7,0.72,0.9),smoothstep(0.1,0.0,abs(frame-0.05)));
 if (frame<0.0) {
  vec2 cc=floor(bp);
  vec3 bg=mix(vec3(0.07,0.05,0.16),vec3(0.13,0.11,0.26),mod(cc.x+cc.y,2.0))+col*0.15;
  float c=cc.x;
  float y=bp.y;
  float yoff=0.0,flash=0.0,kfloor=0.0;
  float jl=jn;
  if (jn-j0<0.5) {
   float f=clamp((e-c*0.05)/0.7,0.0,1.0);
   yoff=9.0*(1.0-f*f);
   jl=jn-1.0;
  } else {
   vec4 ev=gemEvent(jn);
   bool inCol=c>=ev.x&&c<ev.x+ev.z;
   if (e<0.35) {
    jl=jn-1.0;
    if (inCol&&y>=ev.y&&y<ev.y+ev.w) flash=e/0.35;
   } else if (inCol&&y>=ev.y) {
    float f=clamp((e-0.35)/0.35,0.0,1.0);
    yoff=ev.w*(1.0-f*f);
    kfloor=ev.y;
   }
  }
  float yy=y-yoff;
  float r=floor(yy);
  vec3 g=bg;
  if (yy>=0.0&&r>=kfloor) {
   float m=r;
   for (int i=0; i<12; i++) {
    float j=jl-float(i);
    if (j<j0+0.5) break;
    vec4 ev=gemEvent(j);
    if (c>=ev.x&&c<ev.x+ev.z&&m>=ev.y) m+=ev.w;
   }
   float type=floor(hash21(vec2(c*7.31+E*13.7,m*3.17+E*1.91))*7.0);
   g=gem(vec2(fract(bp.x),fract(yy))-0.5,type,bg,flash,max(fp*4.0/cs*1.5,0.02));
  }
  col=g;
 }
 return col;
}
// the starfield screensaver: layers of stars rushing out of the middle,
// each trailing a short streak back toward the centre
vec3 starfield(vec2 s,float t) {
 vec2 c=(s-0.5)*vec2(1.3333,1.0);
 vec3 col=vec3(0.0);
 for (int i=0; i<8; i++) {
  float fi=float(i);
  float z=fract((fi+1.0)/8.0-t*0.12);
  float zz=mix(0.03,1.0,z);
  vec2 off=vec2(fi*7.31,fi*3.17);
  vec2 p=c*zz*7.0+off;
  vec2 cell=floor(p);
  vec2 h=hash22(cell+fi*13.1);
  vec2 sp=cell+0.2+0.6*h;
  vec2 dv=(p-sp)/(zz*7.0);
  vec2 rdir=normalize((sp-off)/(zz*7.0)+1e-4);
  float along=dot(dv,rdir);
  float L=0.06*(1.0-z)*(1.0-z);
  float d=length(vec2(dot(dv,vec2(-rdir.y,rdir.x)),max(abs(along+0.5*L)-0.5*L,0.0)));
  float sz=0.0022+0.004*(1.0-z);
  float b=smoothstep(sz,0.0,d)*smoothstep(1.0,0.8,z);
  col+=b*(0.6+0.9*(1.0-z))*mix(vec3(0.8,0.88,1.0),vec3(1.0,0.95,0.85),h.y);
 }
 return col;
}
vec2 bounce(float t,float seed) {
 vec2 sp=vec2(0.16+0.12*hash11(seed+0.3),0.14+0.1*hash11(seed+0.7));
 return abs(fract(t*sp+hash22(vec2(seed,3.0)))*2.0-1.0)*vec2(1.3333,1.0);
}
float segD(vec2 p,vec2 a,vec2 b) {
 vec2 pa=p-a,ba=b-a;
 return length(pa-ba*clamp(dot(pa,ba)/dot(ba,ba),0.0,1.0));
}
// Mystify: two bouncing quadrilaterals, each dragging a trail of copies
vec3 mystify(vec2 s,float t,float fp) {
 vec2 a=s*vec2(1.3333,1.0);
 vec3 col=vec3(0.0);
 float w=max(0.0035,fp*1.2);
 for (int q=0; q<2; q++) {
  float fq=float(q);
  vec3 qc=0.55+0.45*cos(TAU*(t*0.04+fq*0.45+vec3(0.0,0.33,0.67)));
  for (int k=0; k<6; k++) {
   float tk=t-float(k)*0.11;
   vec2 v0=bounce(tk,fq*4.0+1.0),v1=bounce(tk,fq*4.0+2.0);
   vec2 v2=bounce(tk,fq*4.0+3.0),v3=bounce(tk,fq*4.0+4.0);
   float d=min(min(segD(a,v0,v1),segD(a,v1,v2)),min(segD(a,v2,v3),segD(a,v3,v0)));
   col=max(col,qc*(1.0-float(k)*0.15)*smoothstep(w,w*0.3,d));
  }
 }
 return col*1.3;
}
// The tube: barrel distortion, the selected program, scanlines that fade out
// once they are finer than a pixel, phosphor vignette.
vec3 crt(vec2 s0,float fp) {
 vec2 c=s0*2.0-1.0;
 c*=1.0+0.04*dot(c,c);
 vec2 s=c*0.5+0.5;
 float inside=smoothstep(0.0,0.006,s.x)*smoothstep(1.0,0.994,s.x)*smoothstep(0.0,0.008,s.y)*smoothstep(1.0,0.992,s.y);
 if (inside<=0.0) return vec3(0.0);
 vec3 col;
 if (gMode<0.5) col=desktop(s,gMT,fp);
 else if (gMode<1.5) col=gems(s,gMT,fp);
 else if (gMode<2.5) col=starfield(s,gMT);
 else col=mystify(s,gMT,fp);
 float lp=fp*300.0;
 col*=1.0-SCANLINES*0.45*(0.5+0.5*cos(s.y*300.0*TAU))*smoothstep(0.7,0.3,lp);
 // a soft roll bar
 col*=1.0-0.06*smoothstep(0.08,0.0,abs(fract(s.y-gMT*0.13)-0.5));
 col*=1.0-0.3*pow(dot(c,c)*0.5,2.0);
 return col*inside*gFade;
}
// The little TV on the shelf is playing an 8-bit platformer.
vec3 tvGame(vec2 s,float t) {
 vec2 px=floor(s*vec2(80.0,60.0));
 float x=px.x+floor(t*18.0);
 vec3 col=vec3(0.38,0.55,1.0);
 float cxl=mod(x*0.5,48.0);
 if (abs(px.y-46.0)<3.0&&abs(cxl-20.0)<7.0-abs(px.y-46.0)) col=vec3(1.0);
 float hh=16.0+floor(5.0*sin(x*0.045)+2.5*sin(x*0.13));
 if (px.y<hh) col=vec3(0.12,0.62,0.22)*(0.9+0.1*mod(px.x+px.y,2.0));
 if (px.y<9.0) {
  float by=mod(px.y,4.0);
  float bx=mod(x+4.0*mod(floor(px.y/4.0),2.0),8.0);
  col=(by<0.5||bx<0.5)?vec3(0.25,0.1,0.05):vec3(0.8,0.38,0.15);
 }
 float bk=mod(x,44.0);
 if (px.y>=24.0&&px.y<29.0&&bk>=12.0&&bk<27.0)
  col=mod(bk-12.0,5.0)<0.5||px.y<24.5?vec3(0.3,0.12,0.0):vec3(1.0,0.72,0.1);
 float jy=9.0+floor(max(0.0,sin(t*3.1))*14.0);
 if (px.x>=20.0&&px.x<24.0&&px.y>=jy&&px.y<jy+7.0)
  col=px.y-jy<3.0?vec3(0.2,0.3,0.9):px.y-jy<5.0?vec3(0.95,0.15,0.1):vec3(1.0,0.78,0.6);
 float sc=0.85+0.15*cos(s.y*120.0*TAU);
 vec2 c=s*2.0-1.0;
 return col*sc*(1.0-0.35*pow(dot(c,c)*0.5,2.0))*(0.92+0.08*hash21(px+floor(t*30.0)));
}
// ---- posters ----------------------------------------------------------------
float glyph(vec2 p,float seed) {
 if (p.x<0.1||p.x>0.9||p.y<0.1||p.y>0.9) return 0.0;
 vec2 g=floor((p-0.1)/0.8*vec2(3.0,5.0));
 float h=hash21(vec2(seed,g.x+g.y*3.0));
 return step(g.x==1.0?0.55:0.3,h);
}
vec3 poster(vec2 u,float id,float fp) {
 vec3 col;
 if (id<0.5) {
  // fantasy game poster: dusk sky, big moon, a castle, title in gold
  col=mix(vec3(0.75,0.22,0.12),vec3(0.08,0.05,0.22),smoothstep(0.1,0.85,u.y));
  float md=length((u-vec2(0.6,0.6))*vec2(1.0,1.5))-0.25;
  col=mix(col,vec3(1.0,0.92,0.75)*(0.85+0.15*noise2(u*18.0)),smoothstep(0.01,0.0,md));
  float cx=u.x;
  float h=0.18+0.03*sin(cx*40.0);
  h=max(h,0.34*step(abs(cx-0.62),0.05));
  h=max(h,0.28*step(abs(cx-0.8),0.04));
  h=max(h,0.26*step(abs(cx-0.45),0.12));
  h=max(h,step(abs(cx-0.32),0.06)*(0.42+(0.06-abs(cx-0.32))*2.2));
  h=max(h,step(abs(cx-0.62),0.05)*(0.34+(0.05-abs(cx-0.62))*2.2));
  h+=0.012*step(0.5,fract(cx*50.0))*step(h,0.3);
  col=mix(col,vec3(0.03,0.02,0.05),step(u.y,h));
  col=mix(col,vec3(1.0,0.8,0.2),step(0.3,hash21(floor(u*vec2(60.0,80.0))))*step(u.y,h)*step(0.995,hash21(floor(u*vec2(60.0,80.0))+3.0)));
  vec2 tp=vec2((u.x-0.12)/0.76*7.0,(u.y-0.85)/0.1);
  if (tp.x>0.0&&tp.x<7.0&&tp.y>0.0&&tp.y<1.0)
   col=mix(col,vec3(1.0,0.78,0.2),glyph(fract(tp),floor(tp.x)+10.0));
 } else if (id<1.5) {
  // band poster: halftone on hot pink, three silhouettes, title band
  col=vec3(0.95,0.2,0.45);
  vec2 hp=u*vec2(30.0,45.0);
  float dr=0.18+0.25*u.y;
  col=mix(col,vec3(0.55,0.05,0.3),smoothstep(dr,dr-0.08,length(fract(hp)-0.5))*(1.0-smoothstep(0.3,0.8,fp*45.0)));
  float fig=1e3;
  for (int i=0; i<3; i++) {
   float fi=float(i);
   vec2 c=vec2(0.22+0.28*fi,0.5+0.04*(1.0-abs(fi-1.0)));
   fig=min(fig,length((u-c)*vec2(1.0,0.75))-0.075);
   fig=min(fig,sdBox2(u-vec2(c.x,c.y-0.3),vec2(0.13,0.22))-0.04);
  }
  col=mix(col,vec3(0.04,0.02,0.05),step(fig,0.0));
  if (u.y>0.8) {
   col=vec3(0.05);
   vec2 tp=vec2((u.x-0.1)/0.8*5.0,(u.y-0.83)/0.13);
   if (tp.x>0.0&&tp.x<5.0) col=mix(col,vec3(1.0),glyph(fract(tp),floor(tp.x)+30.0));
  }
 } else if (id<2.5) {
  // space poster: ringed planet over a starfield
  col=mix(vec3(0.02,0.03,0.12),vec3(0.15,0.05,0.25),u.y);
  col+=vec3(1.0)*step(0.985,hash21(floor(u*vec2(120.0,90.0))))*0.8;
  vec2 pc=(u-vec2(0.55,0.52))*vec2(1.4,1.0);
  float pd=length(pc)-0.26;
  vec3 pcol=mix(vec3(0.95,0.6,0.25),vec3(0.6,0.3,0.15),0.5+0.5*sin(pc.y*40.0+sin(pc.x*9.0)));
  pcol*=0.5+0.7*clamp(dot(normalize(vec3(pc,sqrt(max(0.0,0.0676-dot(pc,pc))))),normalize(vec3(-0.6,0.5,0.6))),0.0,1.0);
  vec2 rc=rot(-0.35)*pc;
  float ring=abs(length(rc*vec2(1.0,3.6))-0.4);
  float front=step(0.0,-rc.y);
  col=mix(col,pcol,step(pd,0.0));
  col=mix(col,vec3(0.9,0.82,0.65),smoothstep(0.05,0.03,ring)*max(front,step(0.0,pd)));
  if (u.y<0.16) {
   vec2 tp=vec2((u.x-0.15)/0.7*6.0,(u.y-0.04)/0.09);
   if (tp.x>0.0&&tp.x<6.0) col=mix(col,vec3(0.6,0.9,1.0),glyph(fract(tp),floor(tp.x)+50.0));
  }
 } else {
  // racing poster: speed stripes and a low sports car
  col=mix(vec3(1.0,0.85,0.15),vec3(0.95,0.3,0.1),u.y);
  float st=fract((u.x+u.y*0.6)*5.0);
  col=mix(col,vec3(0.1,0.05,0.05),step(0.75,st)*step(u.y,0.75));
  vec2 cp=u-vec2(0.5,0.33);
  float car=sdBox2(cp,vec2(0.38,0.05))-0.03;
  car=min(car,sdBox2(cp-vec2(0.02,0.08),vec2(0.17,0.035))-0.04);
  col=mix(col,vec3(0.8,0.05,0.05),step(car,0.0));
  col=mix(col,vec3(0.4,0.6,0.8),step(sdBox2(cp-vec2(0.03,0.095),vec2(0.13,0.02))-0.02,0.0));
  float wh=min(length(cp-vec2(-0.24,-0.06)),length(cp-vec2(0.25,-0.06)))-0.07;
  col=mix(col,vec3(0.03),step(wh,0.0));
  col=mix(col,vec3(0.7),step(wh,-0.04));
  if (u.y>0.82) {
   vec2 tp=vec2((u.x-0.08)/0.84*6.0,(u.y-0.85)/0.11);
   if (tp.x>0.0&&tp.x<6.0) col=mix(col,vec3(0.05),glyph(fract(tp),floor(tp.x)+70.0));
  }
 }
 float b=min(min(u.x,1.0-u.x),min(u.y,1.0-u.y));
 col=mix(vec3(0.93,0.91,0.86),col,step(0.025,b));
 // pins / tape at the corners
 vec2 cn=abs(u-0.5);
 col=mix(col,vec3(0.85,0.85,0.75),step(length(cn-vec2(0.47,0.47)),0.025)*0.8);
 return col;
}
// ---- surface shading ---------------------------------------------------------
// coverage of a line of half-width w at distance d, box-filtered by the pixel
float ln(float d,float w,float fp) {
 float hw=max(w,fp*0.5);
 return smoothstep(hw+fp*0.5,hw-fp*0.5,abs(d))*min(1.0,w/hw);
}
vec3 spectral(float w) {
 // wavelength in nm -> rough rgb
 float x=(w-400.0)/300.0;
 vec3 c=vec3(max(1.0-pow((x-0.85)/0.32,2.0),0.0)+0.35*max(1.0-pow((x-0.02)/0.1,2.0),0.0),
    1.0-pow((x-0.5)/0.28,2.0),
    1.0-pow((x-0.15)/0.28,2.0));
 return clamp(c,0.0,1.0)*step(0.0,x)*step(x,1.0);
}
vec3 shade(vec3 pw,vec3 n,vec3 rd,float m,float sh,float ao,float t,float pix) {
 vec3 p=pw;
 if (m<1.5||(m>3.5&&m<11.5)||m>16.5) p.x*=gS;
 vec3 alb=vec3(0.5);
 float spec=0.1,gloss=30.0;
 vec3 emit=vec3(0.0);
 bool slatsOn=true;
 float fp=pix*t;
 if (m<1.5) {
  // laminate, grain along the panel
  vec2 gp=abs(n.y)>0.5?p.xz:abs(n.z)>0.5?p.xy:p.yz;
  float g=noise2(vec2(gp.x*3.0,gp.y*70.0)+3.0*fbm2(gp*vec2(2.0,6.0)));
  alb=mix(ramp3(vec3(0.7,0.58,0.42),vec3(0.55,0.38,0.22),vec3(0.2,0.11,0.06),WOOD),
    ramp3(vec3(0.88,0.78,0.6),vec3(0.74,0.56,0.36),vec3(0.34,0.2,0.12),WOOD),g);
  alb*=0.9+0.1*fbm2(gp*20.0);
  if (p.z>0.93) alb*=0.7;
  spec=0.35; gloss=50.0;
 } else if (m<2.5) {
  alb=hueRot(vec3(0.58,0.68,0.72),WALL_HUE*TAU);
  alb*=0.95+0.05*noise2(p.xy*60.0+p.z*40.0);
  if (n.z<-0.5) {
   // posters on the back wall
   vec2 w=p.xy;
   vec4 P0=vec4(-1.2,1.40,0.28,0.42);
   vec4 P1=vec4(1.03,1.30,0.17,0.28);
   vec4 P2=vec4(-0.25,2.16,0.37,0.27);
   vec4 P3=vec4(-1.92,1.58,0.29,0.37);
   vec4 pr=w.x>0.5?P1:w.y>1.84&&w.x>-0.85?P2:w.x>-1.6?P0:P3;
   float pid=mod((w.x>0.5?1.0:w.y>1.84&&w.x>-0.85?2.0:w.x>-1.6?0.0:3.0)+gSeed,4.0);
   vec2 u=(w-pr.xy)/(2.0*pr.zw)+0.5;
   if (u.x>0.0&&u.x<1.0&&u.y>0.0&&u.y<1.0)
    alb=poster(u,pid,fp/(2.0*pr.w))*0.85;
  }
  if (p.y<0.09) alb=vec3(0.86,0.85,0.82);
  spec=0.05;
 } else if (m<3.5) {
  alb=vec3(0.40,0.37,0.34)*(0.8+0.4*noise2(p.xz*300.0)*noise2(p.xz*37.0+5.0));
  spec=0.0;
 } else if (m<4.5) {
  // beige plastic, and the details moulded or printed on it
  alb=ramp3(vec3(0.86,0.86,0.84),vec3(0.80,0.76,0.66),vec3(0.78,0.65,0.42),YELLOWING)*(0.96+0.04*noise2(p.xy*200.0));
  spec=0.25; gloss=25.0;
  if (n.z<-0.6) {
   vec2 pn=rot(0.12)*vec2(p.x-MX-0.168,p.y-1.152);
   if (CLUTTER>0.55&&abs(p.x-MX)<0.21&&p.z<0.5&&max(abs(pn.x),abs(pn.y))<0.036) {
    // a yellow sticky note with a scrawled reminder
    alb=vec3(1.0,0.86,0.3)*(1.0-0.15*step(pn.y,-0.03));
    float lr=fract((0.02-pn.y)*90.0);
    alb=mix(alb,vec3(0.15,0.2,0.55),0.8*ptext(vec2((pn.x+0.03)*220.0,lr*1.7-0.35),floor((0.02-pn.y)*90.0)+60.0,clamp(fp*300.0-0.5,0.0,1.0))*step(abs(pn.x),0.028)*step(pn.y,0.024)*step(-0.028,pn.y));
    spec=0.05;
   } else if (abs(p.x-MX)<0.21&&p.z<0.5&&p.y<0.87) {
    // monitor chin: badge, buttons, power LED
    vec2 c=vec2(p.x-MX,p.y-0.84);
    alb*=1.0-0.12*step(abs(c.y-0.006),0.0015)*step(abs(c.x),0.19);
    float bt=sdBox2(vec2(mod(c.x+0.2,0.025)-0.0125,c.y+0.012),vec2(0.008,0.004));
    alb=mix(alb,alb*0.75,step(bt,0.0)*step(0.06,c.x));
    alb=mix(alb,vec3(0.15),step(sdBox2(c-vec2(0.0,0.022),vec2(0.035,0.007)),0.0));
    float led=length(c-vec2(0.165,-0.012));
    emit+=vec3(0.2,1.0,0.2)*2.0*smoothstep(0.004,0.0015,led);
   } else if (p.x>0.48&&p.z<0.47) {
    // tower front: CD-ROM, blank bay, floppy, power, LEDs, vents, badge
    vec2 c=vec2(p.x-0.58,p.y);
    float cd=sdBox2(c-vec2(0.0,1.105),vec2(0.072,0.019));
    alb=mix(alb,vec3(0.70,0.68,0.61),step(cd,0.0));
    alb*=1.0-0.55*ln(cd,0.0012,fp);
    alb*=1.0-0.6*ln(c.y-1.099,0.0007,fp)*step(abs(c.x),0.066);
    alb=mix(alb,vec3(0.45,0.44,0.4),step(sdBox2(c-vec2(0.052,1.092),vec2(0.008,0.0035)),0.0));
    emit+=vec3(1.0,0.5,0.1)*1.5*step(length(c-vec2(-0.058,1.092)),0.0018)*step(0.5,fract(gMT*0.23));
    float bay=sdBox2(c-vec2(0.0,1.058),vec2(0.072,0.019));
    alb*=1.0-0.45*ln(bay,0.0011,fp);
    float fl=sdBox2(c-vec2(0.0,1.01),vec2(0.05,0.01));
    alb=mix(alb,vec3(0.70,0.68,0.61),step(fl,0.0));
    alb*=1.0-0.55*ln(fl,0.0011,fp);
    alb=mix(alb,vec3(0.08),ln(c.y-1.013,0.0016,fp)*step(abs(c.x+0.006),0.036));
    alb=mix(alb,vec3(0.5,0.48,0.44),step(sdBox2(c-vec2(0.036,1.004),vec2(0.007,0.0028)),0.0));
    float pb=length(c-vec2(0.0,0.87))-0.016;
    alb=mix(alb,alb*0.82,step(pb,0.0));
    alb*=1.0-0.45*ln(pb,0.0012,fp);
    emit+=vec3(0.2,1.0,0.25)*2.0*smoothstep(0.0035,0.0018,length(c-vec2(-0.03,0.83)));
    float hdd=step(0.55,hash11(floor(gMT*14.0)))*step(0.3,fract(gMT*0.37));
    emit+=vec3(1.0,0.6,0.1)*2.5*hdd*smoothstep(0.0035,0.0018,length(c-vec2(0.03,0.83)));
    // badge
    float bg=sdBox2(c-vec2(0.0,0.94),vec2(0.03,0.009));
    alb=mix(alb,vec3(0.55,0.56,0.6),step(bg,0.0));
    alb=mix(alb,vec3(0.1,0.25,0.6),step(bg,0.0)*step(abs(c.y-0.94),0.003));
    // bezel seam and vents
    alb*=1.0-0.4*ln(sdBox2(c-vec2(0.0,0.96),vec2(0.085,0.2)),0.0008,fp);
    alb*=1.0-0.45*step(c.y,0.8)*step(0.765,c.y)*mix(step(0.55,fract(c.y*160.0)),0.45,clamp(fp*320.0-0.5,0.0,1.0))*step(abs(c.x),0.06);
   } else if (p.y>0.75&&p.y<0.985&&p.z<0.535&&(abs(p.x+0.47)<0.06||abs(p.x-0.21)<0.06)) {
    // speaker grilles and the little cone
    vec2 c=vec2(p.x-(p.x<-0.13?-0.47:0.21),p.y-0.865);
    float gr=sdBox2(c-vec2(0.0,0.01),vec2(0.042,0.085))-0.004;
    vec2 gd=fract(c*260.0)-0.5;
    float holes=mix(step(length(gd),0.3),0.3,clamp(fp*260.0-0.3,0.0,1.0));
    alb=mix(alb,mix(vec3(0.22,0.21,0.2),alb*0.6,holes),step(gr,0.0));
    alb=mix(alb,vec3(0.1),step(abs(length(c-vec2(0.0,0.03))-0.03),0.002));
   }
  }
  if (p.y>0.754&&p.y<0.79&&length(p.xz-vec2(0.31,0.30))<0.07) {
   // mouse buttons split
   vec3 mo=p-vec3(0.31,0.753,0.30);
   mo.xz=rot(0.25)*mo.xz;
   alb*=1.0-0.4*step(abs(mo.x),0.0008)*step(0.01,-mo.z);
   alb*=1.0-0.4*step(abs(mo.z+0.012),0.0008);
  }
 } else if (m<5.5) {
  // CRT face: glass over the phosphor
  vec2 su=vec2(gS*(p.x-MX)/0.32+0.5,(p.y-SCY)/0.24+0.5);
  emit=crt(su,fp/0.24)*1.35*GLOW;
  alb=vec3(0.02,0.025,0.022);
  vec3 r=reflect(rd,n);
  float fr=0.04+0.5*pow(1.0-max(dot(-rd,n),0.0),5.0);
  emit+=windowSeen(pw,r)*fr*0.6+vec3(0.05,0.05,0.07)*fr;
  spec=0.6; gloss=120.0;
  slatsOn=true;
 } else if (m<6.5) {
  float jc=floor((p.x+0.77)/0.0104);
  float h=hash11(jc*1.71+0.3+gSeed*0.37);
  vec3 sp=h<0.35?vec3(0.06):0.5+0.5*cos(TAU*(h*2.3+vec3(0.0,0.33,0.67)));
  alb=n.z<-0.5||(p.y<0.8&&n.y>0.5)?sp:vec3(0.15,0.15,0.17);
  if (n.z<-0.5) alb=mix(alb,vec3(0.95),0.6*step(abs(fract((p.x+0.77)/0.0104)-0.5),0.12)*step(abs(p.y-1.37),0.035));
  spec=0.5; gloss=80.0;
 } else if (m<7.5) {
  alb=vec3(0.06,0.06,0.065);
  spec=0.3; gloss=40.0;
  if (p.x>0.33&&p.x<0.39&&p.z<0.67&&p.y>1.36&&p.y<1.56) {
   // TV knobs
   float kb=min(length(vec2(p.x-0.36,p.y-1.52)),length(vec2(p.x-0.36,p.y-1.47)))-0.012;
   alb=mix(alb,vec3(0.3),step(kb,0.0));
   alb=mix(alb,vec3(0.6,0.6,0.55),step(length(vec2(p.x-0.36,p.y-1.40)),0.006));
  }
 } else if (m<8.5) {
  // no-name cola: red can, silver band, a yellow roundel facing the room
  bool c2=length(p.xz-gCan2)<length(p.xz-gCan);
  vec3 cq=p-vec3(c2?gCan2.x:gCan.x,0.8115,c2?gCan2.y:gCan.y);
  float an=atan(cq.z,cq.x);
  alb=c2?vec3(0.08,0.22,0.6):vec3(0.72,0.05,0.06);
  alb=mix(alb,vec3(0.82,0.82,0.86),step(abs(cq.y+0.036),0.007));
  vec2 rq=vec2((an+2.3)*0.033,cq.y-0.012);
  float rd0=length(rq)-0.018;
  alb=mix(alb,vec3(1.0,0.8,0.15),step(rd0,0.0));
  alb=mix(alb,vec3(0.95),step(abs(rd0+0.004),0.0015));
  alb=mix(alb,vec3(0.95),ptext(vec2(an*9.0,(cq.y+0.022)/0.012),21.0,0.0)*step(abs(an+2.3),0.9));
  if (cq.y>0.055||cq.y<-0.058) alb=vec3(0.6,0.6,0.62);
  spec=0.8; gloss=60.0;
 } else if (m<9.5) {
  vec3 uq=p-vec3(gMug.x,0.797,gMug.y);
  alb=vec3(0.92,0.9,0.86);
  if (uq.y<0.033&&length(uq.xz)<0.0355) alb=vec3(0.12,0.06,0.03);
  // "I (heart) my computer"
  float an=atan(uq.x,-uq.z);
  vec2 hu=vec2(an*0.04,uq.y-0.004);
  if (length(uq.xz)>0.038) {
   vec2 hq=(hu-vec2(0.0,-0.011))/0.022;
   hq.x=abs(hq.x);
   float hd=hq.y+hq.x>1.0?length(hq-vec2(0.25,0.75))-0.35355
     :sqrt(min(dot(hq-vec2(0.0,1.0),hq-vec2(0.0,1.0)),
        dot(hq-0.5*max(hq.x+hq.y,0.0),hq-0.5*max(hq.x+hq.y,0.0))))*sign(hq.x-hq.y);
   alb=mix(alb,vec3(0.85,0.05,0.12),step(hd,0.0));
   alb=mix(alb,vec3(0.1),ptext(vec2((hu.x+0.03)*300.0,(hu.y-0.017)/0.008),81.0,0.0)*step(abs(hu.x),0.03)*step(0.017,hu.y));
   alb=mix(alb,vec3(0.1),ptext(vec2((hu.x+0.035)*300.0,(hu.y+0.024)/0.008),82.0,0.0)*step(abs(hu.x),0.035)*step(hu.y,-0.016));
  }
  spec=0.5; gloss=60.0;
 } else if (m<10.5) {
  // lava lamp: wax blobs rising and sinking in red liquid, lit from below
  vec2 lq=vec2(p.x-LX,p.y-0.87);
  float f=0.0;
  for (int i=0; i<5; i++) {
   float fi=float(i);
   float ph=gMT*0.0+TIME*(0.035+0.025*hash11(fi+0.5))+hash11(fi*3.3+1.0);
   vec2 bc=vec2(0.012*sin(TIME*0.3+fi*2.0),0.02+0.19*(0.5-0.5*cos(ph*TAU)));
   float br=0.011+0.008*hash11(fi*7.7+2.0);
   vec2 dd=lq-bc;
   f+=br*br/(dot(dd,dd)+1e-5);
  }
  f+=0.0004/(lq.y*lq.y+0.0001);
  float blob=smoothstep(0.8,1.3,f);
  emit=mix(vec3(0.75,0.1,0.04),vec3(1.0,0.6,0.12)*1.7,blob)*(0.6+0.4*smoothstep(0.25,0.0,lq.y))*LAVA*1.4;
  alb=vec3(0.05);
  spec=0.8; gloss=90.0;
 } else if (m<11.5) {
  alb=vec3(0.28,0.26,0.24);
  spec=1.0; gloss=40.0;
  emit+=vec3(1.0,0.4,0.1)*0.25*LAVA*smoothstep(0.75+0.14,0.75+0.115,p.y)*smoothstep(0.75+0.08,0.75+0.12,p.y);
 } else if (m<12.5) {
  // blinds: lit through from outside
  alb=vec3(0.92,0.88,0.8);
  emit=gSun*0.12*sunMask(p-vec3(0.03,0.0,0.0),false)*(0.5+0.5*max(n.y,0.0));
  slatsOn=false;
  spec=0.1;
 } else if (m<13.5) {
  alb=vec3(0.88,0.87,0.84);
  spec=0.15;
  if (p.x>1.19&&abs(p.y-2.02)<0.023) slatsOn=false;
 } else if (m<14.5) {
  // plaid blanket
  vec2 bp=abs(n.y)>0.5?p.xz:abs(n.x)>0.5?p.zy:p.xy;
  vec2 f=fract(bp/0.14);
  float a1=step(0.5,f.x),a2=step(0.5,f.y);
  alb=mix(vec3(0.84,0.78,0.62),vec3(0.18,0.22,0.36),(a1+a2)*0.5);
  alb=mix(alb,vec3(0.55,0.12,0.1),0.6*max(step(abs(f.x-0.75),0.03),step(abs(f.y-0.75),0.03)));
  alb*=0.85+0.15*noise2(bp*300.0);
  spec=0.0;
 } else if (m<15.5) {
  alb=vec3(0.86,0.84,0.78); spec=0.05;
 } else if (m<16.5) {
  alb=vec3(0.05,0.05,0.055); spec=0.6; gloss=30.0;
 } else if (m<17.5) {
  // big-box game spines
  float gc=floor((p.x+0.255)/0.05);
  float h=hash11(gc*2.9+4.0+gSeed);
  alb=0.5+0.5*cos(TAU*(h+vec3(0.0,0.33,0.67)));
  alb*=h<0.3?0.25:0.9;
  if (n.z<-0.5) {
   float lx=fract((p.x+0.255)/0.05)-0.5;
   alb=mix(alb,vec3(1.0,0.95,0.85),ptext(vec2((p.y-1.33)*70.0,lx*2.4+0.5),gc+90.0,0.0)*step(abs(lx),0.3));
   alb=mix(alb,vec3(0.95,0.7,0.1),step(length(vec2(lx*0.05,p.y-1.5)),0.012));
  }
  if (pw.y<0.6) {
   // the magazine on the bed
   alb=mix(vec3(0.9,0.15,0.1),vec3(0.15,0.2,0.5),step(0.5,fract(pw.x*4.0+pw.z*3.0)));
   alb=mix(alb,vec3(1.0),ptext(vec2(pw.x*40.0,fract(pw.z*14.0)*1.4),floor(pw.z*14.0),0.0)*0.8);
  }
  spec=0.3;
 } else if (m<18.5) {
  alb=vec3(0.1,0.22,0.55);
  vec2 pq=p.xz-vec2(0.30,0.31);
  alb=mix(alb,vec3(0.9,0.85,0.3),step(abs(length(pq-vec2(-0.05,-0.04))-0.035),0.004));
  spec=0.05;
 } else if (m<19.5) {
  // floppies: black on top, blue under, label and shutter
  bool top=p.y>0.7535;
  vec3 fq=top?p-vec3(gFlp.x-0.01,0.7551,gFlp.y+0.01):p-vec3(gFlp.x,0.7517,gFlp.y-0.005);
  fq.xz=rot(top?-0.2:0.35)*fq.xz;
  alb=top?vec3(0.04):vec3(0.08,0.15,0.55);
  if (sdBox2(fq.xz-vec2(0.0,0.012),vec2(0.035,0.03))<0.0) alb=vec3(0.9,0.88,0.82);
  if (sdBox2(fq.xz-vec2(0.008,-0.034),vec2(0.023,0.014))<0.0) alb=vec3(0.65,0.65,0.68);
  spec=0.3;
 } else if (m<20.5) {
  // keycaps: light alphas, darker modifiers
  vec3 kq=p-vec3(-0.10,0.762,0.315);
  kq.y-=0.08*kq.z;
  kq.x*=gS;
  float pt=0.0185;
  vec2 c=vec2(floor((kq.x+0.218)/pt),clamp(floor((kq.z+0.066)/pt),0.0,5.0));
  bool mod1=kq.x>0.062?(c.x>20.5||(c.x<18.5)):(c.x<0.5||c.x>13.5||c.y>4.5||(c.y<0.5&&(c.x<3.5||c.x>10.5)));
  alb=(mod1?vec3(0.63,0.61,0.56):vec3(0.84,0.81,0.73))*ramp3(vec3(1.07,1.13,1.27),vec3(1.0),vec3(0.97,0.86,0.64),YELLOWING);
  // legends, only when they can be seen
  vec2 kl=vec2(fract((kq.x+0.218)/pt),fract((kq.z+0.066)/pt));
  float lg=glyph(vec2(kl.x*3.0-0.45,kl.y*3.0-1.6),c.x+c.y*17.0)*step(0.02,kq.y);
  alb=mix(alb,vec3(0.3,0.28,0.26),lg*0.55*smoothstep(0.0022,0.0011,fp));
  spec=0.3; gloss=30.0;
 } else if (m<21.5) {
  vec2 su=vec2(gS*(p.x-0.215)/0.2+0.5,(p.y-1.47)/0.15+0.5);
  emit=tvGame(su,TIME)*1.1*GLOW;
  alb=vec3(0.02);
  float fr=0.04+0.5*pow(1.0-max(dot(-rd,n),0.0),5.0);
  emit+=windowSeen(pw,reflect(rd,n))*fr*0.5;
  spec=0.6; gloss=120.0;
 } else if (m<22.5) {
  // CD: a dark mirror with a rainbow streak along the diameter that
  // lines up with the light (the tracks are a radial grating)
  vec3 dq=p-vec3(gCD.x,0.7508,gCD.y);
  float rl=length(dq.xz);
  vec2 rr=dq.xz/max(rl,1e-4);
  alb=rl<0.02?vec3(0.55,0.57,0.6):vec3(0.16,0.17,0.19);
  vec3 v=-rd;
  float tgs=dot(normalize(gL.xz+v.xz),vec2(-rr.y,rr.x));
  float tgm=dot(normalize(normalize(vec3(MX,SCY,0.49)-p).xz+v.xz),vec2(-rr.y,rr.x));
  vec3 irid=spectral(400.0+300.0*fract(rl*14.0+0.3*dot(v.xz,rr)));
  float on=step(0.022,rl);
  emit+=irid*on*(gSun*0.2*exp(-tgs*tgs*10.0)*sunMask(pw,true)*sh+gScr*0.5*exp(-tgm*tgm*10.0)+0.04);
  emit+=windowSeen(pw,reflect(rd,n))*0.3*on;
  spec=1.0; gloss=80.0;
 } else {
  // VHS sleeves: black, some white, label stripe
  float vc=floor((p.y-1.312)/0.026);
  float h=hash11(vc*5.1+9.0+gSeed);
  alb=h<0.6?vec3(0.04):h<0.8?vec3(0.85):vec3(0.6,0.15,0.1);
  if (n.z<-0.5) alb=mix(alb,vec3(0.92,0.9,0.8),step(abs(p.x-0.62),0.05)*step(abs(fract((p.y-1.312)/0.026)-0.5),0.25));
  spec=0.4;
 }
 // ---- lighting
 float sm=sunMask(pw,slatsOn)*sh;
 float dif=max(dot(n,gL),0.0);
 vec3 col=alb*gSun*dif*sm;
 // dusk sky fill and warm bounce off the sunlit patches
 vec3 amb=mix(vec3(0.06,0.06,0.1),vec3(0.11,0.11,0.17),0.5+0.5*n.y);
 amb+=gSun*0.04*(0.6+0.4*max(-n.x,0.0))*(0.7+0.3*n.y);
 col+=alb*amb*ao;
 // the monitor screen as a rectangle light
 vec3 sv=vec3(MX*gS,SCY,0.49)-pw;
 float sd2=dot(sv,sv);
 vec3 sl=sv*inversesqrt(sd2);
 float se=max(dot(n,sl),0.0)*max(sl.z,0.0)*0.077/(sd2*PI+0.02);
 col+=alb*gScr*se*1.4*GLOW*mix(1.0,ao,0.5);
 col+=gScr*spec*pow(max(dot(n,normalize(sl-rd)),0.0),gloss)*se*2.0*GLOW;
 // the little TV
 vec3 tvv=vec3(0.215*gS,1.47,0.66)-pw;
 float td2=dot(tvv,tvv);
 vec3 tl=tvv*inversesqrt(td2);
 col+=alb*vec3(0.45,0.55,0.85)*max(dot(n,tl),0.0)*max(tl.z,0.0)*0.03/(td2*PI+0.01)*GLOW;
 // lava lamp glow
 vec3 lv=vec3(LX*gS,0.97,LZ)-pw;
 float ld2=dot(lv,lv);
 float lw=max(dot(n,lv*inversesqrt(ld2))*0.8+0.2,0.0);
 col+=alb*vec3(1.0,0.35,0.1)*LAVA*0.012*lw/(ld2+0.004)*(0.9+0.1*sin(TIME*2.3));
 // sun specular
 vec3 h=normalize(gL-rd);
 col+=gSun*sm*spec*pow(max(dot(n,h),0.0),gloss)*dif;
 return col+emit;
}
vec3 aces(vec3 x) {
 return clamp((x*(2.51*x+0.03))/(x*(2.43*x+0.59)+0.14),0.0,1.0);
}
vec3 tapDir(float k) {
 return k<0.5?vec3(1.0,-1.0,-1.0):k<1.5?vec3(-1.0,-1.0,1.0):k<2.5?vec3(-1.0,1.0,-1.0):vec3(1.0);
}
void main() {
 vec2 frag=gl_FragCoord.xy;
 vec2 uv=(2.0*frag-RENDERSIZE)/RENDERSIZE.y;
 float T=TIME*SPEED;
 // sun, blinds
 gL=normalize(vec3(1.0,mix(0.12,0.85,SUN_HEIGHT),-0.55));
 gSun=mix(vec3(1.0,0.86,0.68),vec3(1.0,0.45,0.16),GOLDEN)*4.2*SUN;
 float sa=mix(0.35,-0.9,BLINDS);
 gCa=cos(sa); gSa=sin(sa);
 // the arrangement: mug, floppies, CD and can dealt to six spots with an
 // affine step (a coprime to 6, so they never share one); ARRANGE 0 is the
 // original layout. CLUTTER decides which of them are out at all.
 gS=MIRROR?-1.0:1.0;
 gSeed=float(ARRANGE);
 float ka=mod(gSeed,2.0)<0.5?1.0:5.0,kb=floor(gSeed/2.0);
 gMug=slot(kb); gFlp=slot(ka+kb); gCD=slot(2.0*ka+kb); gCan=slot(3.0*ka+kb);
 gJc=slot(4.0*ka+kb); gCan2=slot(5.0*ka+kb);
 gOn=step(vec4(0.15,0.3,0.45,0.05),vec4(CLUTTER));
 // what the monitor is showing
 float sel=float(SCREEN);
 gFade=1.0;
 if (sel<0.5) {
  float cyc=mod(TIME+10.0,56.0);
  gMode=floor(cyc/14.0);
  gMT=cyc-gMode*14.0;
  gFade=smoothstep(0.0,0.35,gMT)*smoothstep(14.0,13.75,gMT);
 } else {
  gMode=sel-1.0;
  gMT=TIME;
 }
 gScr=(gMode<0.5?vec3(0.4,0.55,0.8):gMode<1.5?vec3(0.32,0.18,0.45)
  :gMode<2.5?vec3(0.06,0.06,0.08):vec3(0.12,0.1,0.14))*gFade;
 // camera: a slow push in toward the desk and back out, drifting a little
 float k=0.5-0.5*cos(T*0.06);
 vec3 ro=vec3(mix(-0.32,-0.05,k)+0.06*sin(T*0.13),
    mix(1.33,1.18,k)+0.015*sin(T*0.21),
    mix(-1.2,-0.42,k));
 vec3 ta=vec3(mix(0.08,-0.06,k)+0.03*sin(T*0.09),mix(1.1,1.0,k),0.65);
 vec3 fw=normalize(ta-ro);
 vec3 rt=normalize(cross(vec3(0.0,1.0,0.0),fw));
 vec3 up=cross(fw,rt);
 float foc=1.75;
 vec3 rd=normalize(uv.x*rt+uv.y*up+foc*fw);
 float pix=2.0/(foc*RENDERSIZE.y);
 // One loop, one map() call site: march, then 4 normal taps, then 4 AO
 // taps, then the sun shadow ray.
 float t=0.02,m=-1.0,ph=0.0,kk=0.0;
 vec3 hp=ro,nacc=vec3(0.0),n=vec3(0.0,1.0,0.0);
 float ao=1.0,aoAcc=0.0,sh=1.0,st=0.0,stMax=0.0,hh0=0.0;
 bool outside=false;
 for (int i=0; i<240; i++) {
  vec3 sp;
  if (ph<0.5) sp=ro+rd*t;
  else if (ph<1.5) sp=hp+tapDir(kk)*0.0004*max(t,0.5);
  else if (ph<2.5) sp=hp+n*(0.01+0.022*kk);
  else sp=hp+n*0.002+gL*st;
  vec2 h=map(sp);
  if (ph<0.5) {
   if (h.x<0.0006*t) { m=h.y; hp=sp; hh0=h.x; ph=1.0; kk=0.0; }
   else {
    t+=h.x*0.95;
    if (sp.x>1.42) { outside=true; break; }
    if (t>8.0) break;
   }
  } else if (ph<1.5) {
   nacc+=tapDir(kk)*h.x;
   kk+=1.0;
   if (kk>3.5) { n=normalize(nacc); hp-=n*hh0; ph=2.0; kk=0.0; }
  } else if (ph<2.5) {
   aoAcc+=(0.01+0.022*kk-gR)*pow(0.75,kk);
   kk+=1.0;
   if (kk>3.5) {
    ao=clamp(1.0-4.0*aoAcc,0.0,1.0);
    // stop well short of the window: the trim and slats are
    // sunMask()'s job, and at a metre and a half even a 3 cm
    // clearance reads as deep penumbra to k*h/t
    if (dot(n,gL)<=0.0||sunMask(hp,false)<0.001||hp.x>1.0) break;
    ph=3.0; st=0.015; stMax=(1.0-hp.x)/gL.x;
   }
  } else {
   sh=min(sh,14.0*gR/st);
   st+=clamp(h.x,0.003,0.08);
   if (sh<0.005||st>stMax) break;
  }
 }
 sh=clamp(sh,0.0,1.0);
 vec3 col;
 float tEnd=t;
 if (m>0.0) col=shade(hp,n,rd,m,sh,ao,t,pix);
 else if (outside) col=skyOut(rd);
 else col=vec3(0.02);
 // sunbeams through the blinds: in-scattering along the view ray
 float cth=dot(rd,-gL);
 float bl=0.0,tb=min(tEnd,6.0);
 for (int i=0; i<4; i++) bl+=beamLen(ro,rd,tb,0.06*(float(i)-1.5));
 col+=gSun*bl*0.02*BEAMS*(0.5+0.8*pow(max(cth,0.0),3.0));
 // dust motes catching the light, on a few planes in front of the desk
 if (DUST>0.0) {
  for (int l=0; l<4; l++) {
   float fl=float(l);
   float zl=0.42-0.35*fl;
   if (rd.z<=0.0) break;
   float s=(zl-ro.z)/rd.z;
   if (s<0.05||s>tEnd) continue;
   vec3 wp=ro+rd*s;
   float smk=sunMask(wp,true);
   if (smk<0.01) continue;
   vec2 g=wp.xy*45.0+vec2(fl*13.7,0.0);
   g+=vec2(0.6*sin(TIME*0.05+fl),-TIME*0.12*(0.6+0.2*fl));
   vec2 cell=floor(g);
   vec2 hsh=hash22(cell+fl*31.0);
   if (hsh.x>0.09) continue;
   vec2 off=0.25+0.5*hash22(cell+7.7)+0.15*vec2(sin(TIME*0.7+hsh.y*20.0),cos(TIME*0.6+hsh.x*17.0));
   float dd=length(fract(g)-off);
   float px=max(s*pix*45.0,0.02);
   float tw=0.4+0.6*pow(0.5+0.5*sin(TIME*(2.0+3.0*hsh.y)+hsh.x*40.0),3.0);
   col+=gSun*0.22*DUST*smk*tw*exp(-dd*dd/(px*px))*min(1.0,0.025/px);
  }
 }
 // glow round the screen and the window (a cheap bloom)
 vec3 sc=vec3(MX*gS,SCY,0.49)-ro;
 float tc=dot(sc,rd);
 float dp=length(sc-rd*tc);
 col+=gScr*GLOW*0.05*exp(-dp*dp/0.012)/(1.0+tc*0.5);
 vec3 wc=vec3(BLX,1.5,-0.05)-ro;
 float tw2=dot(wc,rd);
 vec3 wd=wc-rd*tw2;
 col+=gSun*0.02*exp(-(wd.y*wd.y/0.3+wd.z*wd.z/0.35))*step(0.0,tw2);
 // ---- grade
 col*=EXPOSURE;
 col=aces(col*0.9);
 col=pow(col,vec3(0.95,1.0,1.06));
 // drugstore print: lifted blacks, warm cast, softer saturation
 vec3 fc=mix(vec3(dot(col,vec3(0.3,0.59,0.11))),col,0.82);
 fc=fc*vec3(1.0,0.96,0.86)*0.9+vec3(0.06,0.05,0.06);
 col=mix(col,fc,FILM);
 col+=(hash21(frag+fract(TIME*7.0)*113.0)-0.5)*(0.02+0.05*FILM);
 vec2 q=frag/RENDERSIZE;
 col*=0.5+0.5*pow(16.0*q.x*q.y*(1.0-q.x)*(1.0-q.y),0.22+0.1*FILM);
 // the orange date stamp in the corner
 if (DATE_STAMP) {
  float hgt=RENDERSIZE.y*0.028;
  vec2 dp2=(frag-vec2(RENDERSIZE.x-hgt*15.0,hgt*1.6))/hgt;
  // cells are 1.7 wide: a digit plus its italic slant spans about 1.4,
  // and the glow needs room to fade before the cell edge
  if (dp2.x>-0.4&&dp2.x<13.6&&dp2.y>-1.4&&dp2.y<1.4) {
   float cell=floor(dp2.x/1.7);
   vec2 lp=vec2(mod(dp2.x,1.7)-0.85,dp2.y);
   float dg=cell<-0.5?-1.0:cell<0.5?9.0:cell<1.5?8.0:cell<3.5?(cell<2.5?-1.0:1.0)
     :cell<4.5?0.0:cell<5.5?-1.0:cell<6.5?0.0:cell<7.5?1.0:-1.0;
   float d=dg<0.0?1e3:seg7(lp,segMask(dg));
   float tick=length((dp2-vec2(0.02,1.0))*vec2(2.0,1.0))-0.12;
   d=min(d,tick);
   vec3 ink=vec3(1.0,0.5,0.12);
   col+=ink*(smoothstep(0.08,0.0,d)*0.9+0.25*exp(-max(d,0.0)*14.0));
  }
 }
 gl_FragColor=vec4(clamp(col,0.0,1.0),1.0);
}