31:Tc18,grid = { //force this block to re-run when randomize button is pressed randomize; return ( d3 .range(zCount) .map((z) => d3.range(yCount).map((y) => { return d3.range(xCount).map((x) => { //create identifiers for this cell and neighboring cells const coordinate = { x, y, z }; const neighbors = neighborModel == "Moore" ? d3 .range(-1, 2, 1) .map((zD) => d3.range(-1, 2, 1).map((yD) => d3.range(-1, 2, 1).map((xD) => ({ x: x + xD, y: y + yD, z: z + zD, id: `${x + xD}-${y + yD}-${z + zD}`, index: (z + zD) * (yCount * xCount) + (y + yD) * xCount + (x + xD) })) ) ) .flat(Infinity) //remove self .filter((d, i) => i != 13) //remove outside of grid .filter((d) => d.x >= 0 && d.y >= 0 && d.z >= 0) .filter((d) => d.x < xCount && d.y < yCount && d.z < zCount) : [ [1, 0, 0], [-1, 0, 0], [0, 1, 0], [0, -1, 0], [0, 0, 1], [0, 0, -1] ] .map((d) => ({ x: x + d[0], y: y + d[1], z: z + d[2], id: `${x + d[0]}-${y + d[1]}-${z + d[2]}`, index: (z + d[2]) * (yCount * xCount) + (y + d[1]) * xCount + (x + d[0]) })) .flat(Infinity) //remove self .filter((d, i) => i != 13) //remove outside of grid .filter((d) => d.x >= 0 && d.y >= 0 && d.z >= 0) .filter( (d) => d.x < xCount && d.y < yCount && d.z < zCount ); //determine if cell should be dead or alive at start const alive = Math.random() < livingChance ? //? Math.floor(Math.random() * (stateCount - 1)) + 1 stateCount - 1 : 0; //build cell return { coordinate, neighbors, alive, id: `${x}-${y}-${z}` }; }); }) ) //remove unneeded data structure .flat(Infinity) //for performance, find the index of all neighbors so we do not need to look them up every frame. .map((cell, i, grid) => { cell.neighborIndices = cell.neighbors.map((nb) => nb.index); return cell; }) ); }32:T156f,simulation = { //check for button press restart; //for managing the calculation of the next simulation step let currentGeneration = grid; let nextGeneration; //loop forever, so long as toggle is true while (true && play) { //wait a bit to give the sim time to calculate and draw await Promises.delay(generationalDelay); //loop through all the cells nextGeneration = currentGeneration.map((currentCell) => { //count up how many neighbors are alive. const neighborSum = currentCell.neighborIndices.reduce((acc, nbIndex) => { if (currentGeneration[nbIndex].alive == stateCount - 1) { return acc + 1; } else { return acc; } }, 0); currentCell.livingNeighbors = neighborSum; //assume the cell does not survive to next generation let survive = 0; if (currentCell.alive == 0) { //if the cell is dead if (birthCounts.includes(neighborSum)) { //but it has the right number of neighbors to become alive survive = 1; } else { //otherwise stay dead survive = 0; } } else if (currentCell.alive == stateCount - 1) { //else the cell is alive... if (surviveCounts.includes(neighborSum)) { //and it has the right number of neighbors to stay alive survive = 1; } else { //otherwise the living cell dies survive = 0; } } else { //else the cell is dying survive = 0; } //store survival state currentCell.survive = survive; return currentCell; }); //update the cells all at once after calculating if they should survive or not //this second step is necessary to avoid updating some cells before others //which would lead to unexpected outcomes nextGeneration = nextGeneration.map((cell) => { let nextGenAlive = 0; //dead and dying if (cell.survive == 0) { nextGenAlive = cell.alive - 1 > 0 ? cell.alive - 1 : 0; } //dead but spawn else if (cell.survive == 1 && cell.alive == 0) { nextGenAlive = stateCount - 1; } //living and survive else if (cell.survive == 1 && cell.alive == stateCount - 1) { nextGenAlive = stateCount - 1; } //decaying else { nextGenAlive = cell.alive - 1 > 0 ? cell.alive - 1 : 0; } return { ...cell, alive: nextGenAlive }; }); //loop one more time to update the cubes in the visualization nextGeneration.map((cell, i) => { //create a matrix to manage 3d placement const matrix = new THREE.Matrix4(); if (cell.alive > 0) { //if cell is alive put it in the right place matrix.setPosition( cell.coordinate.x - xCount / 2, cell.coordinate.y - yCount / 2, cell.coordinate.z - zCount / 2 ); cubes.setMatrixAt(i, matrix); } else { //otherwise put it really far away so the camera can't see it. matrix.setPosition(10000, 10000, 10000); cubes.setMatrixAt(i, matrix); } if (colorLogic == "Decay") { cubes.setColorAt( i, new THREE.Color(colorScale(cell.alive / (stateCount - 1))) ); } else if (colorLogic == "Position") { cubes.setColorAt( i, new THREE.Color( cell.coordinate.x / xCount, cell.coordinate.y / yCount, cell.coordinate.z / zCount ) ); } else { cubes.setColorAt(i, new THREE.Color(colors[0])); } }); //update the simulation state currentGeneration = nextGeneration; //ask THREE.js to update its internal data model cubes.instanceMatrix.needsUpdate = true; cubes.instanceColor.needsUpdate = true; // const iS = isoS.marchingCubes( // [xCount, yCount, zCount], // function (x, y, z) { // return !Number( // currentGeneration[z * (yCount * xCount) + y * xCount + x].alive > 0 // ); // }, // [ // [0, 0, 0], // [xCount, yCount, zCount] // ] // ); // const vertices = new Float32Array(iS.positions.flat(Infinity)); // isoSurface.geometry.setAttribute( // "position", // new THREE.BufferAttribute(vertices, 3) // ); // var indices = new Float32Array(iS.cells.flat(Infinity)); // isoSurface.geometry.setIndex(new THREE.BufferAttribute(indices, 3)); // isoSurface.geometry.attributes.position.needsUpdate = true; ///////// // if ( // currentGeneration // .filter((d) => d.alive > 0) // .map( // (g) => // new THREE.Vector3( // g.coordinate.x - xCount / 2, // g.coordinate.y - yCount / 2, // g.coordinate.z - zCount / 2 // ) // ).length > 3 // ) { // const geom = new THREE.ConvexGeometry( // currentGeneration // .filter((d) => d.alive > 0) // .map( // (g) => // new THREE.Vector3( // g.coordinate.x - xCount / 2, // g.coordinate.y - yCount / 2, // g.coordinate.z - zCount / 2 // ) // ) // ); // convex.geometry = geom; // } //redraw the scene renderer.render(scene, camera); yield nextGeneration; } }33:[] 2f:{"title":"3D Cellular Automata with Persistence","cells":[{"id":0,"value":"# 3D Cellular Automata with Persistence\n-----\nBased on some of the instructions, guidance, and rulesets [here](https://softologyblog.wordpress.com/2019/12/28/3d-cellular-automata-3/) as well as the always useful [tips and tricks](https://discoverthreejs.com/tips-and-tricks/) resource.","pinned":false,"hidden":false,"mode":"md"},{"id":164,"value":"-----\n### 3D Grid and Initial Condition Settings\n","pinned":false,"hidden":false,"mode":"md"},{"id":6,"value":"viewof xCount = Inputs.range([0, 100], { step: 1, label: \"Columns\", value: 40 })","pinned":false,"hidden":false,"mode":"ojs"},{"id":10,"value":"viewof yCount = Inputs.range([0, 100], { step: 1, label: \"Rows\", value: 30 })","pinned":false,"hidden":false,"mode":"ojs"},{"id":313,"value":"viewof zCount = Inputs.range([0, 100], {\n step: 1,\n label: \"Layers\",\n value: 50\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":147,"value":"viewof livingChance = Inputs.range([0, 1], {\n step: 0.001,\n label: \"Living Chance at Start\",\n value: 0.3\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":137,"value":"viewof randomize = Inputs.button(\"Randomize Starting Condition\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":171,"value":"-----\n### Build Cell Grid Data Structure","pinned":false,"hidden":false,"mode":"md"},{"id":940,"value":"livingPercentageAtStart = d3.filter(grid, (d) => d.alive > 0).length /\n grid.length","pinned":false,"hidden":false,"mode":"ojs"},{"id":28,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":175,"value":"-----\n### Simulation Settings","pinned":false,"hidden":false,"mode":"md"},{"id":854,"value":"viewof neighborModel = Inputs.select([\"Moore\", \"Von Neumann\"], {\n label: \"Neighbor Model\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":871,"value":"viewof stateCount = Inputs.range([2, 10], {\n value: 5,\n label: \"State Count\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":168,"value":"viewof generationalDelay = Inputs.range([0, 10000], {\n step: 10,\n label: \"Delay Per Generation\",\n value: 10\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":500,"value":"surviveCounts = d3.range(9, 17)","pinned":true,"hidden":false,"mode":"ojs"},{"id":490,"value":"birthCounts = [5, 7, 12, 13, 15]","pinned":true,"hidden":false,"mode":"ojs"},{"id":1108,"value":"-----\n### Visualization Settings","pinned":false,"hidden":false,"mode":"md"},{"id":790,"value":"viewof bg = Inputs.color({ label: \"Background Color\", value: \"#001b42\" })","pinned":false,"hidden":false,"mode":"ojs"},{"id":1006,"value":"viewof colorLogic = Inputs.select([\"Decay\", \"Single Color\", \"Position\"], {\n label: \"Color Logic\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1051,"value":"viewof colorCount = Inputs.range([0, 10], {\n step: 1,\n label: \"Color Count\",\n value: 3\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1056,"value":"viewof colors = Inputs.form(\n d3.range(colorCount).map((d, i) =>\n Inputs.color({\n label: \"Color \" + i,\n value: d3\n .rgb(d3.interpolateHcl(\"#ff5da4\", \"#FFF500\")(i / (colorCount - 1)))\n .formatHex()\n })\n )\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":978,"value":"colorScale = d3\n .scaleLinear()\n .interpolate(d3.interpolateHcl)\n .domain(\n d3\n .range(colors.length)\n .map((d) => d / colors.length)\n .concat([1])\n )\n .range([\"#000000\"].concat(colors))","pinned":false,"hidden":false,"mode":"ojs"},{"id":86,"value":"$32","pinned":false,"hidden":false,"mode":"ojs"},{"id":193,"value":"-----\n### Visualization","pinned":false,"hidden":false,"mode":"md"},{"id":119,"value":"viewof restart = Inputs.button(\"Restart\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":759,"value":"viewof play = Inputs.toggle({ label: \"Play\", value: true })","pinned":false,"hidden":false,"mode":"ojs"},{"id":418,"value":"renderer.domElement","pinned":false,"hidden":false,"mode":"ojs"},{"id":821,"value":"-----\n### Download the 3D Model\nPress \"Generate Model\" button to capture current simulation state, then click on the vertical ellipsis in the margin below, and finally \"Download Text\". Rename the downloaded file with the \".obj\" extension.","pinned":false,"hidden":false,"mode":"md"},{"id":812,"value":"viewof update = Inputs.button(\"Generate Download Model\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":801,"value":"CAmodel = {\n if (update) {\n var group = new THREE.Group();\n\n var count = cubes.count;\n var geometry = cubes.geometry;\n var material = cubes.material;\n\n for (var i = 0; i < count; i++) {\n var mesh = new THREE.Mesh(geometry, material);\n\n cubes.getMatrixAt(i, mesh.matrix);\n mesh.matrix.decompose(mesh.position, mesh.quaternion, mesh.scale);\n if (mesh.position.x < 500) {\n group.add(mesh);\n }\n }\n\n group.copy(cubes);\n group.updateMatrixWorld(); // ensure correct world matrices of meshes\n\n return exporter.parse(group);\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":740,"value":"-----\n### THREE.js Stuff","pinned":false,"hidden":false,"mode":"md"},{"id":1234,"value":"// exposedCells = simulation.filter((cell) => {\n// const neighborSum = cell.neighborIndices.reduce((acc, val) => {\n// if (simulation[val].alive > 0) {\n// return acc + 1;\n// } else {\n// return acc;\n// }\n// }, 0);\n// if (neighborSum > 0 && neighborSum < 26) {\n// return true;\n// } else {\n// return false;\n// }\n// }).length","pinned":false,"hidden":false,"mode":"ojs"},{"id":1180,"value":"// convex = {\n// const geom = new THREE.ConvexGeometry(\n// grid\n// .filter((d) => d.alive > 0)\n// .map(\n// (g) =>\n// new THREE.Vector3(\n// g.coordinate.x - xCount / 2,\n// g.coordinate.y - yCount / 2,\n// g.coordinate.z - zCount / 2\n// )\n// )\n// );\n\n// const hull = new THREE.Mesh(geom, mat);\n// return hull;\n// }","pinned":false,"hidden":false,"mode":"ojs"},{"id":1196,"value":"mat = {\n const m = new THREE.MeshBasicMaterial();\n m.transparent = true;\n m.opacity = 0.2;\n return m;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":414,"value":"height = 600","pinned":false,"hidden":false,"mode":"ojs"},{"id":536,"value":"geometry = new THREE.BoxGeometry(1, 1, 1)","pinned":false,"hidden":false,"mode":"ojs"},{"id":519,"value":"material = new THREE.MeshPhongMaterial({\n color: 0xffffff,\n flatShading: true,\n transparent: true\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":533,"value":"cubes = {\n const mesh = new THREE.InstancedMesh(\n geometry,\n material,\n xCount * yCount * zCount\n );\n\n let i = 0;\n\n grid.map((cell, i) => {\n const matrix = new THREE.Matrix4();\n matrix.setPosition(\n cell.coordinate.x - xCount / 2,\n cell.coordinate.y - yCount / 2,\n cell.coordinate.z - zCount / 2\n );\n mesh.setColorAt(\n i,\n new THREE.Color(\n cell.coordinate.x / xCount,\n cell.coordinate.y / yCount,\n cell.coordinate.z / zCount\n )\n );\n mesh.setMatrixAt(i, matrix);\n });\n\n return mesh;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":456,"value":"lights = {\n const light1 = new THREE.PointLight(0xffffff, 1.5, 0, 2);\n light1.position.set(xCount * 2, yCount * 2, zCount * 2);\n\n const light2 = new THREE.PointLight(0xffffff, 1.5, 0, 2);\n light2.position.set(-xCount * 2, yCount * 2, -zCount * 2);\n\n return [light1, light2];\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":784,"value":"exporter = new THREE.OBJExporter()","pinned":false,"hidden":false,"mode":"ojs"},{"id":410,"value":"camera = {\n const fov = 45;\n const aspect = width / height;\n const near = 1;\n const far = 1000;\n const camera = new THREE.PerspectiveCamera(fov, aspect, near, far);\n camera.up.set(0, 0, 1);\n camera.position.set(xCount * 1.5, yCount * 1.5, zCount * 1.5);\n camera.lookAt(0, 0, 0);\n return camera;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":412,"value":"scene = {\n const scene = new THREE.Scene();\n scene.background = new THREE.Color(bg);\n scene.add(cubes);\n // scene.add(isoSurface);\n lights.map((light) => scene.add(light));\n\n return scene;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":408,"value":"renderer = {\n const renderer = new THREE.WebGLRenderer({ antialias: true });\n renderer.setSize(width, height);\n renderer.setPixelRatio(devicePixelRatio);\n const controls = new THREE.OrbitControls(camera, renderer.domElement);\n controls.addEventListener(\"change\", () => renderer.render(scene, camera));\n invalidation.then(() => (controls.dispose(), renderer.dispose()));\n return renderer;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":406,"value":"THREE = {\n const THREE = (window.THREE =\n await require(\"three@0.130.0/build/three.min.js\"));\n await require(\"three@0.130.0/examples/js/controls/OrbitControls.js\").catch(\n () => {}\n );\n await require(\"three@0.130.0/examples/js/exporters/OBJExporter.js\").catch(\n () => {}\n );\n await require(\"three@0.130.0/examples/js/geometries/ConvexGeometry.js\").catch(\n () => {}\n );\n await require(\"three@0.130.0/examples/js/math/ConvexHull.js\").catch(() => {});\n\n return THREE;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1390,"value":"// isoSurface = {\n// const buffGeom = new THREE.BufferGeometry();\n\n// const vertices = new Float32Array(iso.positions.flat(Infinity));\n\n// buffGeom.setAttribute(\"position\", new THREE.BufferAttribute(vertices, 3));\n\n// var indices = new Float32Array(iso.cells.flat(Infinity));\n\n// buffGeom.setIndex(new THREE.BufferAttribute(indices, 3));\n\n// return new THREE.Mesh(buffGeom, mat);\n// }","pinned":true,"hidden":false,"mode":"ojs"},{"id":1275,"value":"// iso = {\n// return isoS.marchingCubes(\n// [xCount, yCount, zCount],\n// function (x, y, z) {\n// return !Number(grid[z * (yCount * xCount) + y * xCount + x].alive > 0);\n// },\n// [\n// [0, 0, 0],\n// [xCount, yCount, zCount]\n// ]\n// );\n// }\n","pinned":true,"hidden":false,"mode":"ojs"},{"id":1271,"value":"isoS = import(\"https://cdn.skypack.dev/isosurface@1.0.0?min\")","pinned":true,"hidden":false,"mode":"ojs"}],"files":"$W33","stdlib":"1","titling":"auto"}