diff options
author | Elizabeth Hunt <elizabeth.hunt@simponic.xyz> | 2023-11-27 10:14:50 -0700 |
---|---|---|
committer | Elizabeth Hunt <elizabeth.hunt@simponic.xyz> | 2023-11-27 10:14:50 -0700 |
commit | aa77d733b0cb100b63a824c60548f19526711bbc (patch) | |
tree | 3f65e1e3381afd5dc43dcdc26317cd7085ca4061 /notes/Nov-27.org | |
parent | 4d2d4f5d7a9f584061b85ef707b3cccee814cbc6 (diff) | |
download | cmath-aa77d733b0cb100b63a824c60548f19526711bbc.tar.gz cmath-aa77d733b0cb100b63a824c60548f19526711bbc.zip |
add leslie matrix dominant eigenvalue unit test
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 |