Simulate an approximation of PI value using Python

pi.py

import numpy as np
import multiprocessing as mp
import random
import math

threads = mp.cpu_count()
points = 100_000_000

def simulate_pi(points):
    radius = 1
    points_inside = 0
    for i in range(points):
        x = random.uniform(-radius, radius)
        y = random.uniform(-radius, radius)
        if math.sqrt(x*x + y*y) < radius:
            points_inside += 1

    return 4*points_inside/points

pool = mp.Pool(mp.cpu_count())
results = pool.map(simulate_pi, [round(points/threads) for i in range(threads)])
pool.close()

print(np.mean(results))