diff options
author | Elizabeth Hunt <elizabeth.hunt@simponic.xyz> | 2023-11-15 14:43:22 -0700 |
---|---|---|
committer | Elizabeth Hunt <elizabeth.hunt@simponic.xyz> | 2023-11-15 14:43:22 -0700 |
commit | 4d2d4f5d7a9f584061b85ef707b3cccee814cbc6 (patch) | |
tree | f588a1cf360744b1f319b548f3ff3823aa3894f2 /src/eigen.c | |
parent | c02573cb66847abfcb9c679416cf14d373790792 (diff) | |
download | cmath-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.c | 6 |
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; } |