Elements

Vasilis Bellos
/
Elements
on 8 Nov 2026
4
39
12
1
9979
game.m
function g=game
% Elements
% A one-button material sandbox game
%
% Interactions
%  Heat and cooling
%   Contact                  Result
%   -----------------------  ---------------------------------------------
%   Oil + fire or lava       Oil ignites
%   Wood + fire or lava      Wood ignites; smoke can form above it
%   Glass + fire             Glass becomes sand, which suppresses the fire
%   Ice or snow + fire       Ice/snow melts into water; the fire goes out
%   Sand + fire              Fire goes out
%   Water + fire             Fire becomes steam; water remains
%   Ice meltwater + fire     Fire goes out during the meltwater cooldown
%   Water + lava             Contacting water and lava become stone
%                            Cooling can spread into nearby lava
%   Sand or acid + lava      Sand/acid becomes smoke
%   Glass + lava             Glass disappears; steam may form above it
%   Ice or snow + lava       Ice/snow becomes steam
%   Water + ice              Water freezes after any meltwater cooldown
%
%  Acid
%   Contact                  Result
%   -----------------------  ---------------------------------------------
%   Acid + water             Contacting cells can become steam
%   Acid + ice               Ice and contacting acid become steam
%   Acid + sand, stone,      Affected material dissolves; contacting acid
%     oil, wood, metal,      becomes smoke
%     or snow
%   Acid + fire              Fire becomes smoke
%   Acid + glass/concrete    No corrosion with the default material table
%
%  Flow, settling and lifetime
%   Water seeks level pools; oil floats above water. Denser mobile
%   materials can sink through lighter fluids. Solids form barriers.
%   Resting snow under sufficient load compacts into ice.
%   Fresh concrete absorbs adjacent water and dries sooner. Dry concrete
%   stops absorbing water; it retains its powder transport behavior.
%   Ordinary fire, steam and smoke disappear when their lifetimes end.
%
%  Particle lettering
%   Original letter particles stay fixed without aging or compacting.
%   Contact reactions and brush replacement still affect them. Changed
%   particles lose their anchors; products and fresh paint move normally.
%   Destroyed letter cells are not restored.
%
% Contact refers to an adjacent cell above, below, left or right. Reaction
% rates and material properties are adjustable in the controls below.
%
% Edit these controls, then restart. Click or Space selects the next material.
% Grid dimensions are columns/rows. A finer grid costs more and makes grains
% smaller. Hz sets numerical steps per second, not movement speed.
% FlowRate sets material transport events per second. All clocks use dt.
c=struct( ...
'grid',[128 96],'hz',30,'flowRate',30, ...
'order',[0 2 1 5 4 6 9 13 15 7 14 12 11 3], ... 
'names',["" split("Sand Water Stone Fire Oil Wood Acid Steam Lava Smoke Metal Glass Ice Concrete Snow Boundary")'], ... 
 ... % SOURCE: stamps/second (0 disables), radius in world
 ... % pixels, and replacement of existing material. True overwrites interior cells;
 ... % false preserves solid ground. A wider brush deposits more
'rate',30,'radius',10,'replace',true, ... 
 ... % MOTION: speed in pixels/second; center/span/start in world pixels
 ... % Cruise gives seconds between fresh destinations; reach is arrival radius
 ... % Turn limits are normal/boundary radians/second; response smooths steering
 ... % LookAhead is seconds of travel reserved for a curved boundary turn.
'speed',240,'center',[400 300],'span',[390 290],'start',[155 205], ...
'cruise',[1.8 3.5],'reach',70,'turn',[3 7.5],'response',.06,'lookAhead',.48, ...
'headingGain',[3 10],'exitAngle',.5,'retreat',.8, ... % steering gain; turn completion angle; inward destination span
 ... % LANDSCAPE: terrain x/y profile; lake rows [left right surface y], in pixels
 ... % Spring is [x y attempts/second]; occupied cells block emission
 ... % A zero spring rate disables it. Outlet drains the floor
'terrain',[0 260;60 250;110 355;200 365;265 290;295 300;370 500;480 505;545 410;575 415;650 550;700 545;735 490;755 500;780 610;800 610], ... 
'lakes',[85 265 275;320 545 429;600 755 505], ... 
'spring',[155 275 15],'outlet',true, ... 
 ... % MATERIAL TABLES: rows are IDs 0..15 above, then the protected border (16)
 ... % Kind: 0 fixed/air, 1 powder, 2 liquid, 3 gas. Density orders mobile materials
'kind',[0 1 2 0 3 2 0 2 3 2 3 0 0 0 1 1 0]', ... 
'density',[0 5 3 0 -2 2 0 4 -1 6 -1 0 0 0 5 1 0]', ... 
'mobile',[0 1 1 0 .9 .8 0 1 .8 .35 .65 0 0 0 1 .45 0]', ... 
'side',[0 0 1 0 .4 .8 0 1 .5 .16 .4 0 0 0 0 .2 0]', ... 
'drift',.45,'rest',.04, ... % snow diagonal share; rest before density swaps (s)
 ... % WATER: enable leveling; minimum pool cells, resting time, interval (seconds)
'level',true,'pool',[8 .15 .067], ... 
 ... % REACTIONS: burn/melt/neutral/freeze/corrosion are event rates per second
 ... % Zero disables a rate. Cool is ice-meltwater cooldown in seconds; hydrate is
 ... % seconds removed from concrete's cure time when it consumes adjacent water
'react',true,'burn',[4.9 8.6 30], ... % wood ignition by fire/lava; oil-front events/second
'melt',[4.9 3.2 3.2], ... % lava against glass, ice, snow
'neutral',8.6,'freeze',30,'cool',3.3,'hydrate',2, ... 
'crust',[.9 .7 .4], ... % probabilities of three extra lava-to-stone rings
'corrosion',[0 8.6 0 1.5 0 8.6 8.6 0 0 0 0 3.2 0 Inf 0 8.6 0]', ... 
'life',[0 0 0 0 1.2 0 0 0 1.7 0 2.3 0 0 0 3.5 0 0]', ... 
'lifeRange',[.75 1.25],'woodLife',1.5, ... % lifetime range and burning-wood multiplier
'compact',[2 5 3 .15], ... % snow rest (s), load cells, snow/ice ratio, interval (s)
 ... % LOOK: explicit RGB colors and background, independent of the MATLAB theme
'bg',[.012 .018 .028], ... 
'haze',[.015 .03 .045],'depth',.008, ... % atmospheric light; water attenuation/pixel
'colors',reshape(sscanf('94abc7ebb8461891e65b6f85ff9a26536b2596562b8af32fadc5d9ff48235a667aa5b4c34c89b2a5def7acb0b9ebf4ff000000','%2x'),3,[])'/255, ... 
'texture',.12,'shine',.22,'exposure',1.05,'smooth',"bilinear", ... 
'flicker',[13 .18],'ripple',[.10 1.3 .13 .06 .08],'wood',[.12 1.9 .4], ... 
'strata',[.14 .25 .05], ... % rock layer contrast, vertical/spatial frequencies
'wet',[.30 .27 .23],'bloom',[3 1.5 .7], ... % radius (cells), spread, intensity
'emission',[0 0 .03 0 1 0 0 .12 0 .8 0 0 0 .02 0 0]', ... 
 ... % LABEL/BRUSH: label x/y/font size; footprint uses the source's grid cells
'hud',[770 25 16],'font',"Segoe UI",'help',[8 2 400 580 12], ... 
'marker',[.22 .65 1.2 5 .15], ... % fill/edge opacity, line width, pulse decay/second, pulse opacity
 ... % LETTERS: center x, top y and total width in world pixels
 ... % Actual material cells spell E L E M E N T S. Gap is in letter heights
 ... % Original letter particles stay in place without aging or compacting
 ... % Contact reactions remain active. Products and fresh paint move normally
'heading',[400 95 660],'letterGap',.14, ...
'letterWaterDepth',10, ... % shading depth cap for anchored water, in world pixels
'titleMaterials',[2 1 4 5 6 15 11 12],'title',true);
% Keep grid dimensions integral and at least 4, physics rate positive,
% gas/concrete lifetimes positive, radii nonnegative, bloom spread positive,
% RGB in [0,1], and the motion box inside the canvas. Zero source rate or speed
% is supported. None stops brush deposition; the spring and world keep running.
% Only a click changes material. Keep turn rates positive, response nonnegative
% and retreat in [0,1]. Chemistry is an artistic sandbox model. Disabling reactions
% also stops lifetime decay and snow compaction; sources, drainage and motion continue.
% Compact controls are stored alphabetically; edit named values above.
% State: a=step accumulator, b=brush offsets, d=letter anchors,
% f=destination duration, i=boundary-turn sign, j=destination,
% l=brush graphic, m=hint graphic, o=angular speed, p=brush position,
% r=label graphic, s=texture noise, u=paint accumulator, w=heading,
% z=click pulse, A=destination age, B=selection, C=previous position,
% D=physics tick, q=particle data, t=elapsed simulation time
c=struct2cell(orderfields(c));
g.init=@(ax)I(ax,c);g.update=@U;g.render=@R;
end

function s=I(a,c)
s=cell2struct(cell(31,1),split('c q n x y v h e b s k t a D u B z p w i o A f j g l d C r m Op'));s.c=c;
w=split(char(N('􃆩𒰞𓱝򄇋񫦹𱸟󭕚񰔇񉊭𤺦󺖦󬑑񁞉򦦡󲚴𿜍򮲽򗛉򥽅󳰫񌯊򡲗𾖐𺌗󩥜𹼣񋔃񶤧򹐭𖦠󟵥񉝞󦴰󃑄񔎔򠗴񠘚򫢺󶩂𴮦򂍳󖤲󒖶򠗶򡱗󷙂񲢨򑴼񱃒󝠼🜎𷀿񿜋񕹁𛙕񰏔󚽰󡤟񶄮𪷨񣻫񁦮񡦸񚟜𧎋񙈶󲈔򺌼񫀎񚐒񷳩򓎄򈺉򸗲󶋕𶍂񉴄򍺦򓖧򕞚򃵡򗵎󦃳򗞙񏔰򛘪𰈟𞭊󄽔򙀝󢎇𼖟𯼩𵌩𗥰򜳶󉂱򁑛򌺻򕎍񲺀񽡁򬍽򙙹򾧐񰇪𥛆񸎥󷞃񨊈󍯄󘪃󋉯񏯵񅮎򋷎򻧎񨃆󾒎򱋬򄈇𲽜񇩠񬥙񲿷󫳱󙦳񐙕󡭀󼻀󃺋󽝄񼉑𒗜򄶠󀝝󂌷񱹓񓟬򚪂񡎚񜻗񣬙𸿅򉀅𢍸򇞠𺼆𐾁󗈠𻧐񫛦񀱶򕥏򨪽򏰕󖥉󢈀󱌚􂥵򝰥򽳯񯴴􂚇򉗴󚵂񜨋𛑧񟸍񲳚񼫶򶲊',706,127)),"~");
w(48:end)=num2cell(str2double(w(48:end)));
s.Op={{0,0,@deal,@mtimes,0,@eq,0,@and,@gt,@size,@(p,c,j)p+j*~c,@colon,@times,@horzcat,@ge,@or,@mrdivide,@numel,@H,@transpose,@find,@rand,@L,@sin,@min,@ne,@not,0,0,@max,@vertcat,@logical,@B,@reshape,@uint16,@zeros,@exp,@floor,@mod,@rdivide,@A,@AJ,@(a,b,c)a(b:c),@abs,@any,@atan2,@cos,@cummax,@power,0,@diff,@double,@mpower,@repmat,@round,@S,@convn,@cummin,@eps,@false,@randi,@sign,@sum,@text,@true,@AE,@AG,@AV,0,@accumarray,@cat,@ceil,@cell,@conncomp,@cumsum,@findobj,@graph,@image,@inpolygon,@interp1,@isempty,@ismember,@nnz,@norm,@patch,@unique},split(char(N('򃩌񭀔󑩜󟺳򨓲󻧅𖬩󆈊񰢅󃼷񷽶񚤶󊘂𵃩󯾑򄇻𿀞񲋛𚩒󔘭󗸨󑙖񥩑򣄸𴂗𖳳򽞒󊹭𑡙򡗂񖭶򤶗𕺖򐟩󵉦𐄀񚛽򂴅󛃏񋸇񓻟󤯂󍺂󩻓󰾰􁆑򨨒񋷔򎃾򖊇򐰍򰣱򗧙򠩺򇶚󾿢񻠊𥎃򫏼󹕇򈾵󈦩󂩿𧴼𖰭񼍌𱯑񲥝񝴭󪛞󔖘񍕗󁵉񞗒𫶻񢈝񝼬𽗄󾐙򜳕󂑍񋽲𗶃򾯳򴅐򽚢򟵊𬸑𡑙󍲼󥉜󄎃𘕇󙤩򸓬𻐄񣼋򗡠񼫨񍟌󼨱򥉼󖥎񵗋􁉎𼜯􃾏򐹨񘵃񠭞𿠠𢝊𶨳򪅄𗶬񦓵񍵿񱋒󱸣󰝿􀱶񎉭񏡴񤇲򅺓󍜹󘢀񆬇򑐳򠲇񥯛󒖃󅅙󽩟򧲾񈺼𫳋󫕗𮦱򓈓𚒙󜵀񵰁𘭉񧓱󡘥񁩰󽃖񕬅񓢟򍹻񥚩񙦔󮗚򉁻𵠳򝓐򱁣񺝶񔶣򿜻񦝥񴦶񁨁󌉈򟅢󖠐󚴡󾧱򵏓󻀄񯣪󻊔򲦶򟒷󄶾󆲏󢋶򫚝𨿨𺎉񚑛󴵕򱓿񎈅𢔘񖍍󰗐󏻘󧷁𶆲𐙚󅈑񎀨򹧛񕁡󀨬󝓣򪠷򉎑󸤇񺌝𘾉󩌋󱸸򤻑񭂔󡰕󥐩򎳽񕏎󱼘󮨤񮘹򒹧𐀋󌤶񬗥񵁤𹼊󨩜𕈊񎉈󋦘񫗕򔌅򻥾󻸶򌯸𲪬񌙩񏵣򳅚򌸥񩾥𻨽󁐙󐂇񞵆𦠿𹆆𲕷񖒋󚹃񓒤󩐩򗥆𣼡񁒌򄉝򑴖𦱩񭚲𥸞񂨉񇔆񖱦󟚞񳹩𼨶󢏧󭄸򇶤𗒩򡾔򜅐򴗪򪢼󭓡󔛲񵘒𓉥󎝅𑥮򤂿󯔾񁞮򨝖󏋧򋾎񃍩񍪈𽎾򒍩󘸀񟑼򣾓𖷸󦅯󳥄򖰁򨞝󾨛󽎲񘢎𽽊񸺨𽚫񘮝򡜷󼨍񝈅򁘸󞕵򘉵񇋫񛨅𬚼𼪫󡠔𦧀񪁈󽛟򹟡񼡍򤏒󧅊󾥠񤶚𤡔𥄴񣷜򦼶򽗵򘖶񤘅𱅮򓠇񢶶񍍻󷯯񞰮񵤠񑻤򴉱󖽸󚐰󉪈󓿧󿖄񻢯򰝜󿎾𔘻񋄕󂯫󓊛񶅲񓸁񝓨񻰖𶪓񹅌񻩩򂿪𙹹󨭲𾤤򀮎񯲽񱌀󶪡񤖜򾅔󷒷񏖬𙌔󖮲𓪼򳺘󚸑􁑟񗏇򍋺񃩄񾻤񂬇𸻭񠲽񣡸񄖁򤆓󽢌𱴨󾑎󼐰񳱒𠓒𯔸񑕕򽡻󊱆񞲅񞺗򋪆򁥆񞴞񼥰񓸮񺦪󎼅񊛗󼧻򃈞𡝐󨨽򛈴񱽖𤠒󷠍񛥍񉫍󛑃󌁞𫩿񋮘󺶝򁸍򌡟񹃧񘈦𭌳𕳑򢨵󱁞񕰔󫣒𕈪򾿠񃘺񕨄󩥨񸛱󹱒򍯦򓵸򴀰񰊱񆝢󊏥񇩊󣒤𨺦򥼼򖺍񔷥񴂐󻂸򌎐𡸼𲏡򣓡򧆻򻨺󫍹򱾝񌧧󖡞𝩰𒽫򣸧󖮆𰣙𺎹󫳛򼖀񆑃񧧸𸶬򁙵󶗘񳐩󩊒󔿥𒴁򣵄󉶦𦑕􂴡񒙿񠟸󓫺򊢌򖇖􀬖򒵵􀤛򻠂􁘹󴢉񖈊񊞊񇨊񧖘󑧬𫤑񍌈񢊝𸛾󈘟񞄄񨡆󿣳𝧌𾵨󤨇󎏎󔹫񚭦𼊼򹭉񞘹󻴌򝒭󹠐𔿁󮛿󴹣񽷚򖁅񻟦񟰟𗆿󎷇󈮊󜊥򌯣񽹯񓋫񊤍񗖨󜉱𻕒󹯲󠅃𠚾򕨨򫺴񑮴񁨌󆡍󇼙󼂙񰳵𙸰𼛞򻄖󀁆񡸠𾅡󮪲򆺠𴵜𛮾󾭻𻒕𪄬񲱖񋁎񛷨򤠸򃟚𰇧𺜻󀡾񈲻󍢑񺩜򨀣򹰃󍣍򼐤񲖼򆁠󭏑񔉲򙐏򬊮񨷽🟾󪜥𱚮󣱻񷿦򳴑񓊾񐡋𙷺󆡿𙍛𵯯󁶅򏀦񱇶򞂍󢝵񠏌񓠴󕃑󌅎򌋅􄂷𩦢󺈸򡋅𱍓򎵷󈥨񾸾𕒖󴮶󱍫񲯵󺽤󱛰򰳪򪟛򓌘𼶴򈧦񎴓󬫀񰅥򓇟󫹀񓨴𝛂󀽬𲣰򍦐񹤊󗱮𢁊󵼎񸦥𢀝󗯲񺪪𕨫񇎑𸲖񋰶𖧴󦖸򹳀񦩄𑻅򟬕򬥴񰊃󨊘񅐌󔺺򦎋𩼢𜬺򃞬󖰏򽎙󡸧򂳚󜴤򠷄򤛺񭼀򧚞񁵸򫘷𬥪񹱡򑧤󟒾𪚀򳆋񕁆🡅󕶥񩊣𘮧󞁞𨡺񥥥񊺚󑎕􁙟񖞤񧾠𚯌񐶨󛢟򏌋𥎟򧅷𒗔󈗡񀫠񅷅񦊭񜆧򽇫𧂅󈋺󤽹𖖾򨧮񕏌򍙠򛷖󫊒񘨗򴾺񯌋󮺩󊎡󄻇𶪐𸒆􂤒𒢴󠍑񗢛󕋯𩊯򛌇󤨏𑉪󹇩𱱏񩻱򞳔򧤴񔑲𫛥򜾓󤏗򇗧򙫢򂛹𰪾󪎦󀍻򇶺󬉃򾖑򁏗𴧧򖛴񵱽󌾵𔚩򑏭󳵈򩕯򯲡󂬗򒟸񳷿񪻪𴰡𗙃󑆴񴿡򇝝򲉵򰦩󔹸򽎃񶎽񏈺󮬑𺱻򅸼𣩶񑯩򵆕󝓂򬘇񐩳𥚤󃩄󨋫򽾛󁀐𹄯񨘼󟂂󴔋񸺌񲚆󞲄񯋥𬟩𢿄񄘚𞂅򩯨󁷈𑡒򰤓󿓠󤇄񛁌󓽕򧭸򱶯򖳡𬏐󣪚񓃮󊱫󤘟񩠨󄐗򴝾񖨳𯀆򡜄𺆼󉟵򫙤𴇖򁇑񼲜򖩤󵂨򁇔󼈲񰸔򋬿񐭯񃨱󴪧𛐬𕴹󹾽𼰵񞤋󺧿򛽑򅖅󽧥򨅎𰩾𦘔􀏜񳣯󾈆򅉪񘹐􁹗񅆚𾣦𫖓𚐥󫽊󎲡􀂧󹩻󭎰󶪣𢠲򹟉𒔟𺸴𗫓𤆪󰞩󽖊𪾸󼢡𣄪򪍽󥪹𒝳󺲔򎌋󰊟󽚗𛐸󀞍򳄔🇾𡪟󄝙񄯳񧽉󴹓򤋤󵞕񵴦򟬢󁬸􀦎򱻢󻇏񃽪񼀾𭀘񕳥񳥺𓦺񼍕󷤖󜪬𼋍򕏌񳷛򨚉𒜻󾊴𰈧񤋅񄽕𐥏񗒵󣣾񔲬󷑖񜔞󱼮󅁰򮩏󿝞񔤶󪭫񉣪񻖳󸿫񫦆򪽧򨍵𷠲󑏔𵿩񰰢񋨳𕂵󳈏򬱆򙽉񁰵񡡡𰍩򡉣񗙳񡖒򼺅򂚘񐆙󭌟񺱳򛍡򢠷񁬆򑡟򐒘𞺐󙉕񳕠󋗸򮆎򔌮򭥉󨺲򜾼񰊦񪼶󳖻󷁡󁃧󦊫񟧡󡟭󐦍񉔽󧸵𳖱󹁊񮔾󀼃􀊐񕅔󃛇򌑛󰷔𩸲󷥈𡣴󱴥򊧆𽁥򇡴񞳲񟌥􂞟󯛦񷻹򁄟򜦀𝥩𦓷񠝰򵿇󿳙򋄅󮏑򪰔󈫛򈉦񋢟򕌅𪂵󄸿򿡑񷨊𒁕񗸃󾾛򛟧񬵾񠼤󿊜񍧍񗘐򸅌򾠉𰂾🙻󘽅𶲪󦗮𒫅񺆞󎈂󚄬򴞴񈬮𞶟󴃰򇂰񎤪񇂟񒨝𨈠񠰎򐐅𕵣񿐘񪏕𸅔󹲧񠜨󨻦򞅃󣰈򽭂񢋲񀻜񸘚𜮼򪄕󥽪򽑮񕔨󴐤򉍱񁛰񈺿򦢰򂖃񀫼񙫀񜎘𷒎򵰭򨖪𭪞򮺿񭣼󿙐򣳾𾑕𮩋򓯽򵴤󷊶𠨉򊯱񎭸󆡂񸨴򟭻󧹯񘎂󙏗𾏮𽲟򱫙􂁧򛣜𯖩񑞨𹠫󫑀򛗂񽒶񃟅𧚈󫊌𤓱񎮏󊓢𩼈󉸴񥩉򳌢񰯦𹬮򼘧𸩄򶣭󳵵𢎢󑋏󾵌󿇼󌫡񕅇򵘾򑊈񝮩𙮔𫾯򄏔񇀻񱽚򡚄󫣮󰆞񻖤󖞈򐚼𴻄󿿞󰞠󠷄𐝯󫺠񈩊􃽰񨥁𻼞񖫞񨴓󢡻𧨪󮪏𙔀񼂻󃇧󱢄򬋞򮯵򳺖񬳛򑻈𔪂񰙲񪛆񵖨񓘆󵯪򲖂񦡈񳅇򧄘󾽈򁎟񲨴𘢼󭠯󧵭񼦽񮷞𞛥󭫏󍘱񊨍𾬌􁰶򇘻󳥍񥽣󛥵𖱂򘌛򁛏𔆞󖜹񆉠񠳘񨪰򭇵񦼘򕷍򄙋🲯𩧻󮲚񀚲򎨹񛞥󬅞𭞪𥻧񠡢󛈘򽔝󜑞򚕿򟭪􁨅󣚵󐟜𳓵񔄺򗬥񗉟󎖳񎩧񧇴󨼄򶁄򟄗򖰟򠶺𠵥𩻝񬮕񁝌񮶷񖻝󄗠󿞹򦊛󎤪񵻤󈱯񞕸󎾛򫬜𨜈񛮪񃺸𻴇򙗱󰌤𑺤󺵀񒤶񂾙󱏠񻔁񖭄󔩨񬻯򛏠񙭚򌬫𑏑𵉤󟬦󊫫󍊶򶴯𤐖𢢗򒝳󱆌𤼛򧹯񁡀񥵞󘆰񪤲񠯔󡭼󫶖󙀻򠲨񤟞񯥀𛪅󠃎󍴼󭃅񻘝𝖵𐦲󮾋񱲛򼦾񄋪򕍴񹪷򨯩󞴎񰈮򤦖𬌖򋔙򫎾𨧡񾊓򗬢󿌆񢬜񟜠󖇂󗪨󉔻򀁂𣲌񟖺󐪵𹟲𪚍񞸸𣒯񽰇񬼺𲫸󿱴񮮞񹽗񼀀󲥸󸺂󏮩򨬌񀾋򸩛𼯱󭸏񇎅󣪯򉞋񙔂񒳈􂼟𨆟񩎦󼩇񲃆󆻑񛐥􁱻𘴾𦞳򙗼򗮖󗎲򰦼𢛙񐩱񉾿򀍶𱼞󐧶򷭞󷑹򎽸󹇾󋌾𿰯𗿾󋉎򪃚󊸅󇨉񬶯񻜥򶝊󠛜񚙰𥺪򽒗񯔢󔫂򚮝󪱅𙆄򩼑𱩢򘸇򮈸𴉫򴧢󣐲󀙦󀹓󽱧򊩍򳜻𬄝򔬵󷃅򺽅򄺛񚦜񵉰𥣲򄡌󾌖𾿈𴜄򕩖򇃁򋞳򇜀󚏉򾨿񾚦𙰂򖪑񢨫𥧿򆄣𪘖󾡝󄯻񃹇󧈵򥣃󁋣󀹀󩡙񝩼񙦨񲳄󬞆񏒯🵓󉠉񟺺󣄸𝁡񳣞𐃶񰧡񙱁񜢅򠎕󰑐𼂃󭲤򦁫񠢁󫌤򀪅𾧨񴼧򷲀𜕀򎻼𐾕񗄟񐇁񄐄𙉞򯰜򕣰󼡋񤠒򴱗񘌓򔏏󵸢򧗍󣏗󮅋𨜮񗌟󪃃򕂦򩿅񀁪񫉞򢘙𔧎񴥓񹤫𒯠𦔴񚪁򟽞񁗰𵻮񬍢򸺴󦴟򚎆󔱾񚜒𳶞𵝶򃤛󕆮񂬾𙸂󗘝򛚝񕄃𽽀򽄋򮗷򘢉𻍩񡑤񔇀񎘢󕳢􁥣󈨂󽉿򱛎򻀇󎉞𺏲󉜗𕘺󧤖󩌌񟞓񸘞󡨨󱜉񫧊󖨥񋢱񾾹򮆴򎥘󝌑򼊓򨙟񫄅򝻥򡬙𴝄򘚎񑀡󵄼񯐎󃳔򬔯󾒌󡑨𞡶񾋤񊨪򌶸󴹍򷟰𡑢󌍂򅉴񴜪𞥂󝰤򃤴򪢂򦎂򯻉򅼭󄏷񆶆󢌂𰶫򴲐򌟼𣷨񔜔󥵁򦫖󢂹򻅆񧬮򤐊򗔶򷩰򄂘𒊹񾪕񞧚𲗟񡗞󣁹𬝖򋲌𿚖𳩵𾡵𚤓򾛙򦿌󜏟󌴒󥪍󘳆򇌘񳵔󆊇񷽰򳏄󳨽񙂁𗘴񘜴􀇋񽻶񥈮򅏼򱰛򗟡򅪔𛏶𕵘󎮹󡉿𙸻󤺭򳦌񓩇󘑛𤸂񇟌󚩵🳶񴳄񏧩󾢛򡳎񘄱򿅩񡞕򛓗󒯑򳋑񺳃􁑝󧡯򃸲򗒫񮅥񆣔򊚲񯙳𿳎񶗭󧰻񰋅󥺔򉼃򟖉󂗃𙜜񅎘𮇋򯷚򖜛𸀶󿐿񸗆򣌂򊖾󒇃󻋆񊼳򦡷򣭲񮧳񃺢򇈷𷞇񂮮򍁍򚋏񵸏𫝦񀆼󉢻򘞰񤃹񛛓󅇧򀿰󮪝򑼐򪖹򧠑󧉮񴎪𷻧񻑨򌁑𹺦󞈆򺄠𰐧񂱸𛀹򝦓󇧌󏽟󞗃󫼯񹿟񀚆񏰴򙡙𽗭񒏚𷣣𭫒󜤥󬔵򉾙󨐒󆢍󷳯𙦢򤁻𭾘󢉸󠍊𢰴񺶗񗲽󕹂򣲓󸟰񸭹𽫧󒱈򵳎󤷎񚎾󯶯򟏼𞹅򧡐񲡇񧢝򁼬򼜐򑊉񍮉🴂񍎙𑲊󂜎󣢯񪑄񇟁򐙉󔈠󰩺𢝇򏧗򚢂𥅱󣧾񇨌񖋝򑩡񩗏񤀨𹻗󟨓󇧢󏻠򀶦񈼑򘏓󜓩񚕃𧋏򵱜񼨚󙣣򋼭򟰶򅽏񝸫񼴰𻛥񾙁𯾘񙆟񈔤󥘐𖥍𰙒򠨿󄹠򵀑󭓛򦭃򫄬𮖠􃻄􂘫𦆛𢯲𑼱𸩴󳎋󺏣򧤨𕇑񸌸𯏣򻼳񯁝񯏑񹡃󥏜򷤋𬝽񼍝󲇪򳓢򎢡򸭀󲶨󒨈𑄲񫠯𕏎𣈸𹧨񷹂򯴄𰾒񃾩𮨩󑾣򜳙񋁨񠥺󡿇🬗𺁵񬮙򩕦􃇜񛍅󁁳񼔒󮴕󧖋𽟉𦍆򹆍󊤧񣓘򁴑񩵯򷽍񦉪𓵥񔚝򙭰󤺨򽒊򇱔񇠝򣽪𨃱󦼗󛿁򗴱񍣢𘂣󥱍򳸮󠹪򞲯󃪍񴆳󝰡򂩷𿬓򂜊󙥿򃶈󹅽􂍇󡇽󊤔񡎈򈸞􂈵󌍕񤜕񕣘񷷩𤎑񔮓񃮰𵸁󊩵􃭛󭕈𦢁񺹤󏂧򼏚𮠢򂺥󐣏󹯝򴋺󉠰󵩼򷁊񁓪󐷠񼷞򤴙򃢯򢢺򺾏񊛽𾇹󄧱󖧻񻗎򅍍򸠻𬔹򆶠򙳒󩴙𔭆򋸤򆬤󹯒񜥡򻱄𱆭񆵵򯶃𽶘𩿨񥗝񢳢򜔑򚙖񰚥򬢯񡅏󽝲򄎆򈐡񢦆񂣃𨳫󞚹񄾦񳆋󖮑𴃅󤩤򿧒򛳹󷻤񂀻𩸢򧋼񊐶𯻋𰅮󱥰򥏋𖢶󬷧󯝇𜉛񈉸󚨚񱴀񌤚򳋷򅅝񃫂󫌌󜬙򼺅𯝃򽔃󆎉񚐋𗀄򌿥򹪱򑡌񴇎񧌾󐍆򆯫񦬻򺷦񅽣󉅳򣌆󠗈󒹴񃧏򿥙𚀅𮠪󛊸򪘤󨞒򤢅󡑡󕲕󔘩򾔒𦈊򺰯񆷎𪥙񺝵򠣘򛨼񆨒򃥘𨇠󰼻𡃉񾺐',8388,266)),"~"),[cell(1,146),w']};
s.d=N(['𓪏򦗷𺬹𳜰󏁕򦉰򟂨󥪕󐡺𨱸򧻵򽗥򚧇󹨩󨏝񶷔𵇿񼱞𪢒񏊈󐓖񧯠򥦐񎻰񳗐򼔻򵔞𱯨󠻤𩔗󅂃򭯊񚠐󫜔󣥏򱦱󫅎󥋷󩖒񣤄󽩊򨆖󁨬' char([55422 56676]) '𶪑򝉙󃡡𧔠񯫶򒹁򋫘󡕷󷁀򂅗򎑠򣔝󅇧񈝘󚙅󉩎񸺯򟓏𕾋𹲖󜵛񝈧' char([55422 56527]) '󛤷𧍪򰕳󩤹󙸖𩖣󑦐񯸜񄪩񄗈򦓪򕎁󇗤񔯇񘤌򂋙򪙀񳕃󜫛𡓚񱔐󅻢󐨐񝙜򌄢򹝅򢒗󡓃󭽱𭸾񍡘򼎧򻍁󞒄򭒒񊴻񧵦𪂀󒺌𮤴󷘆񱎚򃗁񉞊򵟿񇬐𶙩񦗴񮶝񿖼񦾒񳯃󶋺𫐿􂋨𞆂󕬅𰌵𿟮񼒛򴹽򅻔𶲖񮩖󹘃񔑝󊤎𾤏򬄗󐨵򎮓𢋞𤙱񙭄󲊀𡱄񸚹񓁮󼃊𛂖򔽌򗊋򋉉񹐃𑒥󰻷󵒡񪜌񖆭󥶤񘴧𨺺򪊿𠂫򆅌󛠣񭔤񾀄񃜬񆒅򌾮󇳟򏴉򯵛𔄁𪑡󺤹򅄯򩇰񰠺򳿡񷁛񶻡򓻐򱞱𚶥񢿩򶀰󩒺񊛮򘖣𛆛򬢔󸇷𤋞'],667,64);
s=T(s,9,0,0,0,0,a);
end

function s=U(s,a,b)
s=T(s,10,0,0,0,0,0,0,a,b);
end
function a=AG(q,z,r,c,d)
% Shift whole falling/rising chains, clearing only their vacated endpoints
m=q(:,:,1);n=numel(m);
k=c{29}(m+1);a=k>0&c{11}(m+1)*d>0;
if d>0,a(m==15)=a(m==15)&rand(nnz(m==15),1)>c{13};end
if ~any(a,'all'),a={q,z};return;end
f=false(size(m));b=2:size(m,1)-1;if d>0,b=flip(b);end
for y=b
 f(y,:)=a(y,:)&~z(y,:)&r(y,:)<c{39}(m(y,:)+1)' & ...
  ((m(y+d,:)==0&~z(y+d,:))|(m(y+d,:)==m(y,:)&f(y+d,:)));
end
u=find(f);v=u+d;e=f;e(v)=false;
q(v+[0 n])=q(u+[0 n]);q(v+2*n)=0;q(find(e)+(0:2)*n)=0;z([u;v])=true;
a={q,z};
end
function a=AV(q,z,u,v,r,c,d,t)
% Disjoint density/diagonal exchanges, including transported lifetimes
m=q(:,:,1);n=numel(m);
b=m(u);e=m(v);h=c{29}(e+1);
a=c{11}(b+1)>0;
k=c{39}(b+1);k(~a)=c{39}(e(~a)+1);
l=r(u);l(~a)=r(v(~a));
f=t(b+17*e+1)&~z(u)&~z(v)&l<k;
f=f&(e==0|h==3|(q(u+2*n)>=c{51}*c{28}&q(v+2*n)>=c{51}*c{28}));
if d~=0
 g=m(u+d);g(~a)=m(u(~a)+1);
 o=g==0|c{29}(g+1)>1;
 f=f&o&b~=14 & ((a&(m(u+1)~=0|b==15))|(~a&m(v-1)~=0));
end
[q,z]=S(q,z,u,v,f);
a={q,z};
end
function [q,z]=S(q,z,u,v,a)
% Exchange disjoint selected pairs without destroying either endpoint
n=size(q,1)*size(q,2);u=u(a);v=v(a);
q([u;v]+[0 n])=q([v;u]+[0 n]);q([u;v]+2*n)=0;z([u;v])=true;
end

function n=H(m)
% Orthogonal adjacency; zero padding never wraps opposite edges
n=conv2(double(m),[0 1 0;1 0 1;0 1 0],'same')>0;
end

function p=B(a,c)
% Convert a per-second reaction rate to a fixed-step probability
p=1-exp(-a/c{28});
end

function l=L(a,n,c)
% Initialize lifetimes in fixed-step ticks, with a bounded random spread
l=uint16(c{28}*c{34}(a+1)*(c{35}(1)+diff(c{35})*rand(n,1)));
end

function b=AE(b)
% Clamp and quantize RGB channels in one native array expression
b=uint8(255*max(0,min(1,b)));
end

function n=AJ(t,r,h)
% Count elapsed periodic events without rounding the requested interval
n=diff(floor([t-1 t]*r/h+1e-9));
end

function R(~,s)
T(s,8);
end

function s=T(s,varargin)
% Execute the fixed operation table; registers hold fields and temporaries
[f,b,r]=s.Op{:};r(2:nargin)=varargin;r(12:42)=struct2cell(s);r(43:111)=s.c;
while r{2}
n=+b{r{2}};r{2}=r{2}+1;d=n(2);
switch n(1)
case 1
r{d}=r{n(3)}(r{n(4:end)});
case 5
a=r{n(3)};a(r{n(5:end)})=r{n(4)};r{d}=a;
case 29
a=r{n(3)};a.(r{n(4)})=r{n(5)};r{d}=a;
case 50
set(r{n(3:end)});r{d}=0;
case 28
r{d}=r{n(3)}{r{n(4)}};
case 69
a=r{n(3)};a=a{r{n(4)}};r{d}=a(r{n(5:end)});
case 2
r{d}=r{n(3)}+r{n(4)};
case 7
r{d}=r{n(3)}-r{n(4)};
otherwise
r{d}=f{n(1)}(r{n(3:end)});
end
end
s=cell2struct(r(12:42),fieldnames(s),2);
end
function a=A(a,v,varargin)
a{varargin{:}}=v;
end
function v=N(t,n,A)
% Decode instruction, constant and outline streams during initialization
[a,v,r,k,z]=deal(zeros(A,16384),[0 0 0],1,1,0);
for i=4:n
if r<1e9,[r,z,k]=deal(r*1e6,z*1e6+t(k)*1024+t(k+1)-56679424,k+2);end
h=mod(cumsum([0 (v(i-(1:3))+1).*A.^(0:2)]),16384)+1;
f=1+a(:,h)*12.^(0:3)';c=cumsum(f);u=floor(r/c(end));j=find(c>z/u,1);
[z,r,v(i),a(j,h)]=deal(z-(c(j)-f(j))*u,f(j)*u,j-1,a(j,h)+1);
end
end
Play Game
Preview

Play Game opens this in MATLAB Drive. Click Run this app there to play.