From aa77d733b0cb100b63a824c60548f19526711bbc Mon Sep 17 00:00:00 2001 From: Elizabeth Hunt Date: Mon, 27 Nov 2023 10:14:50 -0700 Subject: add leslie matrix dominant eigenvalue unit test --- test/eigen.t.c | 24 ++++++++++++++++++++++++ 1 file changed, 24 insertions(+) (limited to 'test') diff --git a/test/eigen.t.c b/test/eigen.t.c index 9b2f6f0..31c3500 100644 --- a/test/eigen.t.c +++ b/test/eigen.t.c @@ -17,6 +17,30 @@ UTEST(eigen, leslie_matrix) { free_matrix(leslie); free_matrix(m); + free_vector(felicity); + free_vector(survivor_ratios); +} + +UTEST(eigen, leslie_matrix_dominant_eigenvalue) { + Array_double *felicity = InitArray(double, {0.0, 1.5, 0.8}); + Array_double *survivor_ratios = InitArray(double, {0.8, 0.55}); + Matrix_double *leslie = leslie_matrix(survivor_ratios, felicity); + Array_double *v_guess = InitArrayWithSize(double, 3, 1.0); + double tolerance = 0.0001; + uint64_t max_iterations = 64; + + double expect_dominant_eigenvalue = 1.22005; + + double approx_dominant_eigenvalue = + dominant_eigenvalue(leslie, v_guess, tolerance, max_iterations); + + EXPECT_NEAR(expect_dominant_eigenvalue, approx_dominant_eigenvalue, + tolerance); + + free_vector(v_guess); + free_vector(survivor_ratios); + free_vector(felicity); + free_matrix(leslie); } UTEST(eigen, dominant_eigenvalue) { -- cgit v1.2.3-70-g09d2