def lambda_handler(event, context): start1 = time.time() integrate_f(0, 100, 1000) ts1 = time.time() - start1 print(ts1) def f(x): return x**2 - x def integrate_f2(a, b, N): s = 0 dx = (b - a) / N for i in range(N): s += f(a + i * dx) return s * dx start2 = time.time() integrate_f2(0, 100, 1000) ts2 = time.time() - start2 print(ts2) print("Cython gets {0:.2f}x speedup than pure Python!".format(ts2 / ts1))
from integrate import integrate_f print integrate_f(0.0,10.0,1000000)
def run_test(msg): print msg with Timer(): result = integrate_f(0.0, 10.0, 10000000) print "result:", result print
from integrate import integrate_f print(integrate_f(1, 2, 6))
from integrate import integrate_f print integrate_f(0.0, 10.0, 1000000)
#!/usr/bin/env python import argparse import sys sys.path.append("..") from integrate import integrate as integrate_f from integrate import f parser = argparse.ArgumentParser(description='integrator') parser.add_argument('a', nargs='?', type=float, default=0.0) parser.add_argument('b', nargs='?', type=float, default=10.0) parser.add_argument('N', nargs='?', type=int, default=int(10e+6)) args = parser.parse_args() a = args.a b = args.b N = args.N print "Numerical solution from (%(a)f to %(b)f with N=%(N)d : \n%(x)f" % \ {'a': a, 'b': b, 'N': N, 'x': integrate_f(f, a, b, N)}
import numpy as np import time import integrate outer_repeats = 2 inner_repeats = 500 n = 10000 dt_list = [] for i in range(outer_repeats): t0 = time.clock() for j in range(inner_repeats): integrate.integrate_f(.1, .3, n) dt = (time.clock() - t0) / inner_repeats # mean dt_list.append(dt) dt_min = np.min(dt_list) # min print '%.2f us' % (dt_min * 1e6) # Uncomment to raise exceptionxs #print integrate.integrate_f(0, 1, 10)
# import pyximport; pyximport.install() from integrate import integrate_f import time start1 = time.time() integrate_f(0, 100, 1000) ts1 = time.time() - start1 print(ts1) def f(x): return x**2 - x def integrate_f2(a, b, N): s = 0 dx = (b - a) / N for i in range(N): s += f(a + i * dx) return s * dx start2 = time.time() integrate_f2(0, 100, 1000) ts2 = time.time() - start2 print(ts2) print("Cython gets {0:.2f}x speedup than pure Python!".format(ts2 / ts1))
import numpy as np import time import integrate outer_repeats = 10 inner_repeats = 500 n = 10000 dt_list = [] for i in range(outer_repeats): t0 = time.clock() for j in range(inner_repeats): integrate.integrate_f(.1, .3, n) dt = (time.clock() - t0) / inner_repeats # mean dt_list.append(dt) dt_min = np.min(dt_list) # min print '%.2f us' % (dt_min * 1e6) print integrate.integrate_f(0, 1, 10)
#!/usr/bin/env python import sys sys.path.append("..") from integrate import integrate, f # integrate_f_with_functional_tools as integrate_f from integrate import integrate_f_with_functional_tools as integrate_f import argparse parser = argparse.ArgumentParser(description='integrator') parser.add_argument('a', nargs='?', type=float, default=0.0) parser.add_argument('b', nargs='?', type=float, default=10.0) parser.add_argument('N', nargs='?', type=int, default=10**6) args = parser.parse_args() a = args.a b = args.b N = args.N print "Numerical solution from (%(a)f to %(b)f with N=%(N)d : \n%(x)f" % \ {'a': a, 'b': b, 'N': N, 'x': integrate_f(a, b, N)} # print "%(x)f" % {'x': integrate(f, a, b, N)}
def run_test(msg): print msg with Timer(): result = integrate_f(0.0, 10.0, 1000000) print "result:", result print
#!/usr/bin/env python import argparse import sys sys.path.append("..") from integrate import integrate_f_with_functional_tools as integrate_f parser = argparse.ArgumentParser(description="integrator") parser.add_argument("a", nargs="?", type=float, default=0.0) parser.add_argument("b", nargs="?", type=float, default=10.0) parser.add_argument("N", nargs="?", type=int, default=10 ** 6) args = parser.parse_args() a = args.a b = args.b N = args.N print "Numerical solution from (%(a)f to %(b)f with N=%(N)d : \n%(x)f" % { "a": a, "b": b, "N": N, "x": integrate_f(a, b, N), }