seb/XM_40017
1
0
Fork 0
mirror of https://github.com/fverdugo/XM_40017.git synced 2025-11-08 22:24:25 +01:00
Code Issues Actions Packages Projects Releases Wiki Activity

Fixing missing loop over k in LEQ

Browse Source
This commit is contained in:
Francesc Verdugo 2024-09-27 13:24:49 +02:00
parent 53ea3f5157
commit 3a82bf996d
1 changed files with 20 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

36
notebooks/LEQ.ipynb
View File

@ -340,15 +340,17 @@
"\n", "\n",
"```julia\n", "```julia\n",
"n,m = size(B)\n", "n,m = size(B)\n",
"for t in (k+1):m\n", "for k in 1:n\n",
" B[k,t] = B[k,t]/B[k,k]\n", " for t in (k+1):m\n",
"end\n", " B[k,t] = B[k,t]/B[k,k]\n",
"B[k,k] = 1\n", " end\n",
"for i in (k+1):n \n", " B[k,k] = 1\n",
" for j in (k+1):m\n", " for i in (k+1):n \n",
" B[i,j] = B[i,j] - B[i,k]*B[k,j]\n", " for j in (k+1):m\n",
" B[i,j] = B[i,j] - B[i,k]*B[k,j]\n",
" end\n",
" B[i,k] = 0\n",
" end\n", " end\n",
" B[i,k] = 0\n",
"end\n", "end\n",
"```" "```"
] ]
@ -527,15 +529,17 @@
"\n", "\n",
"```julia\n", "```julia\n",
"n,m = size(B)\n", "n,m = size(B)\n",
"for t in (k+1):m\n", "for k in 1:n\n",
" B[k,t] = B[k,t]/B[k,k]\n", " for t in (k+1):m\n",
"end\n", " B[k,t] = B[k,t]/B[k,k]\n",
"B[k,k] = 1\n", " end\n",
"for i in (k+1):n \n", " B[k,k] = 1\n",
" for j in (k+1):m\n", " for i in (k+1):n \n",
" B[i,j] = B[i,j] - B[i,k]*B[k,j]\n", " for j in (k+1):m\n",
" B[i,j] = B[i,j] - B[i,k]*B[k,j]\n",
" end\n",
" B[i,k] = 0\n",
" end\n", " end\n",
" B[i,k] = 0\n",
"end\n", "end\n",
"```\n", "```\n",
"\n", "\n",