#!/usr/bin/python3
import numpy as np
from timeit import default_timer as timer

def dot1(A,B):
	dot=0
	for i in range(len(A)):
		dot+=A[i]*B[i]
	return dot

def dot2(a,b):
	return a@b


def main():
#Normat Distribution N(0,1)
	a=np.random.randn(1000)
	b=np.random.randn(1000)
	T=1000
	A=list(a)
	B=list(b)
	start=timer()
	for t in range(T):
		dot1(A,B)
	end=timer()
	t1=end-start
	start=timer()
	for t in range(T):
		dot2(a,b)
	end=timer()
	t2=end-start
	print ("Time with lists:",t1)
	print ("Time with numpy arrays:",t2)
	print ("Ratio:",t1/t2)

if __name__=="__main__":
	main()