diff --git a/Python/micropi.py b/Python/micropi.py new file mode 100755 index 0000000000000000000000000000000000000000..e522e2f24f9ac904e2cd70f0408aed6c77ab9258 --- /dev/null +++ b/Python/micropi.py @@ -0,0 +1,28 @@ +# +# micropi.py +# Neil Gershenfeld 12/26/22 +# calculation of pi by a MicroPython sum +# pi = 3.14159265358979323846 +# + +import time,machine + +NPTS = 1000000 + +@micropython.native +def calc(): + a = 0.5 + b = 0.75 + c = 0.25 + pi = 0 + for i in range(1,(NPTS+1)): + pi += a/((i-b)*(i-c)) + return pi + +machine.freq(250000000) +start_time = time.time() +pi = calc() +end_time = time.time() +mflops = NPTS*5.0/(1.0e6*(end_time-start_time)) +print("NPTS = %d, pi = %f"%(NPTS,pi)) +print("time = %f, estimated MFlops = %f"%(end_time-start_time,mflops)) diff --git a/README.md b/README.md index 58bb2789ee0c1171a058f94b0631215d88fda7ee..c2a03d52f6110833a36dea8ae304b4b788d90f48 100644 --- a/README.md +++ b/README.md @@ -54,7 +54,7 @@ |0.001|prior|Sun SPARCStation 1|C|prior| |0.001|prior|DEC VAX 8650|C|prior| |0.0007|prior|Intel 486|C|prior| -|0.00021|[pi.py](Python/pi.py)|MicroPython, 250 MHz|RP2040|December, 2022| +|0.00031|[micropi.py](Python/micropi.py)|MicroPython, 250 MHz|RP2040|December, 2022| |0.0002|[pi.ino](Arduino/pi.ino)|Arduino, floats|ATSAMD21E|December, 2020| |0.0001|[pi.ino](Arduino/pi.ino)|Arduino, floats|ATtiny1614|December, 2020| |0.00003|prior|Sun 3/60|C|prior|