Merge pull request #2 from alphal00p/float_refactor

Float refactor

Author: GraDje13

Cargo.toml

@@ -23,9 +23,9 @@ lto = "fat"
 smallvec = "1.13.2"
 ahash = "0.8.11"
 num = { version = "0.4.1", features = ["serde"] }
-f128 = { git = "https://github.com/benruijl/f128" }
 statrs = "0.16.0"
 itertools = "0.12.1"
 rand = "0.8.5"
 criterion = "0.5.1"
 serde = { version = "1.0", features = ["derive"] }
+ref-ops = "0.2.5"

benches/gamma_bench.rs

@@ -1,13 +1,13 @@
 use criterion::{criterion_group, criterion_main, Criterion};
-use momtrop::gamma::inverse_gamma_lr;
+use momtrop::gamma::inverse_gamma_lr_impl;
 
 fn criterion_benchmark(c: &mut Criterion) {
     let mut group = c.benchmark_group("gamma sampling benchmarks");
 
     for omega in [0.5, 1.0, 2.0, 10.0] {
         for p in [0.1, 0.3, 0.5, 0.7, 0.9] {
-            group.bench_function(&format!("omega = {omega}, p = {p}"), |b| {
-                b.iter(|| inverse_gamma_lr(omega, p, 50, 5.0))
+            group.bench_function(format!("omega = {omega}, p = {p}"), |b| {
+                b.iter(|| inverse_gamma_lr_impl(omega, p, 50, 5.0))
             });
         }
     }

benches/graph_bench.rs

@@ -34,11 +34,15 @@ fn criterion_benchmark(c: &mut Criterion) {
     let logger = momtrop::log::DummyLogger {};
 
     let loop_signature = vec![vec![1]; 3];
-    let sampler = graph.build_sampler(loop_signature, 3).unwrap();
+    let sampler = graph.build_sampler(loop_signature).unwrap();
     let mut rng = rand::rngs::StdRng::seed_from_u64(69);
     let p1 = Vector::from_array([3.0, 4.0, 5.0]);
     let p2 = Vector::from_array([6.0, 7.0, 8.0]);
-    let edge_data = vec![(None, Vector::new()), (None, p1), (None, (&p1 + &p2))];
+    let edge_data = vec![
+        (None, Vector::new_from_num(&0.0)),
+        (None, p1),
+        (None, (&p1 + &p2)),
+    ];
 
     let x_space_point = vec![rng.r#gen(); sampler.get_dimension()];