30:T4ee,bellman = (from, to) => { const nodesCount = _(graph.nodes).map(n => n['Id epsico']).max() const dist = new Float32Array(nodesCount) const pred = new Uint16Array(nodesCount) const duration = new Float32Array(nodesCount) for (let i = 0; i < nodesCount; i++) { dist[i] = Infinity pred[i] = i duration[i] = Infinity } dist[from] = 0 duration[from] = 0 for (const _node of Object.entries(graph.nodes)) { for (const edge of graph.edges) { const edge_duration = edge.longueur * 60 / edge.vmax if (dist[edge.ido] + edge.longueur < dist[edge.idf]) { pred[edge.idf] = edge.ido dist[edge.idf] = dist[edge.ido] + edge.longueur duration[edge.idf] = duration[edge.ido] + edge_duration } // Edges are symetrical if (dist[edge.idf] + edge.longueur < dist[edge.ido]) { pred[edge.ido] = edge.idf dist[edge.ido] = dist[edge.idf] + edge.longueur duration[edge.ido] = duration[edge.idf] + edge_duration } } } const result = [] let current = to do { result.push(current) current = pred[current] } while (pred[current] != current) return { nodes: result.reverse(), distance: dist[to], duration: duration[to], } }31:[{"id":"e1fb588052b52626891024d699f6c08bae67eabb208f3e3aa61194858681a77abfa314fc1a0107de975bfa6d9c3c726799550bbd4530dfbe7694c310411ea6d6","href":"https://static.observableusercontent.com/files/e1fb588052b52626891024d699f6c08bae67eabb208f3e3aa61194858681a77abfa314fc1a0107de975bfa6d9c3c726799550bbd4530dfbe7694c310411ea6d6","name":"noeuds_alpha_pcc.txt","type":"text/plain","lastModified":1576870205177.217,"size":126185},{"id":"ad26fee976a7f10c4b5c5a6f6fc3b050fe114ec0c4f9326bbd62b09281b975ee8039ff786ea3e1c658546ce8daf58a81bbdda4df7c2dba89d043c2be097e6900","href":"https://static.observableusercontent.com/files/ad26fee976a7f10c4b5c5a6f6fc3b050fe114ec0c4f9326bbd62b09281b975ee8039ff786ea3e1c658546ce8daf58a81bbdda4df7c2dba89d043c2be097e6900","name":"arcs_vit_en.csv","type":"text/csv","lastModified":1576764678417.98,"size":174944}] 2e:{"title":"Freepsico","cells":[{"id":0,"value":"md`# Freepsico\n\nUne tentative de refaire Epsico (application windows uniquement).\n\nCela permet d’estimer les péages.\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":223,"value":"viewof paramètres = {\n const form = html`\n

\n\n `)\n }\n Départ : \n Arrivée :
\n Éviter les LGV :
\n Traction électrique : \n\n`\n form.oninput = () => form.value = {\n departId: form.depart.value,\n arriveId: form.arrivee.value,\n lgv: form.lgv,\n }\n return form\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":302,"value":"md`## Synthèse\n* **Départ** ${parsed.from? graph.nodes[parsed.from]['libelle epsico']: 'Indéfini'} \n* **Arrivée** ${parsed.to? graph.nodes[parsed.to]['libelle epsico']: 'Indéfini'} \n* **Longueur** ${route? Math.ceil(route.route.distance): '?'} km\n* **Durée min théorique** ${route? `${Math.floor(route.route.duration / 60)}:${Math.floor(route.route.duration % 60)}`: ''}\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":298,"value":"{\n if(route) {\n return route.canvas\n }\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":306,"value":"parsed = {\n const from = Number(paramètres.departId)\n const to = Number(paramètres.arriveId)\n return {\n from,\n to,\n }\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":276,"value":"route = {\n const from = Number(paramètres.departId)\n const to = Number(paramètres.arriveId)\n const route = bellman(from, to)\n if (from && to) {\n return {\n canvas: drawPath(route.nodes).canvas,\n route,\n }\n }\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":290,"value":"result = Object()","pinned":true,"hidden":false,"mode":"ojs"},{"id":242,"value":"paramètres","pinned":true,"hidden":false,"mode":"ojs"},{"id":2,"value":"graph = {\n const rawNodes = d3.dsvFormat('\\t').parse(await FileAttachment(\"noeuds_alpha_pcc.txt\").text(), d3.autoType)\n const nodes = _(rawNodes).map(n => [n['Id epsico'], n]).fromPairs().value()\n const edges = d3.dsvFormat(';').parse(await FileAttachment(\"arcs_vit_en.csv\").text(), d3.autoType)\n \n return {\n nodes,\n edges,\n minX: _(nodes).map(n => n.RGI_XX).min(),\n minY: _(nodes).map(n => n.RGI_YY).min(),\n maxX: _(nodes).map(n => n.RGI_XX).max(),\n maxY: _(nodes).map(n => n.RGI_YY).max(),\n proj: function(x, y, width) {\n const pixels_per_meter = width / (this.maxX - this.minX);\n return {\n x: (x - this.minX) * pixels_per_meter,\n y: width - (y - this.minY) * pixels_per_meter\n }\n },\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":9,"value":"backgroundMap = () => {\n const context = DOM.context2d(width, width);\n for (const [_nodeId, node] of Object.entries(graph.nodes)) {\n const coord = graph.proj(node.RGI_XX, node.RGI_YY, width)\n context.beginPath();\n context.arc(coord.x, coord.y, 2, 0, 2 * Math.PI);\n context.stroke();\n }\n \n for (const edge of graph.edges) {\n const start = graph.nodes[edge.ido]\n const end = graph.nodes[edge.idf]\n const startCoord = graph.proj(start.RGI_XX, start.RGI_YY, width)\n const endCoord = graph.proj(end.RGI_XX, end.RGI_YY, width)\n \n context.beginPath();\n if (edge.nrj == 1) {\n context.strokeStyle = '#008000';\n } else {\n context.strokeStyle = '#800000';\n }\n context.moveTo(startCoord.x, startCoord.y);\n context.lineTo(endCoord.x, endCoord.y);\n context.stroke();\n }\n \n return context;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":6,"value":"d3 = require(\"d3-dsv@1\")","pinned":true,"hidden":false,"mode":"ojs"},{"id":125,"value":"$30","pinned":true,"hidden":false,"mode":"ojs"},{"id":133,"value":"drawPath = path => {\n const context = backgroundMap()\n context.beginPath();\n \n let node = graph.nodes[path[0]]\n let coord = graph.proj(node.RGI_XX, node.RGI_YY, width)\n context.moveTo(coord.x, coord.y)\n \n for (const n of path) {\n const node = graph.nodes[n]\n const coord = graph.proj(node.RGI_XX, node.RGI_YY, width)\n context.lineTo(coord.x, coord.y);\n }\n context.strokeStyle = '#000080';\n context.lineWidth = 3;\n context.stroke();\n \n return context\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":30,"value":"_ = require(\"lodash\")","pinned":true,"hidden":false,"mode":"ojs"}],"files":"$W31","stdlib":"1","titling":"auto"}