Ver 0.39.4

- Optimize `integrate`

Author: Axect

Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "peroxide"
-version = "0.39.3"
+version = "0.39.4"
 authors = ["axect <[email protected]>"]
 edition = "2018"
 description = "Rust comprehensive scientific computation library contains linear algebra, numerical analysis, statistics and machine learning tools with farmiliar syntax"
@@ -29,8 +29,8 @@ criterion = { version = "0.5.1", features = ["html_reports"] }
 
 [dependencies]
 csv = { version = "1.3", optional = true, default-features = false }
-rand = { version = "0.8", features = ["small_rng"] }
-rand_distr = "0.4"
+rand = { version = "0.9", features = ["small_rng"] }
+rand_distr = "0.5"
 order-stat = "0.1"
 puruspe = "0.4"
 matrixmultiply = { version = "0.3", features = ["threading"] }

RELEASES.md

@@ -1,3 +1,19 @@
+# Release 0.39.4 (2025-04-11)
+
+## Optimize `integrate`
+
+- Replace the output signature of `gauss_legendre_table` and `kronrod_table` to `&'static [f64]` to avoid unnecessary allocations.
+- Hard code symmetry of weights and nodes into source code to avoid unnecessary allocations.
+- New helper function - `compute_gauss_kronrod_sum_stored`
+  - Reduce the number of function calls (G+K -> K)
+- Change update method of subinterval tolerance (divide by 2 -> divide by sqrt(2))
+- These changes improve the performance of `integrate` by 1.2x - 50x (to integrate highly oscillatory functions)
+
+## Update dependencies
+
+- Update `rand` to `0.9`
+- Update `rand_distr` to `0.5`
+
 # Release 0.39.3 (2025-03-13)
 
 - Update `puruspe` to `0.4.0`

src/fuga/mod.rs

@@ -205,6 +205,7 @@ pub use crate::util::plot::*;
 
 pub use anyhow;
 pub use paste;
+pub use rand::prelude::*;
 
 // =============================================================================
 // Enums

src/macros/matlab_macro.rs

@@ -1,6 +1,3 @@
-extern crate rand;
-pub use self::rand::prelude::*;
-
 /// MATLAB like zeros - zero matrix
 ///
 /// # Examples
@@ -45,18 +42,18 @@ macro_rules! zeros {
 #[macro_export]
 macro_rules! rand {
     () => {{
-        let mut rng = thread_rng();
-        rng.gen_range(0f64..=1f64)
+        let mut rng = rand::rng();
+        rng.random_range(0f64..=1f64)
     }};
 
     ( $m:expr, $n:expr ) => {{
         let r = $m;
         let c = $n;
-        let mut rng = thread_rng();
+        let mut rng = rand::rng();
         let mut m = matrix(vec![0f64; r * c], r, c, Row);
         for i in 0..r {
             for j in 0..c {
-                m[(i, j)] = rng.gen_range(0f64..=1f64);
+                m[(i, j)] = rng.random_range(0f64..=1f64);
             }
         }
         m

src/macros/r_macro.rs

@@ -16,9 +16,6 @@
 //! * `pt`
 //! * `lm`
 
-extern crate rand;
-pub use self::rand::prelude::*;
-
 /// R like concatenate (Type: Vec\<f64\>)
 ///
 /// # Examples
@@ -278,10 +275,10 @@ macro_rules! runif {
         let n: usize = $x0;
         let mut v = vec![0f64; n];
 
-        let mut rng = thread_rng();
+        let mut rng = rand::rng();
 
         for i in 0..n {
-            v[i] = rng.gen_range($start as f64..=$end as f64);
+            v[i] = rng.random_range($start as f64..=$end as f64);
         }
         v
     }};