bump version number

Author: kmolan

BENCHMARKS.md

@@ -86,5 +86,5 @@ Integrand                              | Approximation error |           Notes
 -------------------------------------- | ------------------- | ------------------------------------------------------- |
 $$\int_{-\infty}^\infty x^2 e^{-x^2} \mathrm{d}x$$   | 1e-30              |  Trivial Integration to showcase accuracy levels    |
 $$\int_{-\infty}^\infty (4x^3 - 3x^2)e^{-x^2} \mathrm{d}x$$  | 1e-12      |   High accuracy for more complicated integrands      |
-$$\int_{-\infty}^\infty (Sin(x) - \sqrt{x})e^{-x} \mathrm{d}x$$ | 1e-1 | Poor performance for non-polynomial integrands |
+$$\int_{-\infty}^\infty (Sin(x) - \sqrt{x})e^{-x^2} \mathrm{d}x$$ | 1e-1 | Poor performance for non-polynomial integrands |
 

Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "multicalc"
-version = "0.4.1"
+version = "0.5.0"
 authors = ["akathail <[email protected]>"]
 edition = "2021"
 description = "Rust scientific computing for single and multi-variable calculus"
@@ -20,7 +20,7 @@ num-complex = "0.4.6"
 rand = "0.8"
 
 [features]
-default = ["heap"]
+default = []
 heap = []
 
 # See more keys and their definitions at https://doc.rust-lang.org/cargo/reference/manifest.html

src/numerical_integration/test.rs

@@ -20,7 +20,6 @@ fn test_booles_integration_1()
 
     //simple integration for x, known to be x*x, expect a value of ~4.00
     let val = integrator.get_single(&func, &integration_limit).unwrap();
-    std::println!("{}", val - 4.0);
     assert!(f64::abs(val - 4.0) < 1e-14);
 }
 
@@ -54,7 +53,6 @@ fn test_booles_integration_2()
 
     //partial integration for z, known to be 2.0*x*z + y*z*z/2.0, expect a value of ~15.0 
     let val = integrator.get_single_partial(&func, 2, &integration_limit, &point).unwrap();
-    std::println!("{}", val);
     assert!(f64::abs(val - 15.0) < 0.00001);
 }
 
@@ -76,27 +74,6 @@ fn test_booles_integration_3()
     assert!(f64::abs(val - 24.0) < 0.00001);
 }
 
-#[test]
-fn test_booles_integration_4()
-{
-    //equation is 2.0*x + y*z
-    let func = | args: &[f64; 2] | -> f64 
-    { 
-        return args[0]/(f64::sqrt(args[0]*args[0] + args[1]*args[1]));
-    };
-
-    let integration_limits = [[0.0, 1.0], [0.0, 1.0]];
-    let point = [1.0, 1.0];
-
-    let integrator = iterative_integration::MultiVariableSolver::from_parameters(20, IterativeMethod::Booles);
-
-    //double partial integration for first x then y, expect a value of ~1.50
-    let val = integrator.get(1, [0], &func, &integration_limits, &point).unwrap();
-    let expected_val = 0.414;
-    std::println!("{}", val - expected_val);
-    assert!(f64::abs(val - 1.50) < 0.00001);
-}
-
 #[test]
 fn test_gauss_legendre_quadrature_integration_1()
 {