summaryrefslogtreecommitdiff
path: root/src/eigen.c
diff options
context:
space:
mode:
authorElizabeth Hunt <elizabeth.hunt@simponic.xyz>2023-11-15 14:43:22 -0700
committerElizabeth Hunt <elizabeth.hunt@simponic.xyz>2023-11-15 14:43:22 -0700
commit4d2d4f5d7a9f584061b85ef707b3cccee814cbc6 (patch)
treef588a1cf360744b1f319b548f3ff3823aa3894f2 /src/eigen.c
parentc02573cb66847abfcb9c679416cf14d373790792 (diff)
downloadcmath-4d2d4f5d7a9f584061b85ef707b3cccee814cbc6.tar.gz
cmath-4d2d4f5d7a9f584061b85ef707b3cccee814cbc6.zip
add hw 7 and documentation for q1 and part of q2 in lizfcm api reference
Diffstat (limited to 'src/eigen.c')
-rw-r--r--src/eigen.c6
1 files changed, 2 insertions, 4 deletions
diff --git a/src/eigen.c b/src/eigen.c
index 6ea26c3..40b41b8 100644
--- a/src/eigen.c
+++ b/src/eigen.c
@@ -14,10 +14,8 @@ Matrix_double *leslie_matrix(Array_double *age_class_surivor_ratio,
free_vector(leslie->data[0]);
leslie->data[0] = age_class_offspring;
- for (size_t i = 0; i < age_class_surivor_ratio->size; i++) {
+ for (size_t i = 0; i < age_class_surivor_ratio->size; i++)
leslie->data[i + 1]->data[i] = age_class_surivor_ratio->data[i];
- }
-
return leslie;
}
@@ -37,9 +35,9 @@ double dominant_eigenvalue(Matrix_double *m, Array_double *v, double tolerance,
Array_double *mx = m_dot_v(m, eigenvector_2);
double new_lambda =
v_dot_v(mx, eigenvector_2) / v_dot_v(eigenvector_2, eigenvector_2);
+
error = fabs(new_lambda - lambda);
lambda = new_lambda;
-
free_vector(eigenvector_1);
eigenvector_1 = eigenvector_2;
}