python-grundlage/find-machinell-epsilon.py
@@ -1,13 +1,19 @@ # Credit: https://stackoverflow.com/a/57579379 import sys eps = 1.0 while eps + 1 > 1: eps /= 2 eps *= 2 sys_eps = sys.float_info.epsilon eps = 1.0 while 1.0 + eps > 1.0: eps /= 2.0 eps *= 2.0 print(f"The calculated epsilon: {eps}") print(f" The system epsilon: {sys_eps}") # with assign expression eps = 1.0 while 1.0 + (eps := eps / 2.0) > 1.0: pass eps *= 2.0 print(f"The calculated epsilon: {eps}") print(f" The system epsilon: {sys_eps}")