Rust使用蒙特卡洛方法估算圆周率

使用Rust和蒙特卡洛方法估算圆周率

蒙特卡洛方法简介

蒙特卡洛方法是一种基于概率和随机抽样的算法，用于解决各种数学和物理问题。在估算圆周率的情况下，蒙特卡洛方法通过随机生成点并计算它们落在单位圆内的比例来工作。

Rust编程语言

Rust是一种系统编程语言，它提供了对性能和安全的保证。Rust的内存安全特性使其成为编写高性能并发程序的理想选择。

实现蒙特卡洛方法

rayon是一个用于并行迭代和处理数据的Rust库，它可以帮助你轻松地将迭代任务分配到多个CPU核心上，从而加速计算。
首先，你需要在Cargo.toml中添加依赖：rayon

1 2	[dependencies] rayon = "1.5.1"

然后，你可以使用的并行迭代功能来并行化圆周率的计算。以下是一个示例代码：

use rayon::prelude::*;
use std::time::Instant;

fn main() {
    let now = Instant::now();
    let total_iterations: usize = 1_000_000_000;
    let inside_circle: usize = (0..total_iterations)
        .into_par_iter()
        .filter(|_| {
            let x: f64 = rand::random();
            let y: f64 = rand::random();
            x * x + y * y <= 1.0
        })
        .count();

    let pi_estimate = (inside_circle as f64 / total_iterations as f64) * 4.0;
    println!("Estimated value of Pi: {}", pi_estimate);
    println!("Time taken: {:.2?}", now.elapsed());
}