2c:null 2d:["$","title",null,{"children":"Grid Cells | Matthew Turk | Observable"}] 32:T7d3,kelvinHelmholtzAMR = { var kh_amr = await d3.json("https://gist.githubusercontent.com/matthewturk/660402b5a92aea58b3fc1ea0b59df06b/raw/e18376232b4dbdcec91fc29c7986d9d47115a75e/kh_amr.json"); var complexConfig = defaultConfig([32, 32], [ [ [0, 512.0], [ 20.0 + width/2, 200.0 ] ], [ [0, 512.0], [ 0.0, Math.PI ] ] ] ); var cellData0 = createCells(complexConfig, 0, [32, 32], [0, 0], 0); var cellData1 = createCells(complexConfig, 1, [64, 32], [0, 32], 0); var cellData2 = createCells(complexConfig, 2, [64, 50], [32, 64], 0); cellData0.forEach( d => { d.density = kh[d.j][d.i]; d.level = 0; } ); cellData1.forEach( d => { d.density = kh_amr[0]["density"][d.i][d.j + 32]; d.velocity_x = kh_amr[0]["velocity_x"][d.i][d.j + 32]; d.velocity_y = kh_amr[0]["velocity_y"][d.i][d.j + 32]; d.level = 1; }); cellData2.forEach( d => { d.density = kh_amr[1]["density"][d.i + 32][d.j + 50]; d.velocity_x = kh_amr[1]["velocity_x"][d.i + 32][d.j + 50]; d.velocity_y = kh_amr[1]["velocity_y"][d.i + 32][d.j + 50]; d.level = 2; }); // Remove overlapping data cellData0 = cellData0.filter( c => { var ind = [c.i * 2, c.j * 2]; var overlaps = ((ind[0] >= 0) && (ind[0] < 64) && (ind[1] >= 32) && (ind[1] < 64)); return !overlaps; }); cellData1 = cellData1.filter( c => { var ind = [c.i * 2, c.j * 2]; var overlaps = ((ind[0] >= 32) && (ind[0] < 96) && (ind[1] >= 64) && (ind[1] < 114)); return !overlaps; }); var cellData = cellData0.concat(cellData1).concat(cellData2); var cellPaths = createPaths(complexConfig, cellData, 10); var tree = d3.quadtree().x(d=>d.x[0] + d.width[0]/2).y(d=>d.y[0] + d.height[0]/2).addAll(cellData); return {'kh_amr': kh_amr, 'cellData': cellData, 'cellPaths': cellPaths, 'tree': tree}; }33:[] 30:{"id":"5b56258d81fbd95b","slug":"grid-cells","title":"Grid Cells","thumbnail":"0b4e03e5fe40601163c2e3968bb7f781173c8ccfb3ad5d1b3d436d286ae6187d","access_level":"public","create_time":"$D2019-10-14T16:19:02.054Z","update_time":"$D2019-10-26T18:56:53.239Z","view_time":null,"first_public_time":"$D2019-10-21T17:44:20.517Z","has_unpublished_edits":null,"fork":false,"stars":1,"starred":false,"owner":{"type":"team","id":"1fcd1ee6fffb0cba","login":"matthewturk","name":"Matthew Turk","avatar_url":"https://avatars-next.observableusercontent.com/221c00aeb2f6f114ca9960fbf8f060d9ba33c911e5f64cea2d589ea2d34d7e94","tier":"starter_2024"},"version":697,"publish_version":697,"role":null,"authored":false,"authors":[{"type":"individual","id":"294f57c680999965","login":"matthewturk","name":"Matthew Turk","avatar_url":"https://avatars-next.observableusercontent.com/221c00aeb2f6f114ca9960fbf8f060d9ba33c911e5f64cea2d589ea2d34d7e94","tier":null}],"members":[{"type":"individual","id":"294f57c680999965","login":"matthewturk","name":"Matthew Turk","avatar_url":"https://avatars-next.observableusercontent.com/221c00aeb2f6f114ca9960fbf8f060d9ba33c911e5f64cea2d589ea2d34d7e94","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":1,"collection":null,"type":"notebook","body":{"titling":"auto","title":"Grid Cells","cells":[{"id":0,"value":"md`# Grid Cells`","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":188,"value":"function createPaths(globalConfig, cellData, Nrow = 9) {\n var cellsDisk = cellData.map( (d, flatIndex) => d3rect.rect(\n {x: (d.flatIndex % Nrow) * globalConfig.dds[0][0] + globalConfig.bounds[0][0][0] + d.groupOffset,\n y: Math.floor(d.flatIndex / Nrow) * globalConfig.dds[1][0] + globalConfig.bounds[1][0][0],\n width: globalConfig.dds[0][0],\n height: globalConfig.dds[1][0]}) );\n var cellsMemory = cellData.map( d => d3rect.rect(\n {x: d.x[0], y: d.y[0], width: d.width[0], height: d.height[0]}));\n var cellsCyl = cellData.map( d => d3.arc()({\n innerRadius: d.x[1],\n outerRadius: d.x[1] + d.width[1],\n startAngle: d.y[1],\n endAngle: d.y[1] + d.height[1],\n padAngle: 0.0}) );\n return {'disk': cellsDisk, 'memory': cellsMemory, 'cyl': cellsCyl};\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":459,"value":"function defaultConfig(Ndims, bounds) {\n // first dimension of bounds is the dimension, second is the variable\n // so bounds[0] corresponds to i, bounds[1] to j, bounds[0][0] is x\n // bounds[0][1] is r, bounds[0][0][0] is left, bounds[0][0][1] is right\n var cfg = {Ndims: Ndims, width: [], dds: [], bounds: bounds};\n Ndims.forEach( (nd, i) => {\n cfg.width.push( bounds[i].map( v => (v[1] - v[0]) ) );\n cfg.dds.push( bounds[i].map( v => (v[1] - v[0])/nd ) );\n });\n return cfg;\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":484,"value":"cfg = defaultConfig([6, 8], [ [ [-120.0, 120.0], [ 20.0, 200.0 ] ],\n [ [-160.0, 160.0], [ 0.0, Math.PI ] ] ] );","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":461,"value":"function createCells(globalConfig, level, dims, levelIndex, groupOffset = 0) {\n var dds = globalConfig.dds.map( v1 => v1.map( v2 => (v2 / (2**level)) ) );\n var thisLeft = globalConfig.bounds.map( (v1, i) => v1.map( (v2, j) => (v2[0] + dds[i][j] * levelIndex[i] ) ));\n var data = [];\n d3.range(dims[1]).forEach(j => d3.range(dims[0]).forEach( i=> {\n data.push({\n flatIndex: i * dims[1] + j,\n i: i,\n j: j,\n x: thisLeft[0].map( (l, n) => (i * dds[0][n] + l)),\n y: thisLeft[1].map( (l, n) => (j * dds[1][n] + l)),\n width: dds[0],\n height: dds[1],\n groupOffset: groupOffset\n });\n }));\n return data;\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":417,"value":"rootGrid = createCells(cfg, 0, [6, 8], [0, 0], -220);\n//generateCellData();","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":537,"value":"refinedData = {\n var refinedData = [];\n var refinedData = refinedData.concat(createCells(cfg, 1, [4, 8], [4, 2], 160));\n return refinedData;\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":453,"value":"function ringermacherMead(A, B, N) {\n return function(phi) {\n return A / Math.log(B * Math.tan(phi / (2*N)));\n }\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":692,"value":"kh = d3.json(\"https://gist.githubusercontent.com/matthewturk/660402b5a92aea58b3fc1ea0b59df06b/raw/edbd8f56f68dc55de8c9c1b220a3b8acf2e0935b/kh.json\");","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":688,"value":"$32","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":690,"value":"galaxy0030 = d3.json(\"https://gist.githubusercontent.com/matthewturk/660402b5a92aea58b3fc1ea0b59df06b/raw/9098fe92388acd00c9c4f9c98140f82f4b34b83d/galaxy0030.json\")","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":694,"value":"function findOverlappingNodesPoint(point, biggestCell, fn) {\n function visitTree(node, x0, y0, x1, y1) {\n if ( x0 > point[0] + biggestCell[0]/2 || x1 < point[0] - biggestCell[0]/2\n || y0 > point[1] + biggestCell[1]/2 || y1 < point[1] - biggestCell[1]/2 ) {\n return true;\n }\n if (!node.length) do {\n var d = node.data;\n if ( (point[0] >= d.x[0]) && (point[0] <= d.x[0] + d.width[0])\n && (point[1] >= d.y[0]) && (point[1] <= d.y[0] + d.height[0]) ) {\n fn(d);\n };\n } while (node = node.next);\n }\n return visitTree;\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":696,"value":"function findOverlappingNodesRect(extent, fn) {\n function visitTree(node, x0, y0, x1, y1) {\n if ( x0 > extent[1][0] || x1 < extent[0][0] || y0 > extent[1][1] || y1 < extent[0][1] ) {\n return true;\n }\n if (!node.length) do {\n var d = node.data;\n if ( extent[0][0] < d.x[0] + d.width[0]/2 && extent[0][1] < d.y[0] + d.height[0]/2\n && extent[1][0] > d.x[0] + d.width[0]/2 && extent[1][1] > d.y[0] + d.height[0]/2 ) {\n fn(d);\n };\n } while (node = node.next);\n }\n return visitTree;\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":4,"value":"d3 = require(\"d3@5\");","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":226,"value":"d3rect = require(\"d3-rect\");","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":205,"value":"flubber = require(\"flubber\");","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"}],"files":"$W33","stdlib":"1"}} 31:"$30:body"