From 3a82bf996d31b16cbe6b253d6bde1de97867b148 Mon Sep 17 00:00:00 2001 From: Francesc Verdugo Date: Fri, 27 Sep 2024 13:24:49 +0200 Subject: [PATCH] Fixing missing loop over k in LEQ --- notebooks/LEQ.ipynb | 36 ++++++++++++++++++++---------------- 1 file changed, 20 insertions(+), 16 deletions(-) diff --git a/notebooks/LEQ.ipynb b/notebooks/LEQ.ipynb index b2e951c..a393385 100644 --- a/notebooks/LEQ.ipynb +++ b/notebooks/LEQ.ipynb @@ -340,15 +340,17 @@ "\n", "```julia\n", "n,m = size(B)\n", - "for t in (k+1):m\n", - " B[k,t] = B[k,t]/B[k,k]\n", - "end\n", - "B[k,k] = 1\n", - "for i in (k+1):n \n", - " for j in (k+1):m\n", - " B[i,j] = B[i,j] - B[i,k]*B[k,j]\n", + "for k in 1:n\n", + " for t in (k+1):m\n", + " B[k,t] = B[k,t]/B[k,k]\n", + " end\n", + " B[k,k] = 1\n", + " for i in (k+1):n \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", - " B[i,k] = 0\n", "end\n", "```" ] @@ -527,15 +529,17 @@ "\n", "```julia\n", "n,m = size(B)\n", - "for t in (k+1):m\n", - " B[k,t] = B[k,t]/B[k,k]\n", - "end\n", - "B[k,k] = 1\n", - "for i in (k+1):n \n", - " for j in (k+1):m\n", - " B[i,j] = B[i,j] - B[i,k]*B[k,j]\n", + "for k in 1:n\n", + " for t in (k+1):m\n", + " B[k,t] = B[k,t]/B[k,k]\n", + " end\n", + " B[k,k] = 1\n", + " for i in (k+1):n \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", - " B[i,k] = 0\n", "end\n", "```\n", "\n",