From 4d2d4f5d7a9f584061b85ef707b3cccee814cbc6 Mon Sep 17 00:00:00 2001 From: Elizabeth Hunt Date: Wed, 15 Nov 2023 14:43:22 -0700 Subject: add hw 7 and documentation for q1 and part of q2 in lizfcm api reference --- src/eigen.c | 6 ++---- 1 file changed, 2 insertions(+), 4 deletions(-) (limited to 'src/eigen.c') 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; } -- cgit v1.2.3-70-g09d2