approximate_pi.cpp

Go to the documentation of this file.

 
#include <cassert>   
#include <cstdlib>   
#include <iostream>  
#include <vector>    
 
namespace math {
 
using Point = struct {
    double x;
    double y;
};
 
double approximate_pi(const std::vector<Point> &pts) {
    double count = 0;  // Points in circle
    for (Point p : pts) {
        if ((p.x * p.x) + (p.y * p.y) <= 1) {
            count++;
        }
    }
    return 4.0 * count / static_cast<double>(pts.size());
}
}  // namespace math
 
static void tests() {
    std::vector<math::Point> rands;
    for (std::size_t i = 0; i < 100000; i++) {
        math::Point p;
        p.x = rand() / static_cast<double>(RAND_MAX);  // 0 <= x <= 1
        p.y = rand() / static_cast<double>(RAND_MAX);  // 0 <= y <= 1
        rands.push_back(p);
    }
    assert(math::approximate_pi(rands) > 3.135);
    assert(math::approximate_pi(rands) < 3.145);
 
    std::cout << "All tests have successfully passed!" << std::endl;
}
 
int main() {
    tests();  // run self-test implementations
    return 0;
}