epl692/pi
1
0
Fork 0
mirror of https://git.sdf.org/epl692/pi.git synced 2026-02-10 19:54:09 -05:00
Code Issues Packages Projects Releases Wiki Activity

Parallelize binary-splitting using rayon::join

Browse Source
This commit is contained in:
Eric Lynema
2025-12-20 23:23:56 -05:00
parent ed5c349796
commit 73d5325203
1 changed files with 4 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
src/main.rs
View File

@@ -1,5 +1,6 @@
use clap::Parser; use clap::Parser;
use rug::{Float, Integer, ops::Pow}; use rug::{Float, Integer, ops::Pow};
use rayon::join;
use std::fs::File; use std::fs::File;
use std::io::Write; use std::io::Write;
use std::path::PathBuf; use std::path::PathBuf;
@@ -40,8 +41,9 @@ fn bs(a: u64, b: u64) -> (Integer, Integer, Integer) {
return (p, q, t); return (p, q, t);
} }
let m = (a + b) / 2; let m = (a + b) / 2;
let (p1, q1, t1) = bs(a, m); let (left, right) = join(|| bs(a, m), || bs(m, b));
let (p2, q2, t2) = bs(m, b); let (p1, q1, t1) = left;
let (p2, q2, t2) = right;
let p = (&p1 * &p2).into(); let p = (&p1 * &p2).into();
let q = (&q1 * &q2).into(); let q = (&q1 * &q2).into();
let t1q2: Integer = (&t1 * &q2).into(); let t1q2: Integer = (&t1 * &q2).into();