fix: Use reproducible_dmatrix() for Cholesky benches

Random matrices may be not positive-definite and Cholesky decomposition
benchmarks panic because of that:

	Benchmarking cholesky_decompose_unpack_100x100: Warming up for 3.0000 s
	thread 'main' panicked at benches/linalg/cholesky.rs:38:45:
	called `Option::unwrap()` on a `None` value

Author: im-0

benches/linalg/cholesky.rs

@@ -1,10 +1,10 @@
-use na::{Cholesky, DMatrix, DVector};
+use na::{Cholesky, DVector};
 
 fn cholesky_100x100(bh: &mut criterion::Criterion) {
     bh.bench_function("cholesky_100x100", |bh| {
         bh.iter_batched(
             || {
-                let m = DMatrix::<f64>::new_random(100, 100);
+                let m = crate::reproducible_dmatrix(100, 100);
                 &m * m.transpose()
             },
             |m| Cholesky::new(m),
@@ -17,7 +17,7 @@ fn cholesky_500x500(bh: &mut criterion::Criterion) {
     bh.bench_function("cholesky_500x500", |bh| {
         bh.iter_batched(
             || {
-                let m = DMatrix::<f64>::new_random(500, 500);
+                let m = crate::reproducible_dmatrix(500, 500);
                 &m * m.transpose()
             },
             |m| Cholesky::new(m),
@@ -31,7 +31,7 @@ fn cholesky_decompose_unpack_100x100(bh: &mut criterion::Criterion) {
     bh.bench_function("cholesky_decompose_unpack_100x100", |bh| {
         bh.iter_batched(
             || {
-                let m = DMatrix::<f64>::new_random(100, 100);
+                let m = crate::reproducible_dmatrix(100, 100);
                 &m * m.transpose()
             },
             |m| {
@@ -46,7 +46,7 @@ fn cholesky_decompose_unpack_500x500(bh: &mut criterion::Criterion) {
     bh.bench_function("cholesky_decompose_unpack_500x500", |bh| {
         bh.iter_batched(
             || {
-                let m = DMatrix::<f64>::new_random(500, 500);
+                let m = crate::reproducible_dmatrix(500, 500);
                 &m * m.transpose()
             },
             |m| {
@@ -62,7 +62,7 @@ fn cholesky_solve_10x10(bh: &mut criterion::Criterion) {
     bh.bench_function("cholesky_solve_10x10", |bh| {
         bh.iter_batched_ref(
             || {
-                let m = DMatrix::<f64>::new_random(10, 10);
+                let m = crate::reproducible_dmatrix(10, 10);
                 let m = &m * m.transpose();
                 let v = DVector::<f64>::new_random(10);
                 let chol = Cholesky::new(m).unwrap();
@@ -78,7 +78,7 @@ fn cholesky_solve_100x100(bh: &mut criterion::Criterion) {
     bh.bench_function("cholesky_solve_100x100", |bh| {
         bh.iter_batched_ref(
             || {
-                let m = DMatrix::<f64>::new_random(100, 100);
+                let m = crate::reproducible_dmatrix(100, 100);
                 let m = &m * m.transpose();
                 let v = DVector::<f64>::new_random(100);
                 let chol = Cholesky::new(m).unwrap();
@@ -94,7 +94,7 @@ fn cholesky_solve_500x500(bh: &mut criterion::Criterion) {
     bh.bench_function("cholesky_solve_500x500", |bh| {
         bh.iter_batched_ref(
             || {
-                let m = DMatrix::<f64>::new_random(500, 500);
+                let m = crate::reproducible_dmatrix(500, 500);
                 let m = &m * m.transpose();
                 let v = DVector::<f64>::new_random(500);
                 let chol = Cholesky::new(m).unwrap();
@@ -110,7 +110,7 @@ fn cholesky_inverse_10x10(bh: &mut criterion::Criterion) {
     bh.bench_function("cholesky_inverse_10x10", |bh| {
         bh.iter_batched_ref(
             || {
-                let m = DMatrix::<f64>::new_random(10, 10);
+                let m = crate::reproducible_dmatrix(10, 10);
                 let m = &m * m.transpose();
                 Cholesky::new(m).unwrap()
             },
@@ -124,7 +124,7 @@ fn cholesky_inverse_100x100(bh: &mut criterion::Criterion) {
     bh.bench_function("cholesky_inverse_100x100", |bh| {
         bh.iter_batched_ref(
             || {
-                let m = DMatrix::<f64>::new_random(100, 100);
+                let m = crate::reproducible_dmatrix(100, 100);
                 let m = &m * m.transpose();
                 Cholesky::new(m).unwrap()
             },
@@ -138,7 +138,7 @@ fn cholesky_inverse_500x500(bh: &mut criterion::Criterion) {
     bh.bench_function("cholesky_inverse_500x500", |bh| {
         bh.iter_batched_ref(
             || {
-                let m = DMatrix::<f64>::new_random(500, 500);
+                let m = crate::reproducible_dmatrix(500, 500);
                 let m = &m * m.transpose();
                 Cholesky::new(m).unwrap()
             },