summaryrefslogtreecommitdiff
path: root/notes/Nov-27.org
diff options
context:
space:
mode:
authorElizabeth Hunt <elizabeth.hunt@simponic.xyz>2023-11-27 10:14:50 -0700
committerElizabeth Hunt <elizabeth.hunt@simponic.xyz>2023-11-27 10:14:50 -0700
commitaa77d733b0cb100b63a824c60548f19526711bbc (patch)
tree3f65e1e3381afd5dc43dcdc26317cd7085ca4061 /notes/Nov-27.org
parent4d2d4f5d7a9f584061b85ef707b3cccee814cbc6 (diff)
downloadcmath-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.org18
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