diff options
Diffstat (limited to 'notes/Nov-27.org')
-rw-r--r-- | notes/Nov-27.org | 18 |
1 files changed, 18 insertions, 0 deletions
diff --git a/notes/Nov-27.org b/notes/Nov-27.org new file mode 100644 index 0000000..c4d38b0 --- /dev/null +++ b/notes/Nov-27.org @@ -0,0 +1,18 @@ +x^{k+1} = D^{-1}(b - (L + U) x^k) +x^{k + 1} \rightarrow Ax^k + + +#+BEGIN_SRC c + loop while (err > tol && iter < maxiter) { + for (int i = 0; i < n; i++) { + sum = b[i]; + for (int j = 0; j < i; j++) { + sum = sum - a[i][x] * x_0[i]; + } + for (int j = i; j < n; j++) { + sum = sum + a[i][j] * x_0[j]; + } + x_1[i] = sum / a[i][i]; + } + } +#+END_SRC |