From aa77d733b0cb100b63a824c60548f19526711bbc Mon Sep 17 00:00:00 2001 From: Elizabeth Hunt Date: Mon, 27 Nov 2023 10:14:50 -0700 Subject: add leslie matrix dominant eigenvalue unit test --- notes/Nov-27.org | 18 ++++++++++++++++++ 1 file changed, 18 insertions(+) create mode 100644 notes/Nov-27.org (limited to 'notes/Nov-27.org') 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 -- cgit v1.2.3-70-g09d2