Question

I have tried many algorithms for finding π using Monte Carlo.
One of the solutions (in Python) is this:

def calc_PI():
    n_points = 1000000
    hits = 0

    for i in range(1, n_points):
        x, y = uniform(0.0, 1.0), uniform(0.0, 1.0)

        if (x**2 + y**2) <= 1.0:
            hits += 1

    print "Calc2: PI result", 4.0 * float(hits) / n_points

The sad part is that even with 1000000000 the precision is VERY bad (3.141…).

Is this the maximum precision this method can offer?
The reason I choose Monte Carlo was that it’s very easy to break it in parallel parts.
Is there another algorithm for π that is easy to break into pieces and calculate?

Answer 1

This is a classic example of Monte Carlo. But if you’re trying to break the calculation of pi into parallel parts, why not just use an infinite series and let each core take a range, then sum the results as you go?

http://mathworld.wolfram.com/PiFormulas.html