Fix upper bound for diagonal elements in Jacobi method

Andy Oertel · Andy Oertel · commit 7893f6da57ef · 2021-12-29T13:12:35.000+01:00
diff --git a/benches/jacobi_method.rs b/benches/jacobi_method.rs
@@ -67,7 +67,7 @@ fn benchmark_jacobi(_c: &mut Criterion) {
                         diag_elem += matrix[(i, j)];
                     }
                 }
-                diag_elem = rng.gen_range((1.0 + diag_elem)..(diag_elem * diag_elem));
+                diag_elem = rng.gen_range((1.0 + diag_elem)..((diag_elem * diag_elem) + 1.0));
                 diag[i] = diag_elem;
                 matrix[(i, i)] = diag_elem;
                 off_diag[(i, i)] = 0.;

Original file line number	Diff line number	Diff line change
`@@ -67,7 +67,7 @@ fn benchmark_jacobi(_c: &mut Criterion) {`
`67`	`67`	`diag_elem += matrix[(i, j)];`
`68`	`68`	`}`
`69`	`69`	`}`
`70`		`- diag_elem = rng.gen_range((1.0 + diag_elem)..(diag_elem * diag_elem));`
	`70`	`+ diag_elem = rng.gen_range((1.0 + diag_elem)..((diag_elem * diag_elem) + 1.0));`
`71`	`71`	`diag[i] = diag_elem;`
`72`	`72`	`matrix[(i, i)] = diag_elem;`
`73`	`73`	`off_diag[(i, i)] = 0.;`