2d:null 2e:["$","title",null,{"children":"RREF | Vinícius Mello | Observable"}] 33:Tfa3,{ let mcell = []; let m = M.length; let n = M[0].length; let sz = 60; let matrix = html`

`; let { $, parse, Eval, equal, Integer, Plus, Times, Divide, copy } = Guacyra.Kernel; let { latex } = Guacyra.Formatting; let A = copy($`Table(0,[i,1,${m}],[j,1,${n}])`); for (let i = 1; i <= m; ++i) { for (let j = 1; j <= n; ++j) { A[i][j] = Eval(parse(`${M[i - 1][j - 1]}`)); } } for (let i = 1; i <= m; ++i) { let mrow = []; mcell.push(mrow); for (let j = 1; j <= n; ++j) { let cell = html`

`; cell.appendChild(tex`\sf ${latex(A[i][j])}`); mrow.push(cell); matrix.appendChild(cell); } } const log = html`

`; matrix.appendChild(log); const appendLog = text => { log.appendChild(html`

${text}

`); log.scrollTop = log.scrollHeight; }; const setClass = (el, cl) => { el.classList.remove(el.classList.value.split(' ').pop()); el.classList.add(cl); }; const getClass = el => { return el.classList.value.split(' ').pop(); }; let pivot = mcell[0][0]; let curr = mcell[0][0]; yield matrix; await Promises.delay(mutable delay); let zero = (i, j) => equal(A[i][j], Integer(0)); for (let s of RowReduce(m, n, zero)) { if (s.action === 'swap') { let i = s.i; let ip = s.ip; if (i == ip) continue; appendLog(tex`\sf l_{${ip}}\leftrightarrow l_{${i}}`); await Promises.delay(mutable delay / 2); anime({ targets: mcell[i - 1], top: sz * (ip - 1) }); anime({ targets: mcell[ip - 1], top: sz * (i - 1) }); [A[i], A[ip]] = [A[ip], A[i]]; [mcell[i - 1], mcell[ip - 1]] = [mcell[ip - 1], mcell[i - 1]]; } if (s.action === 'pivot') { let j = s.j; let i = s.i; pivot = mcell[i - 1][j - 1]; setClass(pivot, 'pivotcell'); } if (s.action === 'search') { let j = s.j; let i = s.i; let cl = getClass(curr); if (cl === 'searchcell') setClass(curr, 'normalcell'); curr = mcell[i - 1][j - 1]; setClass(curr, 'searchcell'); } if (s.action === 'scale') { let jj = s.j; let i = s.i; const k = Eval(Divide(1, A[i][jj])); if (equal(k, Integer(1))) continue; const tt = latex($`${k}*aaa`).replace('aaa', `l_{${i}}`); appendLog(tex`\sf l_{${i}}\rightarrow ${tt}`); for (let j = jj; j <= n; ++j) { A[i][j] = Eval(Times(A[i][j], k)); let cell = mcell[i - 1][j - 1]; let cl = getClass(cell); setClass(cell, 'compcell'); await Promises.delay(mutable delay / 2); cell.innerHTML = ''; cell.appendChild(tex`\sf ${latex(A[i][j])}`); await Promises.delay(mutable delay / 2); setClass(cell, cl); } } if (s.action === 'eliminate') { let jp = s.j; let i = s.i; let ip = s.ip; if (i == ip) continue; const k = Eval(Times(-1, Divide(A[i][jp], A[ip][jp]))); if (equal(k, Integer(0))) continue; const tt = latex($`aaa+${k}*bbb`) .replace('aaa', `l_{${i}}`) .replace('bbb', `l_{${ip}}`); appendLog(tex`\sf l_{${i}}\rightarrow ${tt}`); for (let j = jp; j <= n; ++j) { A[i][j] = Eval(Plus(A[i][j], Times(k, A[ip][j]))); let cell = mcell[i - 1][j - 1]; setClass(cell, 'compcell'); await Promises.delay(mutable delay / 2); cell.innerHTML = ''; cell.appendChild(tex`\sf ${latex(A[i][j])}`); await Promises.delay(mutable delay / 2); setClass(cell, 'normalcell'); } } await Promises.delay(mutable delay); yield matrix; } await Promises.delay(mutable delay); yield matrix; }34:[] 31:{"id":"31212eba05bbaacf","slug":"rref","title":"RREF","thumbnail":"5a95a3226b2a17d239e86beca4ca4126febe9de91973b864bdf478c5c0c7d08c","access_level":"public","create_time":"$D2021-08-18T20:58:27.958Z","update_time":"$D2022-11-09T12:10:34.396Z","view_time":null,"first_public_time":"$D2021-08-18T21:53:51.174Z","has_unpublished_edits":null,"fork":{"id":"2706323c94052c67","slug":"gaussian-elimination","title":"Gaussian Elimination","owner":{"id":"7a4ecd430a93cd6f","login":"vinicius-mello","name":"Vinícius Mello"}},"stars":0,"starred":false,"owner":{"type":"team","id":"7a4ecd430a93cd6f","login":"vinicius-mello","name":"Vinícius Mello","avatar_url":"https://avatars-next.observableusercontent.com/b5d4671aa488d8bcb984f575f985332403cde5f63942ea883a6a1a630941ec05","tier":"starter_2024"},"version":580,"publish_version":null,"role":null,"authored":false,"authors":[{"type":"individual","id":"ecc4854ab7ad95b6","login":"vinicius-mello","name":"Vinícius Mello","avatar_url":"https://avatars-next.observableusercontent.com/b5d4671aa488d8bcb984f575f985332403cde5f63942ea883a6a1a630941ec05","tier":null}],"members":[{"type":"individual","id":"ecc4854ab7ad95b6","login":"vinicius-mello","name":"Vinícius Mello","avatar_url":"https://avatars-next.observableusercontent.com/b5d4671aa488d8bcb984f575f985332403cde5f63942ea883a6a1a630941ec05","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":2,"collection":null,"type":"notebook","body":{"title":"RREF","cells":[{"id":0,"value":"md`# RREF`","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":14,"value":"M = [\n [1, 0, 0, 2],\n [0, 1, 1, 1],\n [1, 1, -1, 1]\n]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":141,"value":"$33","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":169,"value":"mutable delay = 200","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":329,"value":"speedControl = {\n var input = DOM.range(50, 2000, 50);\n input.value = mutable delay;\n input.oninput = () => (mutable delay = input.valueAsNumber);\n return input;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":128,"value":"function* RowReduce(m, n, zero) {\n let i;\n let j;\n let ip;\n let sm = {\n init: function() {\n i = 1;\n j = 1;\n ip = 1;\n return 'search';\n },\n search: function() {\n if (zero(i, j)) {\n i = i + 1;\n if (i > m) {\n i = ip;\n j = j + 1;\n }\n if (j > n) return 'stop';\n return 'search';\n } else {\n return 'pivot';\n }\n },\n pivot: function() {\n return 'scale';\n },\n scale: function() {\n return 'swap';\n },\n swap: function() {\n i = 1;\n return 'eliminate';\n },\n eliminate: function() {\n i = i + 1;\n if (i > m) {\n j = j + 1;\n if (j > n) return 'stop';\n ip = ip + 1;\n i = ip;\n if (i > m) return 'stop';\n return 'search';\n } else return 'eliminate';\n }\n };\n let st = 'init';\n do {\n st = sm[st]();\n yield { action: st, i, j, ip };\n } while (st !== 'stop');\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":254,"value":"html``","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":40,"value":"anime = require('https://unpkg.com/animejs/lib/anime.min.js')","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":8,"value":"Guacyra = require('guacyra')","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"}],"files":"$W34","stdlib":"1","titling":"auto"}} 32:"$31:body"